# Energyscapes and prey fields shape a North Atlantic seabird wintering hotspot under climate change

F. Amélineau, Jérôme Fort, D. Mathewson, D. Speirs, N. Courbin, Samuel

Perret, W. Porter, R. Wilson, David Grémillet

► To cite this version:

F. Amélineau, Jérôme Fort, D. Mathewson, D. Speirs, N. Courbin, et al. . Energyscapes and prey fields shape a North Atlantic seabird wintering hotspot under climate change. Royal Society Open Science, 2018, 5 (1), 10.1098/rsos.171883. hal-01774954

HAL Id: hal-01774954

https://hal.science/hal-01774954

Submitted on 31 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

Check for
updates

**Cite this article:** Amélineau F, Fort J, Mathewson PD, Speirs DC, Courbin N, Perret S, Porter WP, Wilson RJ, Grémillet D. 2018 Energyscapes and prey fields shape a North Atlantic seabird wintering hotspot under climate change. *R. Soc. open sci.* **5**: 171883.
http://dx.doi.org/10.1098/rsos.171883

Received: 14 November 2017

Accepted: 30 November 2017

**Subject Category:**

Biology (whole organism)

**Subject Areas:**

ecology/environmental science/physiology

**Keywords:**

bioenergetics, biologging, habitat modelling, little auk (*Alle alle*), migration, spatial ecology

**Author for correspondence:**

F. Amélineau

e-mail: francoise.amelineau@gmail.com

Electronic supplementary material is available
online at https://dx.doi.org/10.6084/m9.
figshare.c.3965346.

# Energyscapes and prey fields shape a North Atlantic seabird wintering hotspot under climate change

F. Amélineau<fnref n="1" />,<fnref n="2" />, J. Fort<fnref n="2" />, P. D. Mathewson<fnref n="3" />,
D. C. Speirs<fnref n="4" />, N. Courbin<fnref n="1" />, S. Perret<fnref n="1" />, W. P. Porter<fnref n="3" />,
R. J. Wilson<fnref n="4" /> and D. Grémillet<fnref n="1" />,<fnref n="5" />

<figure><img src="image_4.png" /></figure>

FA, 0000-0001-9723-5858; NC, 0000-0002-6428-5012;
DG, 0000-0002-7711-9398

There is an urgent need for a better understanding of animal
migratory ecology under the influence of climate change. Most
current analyses require long-term monitoring of populations
on the move, and shorter-term approaches are needed. Here,
we analysed the ecological drivers of seabird migration within
the framework of the energyscape concept, which we defined
as the variations in the energy requirements of an organism
across geographical space as a function of environmental
conditions. We compared the winter location of seabirds with
their modelled energy requirements and prey fields throughout
the North Atlantic. Across six winters, we tracked the migration
of 94 little auks (*Alle alle*), a key sentinel Arctic species, between
their East Greenland breeding site and wintering areas off
Newfoundland. Winter energyscapes were modelled with
Niche Mapper™, a mechanistic tool which takes into account
local climate and bird ecophysiology. Subsequently, we used
a resource selection function to explain seabird distributions
through modelled energyscapes and winter surface
distribution of one of their main prey, *Calanus finmarchicus*.
Finally, future energyscapes were calculated according to
IPCC climate change scenarios. We found that little auks
targeted areas with high prey densities and moderately
elevated energyscapes. Predicted energyscapes for 2050 and

<footnote n="1">¹CEFE UMR 5175, CNRS – Université de Montpellier – Université Paul-Valéry Montpellier – EPHE, Montpellier, France</footnote>

<footnote n="2">²Littoral Environnement et Sociétés (LIENSs), UMR 7266 CNRS-Université de La Rochelle, La Rochelle, France</footnote>

<footnote n="3">³Department of Integrative Biology, University of Wisconsin, Madison, WI, USA</footnote>

<footnote n="4">⁴Department of Mathematics and Statistics, University of Strathclyde, Livingstone Tower, 26 Richmond Street, Glasgow G1 1XQ, Scotland, UK</footnote>

<footnote n="5">⁵Percy FitzPatrick Institute, DST/NRF Centre of Excellence, University of Cape Town, Rondebosch, South Africa</footnote>
2095 showed a decrease in winter energy requirements under the high emission scenario, which may be beneficial if prey availability is maintained. Overall, our study demonstrates the great potential of the energyscape concept for the study of animal spatial ecology, in particular in the context of global change.

# 1. Introduction

Theoretical and empirical studies have demonstrated that migration evolves to maximize fitness in a seasonal environment [1,2]. Migratory animals thereby target the most profitable areas [3], balancing their energy requirements with available resources, within ever-changing ecological landscapes. Such complex spatio-temporal match/mismatch of migratory species with their biotic and abiotic environments shape migratory dynamics and the fate of populations on the move [3,4]. Migrating animals are thus particularly vulnerable to climate change and resulting environmental modifications [3,5,6]. Among birds, most long-distance migrants breed at high latitudes, where climatic changes have the strongest amplitude and seasonally travel to more favourable wintering grounds [7,8]. Overall, this puts a strong emphasis on the migratory ecology of birds in a warming Arctic [9].

Bird sensitivity to climate change has been mainly studied during the breeding season. In particular, drastic changes in phenology have been noticed worldwide, triggered by shorter winters in polar and temperate regions [9–11]. Such phenological changes, that differ between species, can lead to a mismatch between food availability and demands for reproduction with strong impacts on breeding success, particularly for migrating species [5,10]. However, the impacts of climatic changes during winter have been the focus of fewer studies [6,12], probably because it is far more challenging to monitor individuals outside the breeding season [3].

A prerequisite to understanding how wintering animals are affected by climate change is to quantify the ecological benefits of seasonal habitat choice. For this purpose, it is essential to define and use ecological metrics that will allow researchers to rate and compare the profitability of wintering areas. According to evolutionary theory, these fitness proxies should show functional relationships with the capacity of each individual to survive and reproduce [13]. Indeed, a series of studies have identified the impact of wintering conditions on adult survival probabilities [14,15], while others have identified significant carry-over effects of such conditions on survival and reproduction in subsequent months [16], or even years [17].

Despite the great success and the necessity of such studies, they require long-term population monitoring and individual phenotyping, which are both extremely time-consuming and costly, particularly for long-lived species, like many migratory birds. Animal energetics offer a powerful, short-term alternative [18]. In particular, within the rapidly emerging field of movement ecology [19], the concept of energy landscapes (hereafter ‘energyscapes’) seems extremely relevant and attractive for mechanistic explanations of animals’ geographical distribution. Specifically, Wilson, Shepard and collaborators defined energyscapes as ‘environmentally dependent variation in the cost of transport, driven by variation in parameters such as incline, substrate type, vegetation, current speed, or direction’ [20,21]. Here, we propose to broaden the concept of energyscape and define it as variation in the energy requirements of an organism across geographical space as a function of environmental conditions. In addition to the cost of transport, our definition thereby includes all costs associated with body maintenance and thermoregulation. Energyscapes are, consequently, highly sensitive to climatic conditions.

The recent development of mechanistic energetic models allows energyscapes to be estimated for any animal at any location around the globe [22–26]. This methodology, which in birds compares well with the accuracy of empirical measurements [27,28], provides a unique opportunity to develop and expand the concept of energyscapes for the study of animal migration.

In this study, we determined monthly energyscapes for an Arctic migrant, the little auk (Alle alle), and used them to test hypotheses related to wintering ecology in a climate change context. Little auks are the most numerous seabird of the Arctic (with an overall population estimated between 40 and 80 million individuals) and they are key components of Arctic marine food webs. Their winter migration across the North Atlantic has only recently been depicted through the use of miniaturized electronic geolocators [24]. Because the little auk is the smallest of all seabirds in the Atlantic Arctic, with a body mass of 150 g, its morphology makes it a particularly attractive model for the study of avian energetics in a migratory and wintering context. Indeed, its mass-specific resting metabolic rate is approximately
seven times that of Emperor penguins (*Aptenodytes forsteri*) [29], and therefore energetic constraints
acting upon little auk winter migratory decisions are predicted to be drastic. Finally, little auks feed
on zooplankton, especially Calanoid copepods, which show a strong response to the climate-induced
increase of North Atlantic surface water temperatures [30]. It is therefore an excellent model species for
the study and forecast of the marine ecological consequences of rapid warming in the North Atlantic.

We took advantage of a large dataset on little auk migratory movements in the North Atlantic. This information was collected over 6 years using geolocators for a little auk population which breeds in East Greenland and predominantly overwinters off Newfoundland. Once the location of this wintering hotspot had been identified, we used the mechanistic model Niche Mapper™ and remote-sensing information to assess the energyscapes of birds within the North Atlantic. Then, we used resource selection functions (RSFs) to compare little auk spatial distribution with the modelled energyscapes and abundance of zooplankton prey. This approach allowed us to characterize the energetic strategy of wintering little auks, and to go further than the correlative approach of [24] in the definition of ecological drivers of their distribution. Finally, using climate forecasting models, we modelled the future energyscapes of little auks in their wintering region, to assess the potential impacts of forthcoming North Atlantic climate warming. As temperature is predicted to be the main driver of little auk energy requirements during winter according to [22], we expected that little auk energy demands would consequently decrease due to this warming.

We tested the hypotheses that (i) little auks optimize the position of their wintering location, so as to minimize their winter daily energy requirements, and maximize their use of zooplankton density gradients; and (ii) climate change will significantly modify little auk energy requirements at their current wintering location.

## 2. Methods

### 2.1. Seabird winter geolocation

Fieldwork took place at Ukaleqarteq (Kap Höegh, Liverpool Land, 70°44′ N, 21°35′ W), East Greenland. Breeding adults were equipped with light-level archival tags (GLS) each summer from 2009 to 2014 and recaptured the following years. All GLSs weighed between 0.8 and 1.5 g (0.6–1.1% of the lowest equipped bird weight). GLS types, technical characteristics and deployment/recapture details are available in electronic supplementary material, file S1. In total 244 GLSs were deployed, 102 were retrieved and 94 data files were exploitable. Three birds were equipped for 2 years, and one bird for 3 years. Birds were captured either in their underground nests or with a lasso placed on the rocks surrounding their nests. They were weighed and fitted with a metal ring on which a GLS was attached using a cable tie. Previous investigations showed that GLS deployments had no measurable impacts on little auk body condition [24].

GLS data were analysed with British Antarctic Survey (BAS) softwares (TransEdit and Locator) for BAS and Biotrack devices, and with Intiproc® for Migrate Technology devices. In both cases, the light intensity threshold was set to 10 lux, and the sun elevation angle to -3°. Two positions per day were obtained for each bird, one at local noon and one at local midnight. Points on land and outside the study area (70° W–15° E, 30° N–80° N) were removed. All bird handling procedures were approved by the Government of Greenland (Permits nos. 66.01.13, 2011–047447, 2012-065815, 2013-083634 and 2014-098814) and validated by the ethics committee of the French Polar Institute (Permit no. MP/53/06/12). Tracking data are available on Movebank (https://www.movebank.org/) under the study named 'Adaclim'.

### 2.2. Modelling seabird energyscapes

To model seabird energy requirements at a point in time and space, we used the mechanistic model Niche Mapper™, which evaluates the daily energy requirements of an individual using the biophysical properties of seabird bodies exposed to specific microclimatic conditions [22,31]. Niche Mapper consists of two submodels: a microclimate model and an animal model. The microclimate model uses macroclimate data (sea surface temperature (SST), air temperature, cloud cover, relative humidity and wind speed), substrate properties, geographic location and time of year to calculate hourly environmental conditions at the animal's height. Microclimate model calculations are detailed in [32]. The animal model uses the outputs from the microclimate model to iteratively solve a heat balance equation to find the metabolic rate needed for the animal to maintain its body temperature, accounting
for convective, radiative, evaporative and solar heat fluxes with its microenvironment. Heat balance calculations are detailed in [33]. Little auks were modelled as a series of simple shapes that have well-understood heat transfer properties that enable surface temperature calculations—and thus heat flux calculations—given a certain core temperature to maintain: two ellipsoids for the head and torso and two cylinders for the featherless parts of the legs. Legs were modelled to allow heat loss, with possible changes in peripheral blood flow as a function of bird heat balance and external temperature. For time not spent diving or flying, little auks were modelled as floating on the ocean surface with legs and 25% of torso submerged in the water (electronic supplementary material, file S2). Auks were modelled as diving 24% of each day and flying 9% of each day (electronic supplementary material, file S2, [22]). All model input values are detailed in electronic supplementary material, file S2. Monthly average environmental input data were downloaded from the International Comprehensive Ocean-Atmosphere Data Set (ICOADS, http://icoads.noaa.gov/, observed data, 1° × 1° resolution). The model was run to predict little auk daily energy requirements during the winter months (November, December, January and February) between 30°N to 80°N and 70°W to 15°E with a 1° × 1° grid size. A sensitivity analysis was performed to identify input variables which had the strongest influence upon modelled energy requirements (electronic supplementary material, file S3).

## 2.3. Zooplankton winter abundance

We used the distribution and abundance of the copepod *Calanus finmarchicus*, one of little auk main prey, as a proxy for little auk prey availability [24,34]. Since temporally continuous observational data on *C. finmarchicus* are not available on the spatial scale required for our study, we used abundance estimates obtained from the ocean-scale population model of [35]. This is an explicitly spatial model (at 0.5° × 0.25° grid size) where at each discrete location the population is divided into surface (0–100 m depth) individuals and deep diapausing individuals. The surface population is further divided into discrete development classes that map on to the naupliar and copepodite stages. Development rate and egg production depend on the local temperature and food availability. A proportion of individuals entering the pre-adult (CV) copepodite stage join the diapausing population and emerge as surface adults (CVI) in the spring. Temperature and spatial transport for the model was determined using outputs from the Ocean Circulation and Climate Advanced Modelling Project [36] ocean circulation model, while phytoplankton food was estimated from SeaWIFS satellite observations (https://oceancolor.gsfc.nasa.gov/data/seawifs/). Model outputs have been successfully compared with field data from the continuous plankton recorder, *in situ* winter distributions of diapausers and, at some locations, copepod time series. A full formal description of the approach is given in [35]. The spatial extent of the model encompasses our entire study area, as defined in the previous section, and we used the modelled abundances of the surface CV and CVI copepodite stages as our proxy for little auk prey. The current version of this model does not yet allow forward projections of *C. finmarchicus* abundance due to a lack of reliable projections of phytoplankton abundances.

## 2.4. Habitat selection of little auk

We determined how little auks balanced the cost of their daily energy requirements with the benefits of searching for *C. finmarchicus* for each winter month (November to February). We used RSFs [37] that assessed whether a given habitat feature is used disproportionately relative to availability (i.e. selection or avoidance). RSFs compared the daily energy requirements of little auks and the density of *C. finmarchicus* at GLS locations with those expected at an equal set of random locations generated within the 95% kernel density of observed locations for each month (figure 1). This kernel density contour was chosen to represent available locations (both favourable and less favourable wintering locations). RSFs were fitted using generalized linear mixed models with a binomial distribution for errors and a logit link. We added a random intercept to account for the unbalanced number of locations collected across individuals [38]. We assessed empirical standard errors that are robust to both among- and within-individual correlations (i.e. serial correlation) and that provide robust estimates of significance [39]. The RSF took the form $w(\mathbf{x}_{ij}) = \exp(\beta_0 + \beta_1 x_{1ij} + \beta_2 x_{2ij} + \dots + \beta_n x_{nij} + \gamma_{0j})$, where $w(\mathbf{x}_{ij})$ is the relative probability of selection for little auk, $\mathbf{x}$ is a vector, $\beta_0$ is the mean intercept, $\beta$ is the estimated fixed regression coefficient for continuous covariate $x$, $i$ represents the $i$th observation, $j$ represents the $j$th individual and $\gamma_{0j}$ is the random effect on the intercept $\beta_0$ for animal $j$. We included the squared energy requirements and the squared density of *C. finmarchicus* to allow for quadratic effects, and an interaction
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Little auk winter distribution (number of individuals per 1° × 1° cell, green), little auk energyscape (kJ d⁻¹, red), little auk prey distribution (log₁₀ density (m⁻²)) and relative probability of habitat selection by little auks for each month in the North Atlantic (brown). Monthly 95% kernel contours of little auk positions are presented in red. The position of the colony (red plus symbol) and the 50% kernel of winter positions (green) are presented in the top-left panel. White areas represent an absence of data.</figcaption></figure>

term between prey density and energy requirements. Both the daily energy requirements and density of *C. finmarchicus* were centred to avoid collinearity issue and multicollinearity was low in all RSFs (variance inflation factors ≤6.6 in all models, [40]). We evaluated model robustness using k-fold cross-validation, by developing RSFs with 80% of the locations (training set), and then by testing the predictive power of these RSFs with the 20% withheld locations (testing set) [39,41]. RSFs were performed with the GLIMMIX procedure of SAS v. 9.2 software (SAS Inst.).

## 2.5. Energyscape projections

We investigated changes in energy requirements within the core wintering area of little auks (defined as the 50% kernel density contour of the positions from 1 November to 28 February) for three decades centred in 2010, 2050 and 2095, using climatic projections based on two IPCC (Intergovernmental Panel on Climate Change) scenarios: one reflecting a low greenhouse gas concentration trajectory (Representative Concentration Pathway, RCP 2.6) and one reflecting a high concentration trajectory (RCP 8.5). Based on model comparisons by [42] for our studied area, we chose the Canadian Earth System Model v. 2 (CanESM2), which reproduced the best observed climatic trends off Newfoundland. Climatic predictions for sea-surface temperatures, air temperatures, cloud cover and relative humidity were downloaded from http://www.cccma.ec.gc.ca/data/cgcm4/index.shtml. As Niche Mapper™ requires a minimum and a maximum value for each input variable, and only mean values were available from CanESM2, we calculated min/max from the observed amplitude within the same variables from ICOADS data. Niche Mapper™ was run for each year of the three decades using both scenarios (RCP 2.6 and 8.5). Then, for each month and each run, we calculated the mean daily energy requirements of little auks within their core wintering area (figure 4, core wintering area presented in the top-left panel of figure 1, [43]).
# 3. Results

## 3.1. Little auk winter distribution and timing of migration

Figure 1 presents the density of bird positions per 1° × 1° cell for each winter month (wintering period). On figure 2, bird positions were summed by longitude and latitude for each non-breeding month (migration + wintering periods). Over the study period (2009–2015), GLS recordings showed that little auks migrated to their wintering grounds in October (figure 2). Birds wintered off Newfoundland, between 40° N and 55° N and 35° W and 55° W, at nearly 3000 km from their colony. They remained within this area until February, and migrated towards their breeding grounds in March (figure 2). Figure 2 shows a clear overlap between little auk distribution and peaks of *C. finmarchicus* abundance over the winter months. Hence, the timing of migration matched the narrowing (October) and broadening (March) of the *C. finmarchicus* distribution in the North Atlantic.

## 3.2. Little auk energyscape

The sensitivity analysis revealed that little auk daily energy requirements are mainly driven by SSTs and air temperatures (electronic supplementary material, file S3). Little auk daily energy requirements in the North Atlantic followed a latitudinal gradient, with higher energy requirements at higher latitudes all year round (figures 1 and 2b). During winter, daily energy requirements were highest around Newfoundland, within the cold Labrador Current, and increased gradually from December to March. Visually, this peak of seabird energyscape overlapped with the highest predicted densities of *C. finmarchicus*, especially in January and February, where the prey spatial distributions were less widespread (figure 1). Consequently, in January and February, the little auk distribution became slightly uncoupled from the prey distributions, birds remaining in areas with lower energy requirements in the southeast of prey distribution (figures 1 and 2b).

## 3.3. Little auk habitat selection

All RSFs were robust to k-fold cross validation and had a high power to predict spatial distribution of little auk (table 1). From November to January, little auks had a higher likelihood of selecting areas with high *C. finmarchicus* densities and moderate levels of daily energy requirements (figure 3a-c). In December, little auks were more likely to occur where their energy requirements were around 260 kJ d⁻¹. In November and January, they experienced increased daily energy requirements, around 340 kJ d⁻¹. As the daily energy requirements progressively increased, and the prey spatial distribution strongly overlapped high energyscapes around Newfoundland in January and February (figures 1 and 2), little auk made a trade-off between the benefit of food resource and energy requirements. In January, little auks also selected in a lesser extent areas with low density of *C. finmarchicus* (figure 3c). The trade-off was still stronger in February, little auks being more likely to select areas with high daily energy requirements (figures 1–3d).

## 3.4. Predictions

Under the low emission scenario, projections of little auk energy requirements did not differ with time, except in January where a slight difference between 2010 and 2095 was found (figure 4, table 2, November: *F*<sub>2,27</sub> = 2.69, *p* = 0.086, December: *F*<sub>2,27</sub> = 0.49, *p* = 0.61, January: *F*<sub>2,27</sub> = 3.7, *p* = 0.038, February: *F*<sub>2,27</sub> = 2.88, *p* = 0.074). Under the high emission scenario, there was a significant decrease in energy requirements with time, for each month (figure 4, table 2, November: *F*<sub>2,27</sub> = 63.5, December: *F*<sub>2,27</sub> = 99.1, January: *F*<sub>2,27</sub> = 42.8, February: *F*<sub>2,27</sub> = 46.3, all *p* < 0.0001).

# 4. Discussion

Our detailed analysis of little auk wintering energyscapes and prey fields in a climate change context, support our first hypothesis: the wintering grounds of little auks were shaped by both prey availability and energyscapes. Yet, prey fields seemed more important than energy requirements to explain bird winter locations. Our results also partially support our second hypothesis: forthcoming climate warming will substantially modify energyscapes for little auks wintering off Newfoundland under a high emission scenario only. Overall, our study highlights the relevance of the emerging energyscape concept as a
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Number of little auk GLS positions (black), mean *Calanus finmarchicus* density in log₁₀ density m⁻² (red, dashed line) and little auk mean daily energy requirements in kJ d⁻¹ (blue, dotted line) per longitude (a) and latitude (b) and per month. The numbers of positions are low in August and May because most of the birds are already above the polar circle and permanent daylight does not allow positioning using GLS recorders. Similarly, bird position is not indicated during the summer (June and July), but their position was known as they were breeding at our study site in East Greenland.</figcaption></figure>
<table><caption>**Table 1.** Generalized linear mixed effect models for habitat selection for little auk (n = 94), with parameter estimates (β) and their 95% confidence interval (CI) for each winter month. *Calanus finmarchicus* density (Cfin, log₁₀ density m⁻²) and energy requirements (ER) were centred. 95% CI exclude 0 at α < 0.05 (*), < 0.01 (**), < 0.001 (****).</caption><thead><tr><th></th><th>November β 95% CI</th><th>December β 95% CI</th><th>January β 95% CI</th><th>February β 95% CI</th></tr></thead><tbody><tr><td>intercept</td><td>0.073 ± 0.088</td><td>0.31 ± 0.077***</td><td>0.32 ± 0.11**</td><td>0.36 ± 0.11**</td></tr><tr><td>Cfin</td><td>0.19 ± 0.031***</td><td>0.20 ± 0.030***</td><td>0.081 ± 0.029**</td><td>0.10 ± 0.046*</td></tr><tr><td>Cfin²</td><td>-5.8 × 10⁻³ ± 1.8 × 10⁻³**</td><td>-7.7 × 10⁻³ ± 2.1 × 10⁻³***</td><td>-5.3 × 10⁻³ ± 1.4 × 10⁻³***</td><td>-5.7 × 10⁻³ ± 2.3 × 10⁻³*</td></tr><tr><td>ER</td><td>0.017 ± 7.5 × 10⁻³*</td><td>-4.3 × 10⁻³ ± 6.1 × 10⁻³</td><td>3.0 × 10⁻³ ± 4.2 × 10⁻³</td><td>0.016 ± 3.9 × 10⁻³***</td></tr><tr><td>ER²</td><td>-2.0 × 10⁻³ ± 8.7 × 10⁻⁵*</td><td>-4.0 × 10⁻⁵ ± 4.1 × 10⁻⁵</td><td>-7.0 × 10⁻³ ± 4.3 × 10⁻⁵</td><td>-9.0 × 10⁻³ ± 3.3 × 10⁻³**</td></tr><tr><td>Cfin * ER</td><td>-2.0 × 10⁻⁴ ± 1.4 × 10⁻³</td><td>-1.8 × 10⁻³ ± 6.4 × 10⁻⁴**</td><td>-1.1 × 10⁻⁴ ± 2.6 × 10⁻⁴</td><td>-1.0 × 10⁻³ ± 2.7 × 10⁻⁴***</td></tr><tr><td>k-fold (r<sub>s</sub>)</td><td>0.93 ± 0.039</td><td>0.87 ± 0.058</td><td>0.90 ± 0.046</td><td>0.91 ± 0.041</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Relative probability of selection for little auk as a function of their daily energy requirements, and *C. finmarchicus* density modelled for each during the wintering period. Little auk daily energy requirements are only represented in the range of observed values for each month.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** 10-year average little auk energetic requirements in 2010, 2050 and 2095 for two emission scenarios. Values are means ± s.d. Energetic requirements are calculated for little auk core wintering areas, defined as the 50% kernel density contour of the GLS positions obtained each year (2009–2015) from 1 November to 28 February (represented in top-left panel of figure 1). Low emission = RCP 2.6, black, circles. High emission = RCP 8.5, red, triangles. For each scenario and month, means with different letters are significantly different (Tukey post hoc test, table 2). No difference was found for RCP 2.6 in November, December and February (table 2).</figcaption></figure>

methodological framework for the study of animal migration ecology and evolution in the context of global change [44]. Nevertheless, our results also demonstrate the need to include prey fields and their energetic profitability in future energyscape studies.

## 4.1. Seabird wintering in the Northwest Atlantic

In this study, we took advantage of existing, detailed knowledge of Calanoid copepod distribution and habitat models in the North Atlantic [35,45]. While the importance of *Calanus finmarchicus* in the winter diet of little auks is still a matter of debate, we considered that *C. finmarchicus* density and distribution was a good proxy for overall little auk prey availability, as they fed either on those prey, or on species from the next trophic level, such as amphipods and krill when they wintered off Newfoundland [31,46,47]. Interestingly, little auk timing of migration correlated closely with *C. finmarchicus* range restriction in autumn, and range expansion in spring (figure 2). This observation also strongly suggests a tight link between little auks and *C. finmarchicus* off Newfoundland in winter.

Little auks did not winter in areas where their energyscape was the most favourable, i.e. the East Atlantic (figure 1). In fact, they were constrained because highly productive areas were also located where their energy requirements are higher, within and around the cold Labrador Current off Newfoundland (figures 1 and 2). Within this area, they fine-tuned their spatial distribution in respect to the energyscape, targeting the southeastern range of *C. finmarchicus* during winter, where prey items were still abundant and bird energyscapes were more advantageous (figures 1 and 2). They seemed to avoid shelf waters,
<table><caption>**Table 2.** One-way ANOVAs and Tukey post hoc tests comparing mean energy requirements per decade, for each scenario and each month. d.f., degrees of freedom; SS, sum of squares; MS, mean square.</caption><thead><tr><th rowspan="2">scenario</th><th rowspan="2">month</th><th colspan="6">ANOVA</th><th colspan="3">Tukey post hoc test, adjusted p-values</th></tr><tr><th>factor</th><th>d.f.</th><th>SS</th><th>MS</th><th>F</th><th>p</th><th>2010-2050</th><th>2050-2095</th><th>2010-2095</th></tr></thead><tbody><tr><td rowspan="6">low emission<br>RCP 2.6</td><td rowspan="2">November</td><td>decade</td><td>2</td><td>26.91</td><td>13.46</td><td>2.69</td><td>0.086</td><td></td><td></td><td></td></tr><tr><td>residuals</td><td>27</td><td>134.94</td><td>5</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">December</td><td>decade</td><td>2</td><td>6.57</td><td>3.28</td><td>0.49</td><td>0.61</td><td></td><td></td><td></td></tr><tr><td>residuals</td><td>27</td><td>180.07</td><td>6.67</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">January</td><td>decade</td><td>2</td><td>189.4</td><td>94.68</td><td>3.7</td><td>0.038</td><td>0.2</td><td>0.64</td><td>0.033</td></tr><tr><td>residuals</td><td>27</td><td>691.3</td><td>25.6</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">February</td><td>decade</td><td>2</td><td>531.7</td><td>265.86</td><td>2.88</td><td>0.074</td><td></td><td></td><td></td></tr><tr><td>residuals</td><td>27</td><td>2497</td><td>92.48</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="6">high emission<br>RCP 8.5</td><td rowspan="2">November</td><td>decade</td><td>2</td><td>602.6</td><td>301.3</td><td>63.5</td><td>&lt;0.0001</td><td>0.0004</td><td>&lt;0.0001</td><td>&lt;0.0001</td></tr><tr><td>residuals</td><td>27</td><td>128.1</td><td>4.74</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">December</td><td>decade</td><td>2</td><td>709.2</td><td>354.6</td><td>99.05</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td></tr><tr><td>residuals</td><td>27</td><td>96.7</td><td>3.6</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">January</td><td>decade</td><td>2</td><td>1336.8</td><td>668.4</td><td>42.84</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0.0008</td><td>&lt;0.0001</td></tr><tr><td>residuals</td><td>27</td><td>421.3</td><td>15.6</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">February</td><td>decade</td><td>2</td><td>3806</td><td>1903.2</td><td>46.25</td><td>&lt;0.0001</td><td>0.0001</td><td>0.0009</td><td>&lt;0.0001</td></tr><tr><td>residuals</td><td>27</td><td>1111</td><td>41.1</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

potentially because at these locations cold water stemming from the Labrador Current [48] increased their energy requirements. As in summer, they also targeted the edges of the continental shelf such as around the Grand Banks, where upwelling concentrates prey [49]. A second wintering strategy was also observed: according to the RSF, there was a medium (January) and high (February) probability of little auk occurrence with medium energy requirements when there was little or no *C. finmarchicus* (figure 3c,d), which suggests that some of the little auks from East Greenland could rely on different prey/food web at the end of winter.

Little auk wintering areas are also crucial for many other seabird species [50–57]. Indeed, the Northwest Atlantic is a major wintering hotspot for seabirds from various breeding grounds [50–53], and even for birds from the Southern Hemisphere during the boreal summer [57]. It is also a migration stopover for some long-distance migrants [55,56], and Southwest Greenland itself hosts seabirds from the whole Arctic [54]. Therefore, the energyscape concept illustrated here, as well as associated analytical tools, will also be of great use for the general ecological understanding of the evolution of aquatic bird migration in the North Atlantic, in a global change context.

## 4.2. Wintering strategies under forecasted climate change

Under the high emission scenario (RCP 8.5), little auk energy requirements should decrease during the twenty-first century within their current core wintering areas (figure 4, table 2). Under the low emission scenario (RCP 2.6), their energy requirements are predicted to decrease only in January (table 2, figure 4). The RCP 2.6 represents a high mitigation scenario that aims to keep global warming to less than 2°C above pre-industrial temperatures [7]. A decrease in energy requirements during winter is therefore likely for little auks, and could be beneficial for them as well as for other wintering seabirds. In order to anticipate the future distributions of wintering migrants, it is also necessary to know how their prey will react to ongoing environmental changes, and two scenarios can be considered: (1) if prey biomass remains constant, winter visitors could stay in their current wintering areas, especially if those become energetically more profitable, due to climate warming. However, (2) if prey abundance decreases in this area, they may have to move to follow their current prey distribution and/or target different prey species. It is therefore imperative to understand whether little auk wintering areas off Newfoundland
will remain highly productive during winter, thereby supporting scenario (1), or if changes in ocean
circulation may modify primary productivity and associated food webs. Higher temperatures are known
to favour smaller zooplankton species [58,59]. Therefore, prey biomass has to be maintained to sustain
the rest of the food web, including migratory bird populations. Yet, over the past 40 years, a decrease
in total zooplankton biomass has been observed in the Northeast Atlantic and the same phenomenon is
predicted to occur in the Northwest Atlantic [60].

Thanks to the continuous plankton recorder survey [61], the spatio-temporal evolution of zooplankton species has been well studied in the North Atlantic during the last 50 years. A northward shift in zooplankton assemblages and in the copepod *Calanus finmarchicus* has been demonstrated for the Northeast Atlantic [45,62] as well as for the Northwest Atlantic [63], whereas no changes have been found so far for the Arctic species *C. glacialis* and *C. hyperboreus* that form an important food link at higher latitudes [62]. The authors speculate that this may be due to the southward penetration of water from the Labrador current into the Northwest Atlantic [62]. These studies suggest that seabird visitors may have to move northwards to target their preferred prey during winter. It is, however, not clear whether prey will move further north in the Labrador Sea. Recent climatic models suggest a cooling of the Labrador Sea due to a local collapse of deep-ocean convection [64] that would prevent *Calanus finmarchicus* range expansion.

Currently, seabirds (including little auks) that winter at higher latitudes in the Labrador Sea may rely on different prey communities [31,47,65,66]. Indeed, energy requirements are higher further north [22], but prey of these areas are richer in lipids [67,68]. In particular, copepods from colder waters are bigger and richer in lipids [63]. Changes in winter distribution will be facilitated in species that show a low individual consistency in their wintering areas. While seabirds were until recently seen as highly consistent in their wintering locations [69], this understanding is changing [70], and among Alcids in particular, individual consistency varies strongly between species, suggesting some plasticity in wintering ecology [71,72]. In this study, preliminary results on little auks equipped with GLSs during two or three successive years confirm that they are not always wintering at the same place (electronic supplementary material, file S4).

Finally, a northward shift of little auk wintering grounds could be beneficial because it reduces overall migration distance. In the case of little auks from East Greenland, wintering in the Labrador Sea instead of the area around the Grand Banks of Newfoundland would reduce migration distance by about 1000 km, corresponding to 33% of their current migration distance. Such a shortening in migration distance has been observed in some terrestrial bird species [73,74], and an experimental study has even shown that residency can evolve rapidly in a passerine population if selection pressure towards shorter migration is maintained [75]. Migrating closer to the breeding area is not only beneficial because of a decrease in travel costs, but also because it allows a better detection of the environmental conditions occurring at the breeding grounds, so as to match resource phenology at the breeding site [74].

However, one potential barrier to the poleward shift of little auk wintering areas is the decrease
in daylight duration and the polar night, but this does not seem to directly impact Arctic seabirds, as
individual seabirds, including little auks, have been found to actively forage during the polar night in
Greenland, Spitsbergen or the Barents Sea ([76,77], Fort J et al. unpublished). Yet the polar night also limits
winter primary productivity at high latitudes, as well as the carrying capacity of the coastal ecosystem
with respect to apex predators. In this context, poleward range expansions may also enhance inter- and
intra-specific competition [78,79]. In our study system, little auks from other populations are already
wintering further north in the Labrador Sea, and could become competitors if the whole wintering area
of this species is narrowed [80]. Lastly, a poleward shift may result in a mismatch between prey blooms
and the timing of migration, as phytoplankton blooms occur later at higher latitudes because of light
restriction during winter [81].

4.3. Outlook

With the present study, we have broadened the concept of energy landscapes as defined by Wilson, Shepard and collaborators [20,21], and redefined it as the spatial variations in energy requirements of an animal at a specific moment in time. Therefore, in contrast to [21], this approach does not require recording detailed accelerometry data to estimate transport cost for each study individual, but rather some knowledge of environmental conditions encountered through time, and on species time-budgets and metabolism. Energyscapes calculated with Niche Mapper™ can consequently be used in a much wider range of species and ecological contexts [82,83]. These two approaches could nonetheless be used complementarily if accelerometry data help fine-tune NicheMapper™ input parameters.
Our study of little auk wintering ecology has also demonstrated the overarching importance of prey fields, for a thorough understanding of individual strategies and population biogeography under climate change. Therefore, once functional relationships between prey availability and predatory performance become known for little augs [84], we propose further expanding the energyscape concept to include this information. Energyscapes would then be defined as the energetic profitability, for a given species at a given time, thereby setting an exciting target for future investigations of species biogeography in a changing world.

Ethics. All bird handling procedures were approved by the Government of Greenland (Permit nos. 66.01.13, 2011-047447, 2012-065815, 2013-083634 and 2014-098814) and validated by the ethics committee of the French Polar Institute (Permit no. MP/53/06/12).

Data accessibility. Our data are deposited on Movebank (https://www.movebank.org/) under the study named 'Adaclim'.

Authors' contributions. D.C., J.F. and F.A. designed the study, collected fieldwork data on birds and wrote the first draft of the manuscript. P.D.M., D.C.S., N.C., W.P.P., R.J.W., J.F. and F.A. performed the analyses. All authors contributed to interpreting the data and revising the manuscript.

Competing Interests. We have no competing interests.

Funding. This study was funded by the French Polar Institute Paul-Emile Victor, through the ADACLIM programme (no. 388) led by D.G. and J.F. This study also represents a contribution to the ARCTOX research project (European commission Marie Curie CIG, Project 631203 to J.F.).

Acknowledgements. We thank Nanu Travel for the logistical support, Cyril Bernard and Xavier Morin for their help in the analyses and all the fieldworkers involved in this study.

## References

1. Alerstam T, Hedenström A, Åkesson S. 2003 Long-distance migration: evolution and determinants. *Oikos* **103**, 247–260. (doi:10.1034/j.1600-0706.2003.12559.x)

2. Somveille M, Rodrigues AS, Manica A. 2015 Why do birds migrate? A macroecological perspective. *Glob. Ecol. Biogeogr.* **24**, 664–674. (doi:10.1111/geb.12298)

3. Newton I. 2010 *The migration ecology of birds*. New York, NY: Academic Press. See https://books.google.fr/books?hl=f&lr=&id=BnldlbshDWtgC&oi=fnd&pg=PP1&dq=newton-the-migration+ecology+of+birds&ts=JX90_YLH0&sig=TlyxMJmWSpexgWF96QpMF5mbsE.

4. Somveille M, Manica A, Butchart SH, Rodrigues AS. 2013 Mapping global diversity patterns for migratory birds. *PLoS ONE* **8**, e70907. (doi:10.1371/journal.pone.0070907)

5. Robinson RA et al. 2009 Travelling through a warming world: climate change and migratory species. *Endanger. Species Res.* **7**, 87–99. (doi:10.3354/esr00095)

6. van Gils JA et al. 2016 Body shrinkage due to Arctic warming reduces red knot fitness in tropical wintering range. *Science* **352**, 819–821. (doi:10.1126/science.aad6351)

7. IPCC. 2013 Climate change 2013: the physical science basis. In *contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change*, (eds TF Stocker, D Qin, GK Plattner, M Tignor, SK Allen, J Boschung, A Nauels, Y Xia, V Bex, PM Midgley), p. 1535. Cambridge, UK: Cambridge University Press.

8. Meltofte H et al. 2013 Arctic biodiversity assessment. Synthesis. Conservation of Arctic Flora and Fauna (CAFF). See https://oaarchive.arctic-council.org/handle/11374/232.

9. Ward DH, Helmericks J, Hupp JW, McManus L, Budde M, Douglas DC, Tape KD. 2016 Multi-decadal trends in spring arrival of avian migrants to the central Arctic coast of Alaska: effects of

environmental and ecological factors. *J. Avian Biol.* **47**, 197–207. (doi:10.1111/jav.00074)

10. Parmesan C. 2006 Ecological and evolutionary responses to recent climate change. *Annu. Rev. Ecol. Evol. Syst.* **37**, 637–669. (doi:10.1146/annurev.ecolsys.37.091305.110100)

11. Hoye TT, Post E, Meltofte H, Schmidt NM, Forchhammer MC. 2007 Rapid advancement of spring in the High Arctic. *Curr. Biol.* **17**, R449–R451. (doi:10.1016/j.cub.2007.04.047)

12. Louzao M, Afán I, Santos M, Brereton T. 2015 The role of climate and food availability on driving decadal abundance patterns of highly migratory pelagic predators in the Bay of Biscay. *Front. Ecol. Evol.* **3**. (doi:10.3389/fevo.2015.00090)

13. Alves JA, Gunnarsson TG, Hayhow DB, Appleton GF, Potts PM, Sutherland WJ, Gill JA. 2013 Costs, benefits, and fitness consequences of different migratory strategies. *Ecology* **94**, 11–17. (doi:10.1890/12-0737.1)

14. Sandvik H, Erikstad KE, Barrett RT, Yoccoz NG. 2005 The effect of climate on adult survival in five species of North Atlantic seabirds. *J. Anim. Ecol.* **74**, 817–831. (doi:10.1111/j.1365-2656.2005.00981.x)

15. Reiertsen TK et al. 2014 Prey density in non-breeding areas affects adult survival of black-legged kittiwakes *Rissa tridactyla*. *Mar. Ecol. Prog. Ser.* **509**, 289–302. (doi:10.3354/meps10825)

16. Bogdanova MI, Daunt F, Newell M, Phillips RA, Harris MP, Wanless S. 2011 Seasonal interactions in the black-legged kittiwake, *Rissa tridactyla*: links between breeding performance and winter distribution. *Proc. R. Soc. B* **278**, 2412–2418. (doi:10.1098/rspb.2010.2601)

17. O'Connor CM, Norris DR, Crossin GT, Cooke SJ. 2014 Biological carryover effects: linking common concepts and mechanisms in ecology and evolution. *Ecosphere* **5**, 1–11. (doi:10.1890/ES13-00388.1)

18. Thomas DW, Blondel J, Perret P, Lambrechts MM, Speakman JR. 2001 Energetic and fitness costs of

mismatching resource supply and demand in seasonally breeding birds. *Science* **291**, 2598–2600. (doi:10.1126/science.1057487)

19. Patterson TA, Thomas L, Wilcox C, Ovaskainen O, Matthiopoulos J. 2008 State-space models of individual animal movement. *Trends Ecol. Evol.* **23**, 87–94. (doi:10.1016/j.tree.2007.10.009)

20. Wilson RP, Quintana F, Hobson VJ. 2012 Construction of energy landscapes can clarify the movement and distribution of foraging animals. *Proc. R. Soc. B* **279**, 975–980. (doi:10.1098/rspb.2011.1544)

21. Shepard ELC, Wilson RP, Rees WG, Grundy E, Lambertucci VS. 2013 Energy landscapes shape animal movement ecology. *Am. Nat.* **182**, 298–312. (doi:10.1086/671257)

22. Fort J, Porter WP, Grémillet D. 2009 Thermodynamic modelling predicts energetic bottleneck for seabirds wintering in the Northwest Atlantic. *J. Exp. Biol.* **212**, 2483–2490. (doi:10.1242/jeb.032300)

23. Kearney M, Porter W. 2009 Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges. *Ecol. Lett.* **12**, 334–350. (doi:10.1111/j.1461-0248.2008.01277.x)

24. Fort J, Beaugrand G, Grémillet D, Phillips RA. 2012 Biologging, remotely-sensed oceanography and the continuous plankton recorder reveal the environmental determinants of a seabird wintering hotspot. *PLoS ONE* **7**, e41194. (doi:10.1371/journal.pone.0041194)

25. Deville A-S, Labaude S, Robin J-P, Béchet A, Gauthier-Clerc M, Porter W, Fitzpatrick M, Mathewson P, Grémillet D. 2014 Impacts of extreme climatic events on the energetics of long-lived vertebrates: the case of the greater flamingo facing cold spells in the Camargue. *J. Exp. Biol.* **217**, 3700–3707. (doi:10.1242/jeb.106344)

26. Mathewson PD, Moyer-Horner L, Beever EA, Briscoe NJ, Kearney M, Yahn JM, Porter WP. 2017 Mechanistic variables can enhance predictive models of endotherm distributions: the American
pika under current, past, and future climates. *Glob. Change Biol.* **23**, 1048–1064. (doi:10.1111/gcb.13454)

27. Fort J, Porter WP, Grémillet D. 2011 Energetic modelling: a comparison of the different approaches used in seabirds. *Comp. Biochem. Physiol. A. Mol. Integr. Physiol.* **158**, 358–365. (doi:10.1016/j.cbpa.2010.05.004)

28. Fitzpatrick MJ, Mathewson PD, Porter WP. 2015 Validation of a mechanistic model for non-invasive study of ecological energetics in an endangered wading bird with counter-current heat exchange in its legs. *PLoS ONE* **10**, e136677. (doi:10.1371/journal.pone.0136677)

29. Ellis HI, Gabrielsen GW. 2002 Energetics of free-ranging seabirds. In *Biology of marine birds* (eds EA Schreiber, J Burger), pp. 359–407. Boca Raton, FL: CRC Press.

30. Beaugrand G, Luczak C, Edwards M. 2009 Rapid biogeographical plankton shifts in the North Atlantic ocean. *Glob. Change Biol.* **15**, 1790–1803. (doi:10.1111/j.1365-2486.2009.01848.x)

31. Fort J, Cheral Y, Harding AM, Egevang D, Steen H, Kuntz G, Porter WP, Grémillet D. 2010 The feeding ecology of little aiks raises questions about winter zooplankton stocks in North Atlantic surface waters. *Biol. Lett.* **6**, 682–684. (doi:10.1098/rslb.2010.0082)

32. Fuentes M, Porter WP. 2013 Using a microclimate model to evaluate impacts of climate change on sea turtles. *Ecol. Model.* **251**, 150–157. (doi:10.1016/j.ecolmodel.2012.12.020)

33. Mathewson PD, Porter WP. 2013 Simulating polar bear energetics during a seasonal fast using a mechanistic model. *PLoS ONE* **8**, e72863. (doi:10.1371/journal.pone.0072863)

34. Karnovsky N et al. 2010 Foraging distributions of little aiks *Alle alle* across the Greenland Sea: implications of present and future Arctic climate change. *Mar. Ecol. Prog. Ser.* **415**, 283–293. (doi:10.3354/meps08749)

35. Speirs DC, Gurney WS, Heath MR, Horbelt W, Wood SN, De Cuevas BA. 2006 Ocean-scale modelling of the distribution, abundance, and seasonal dynamics of the copepod *Calanus finmarchicus*. *Mar. Ecol. Prog. Ser.* **313**, 173–192. (doi:10.3354/meps313173)

36. Coward AC, de Cuevas BA. 2005 The OCCAM 66 level model: model description, physics, initial conditions and external forcing. Internal Document 99. Southampton, UK: Southampton Oceanography Centre.

37. Manly BFJ, McDonald LL, Thomas DL, McDonald TL, Erickson WP. 2002 *Resource selection by animals: statistical analysis and design for field studies*. Nordrecht, The Netherlands: Kluwer

38. Gillies CS, Hebblewhite M, Nielsen SE, Krawchuk MA, Aldridge CL, Frail JL, Saher DJ, Stevens CE, Jerde CL. 2006 Application of random effects to the study of resource selection by animals. *J. Anim. Ecol.* **75**, 887–898. (doi:10.1111/j.1365-2656.2006.01106.x)

39. Koper N, Manseau M. 2009 Generalized estimating equations and generalized linear mixed-effects models for modelling resource selection. *J. Appl. Ecol.* **46**, 590–599. (doi:10.1111/j.1365-2664. 2009.01642.x)

40. Dormann CF et al. 2013 Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. *Ecography* **36**, 27–46. (doi:10.1111/j.1600-0587.2012.07348.x)

41. Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FK. 2002 Evaluating resource selection functions. *Ecol.*

Model. **157**, 281–300. (doi:10.1016/S0304-3800 (02)00200-4)

42. Han G, Colbourne E, Pepin P, Xie Y. 2015 Statistical projections of ocean climate indices off Newfoundland and Labrador. *Atmos. Ocean* **53**, 556–570. (doi:10.1080/07055900.2015.1047732)

43. Phillips RA, Silk JRD, Croxall JP, Afanaseyv V, Briggs DR. 2004 Accuracy of geolocation estimates for flying seabirds. *Mar. Ecol. Prog. Ser.* **266**, 265–272. (doi:10.3354/meps26625)

44. Masello JF, Kato A, Sommerfeld J, Mattern T, Quillfeldt P. 2017 How animals distribute themselves in space: variable energy landscapes. *Front. Zool.* **14**, 33. (doi:10.1186/s12983-017-0219-8)

45. Reygondeau G, Beaugrand G. 2011 Future climate-driven shifts in distribution of *Calanus finmarchicus*. *Glob. Change Biol.* **17**, 756–766. (doi:10.1111/j.1365-2486.2010.02310.x)

46. Fort J, Cheral Y, Harding A, Welcker J, Jakubas D, Steen H, Karnovsky NJ, Grémillet D. 2010 Geographic and seasonal variability in the isotopic niche of little aiks. *Mar. Ecol. Prog. Ser.* **414**, 293–302. (doi:10.3354/meps08721)

47. Rosing-Asvid A, Hedeholm R, Arendt KE, Fort J, Robertson GJ. 2013 Winter diet of the little auk (*Alle alle*) in the Northwest Atlantic. *Polar Biol.* **36**, 1601–1608. (doi:10.1007/s00300-013-1379-4)

48. Han G, Lu Z, Wang Z, Helbig J, Chen N, de Young B. 2008 Seasonal variability of the Labrador Current and shelf circulation off Newfoundland. *J. Geophys. Res. Oceans* **113**, C10013. (doi:10.1029/2007JC004376)

49. Amélineau F, Grémillet D, Bonnet D, Bot TL, Fort J. 2016 Where to forage in the absence of sea ice? Bathymetry as a key factor for an Arctic seabird. *PLoS ONE* **11**, e0157764. (doi:10.1371/journal.pone.0157764)

50. Frederiksen M et al. 2012 Multicolony tracking reveals the winter distribution of a pelagic seabird on an ocean basin scale. *Divers. Distrib.* **18**, 530–542. (doi:10.1111/j.1472-4642.2011.00864.x)

51. Frederiksen M et al. 2016 Migration and wintering of a declining seabird, the thick-billed murre *Uria lomvia*, on an ocean basin scale: conservation implications. *Biol. Conserv.* **200**, 26–35. (doi:10.1016/j.biocon.2016.05.011)

52. McFarlane Tranquilla LA et al. 2013 Multiple-colony winter habitat use by murres *Uria* spp. in the Northwest Atlantic Ocean: implications for marine risk assessment. *Mar. Ecol. Prog. Ser.* **472**, 287–303. (doi:10.3354/meps10053)

53. Linneberg JF, Fort J, Guilford T, Reuleaux A, Mosbech A, Frederiksen M. 2013 Sympatric breeding aiks shift between dietary and spatial resource partitioning across the annual cycle. *PLoS ONE* **8**, e72987. (doi:10.1371/journal.pone.0072987)

54. Boertmann D, Lyngs P, Merkel FR, Mosbech A. 2004 The significance of Southwest Greenland as winter quarters for seabirds. *Bird Conserv. Int.* **14**, 87–112. (doi:10.1017/S0959270904000127)

55. Egevang C, Stenhouse LJ, Phillips RA, Petersen A, Fox JW, Silk JRD. 2010 Tracking of Arctic terns *Sterna paradisaea* reveals longest animal migration. *Proc. Natl Acad. Sci. USA* **107**, 2078–2081. (doi:10.1073/pnas.0904993107)

56. Gil O et al. 2013 Trans-equatorial migration routes, staging sites and wintering areas of a high-Arctic avian predator: the long-tailed skua (*Stercorarius longicaudus*). *PLoS ONE* **8**, e64614. (doi:10.1371/journal.pone.0064614)

57. Hedd A, Montevecchi W, Otley H, Phillips R, Fifield D. 2012 Trans-equatorial migration and habitat use by sooty shearwaters *Puffinus griseus* from the South Atlantic during the nonbreeding season. *Mar. Ecol. Prog. Ser.* **449**, 277–290. (doi:10.3354/meps09538)

58. Daufresne M, Lengfellner K, Sommer U. 2009 Global warming benefits the small in aquatic ecosystems. *Proc. Natl Acad. Sci. USA* **106**, 12 788–12 793. (doi:10.1073/pnas.0902080106)

59. Garzke J, Ismar SMH, Sommer U. 2014 Climate change affects low trophic level marine consumers: warming decreases copepod size and abundance. *Oecologia* **177**, 849–860. (doi:10.1007/s00442-014-3130-4)

60. Pitois SG, Fox CJ. 2006 Long-term changes in zooplankton biomass concentration and mean size over the Northwest European shelf inferred from continuous plankton recorder data. *ICES J. Mar. Sci. J. Cons.* **63**, 785–798. (doi:10.1016/j.icesjms. 2006.03.009)

61. Reid PC, Colebrook JM, Matthews JBL, Aiken J, Team CPR. 2003 The continuous plankton recorder: concepts and history, from plankton indicator to undulating recorders. *Prog. Oceanogr.* **58**, 117–173. (doi:10.1016/j.pcean.2003.08.002)

62. Chust G, Castellani C, Licandro P, Ibaibarriaga L, Sagarminaga Y, Irigoien X. 2014 Are *Calanus* spp. shifting poleward in the North Atlantic? A habitat modelling approach. *ICES J. Mar. Sci. J. Cons.* **71**, 241–253. (doi:10.1093/icesjms/fst147)

63. Wilson RJ, Speirs DC, Heath MR. 2015 On the surprising lack of differences between two congeneric calanoid copepod species, *Calanus finmarchicus* and *C. helgolandicus*. *Prog. Oceanogr.* **134**, 413–431. (doi:10.1016/j.pcean.2014.12.008)

64. Sugin G, Swingedouw D, Drijfhouf S, Mary Y, Bennabi A. 2017 Abrupt cooling over the North Atlantic in modern climate models. *Nat. Commun.* **8**, ncomms14375. (doi:10.1038/ncomms14375)

65. Falk K, Durinck J. 1993 The winter diet of thick-billed murres, *Uria lomvia*, in western Greenland, 1988–1989. *Can. J. Zool.* **71**, 264–272. (doi:10.1139/z93-038)

66. Elliot RD, Ryan PC, Lidster WW. 1990 The winter diet of thick-billed murres in coastal Newfoundland waters. *Stud. Avian. Biol.* **14**, 125–138.

67. Pedersen CE, Falk K. 2001 Chick diet of dovekies *Alle alle* in Northwest Greenland. *Polar Biol.* **24**, 53–58. (doi:10.1007/s003000000173)

68. Lawson JW, Magalhães AM, Miller EH. 1998 Important prey species of marine vertebrate predators in the Northwest Atlantic: proximate composition and energy density. *Mar. Ecol. Prog. Ser.* **164**, 13–20. (doi:10.3354/meps164013)

69. Phillips RA, Silk JRD, Croxall JP, Afanaseyv V, Bennett VJ. 2005 Summer distribution and migration of nonbreeding albatrosses: individual consistencies and implications for conservation. *Ecology* **86**, 2386–2396. (doi:10.1890/04-1885)

70. Dias MP, Granadeiro JP, Phillips RA, Alonso H, Catry P. 2011 Breaking the routine: individual Cory's shearwaters shift winter destinations between hemispheres and across ocean basins. *Proc. R. Soc. B* **278**, 1786–1793. (doi:10.1098/rspb.2010. 2114)

71. Guilford T, Freeman R, Boyle D, Dean B, Kirk H, Phillips R, Perrins C. 2011 A dispersive migration in
the Atlantic puffin and its implications for
migratory navigation. *PLoS ONE* 6, e21336.
(doi:10.1371/journal.pone.0021336)

72. Tranquilla LAM, Montevecchi WA, Fifield DA, Hedd A, Gaston AJ, Robertson GJ, Phillips RA. 2014 Individual winter movement strategies in two species of murre (*Uria* spp.) in the Northwest Atlantic. *PLoS ONE* 9, e90583. (doi:10.1371/journal.pone.0090583)

73. Visser ME, Perdeck AC, Van Balen JH, Both C. 2009 Climate change leads to decreasing bird migration distances. *Glob. Change Biol.* **15**, 1859–1865. (doi:10.1111/j.1365-2486.2009.01865.x)

74. Heath JA, Steenhof K, Foster MA. 2012 Shorter migration distances associated with higher winter temperatures suggest a mechanism for advancing nesting phenology of American kestrels *Falco sparverius*. *J. Avian Biol.* **43**, 376–384. (doi:10.1111/j.1600-048X.2012.05595.x)

75. Pulido F, Berthold P. 2010 Current selection for lower migratory activity will drive the evolution of

residency in a migratory bird population. *Proc. Natl Acad. Sci. USA* **107**, 7341–7346. (doi:10.1073/pnas.0910361107)

76. Berge J *et al*. 2015 In the dark: a review of ecosystem processes during the Arctic polar night. *Prog. Oceanogr.* **139**, 258–271. (doi:10.1016/j.pocean.2015.08.005)

77. Grémillet D, Kuntz G, Gilbert C, Woakes AJ, Butler PJ, le Maho Y. 2005 Cormorants dive through the Polar night. *Biol. Lett.* **1**, 469–471. (doi:10.1098/rsbl.2005.0356)

78. Urban MC, Zarnetske PL, Skelly DK. 2013 Moving forward: dispersal and species interactions determine biotic responses to climate change. *Ann. N. Y. Acad. Sci.* **1297**, 44–60. (doi:10.1111/nyas.12184)

79. Alexander JM, Diez JM, Levine JM. 2015 Novel competitors shape species' responses to climate change. *Nature* **525**, 515. (doi:10.1038/nature14952)

80. Fort J *et al*. 2013 Multicolony tracking reveals potential threats to little auks wintering in the North Atlantic from marine pollution and shrinking

sea ice cover. *Divers. Distrib.* **19**, 1322–1332.
(doi:10.1111/ddi.12105)

81. Siegel DA, Doney SC, Yoder JA. 2002 The North Atlantic spring phytoplankton bloom and Sverdrup's critical depth hypothesis. *Science* **296**, 730–733. (doi:10.1126/science.1069174)

82. Long RA, Bowyer RT, Porter WP, Mathewson P, Monteith KL, Findholt SL, Dick BL, Kie JG. 2016 Linking habitat selection to fitness-related traits in herbivores: the role of the energy landscape. *Oecologia* **181**, 709–720.

83. Dudley PN, Bonazza R, Porter WP. 2016 Climate change impacts on nesting and internesting leatherback sea turtles using 3D animated computational fluid dynamics and finite volume heat transfer. *Ecol. Model.* **320**, 231–240. (doi:10.1016/j.ecolmodel.2015.10.012)

84. Enstipp MR, Grémillet D, Jones DR. 2007 Investigating the functional link between prey abundance and seabird predatory performance. *Mar. Ecol. Prog. Ser.* **331**, 267–279. (doi:10.3354/meps331267)