# Reproductive constraints influence habitat accessibility, segregation, and preference of sympatric albatross species

Michelle Kappes, Scott Shaffer, Yann Tremblay, David Foley, Daniel Palacios, Steven Bograd, Daniel Costa

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

To cite this version:

Michelle Kappes, Scott Shaffer, Yann Tremblay, David Foley, Daniel Palacios, et al.. Reproductive constraints influence habitat accessibility, segregation, and preference of sympatric albatross species. Movement Ecology, 2015, 3, pp.34. 10.1186/s40462-015-0063-4. hal-01922016

HAL Id: hal-01922016

https://hal.umontpellier.fr/hal-01922016

Submitted on 7 Jun 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

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

# Reproductive constraints influence habitat accessibility, segregation, and preference of sympatric albatross species

Michelle A. Kappes<fnref n="1" />,<fnref n="4" />*, Scott A. Shaffer<fnref n="1" />,⁵, Yann Tremblay<fnref n="1" />,⁶, David G. Foley²,³, Daniel M. Palacios²,³, Steven J. Bograd² and Daniel P. Costa<fnref n="1" />

## Abstract

**Background:** The spatiotemporal distribution of animals is dependent on a suite of factors, including the distribution of resources, interactions within and between species, physiological limitations, and requirements for reproduction, dispersal, or migration. During breeding, reproductive constraints play a major role in the distribution and behavior of central place foragers, such as pelagic seabirds. We examined the foraging behavior and marine habitat selection of Laysan (*Phoebastria immutabilis*) and black-footed (*P. nigripes*) albatrosses throughout their eight month breeding cycle at Tern Island, Northwest Hawaiian Islands to evaluate how variable constraints of breeding influenced habitat availability and foraging decisions. We used satellite tracking and light-based geolocation to determine foraging locations of individuals, and applied a biologically realistic null usage model to generate control locations and model habitat preference under a case-control design. Remotely sensed oceanographic data were used to characterize albatross habitats in the North Pacific.

**Results:** Individuals of both species ranged significantly farther and for longer durations during incubation and chick-rearing compared to the brooding period. Interspecific segregation of core foraging areas was observed during incubation and chick-rearing, but not during brooding. At-sea activity patterns were most similar between species during brooding; neither species altered foraging effort to compensate for presumed low prey availability and high energy demands during this stage. Habitat selection during long-ranging movements was most strongly associated with sea surface temperature for both species, with a preference for cooler ocean temperatures compared to overall availability. During brooding, lower explanatory power of habitat models was likely related to the narrow range of ocean temperatures available for selection.

**Conclusions:** Laysan and black-footed albatrosses differ from other albatross species in that they breed in an oligotrophic marine environment. During incubation and chick-rearing, they travel to cooler, more productive waters, but are restricted to the low-productivity environment near the colony during brooding, when energy requirements are greatest. Compared to other albatross species, Laysan and black-footed albatrosses spend a greater proportion of time in flight when foraging, especially during the brooding period; this strategy may be adaptive for locating dispersed prey in an oligotrophic environment.

**Keywords:** Laysan albatross, Black-footed albatross, Foraging behavior, Activity patterns, Spatial segregation, Habitat selection, Satellite tracking, Geolocation, Spatial modeling

<footnote>* Correspondence: michelle.kappes@oregonstate.edu</footnote>

<footnote n="1">¹Department of Ecology and Evolutionary Biology, University of California
Santa Cruz, 100 Shaffer Road, Santa Cruz, California 95060, USA</footnote>

<footnote n="4">⁴Present address: Department of Fisheries and Wildlife, Oregon State
University, 104 Nash Hall, Corvallis, Oregon 97331, USA</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

To maximize fitness, animals can optimize energy acquisition through the selection of favorable habitats [1–4]. Foraging habitat use is constrained not only by the distribution of resources, but also by the physiological capabilities of the animal, memory and learned behaviors, intra- and interspecific interactions, and requirements for activities other than foraging, such as reproduction or migration [5–11].

Central place foragers, animals constrained by the need to return to a particular place (e.g., breeding colony) after a foraging trip, have a limited window of time, and therefore space, in which to search for and obtain food [12, 13]. Marine species that forage at sea while conducting breeding activities on land exemplify this behavior. The proximity of suitable foraging habitat to breeding sites has profound implications for behavior, habitat use, and energetics of central place foragers [14–20]. The extent to which marine animals are constrained in this context can depend on the particular stage of reproduction, especially when coupled with temporally changing energetic demands of offspring [15, 18, 21–24].

Albatrosses are pelagic seabirds that exhibit three distinct reproductive stages: (1) the incubation period (parents alternate between fasting at the nest to incubate the egg, and foraging at sea); (2) the brooding period (breeding pairs alternate between fasting at the nest to brood and provision the chick, and foraging at sea); and (3) the chick-rearing stage (breeding pairs forage independently at sea, returning to the nest periodically to quickly provision the chick). Differing energetic demands during these reproductive stages [18, 21, 25] constrain the duration and range of foraging movements, and ultimately influence the accessibility of foraging habitats to breeding adults [20, 26–29]. Habitat accessibility is additionally influenced by spatiotemporal variation in the distribution of preferred marine habitats e.g., [30–32].

Laysan (*Phoebastria immutabilis*) and black-footed (*P. nigripes*) albatrosses are long-ranging pelagic seabirds of the North Pacific that breed primarily in the warm, oligotrophic marine environment of the Northwest Hawaiian Islands [33, 34]. This is in contrast to other albatross species that nest in closer proximity to more productive marine habitats, such as subtropical and polar convergences, continental shelf-breaks and slopes, and coastal upwelling zones [35, 36]. Albatrosses use energetically efficient gliding flight to make long distance movements, allowing Laysan and black-footed albatrosses to travel to more productive habitats of the North Pacific during the incubation and chick-rearing stages [28, 29, 32]. During brooding, however, they are restricted to the oligotrophic waters close to their breeding colonies, where prey abundance is likely

lower [37, 38]. Even the equatorially nesting waved albatross (*Phoebastria irrorata*) has access to higher productivity waters during brooding due to the close proximity of the Humboldt Current, equatorial fronts, and localized upwelling [39].

According to a model developed by Ricklefs [21], pelagic seabirds expend more energy per day during brooding than incubation or chick-rearing. This is due to the fact that brooding adults fast while providing the daily needs of the rapidly-growing chick, whereas chick-rearing adults provide half of the chick's daily energy requirements and do not need to fast [40]. In the wandering albatross (*Diomedea exulans*), individuals respond by expending more energy while foraging during brooding compared to incubation [18], which allows them to maximize the rate of energy delivery to chicks during this time-limited breeding stage [18, 41]. It has been suggested that some albatrosses may couple the high energetic demands of the brooding period with seasonal increases in prey abundance [15], however, it is unknown how Hawaiian albatrosses locate sufficient prey resources to meet these high energetic demands when foraging in a low-productivity environment. Other top predators that make use of oligotrophic environments may exhibit lower abundances in the least productive waters or use low-productivity waters primarily for breeding, rather than feeding, activities [42–44].

Both Laysan and black-footed albatrosses are known to extensively use an oceanic region known as the North Pacific Transition Zone (NPTZ) when making long-ranging movements [28, 29, 32]. The NPTZ is bounded by the eastward flowing currents of the subtropical and subarctic gyres [45], and exhibits a sharp meridional transition in surface phytoplankton chlorophyll-a concentration, a feature known as the Transition Zone Chlorophyll Front [TZCF; [46]. The dynamics of the front act to aggregate phytoplankton and particulate matter, attracting mobile organisms; because zooplankton and other actively swimming or buoyant organisms can maintain their position in the front, the resulting prey aggregation serves to attract higher-trophic-level predators [46–49]. The latitudinal positional of the TZCF fluctuates on seasonal, interannual, and decadal time scales [50]; this spatiotemporal variation alters the distribution of marine habitats used by Hawaiian albatrosses and their proximity to breeding colonies. Expansion of oligotrophic habitat in the subtropical gyre of the North Pacific [51], as well as changes to the positioning of the North Pacific Current [52] due to climate change, could have negative effects on Laysan and black-footed albatross populations if preferred habitats become more distant from the colony during critical portions of the breeding season.

Here we examine the foraging behavior and habitat preference of sympatrically nesting Hawaiian albatrosses
throughout the breeding cycle in order to evaluate how reproductive constraints impact foraging distribution, habitat segregation, at-sea activity patterns, and habitat preference in the context of their oligotrophic breeding environment. Based on previous studies, we expected that foraging movements would be most constrained during the brooding period for both species and least constrained during the incubation and chick-rearing periods [28, 29, 32]. We also expected to see the greatest overlap in foraging distributions during brooding when adults take shorter trips to frequently provision young chicks [29]. We therefore hypothesized that activity patterns of Hawaiian albatrosses during brooding would differ from the incubation and chick-rearing periods in response to (1) the need to maximize energy delivery to rapidly-growing chicks, (2) the use of a low-productivity foraging environment, and (3) the greater potential for inter- and intraspecific competition due to contraction of foraging ranges. We also hypothesized that greater interspecific differences in activity patterns would be observed during brooding; divergent foraging strategies could reduce competitive interactions between species despite high overlap in foraging distributions. Finally, we hypothesized that each species would display consistent environmental associations in selecting habitats throughout the breeding season, but that reproductive constraints and seasonal variation in marine habitats would influence the composition of utilized habitats during the different breeding stages.

# Methods

## Study area and tracking activities

Breeding Laysan and black-footed albatrosses were studied at Tern Island (23.87° N, 166.28° W), French Frigate Shoals, Northwest Hawaiian Islands, during the incubation, brooding, and chick-rearing periods. The

incubation period for both species lasts approximately 60 days, beginning in mid-late November with the laying of a single egg and ending with the onset of hatching in late January. The brooding period typically lasts 2–3 weeks, until chicks are left alone at the nest in mid-February. The chick-rearing period lasts approximately 4–5 months; after a period of fasting, chicks fledge independently in June and July [33, 34].

We used a combination of satellite tracking and light-level based geolocation to determine at-sea locations of adult Hawaiian albatrosses throughout the reproductive period. Satellite tags were used for short-term deployments during incubation and brooding, whereas geolocation tags were used to obtain foraging positions during chick-rearing, and to supplement data from the incubation period. While geolocation tag deployments also spanned the brooding period, the single position derived per day by light-level based geolocation did not have sufficient resolution to accurately capture movements during brooding when average trip durations were 2–3 days (Table 1; satellite tags are accurate to <10 km [53]; geolocation tags are accurate to ~200 km [54]). Tracking activities were conducted during five consecutive breeding seasons, from 2002–03 through 2006–07; satellite tracks were obtained during the 2002–03 to 2005–06 seasons and geolocation data were obtained during the 2003–04 to 2006–07 seasons. Sex of tracked individuals was determined from blood samples [55]; for six individuals for which we did not obtain blood samples, sex was determined by comparison of culmen lengths [28].

## Satellite tracking

One hundred and forty-seven adult albatrosses (76 Laysan and 71 black-footed) were equipped with satellite platform terminal transmitters (30 g Pico-100,

<table><caption>**Table 1** Summary characteristics (Mean ± SD) of Laysan and black-footed albatross foraging trips. To reduce the influence of individuals tracked for multiple foraging trips, a single trip per individual was randomly-selected for each reproductive stage to include in the calculation of mean values</caption><thead><tr><th>Species</th><th></th><th>Incubation</th><th>Brooding</th><th>Chick-Rearing</th></tr></thead><tbody><tr><td rowspan="5">Laysan albatross</td><td>Number of individuals tracked</td><td>58</td><td>38</td><td>26</td></tr><tr><td>Total number of foraging trips</td><td>72</td><td>40</td><td>153</td></tr><tr><td>Trip duration (days)</td><td>17.6 ± 7.45<sup>a,*</sup></td><td>2.62 ± 0.76<sup>b</sup></td><td>14.5 ± 3.91<sup>c</sup></td></tr><tr><td>Maximum distance from colony (km)</td><td>2433 ± 837<sup>a</sup></td><td>420 ± 266<sup>b</sup></td><td>2489 ± 676<sup>a,*</sup></td></tr><tr><td>Azimuth to most distant point from colony (°)</td><td>350 ± 26.5<sup>a,*</sup></td><td>35.4 ± 64.5<sup>b</sup></td><td>358 ± 32.4<sup>a,*</sup></td></tr><tr><td rowspan="5">Black-footed albatross</td><td>Number of individuals tracked</td><td>51</td><td>35</td><td>23</td></tr><tr><td>Total number of foraging trips</td><td>60</td><td>37</td><td>83</td></tr><tr><td>Trip duration (days)</td><td>13.3 ± 5.07<sup>a,*</sup></td><td>2.67 ± 0.87<sup>b</sup></td><td>14.4 ± 5.14<sup>a</sup></td></tr><tr><td>Maximum distance from colony (km)</td><td>2045 ± 1035<sup>a</sup></td><td>313 ± 144<sup>b</sup></td><td>2883 ± 998<sup>c,*</sup></td></tr><tr><td>Azimuth to most distant point from colony (°)</td><td>22.8 ± 32.5<sup>*</sup></td><td>47.4 ± 65.7</td><td>31.9 ± 35<sup>*</sup></td></tr></tbody></table>

Different lowercase letters indicate significant differences (P < 0.05) among reproductive stages; asterisks indicate significant differences between species
Microwave Telemetry, Columbia, MD; or 35 g SPOT4, Wildlife Computers, Redmond, WA); tags were attached to dorsal feathers with adhesive tape (tesa®, Hamburg, Germany), and transmissions were downloaded via the Argos satellite system (Service Argos, Inc., Largo, MD). Satellite tags were programmed to transmit continuously every 90 s, with the exception of 19 foraging trips, when tags were programmed to use a 6:18 h on:off duty cycle (11 trips), 20:4 h on:off duty cycle (6 trips), or 9:15 h on:off duty cycle (2 trips) to conserve battery life on anticipated longer trips. In three instances, the battery on the satellite transmitter failed or the satellite transmitter fell off the bird before completion of the foraging trip. During 2002–03, 2004–05, and 2005–06, satellite-tracked individuals were also equipped with archival tags (10 g Lotek LTD 2400, Lotek Wireless, St. John’s, Newfoundland) attached to a plastic leg band so that temperature recordings (±0.05 °C) every 12–40 s could be used to characterize foraging activity while at sea [56]. The combined mass of devices deployed on individuals was less than 2 % of total bird body mass in all cases, below the recommended limit for studies involving albatrosses [57].

Before calculating trip characteristics, satellite loca-
tions were first delimited by observations of departure
and arrival times at the breeding colony or by visual
inspection of the tracks. To remove unlikely locations
from the data set, we applied an iterative forward/
backward averaging speed filter [58] implemented in
Matlab (The MathWorks, Natick, MA) with a max-
imum speed limit of 80 km h⁻¹ (following [29, 59])
to remove unrealistic flight speeds [60]. Satellite
tracks were then interpolated to every 10 min using a
Bézier curve with µ = 0.3 [µ controls the elasticity of
the curve; following [61], and subsampled to two lo-
cations per day to match the temporal scale of geolo-
cation positions (see below).

Tracking with geolocation tags

Archival geolocation tags (10 g Lotek LTD 2400, Lotek Wireless, St. John's, Newfoundland) were deployed and successfully recovered from 34 Laysan and 26 black-footed albatrosses. Each tag recorded ambient light intensity and temperature every 480 s or 540 s to determine a single daily location: longitude based on the establishment of local noon in comparison to Universal Time, and latitude based on day length for the established longitude [62, 63]. Temperature sensors on the geolocation tags allowed a refinement of the location data based on sea surface temperature [SST; [54, 64], as well as providing a record of foraging activity while at sea.

Geolocation tags were deployed for up to one year,
only a portion of which overlapped our on-colony re-
search activities, therefore we had limited information

on the presence or absence of geolocation-tagged indi-
viduals at the colony. Due to the nature of the SST-
processing algorithm [54, 64], it was necessary to first
delimit individual foraging trips so that on-colony loca-
tions were removed from the dataset. We calculated
minimum daily temperatures recorded by geolocation
tags to determine days when albatrosses were likely at
sea; evaporative cooling due to immersion in water
would lead to cooler minimum temperatures on days
when albatrosses were off-colony. In order to detect
shifts in minimum daily temperatures indicative of on-
or off-colony periods, we used an algorithm developed
for detecting climate regime shifts [65], which imple-
ments a sequential version of the partial CUSUM
method combined with the t-test, and is available for
download at: http://www.beringclimate.noaa.gov/regimes/.
For those individuals for which exact colony attendance
patterns were known (25 individuals during incubation),
this detection algorithm correctly classified 99.0 % of days
when birds were known to be on-colony and 83.6 % of
days when birds were known to be at-sea (n = 548 bird
days). This method was therefore effective in delimiting
foraging departures, but some at-sea locations may have
been excluded. Estimated departure and arrival times for
geolocation-derived foraging trips were calculated based
on Argos-derived transit rates during the incubation
period [32], and distance to the colony from the first and
last off-colony location.

To remove unlikely light-based positions, we applied
the same speed filter as above, with a more conserva-
tive speed limit of 50 km h⁻¹ to account for greater
error in geolocation position estimates [54], and inter-
polated to two positions per day using a Hermite
spline [following [61].

Foraging movements

Descriptive characteristics of all Argos- and geolocation-
derived tracks were calculated to compare foraging
behavior during each reproductive stage. Maximum dis-
tance traveled from the colony was calculated using
great-circle distances to account for the earth's curva-
ture. Destination bearings, as defined by the azimuth to
the most distant point, were calculated to describe over-
all direction of foraging trips. Trip duration was defined
as either the time elapsed between the observed depart-
ure and arrival of the bird, or the time elapsed between
the estimated departure and arrival times.

We used linear mixed-effects models [66], with indi-
vidual as a random effect, to compare maximum ranges
and trip durations between species and reproductive
stages. Mixed-effects models were followed by contrast
analysis with the multcomp package [67] in the program
R [68]. To compare destination bearings between species
and breeding stages, we used Watson-Williams tests for
circular data in the CircStat toolbox in Matlab [69] after ensuring data followed the von Mises distribution (equivalent to normal distribution for circular data). Accounting for random effects has not been developed in the circular statistics modeling framework, therefore we randomly selected a single trip for each individual during each reproductive stage for this analysis.

## Foraging distribution

To determine patterns of interspecific habitat segregation during each reproductive stage, we used kernel estimation [70] to determine utilization distributions (UD) for each species. Because individuals contributed varying numbers of foraging trips to the overall data set, a single trip was randomly selected for each individual during each reproductive stage so that the influence of individual behaviors on estimation of kernel densities was reduced. Geographic coordinates of interpolated albatross locations were transformed using a Lambert Cylindrical Equal Area projection [71], and UD s were computed on a 50-km grid using the R package adehabitat [72]. To allow comparisons between species, the smoothing factor (h) was set to the mean of the h values calculated from each species, as determined using least-squares cross-validation [73]. We then employed a randomization analysis to test the null hypotheses that there was no spatial segregation in foraging distributions between species during each reproductive stage [74]. For each comparison, species was randomly assigned to tracks using the same species ratio observed, and kernel analysis applied. The area of overlap divided by the area of the larger of the two UD polygons was used as the test statistic following Breed et al. [74], for the 25 % (core area), 50 % (focal region), and 95 % (foraging range) UD s [29]. Each test was iterated 500 times, and the P-value was calculated as the proportion of random overlaps smaller than the observed overlap [74].

## Activity patterns at sea

We determined the proportion of time albatrosses spent in flight and the frequency of landings on the sea surface during daylight and nighttime hours to characterize foraging activity patterns during each reproductive stage using geolocation tags equipped with temperature sensors. Albatrosses are surface-feeders and, by necessity, must land on the sea surface in order to consume prey. Therefore, landing rates are indicative of the level of feeding effort. Flight costs measured for albatrosses have demonstrated that the most energetically demanding activities albatrosses engage in at sea are take-offs and landings [75]; landing rates were also highly correlated with field metabolic rates in a study of wandering albatrosses [76].

For albatrosses equipped with temperature recorders, we implemented an algorithm (Iknos toolkit for Matlab; Y. Tremblay, unpublished) designed to identify landings based on rapid changes in temperature, and stable periods associated with sitting on the sea surface [56]. We used civil twilight (sun no more than 6° below the horizon) to define daylight hours, based on temporally-matched tracking locations and NOAA's solar calculator, as implemented in the maptools package in R [77].

To test whether percent time in flight and landing rates differed between high-resolution temperature records (12–40 s) and low-resolution records (480 s and 540 s), we subsampled high resolution records to the lower resolution and compared these metrics. While percent time in flight did not differ between high resolution and subsampled records (paired t-test: $t_{140} = 1.11$, P = 0.27), landing rates were significantly lower in the subsampled group (paired t-test: $t_{140} = 20.9$, P < 0.0001). To compare all temperature records in subsequent analyses (only low resolution records were available during the chick-rearing period), we used histogram matching to rescale landing rates from low resolution records [78]. Histogram matching is an image processing technique used to rescale lower resolution data so that the histogram of the data after transformation matches that of the reference data [79]. This allowed us to present landing rates comparable to values in the literature determined from high resolution temperature loggers; significance tests based on histogram-matched landing rates yielded the same conclusions as tests based on landing rates from all the data subsampled to the lowest resolution available.

Percent time in flight was arcsine transformed and landing rate was log transformed prior to analysis to meet normality assumptions. We used linear mixed-effects models [66] with individual as a random effect to compare percent time in flight and landing rates between species and reproductive stages, followed by contrast analysis with the multcomp package in R [67].

## Habitat preference

We followed the analytical framework of Aarts et al. [80] to model habitat preference of Laysan and black-footed albatrosses during each stage of breeding. We adopted a case-control design such that each tracking location was temporally matched with three randomly generated control locations. Because breeding albatrosses are central place foragers, it is unrealistic to assume that all points within the study area are equally accessible [80, 81]. Therefore, we adopted a simple null model of usage that assumes that the accessibility $\alpha$ [81] of a point in space is inversely related to the distance from the colony ($d_c$) [82]. We assumed that locations were not accessible beyond the maximum range observed for each breeding
stage for each species; control locations were then quasi-randomly selected at a rate proportional to $\alpha$ within this range for each species-stage combination. Because this null model may over- or under-predict true accessibility, we also included $d_c$ as a candidate covariate in our habitat preference models [83] as suggested by Aarts et al. [80].

Laysan and black-footed albatrosses breeding at Tern Island rarely made southward departures from the colony (Fig. 1); therefore we also modeled habitat preference using a more restrictive null model of usage to reflect the northern bias of tracking locations (see Appendix A in Additional file 1) and to ensure conclusions were not sensitive to the choice of null usage model. Final habitat models and response curves of selected covariates were generally similar irrespective of the null usage model (see Results and Appendix A in Additional file 1), therefore we present only the results based on the simpler null usage

model, which assumed that locations were accessible within the maximum range observed for each species-stage, regardless of the direction from the colony. This required fewer assumptions regarding accessibility of habitats, and followed the implementation of Wakefield et al. [82].

Habitat preference models were implemented using a binomial distribution and logit link (the inverse of the logistic function) to relate the response (tracking locations, value of 1 (Fig. 1) and control locations, value of 0 (Fig. 2)) to environmental covariates. Geolocation tracks (incubation and chick-rearing) were only included in habitat analyses if they included at least five filtered off-colony locations, to ensure adequate coverage of sampled habitats along each track. We applied generalized additive mixed models (GAMMs) to allow for the possibility of a nonlinear response to environmental covariates [84], and to account for non-independence of

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /></figure>

points within trips and the variable number of trips contributed by each individual [85].

## Environmental covariates

We selected marine habitat variables based on their potential to characterize physical features that may stimulate or aggregate albatross prey resources. Because Hawaiian albatrosses are surface-feeders foraging within the top meter of the water column, they rely on neustonic or vertically migrating prey, as well as on carrion and fisheries discards [33, 34, 86, 87]. Hawaiian albatrosses consume a diverse array of prey items, with Ommastrephid squid and flying fish (Exocoetidae) eggs comprising the largest proportion of the Laysan and black-footed albatross diet, respectively [86]. The diversity of prey types consumed by Hawaiian albatrosses reflects the variety of marine habitats encountered during their long-ranging movements; for instance, flying fish eggs are likely unavailable to black-footed albatrosses when foraging in the temperate waters of the California Current, where squid may be a larger component of their diet.

Previous studies have demonstrated a relationship between catches of Ommastrephid squid and SST [88–90], productivity [91], and position relative to frontal features of the North Pacific Transition Zone (NPTZ) [90]. Areas of surface convergence, such as those found near fronts and eddies, can also aggregate algae or floating material which many species of flying fish use to attach their non-buoyant eggs with sticky filaments [92, 93]. Adult flying fish, also prey of Hawaiian albatrosses [86], can be found at the outer edges of rapidly rotating eddies [93, 94].

The habitat relationships of these prey items informed our selection of environmental parameters with which to

investigate albatross habitat use. We used SST and primary productivity (PP) to characterize the regional thermal and phytoplankton production regimes; latitudinal gradients in SST and proximity to the Transition Zone Chlorophyll Front (TZCF) to describe large-scale frontal characteristics; wind stress curl to describe wind-driven oceanic upwelling or downwelling; and sea surface height anomaly (SSHa) and eddy kinetic energy (EKE) to investigate overall intensity of eddy activity [95–97]. As albatrosses use near-surface winds to engage in gliding flight [98], we also characterized wind speed and direction for all locations used in the habitat analysis. Because Hawaiian albatrosses periodically visit the continental margins of the western coast of North America and the Aleutian Islands during breeding [28, 29, 32], and bathymetric features can stimulate local production and aggregate prey [99, 100], we also used sea floor depth to characterize marine habitats; bathymetry data was extracted from NOAA's ETOPO2 data set (http://www.ngdc.noaa.gov/mgg/global/etopo2.html).

Environmental data were obtained by querying the NOAA OceanWatch Live Access Server using Matlab and ERDDAP (http://coastwatch.pfel.noaa.gov/erddap/). Where possible, we used satellite products with the finest temporal resolution available (daily or weekly); we used monthly composites when necessary to compensate for data gaps due to cloud cover.

We used a blended product of SST derived from both microwave and infrared sensors carried on multiple platforms at a spatial resolution of 0.1 degrees and as a 5-day temporal composite [101]. As SSTs in the central North Pacific are distributed in broad latitudinal bands, latitudinal gradients in SST (dySST) were used to describe frontal structure [102], by computing the local derivative of adjacent 0.1-degree pixels in the north–south direction.
Monthly PP estimates [103] were derived at a spatial resolution of 0.1 degrees from monthly chlorophyll-a values and photosynthetically available radiation obtained from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) on the Orbview-2 satellite, and SST values from the AVHRR Pathfinder Oceans Project.

Monthly chlorophyll-a values obtained from SeaWiFS were also used to calculate distance to the TZCF. Following Polovina et al. [46] and Bograd et al. [50], the TZCF was defined as the 0.2 mg m⁻³ chlorophyll-a contour.

Daily composites of wind velocity fields were gathered by the SeaWinds scatterometer aboard NASA's QuikSCAT satellite at a resolution of 0.25 degrees and a reference height of 10 m, and used to assign daily values of wind speed and direction to each location.

The NOAA CoastWatch West Coast Node also processes wind vector fields to calculate wind stress curl from these data. We used wind stress curl as a metric of surface convergence/divergence [50], as a proxy for local aggregation of neustonic or buoyant prey (e.g., fish eggs, fish and squid larvae, zooplankton, fishing discards). For each location of interest, we calculated a mean wind stress curl from the previous 30 days to get a representation of past physical forcing in the area of interest.

EKE was calculated from surface ocean currents derived from four satellite altimeters (Jason-1, ENVISAT, ERS-1 and 2, and TOPEX/Poseidon) provided by the AVISO program at a resolution of 0.25 degrees. For each location of interest, a 10 x 10-degree spatial mean was removed from each 7-day composite of zonal and meridional surface ocean currents to generate zonal and meridional velocity anomalies (u' and v'). The EKE (per unit mass) was then calculated as

$$EKE = \frac{1}{2} (u'^2 + v'^2) \qquad (1)$$

SSHa was derived from AVISO sea level data compared to the mean geoid as measured from 1993–1995.

Oceanographic data were extracted so that composite data matched closest in time to tag-derived locations (and matching control locations). The median value of each oceanographic variable was calculated for grid cells falling within the approximate error of Argos or geolocation positions. To estimate error of Argos locations, we conducted trials on seven satellite tags affixed to permanent structures at the University of California Santa Cruz, Long Marine Lab, for a period of three weeks. Argos locations were then compared to GPS-derived locations; mean error for all location qualities was 0.06° (6.7 km) from a total of 2,585 location fixes. To be conservative, satellite tag-derived locations (and matching control locations) were extracted within a 0.15° (16.7 km) longitude by 0.15° (16.7 km) latitude grid centered on each at-sea location. For geolocation-derived positions, data were extracted within a 1° (111 km) longitude by 2° (222 km) latitude grid centered on each location [the approximate error of the geolocation method; [54].

## Model fitting and selection

We modeled habitat preference for each species and breeding stage considering the following candidate covariates: SST (°C), PP (mg C m⁻²day⁻¹), dySST (°C km⁻¹), minimum distance to the TZCF (dTZCF; km), wind stress curl (Pa m⁻¹), SSHa (cm), EKE (cm² s⁻²), wind speed (m s⁻¹), wind direction (°), sea floor depth (m), and d_c (km). EKE and PP were log-transformed to achieve an even spread of covariate values, thus avoiding undue leverage of a few high values [104]. To ensure that there was not strong collinearity among parameters, we used generalized linear models with a binomial distribution and logit link to relate the response variable (presence/control) to all covariates, and calculated tolerance values (inverse of the variance inflation factor) for each species-stage; all tolerance values were greater than 0.1, the approximate guide suggested by Quinn and Keough [105].

GAMMs were implemented within the gamm4 package in R [106]. All candidate covariates were fitted using cubic regression splines; a cyclic spline was used to model wind direction. We first investigated the potential effects of transmitter duty-cycle and tag type (geolocation versus satellite tracking), as well as year and sex effects, by comparing intercept-only GAMMs for each species and reproductive stage to models including each of the above terms. In all cases, Akaike's information criterion (AIC) values were not improved by the addition of these terms, therefore they were not subsequently included in models of habitat preference.

To arrive at an inferential model for each species and breeding stage, we used AIC as a guide, and relied on cluster-level cross-validation for final model selection [80]. This approach was taken due to the fact that model selection based on AIC alone can lead to overparameterized models when data are spatially or temporally autocorrelated, as is expected with both tracking and oceanographic data [80]. Therefore, we randomly selected two-thirds of the individuals tracked during each breeding stage for each species to use in the first stage of model selection. We started with intercept-only models and used forward model selection with AIC to arrive at a suite of candidate models for each species-stage. We then fit data from the remaining individuals to the sequence of models obtained from forward selection. The model with the lowest AIC from the cross-validation step was retained as the inferential model. To reduce the chance of over-fitting, we then replaced each
spline in turn by a linear term and selected a final inferential model for each species-stage based on AIC values. Candidate GAMMs were fit using maximum likelihood to allow comparison of models with different fixed effect structures [107]; final inferential models were fit using restricted maximum likelihood estimation [108]. We used the proportion of the total deviance explained by the model to assess goodness of fit [107], and the percent deviance contributed by each coefficient to assess how much variability in the response could be explained by each main effect [109].

## Results

### Foraging movements

Both Laysan and black-footed albatrosses ranged significantly farther and were at sea for longer durations during incubation and chick-rearing, as compared to the brooding period ($P < 0.001$ for all pair-wise tests; Table 1). Although the majority of foraging trips occurred in

pelagic waters of the North Pacific, Laysan and black-footed albatrosses regularly traveled to neritic environments of the Aleutian islands and the western coast of North America, respectively, especially during the chick-rearing period (Fig. 1). Black-footed albatrosses ranged farther during chick-rearing, as compared to the incubation period, while Laysan albatrosses had significantly longer trip durations during incubation as compared to chick-rearing (Table 1). In both species there was a greater spread in the distribution of trip lengths and foraging ranges during incubation and chick-rearing compared to brooding (Fig. 3a and b). Some overlap in foraging ranges and trip durations occurred during all three stages, when individuals took shorter trips during incubation and chick-rearing (Fig. 3a and b). Laysan albatrosses changed destination bearings between reproductive stages, whereas black-footed albatrosses demonstrated consistency in direction of travel between stages (Table 1).

<figure><img src="image_7.png" /></figure>
Laysan albatrosses traveled for longer durations than black-footed albatrosses during the incubation period, however maximum distances did not differ between species during this stage (Table 1). Black-footed albatrosses traveled farther than Laysan albatrosses during chick-rearing, but species traveled similar distances during brooding. Azimuth to most distant point was significantly different between species during incubation and chick-rearing, with Laysan albatrosses demonstrating a more northerly distribution during incubation, and black-footed albatrosses demonstrating a more easterly distribution, especially during chick-rearing (Fig. 1; Table 1).

## Foraging distribution

Based on randomization tests of spatial overlap, Laysan and black-footed albatrosses demonstrated significant spatial segregation of focal (50 % UD) and core (25 % UD) foraging areas during incubation (focal: 64.7 % overlap, *P* = 0.014; core: 50.2 % overlap, *P* = 0.006) and chick-rearing (focal: 42.9 % overlap, *P* = 0.002; core: 24.8 % overlap, *P* = 0.004), but not during brooding (focal: 74.2 % overlap, *P* = 0.32; core: 72.0 % overlap, *P* = 0.35; Fig. 4). Overall foraging ranges (95 % UD) demonstrated a different pattern; during incubation as well as brooding, there was significant interspecific segregation of foraging ranges (incubation: 54.3 % overlap, *P* = 0.002; brooding: 51.0 % overlap, *P* = 0.012), but no spatial segregation of foraging ranges during chick-rearing (68.6 % overlap, *P* = 0.14; Fig. 4).

## Activity patterns at sea

Laysan albatrosses spent a greater proportion of time in flight during brooding (85.2 %), as compared to the incubation (74.0 %) and chick-rearing periods (78.5 %; Table 2). This difference was reflected in both the percent time spent in flight during the day and at night (Table 2). Black-footed albatrosses did not vary in terms of the proportion of time in flight between reproductive stages (82.6 % overall); however, landing rates were higher during chick-rearing (0.65 landings h⁻¹) as compared to the incubation period (0.52 landings h⁻¹; Table 2). This difference was accounted for by the greater frequency of landings during the day during chick-rearing (0.84 landings h⁻¹) versus incubation (0.62 landings h⁻¹); there was no significant difference in landing rate at night between reproductive stages (0.46 landings h⁻¹ overall). Laysan albatrosses did not demonstrate differences in landing rate as a function of breeding stage (0.84 landings h⁻¹ overall), either during the day (0.96 landings h⁻¹ overall) or at night (0.72 landings h⁻¹ overall).

During incubation, black-footed albatrosses spent a greater proportion of time in flight compared to Laysan

albatrosses (Table 2). Percent time in flight did not differ between species during daylight hours; however, black-footed albatrosses spent more time in flight at night during incubation compared to Laysan albatrosses (Table 2). Laysan albatrosses landed on the sea surface more frequently on average than black-footed albatrosses during incubation and chick-rearing. Laysan albatross landing rates during the day were more frequent than those of black-footed albatrosses during the incubation period only; landing rates at night were more frequent for

<figure><img src="image_8.png" /></figure>
<table><caption>**Table 2** Summary characteristics (Mean ± SD) of at-sea activity patterns of Laysan and black-footed albatrosses. To reduce the influence of individuals tracked for multiple foraging trips, a single trip per individual was randomly-selected for each reproductive stage to include in the calculation of mean values</caption><thead><tr><th>Species</th><th></th><th>Incubation</th><th>Brooding</th><th>Chick-Rearing</th></tr></thead><tbody><tr><td rowspan="6">Laysan albatross</td><td>Time in flight (%)</td><td>74.0 ± 10.9<sup>a,*</sup></td><td>85.2 ± 10.4<sup>b</sup></td><td>78.5 ± 12.3<sup>a</sup></td></tr><tr><td>Day time in flight (%)</td><td>83.6 ± 8.68<sup>a</sup></td><td>90.1 ± 7.6<sup>b</sup></td><td>86.8 ± 9.08<sup>a</sup></td></tr><tr><td>Night time in flight (%)</td><td>66.0 ± 15.9<sup>a,*</sup></td><td>80.4 ± 15.8<sup>b</sup></td><td>66.8 ± 19.0<sup>a</sup></td></tr><tr><td>Landings per hour</td><td>0.85 ± 0.37<sup>*</sup></td><td>0.85 ± 0.44</td><td>0.82 ± 0.44<sup>*</sup></td></tr><tr><td>Day landings per hour</td><td>0.93 ± 0.33<sup>*</sup></td><td>0.98 ± 0.49</td><td>1.01 ± 0.49</td></tr><tr><td>Night landings per hour</td><td>0.77 ± 0.49<sup>*</sup></td><td>0.71 ± 0.51<sup>*</sup></td><td>0.64 ± 0.45<sup>*</sup></td></tr><tr><td rowspan="6">Black-footed albatross</td><td>Time in flight (%)</td><td>81.5 ± 8.26<sup>*</sup></td><td>84.8 ± 10.2</td><td>81.3 ± 8.05</td></tr><tr><td>Day time in flight (%)</td><td>84.9 ± 8.43</td><td>89.1 ± 7.72</td><td>86.9 ± 7.28</td></tr><tr><td>Night time in flight (%)</td><td>78.6 ± 10.1<sup>*</sup></td><td>80.5 ± 15.7</td><td>73.0 ± 12.2</td></tr><tr><td>Landings per hour</td><td>0.52 ± 0.13<sup>a,*</sup></td><td>0.71 ± 0.49<sup>ab</sup></td><td>0.65 ± 0.17<sup>b,*</sup></td></tr><tr><td>Day landings per hour</td><td>0.62 ± 0.15<sup>a,*</sup></td><td>0.92 ± 0.53<sup>ab</sup></td><td>0.84 ± 0.25<sup>b</sup></td></tr><tr><td>Night landings per hour</td><td>0.41 ± 0.16<sup>*</sup></td><td>0.52 ± 0.56<sup>*</sup></td><td>0.46 ± 0.17<sup>*</sup></td></tr></tbody></table>

Different lowercase letters indicate significant differences (P < 0.05) among reproductive stages; asterisks indicate significant differences between species

Laysan albatrosses during all three breeding stages (Table 2).

Both species demonstrated diel patterns in foraging activity during all reproductive stages (Table 2). Laysan albatrosses spent a higher proportion of time in flight during the day than at night (incubation: $t_{100} = 6.70$, $P < 0.0001$; brooding: $t_{74} = 3.51$, $P = 0.0008$; chick-rearing: $t_{50} = 5.10$, $P < 0.0001$) and landed more frequently during the day (incubation: $t_{100} = 3.08$, $P = 0.003$; brooding: $t_{74} = 2.97$, $P = 0.004$; chick-rearing: $t_{50} = 3.83$, $P = 0.0004$). Black-footed albatrosses demonstrated the same pattern; they spent more time in flight during the day (incubation:

$t_{86} = 3.19$, $P = 0.002$; brooding: $t_{68} = 2.40$, $P = 0.02$; chick-rearing: $t_{44} = 4.43$, $P < 0.0001$) and landed more frequently during daylight than at night (incubation: $t_{86} = 6.51$, $P < 0.0001$; brooding: $t_{68} = 4.90$, $P = 0.0008$; chick-rearing: $t_{44} = 6.06$, $P < 0.0001$).

## Habitat preference

Final inferential habitat preference models included environmental covariates descriptive of regional oceanic production (SST and PP), large-scale fronts ($d_{TZCF}$), mesoscale activity (EKE and SSHa), and sea floor depth (Table 3). Environmental covariates not retained in final

<table><caption>**Table 3** Model selection results from generalized additive mixed modeling of Laysan and black-footed albatross habitat preference. Final models were arrived at by forward selection based on AIC using a subset of individuals, followed by cross-validation using the remaining individuals</caption><thead><tr><th></th><th></th><th colspan="6">Model terms [% deviance explained by coefficient]</th><th>% deviance explained by model</th></tr></thead><tbody><tr><td rowspan="2">Laysan albatross</td><td rowspan="2">Incubation</td><td>s(SST)</td><td>s(d<sub>c</sub>)</td><td>s(deepth)</td><td>SSHa</td><td>s(d<sub>TZCF</sub>)</td><td>EKE</td><td>39.8 %</td></tr><tr><td>[60.6]</td><td>[30.6]</td><td>[5.34]</td><td>[1.79]</td><td>[0.90]</td><td>[0.75]</td><td></td></tr><tr><td rowspan="2">Brooding</td><td rowspan="2"></td><td>d<sub>c</sub></td><td>PP</td><td></td><td></td><td></td><td></td><td>26.3 %</td></tr><tr><td>[82.0]</td><td>[18.0]</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Chick-rearing</td><td rowspan="2"></td><td>s(SST)</td><td>s(d<sub>c</sub>)</td><td>s(deepth)</td><td>s(PP)</td><td>s(d<sub>TZCF</sub>)</td><td></td><td>39.2 %</td></tr><tr><td>[45.7]</td><td>[45.1]</td><td>[4.08]</td><td>[2.64]</td><td>[2.52]</td><td></td><td></td></tr><tr><td rowspan="2">Black-footed albatross</td><td rowspan="2"></td><td rowspan="2">Incubation</td><td>s(SST)</td><td>s(d<sub>c</sub>)</td><td>SSHa</td><td></td><td></td><td>38.8 %</td></tr><tr><td>[75.9]</td><td>[21.9]</td><td>[2.18]</td><td></td><td></td><td></td></tr><tr><td rowspan="2">Brooding</td><td rowspan="2"></td><td>s(SST)</td><td>d<sub>c</sub></td><td></td><td></td><td></td><td></td><td>11.9 %</td></tr><tr><td>[60.5]</td><td>[39.5]</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Chick-rearing</td><td rowspan="2"></td><td>s(SST)</td><td>s(d<sub>c</sub>)</td><td>s(d<sub>TZCF</sub>)</td><td>s(EKE)</td><td>s(PP)</td><td>depth</td><td>28.3 %</td></tr><tr><td>[45.5]</td><td>[19.8]</td><td>[12.2]</td><td>[11.2]</td><td>[6.44]</td><td>[4.84]</td><td></td></tr></tbody></table>

Cubic regression splines retained in the final models are represented by s() [84]. SST: sea surface temperature. d<sub>c</sub>: distance to the breeding colony. SSHa: sea surface height anomaly; d<sub>TZCF</sub>: distance to the Transition Zone Chlorophyll Front; EKE: eddy kinetic energy; PP: primary productivity
models were: dySST, wind stress curl, and wind speed or direction. SST was retained as a smooth term, and explained the greatest proportion of explained deviance in final habitat models selected for all reproductive stages for black-footed albatrosses, and during incubation and chick-rearing for Laysan albatrosses (46-76%; Table 3).

The range of accessible thermal habitats was considerably wider during incubation and chick-rearing (~0-31 °C), as compared to brooding for both species (~15-26 °C; Fig. 5a and 6a). Laysan and black-footed albatrosses generally demonstrated a preference for cooler water temperatures compared to availability (Fig. 5a and 6a), especially during incubation and chick-rearing. During these two breeding stages, Laysan and black-footed albatrosses demonstrated differing response curves between SST and habitat preference, after accounting for the effects of other covariates in the final habitat models (Fig. 7a and c; Fig. 8a and c). Laysan albatross preference was highest at a broad range of cool water temperatures (~0-12 °C), decreasing steeply at temperatures greater than ~20 °C (Fig. 7a and c). Black-footed albatross habitat preference was highest at intermediate temperatures, peaking at ~14 °C during incubation and ~7 °C during chick-rearing (Fig. 8a and c).

All final habitat preference models included $d_c$, which was retained as a smooth term during incubation and chick-rearing, but as a linear term during brooding for both species (Table 3). In all cases, $d_c$ had a negative slope (Fig. 7 and 8), indicating that spatial use of foraging locations decreased at a greater rate than $1/d_c$. Distance to the colony accounted for the greatest proportion of explained deviance in the final habitat model for Laysan albatrosses during brooding, and the second greatest proportion of explained deviance in all other species-stages (20-82%; Table 3). During incubation and chick-rearing, both Laysan and black-footed albatrosses demonstrated a preference for habitats within ~2000 km of the Tern Island colony (Fig. 7 and 8).

For Laysan albatrosses during brooding, PP was the only environmental covariate retained in the final GAMM; Laysan albatrosses during this stage demonstrated a preference for more productive habitats (Fig. 7b), however, the slope of this relationship was fairly shallow. Generally, a narrower range of primary production was available to foraging albatrosses during brooding (~150-600 $mg C m^{-2} day^{-1}$) as compared to incubation (~30-1000 $mg C m^{-2} day^{-1}$) and chick-rearing; the most productive habitats were available during chick-rearing (~50-2000 $mg C m^{-2} day^{-1}$; Fig. 5b and 6b). During chick-rearing, PP was retained as a smooth term in the final habitat models for both species; preference was generally highest at intermediate levels (~350-1000 $mg C m^{-2} day^{-1}$; Fig. 7c and Fig. 8c), although black-footed albatrosses preference also

increased at the highest productivity values (>1350 $mg C m^{-2} day^{-1}$; Fig. 8c).

Laysan albatrosses during incubation and chick-rearing, and black-footed albatrosses during chick-rearing, demonstrated a relationship between $d_{TZCF}$ and habitat preference (Table 3). Preference generally decreased after distances of ~2000-3000 km (Fig. 7a and c; Fig. 8c), indicating that this response may be related to the northerly bias in albatross tracking locations (Fig. 1) rather than a tight association with the TZCF itself. The percent explained deviance for this variable was generally low for Laysan albatrosses (<3%), but explained 12.2% of the deviance in the final chick-rearing habitat model for black-footed albatrosses (Table 3).

EKE accounted for the next highest proportion of explained deviance in the black-footed albatross habitat model during chick-rearing (Table 3). Preference was for lower values of EKE (Fig. 8c), contrary to our expectation. EKE was also retained in the final model for Laysan albatrosses during incubation, and demonstrated a positive slope (Fig. 7a), however, this term accounted for a very small proportion of explained deviance (0.75%; Table 3). SSHa, another indicator of mesoscale activity, was retained in the final incubation models for both species, but accounted for a relatively small proportion of explained deviance (~2%; Table 3); in both cases preference had a positive relationship with SSHa (Fig. 7a and 8a).

Sea floor depth was retained as a smooth term in the incubation and chick-rearing models for Laysan albatrosses, and as a linear term in the chick-rearing model for black-footed albatrosses (Table 3). Black-footed albatrosses demonstrated a general preference for shallower waters (Fig. 8c), whereas Laysan albatrosses demonstrated a U-shaped response, with a preference for both the shallowest (<2000 m) and deepest habitats (>5000 m; Fig. 7a and c). Sea floor depth accounted for 4-5% of the explained deviance in final habitat models for these species-stages (Table 3).

# Discussion

## Foraging movements

Temporal changes in foraging behavior have been linked to changing energetic demands of offspring in a variety of central place foragers [23,24,110-112]. Consistent with previous research, Laysan and black-footed albatrosses spent more time at sea and ranged farther from the colony during incubation and chick-rearing than during brooding [26, 28, 29, 32]. Albatrosses have the ability to make long-distance movements during these breeding stages due to the fasting capabilities of adults [incubation spans can last as long as 50 days; [113] and post-guard chicks, which fast between intermittent feeds
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Kernel density of environmental covariates (a-f) at Laysan albatross tracking and control locations. Only environmental covariates retained in final habitat preference models for the species (Table 3) are provided. SST: sea surface temperature; PP: primary productivity; TZCF: Transition Zone Chlorophyll Front; SSHa: sea surface height anomaly; and EKE: eddy kinetic energy</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 6** Kernel density of environmental covariates (a-f) at black-footed albatross tracking and control locations. Only environmental covariates retained in final habitat preference models for the species (Table 3) are provided. SST: sea surface temperature; PP: primary productivity; TZCF: Transition Zone Chlorophyll Front; SSHa: sea surface height anomaly; and EKE: eddy kinetic energy</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Fig. 7** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Fig. 7** Effects of covariates in final GAMMs for Laysan albatrosses during incubation (a), brooding (b), and chick-rearing (c). The contribution of each retained covariate to the linear predictor is plotted on the scale of the link function (y-axes); the plots can therefore be interpreted as population-level habitat preferences [80]. Dashed lines indicate approximate 95 % confidence intervals. GAMM: generalized additive mixed model; SST: sea surface temperature; SSHa: sea surface height anomaly; TZCF: Transition Zone Chlorophyll Front; EKE: eddy kinetic energy; and PP: primary productivity

[40]. During brooding, adults are limited to short trips when young chicks require frequent meals [21].

When conducting long-ranging trips, Laysan albatrosses generally foraged at more northerly latitudes compared to black-footed albatrosses. Ommastrephid squid, which comprise the largest proportion of the Laysan albatross diet, also have a more northerly distribution in the NPTZ compared to flying fishes [90], the eggs of which comprise the largest proportion of the black-footed albatross diet [86]. During chick-rearing, black-footed albatrosses traveled farther north and east than during incubation, reaching greater distances than Laysan albatrosses. The northward shift may be a response to seasonal changes in the distribution of prey associated with the NPTZ [50]; the eastward shift is likely a response to prey aggregations associated with productive coastal upwelling along the western coast of North America during this time of year [114, 115]. Although flying fish eggs dominate the diet of black-footed albatrosses during chick-rearing [86], we expect that this reflects foraging events in subtropical waters, and that squid and other prey types would be consumed by black-footed albatrosses foraging in the California Current.

## Interspecific habitat segregation

We found that Laysan and black-footed albatrosses demonstrated spatial segregation of core and focal foraging areas during the incubation and chick-rearing periods, which agrees with previous studies conducted during both the breeding [28, 29, 32] and non-breeding periods [116]; Shaffer et al., unpublished data]. However, habitat segregation between sympatric species can vary depending on reproductive stage, and the strength of these differences may be affected by prey availability, energetic constraints, or changes in interspecific encounters [117–120]. During brooding, when reproductive constraints restrict movements, overlap of core and focal foraging areas was greatest [see also [29]. Although our results showed that overall foraging ranges (95 % UD) were spatially segregated during brooding, this was likely due to a few longer-distance trips made by Laysan albatrosses, reflected in the high variability of maximum distances reached during brooding (Table 1).

Greater dietary partitioning during brooding could provide a means of reducing competitive interactions due to shared feeding areas. There is limited information,

however, on changes in diet throughout the reproductive period in these species; in other albatross species, diet can change considerably over the course of the breeding season [121]. The most comprehensive dietary dataset for Hawaiian albatrosses are mainly representative of the chick-rearing period and are not separated for each stage [86]. Other diet studies have focused on by-caught birds in fisheries [87] or analysis of regurgitated pellets from chicks during the rearing period only [122]. Comparison of a limited number of chick diets from the brooding period suggests some dietary segregation during this stage; squid and fish were found more frequently in the diet of Laysan albatross chicks, and flying fish eggs were found more frequently in the diet of black-footed albatross chicks (M.A. Kappes, unpublished data).

## Activity patterns at sea

Contrary to our expectation, neither species changed the frequency of landings on the sea surface during the brooding period to meet higher energy requirements during this stage. Although black-footed albatrosses landed more frequently on average during brooding, differences were not statistically significant due to high variability in landing rate during this stage (Table 2). Laysan albatrosses did, however, spend a greater proportion of time in flight during brooding as compared to the incubation and chick-rearing periods. Weimerskirch et al. [123] suggests that periods of sustained movement in wandering albatrosses is an adaptation to locating dispersed and unpredictable prey, a strategy which may be relevant to Hawaiian albatrosses foraging in a warm, oligotrophic environment during brooding.

The brooding period is also when core foraging areas of Laysan and black-footed albatrosses overlap the most, so species differences in timing of foraging activities could minimize competition for prey resources during this stage. It has been proposed that higher levels of rhodopsin in the eyes of Laysan albatrosses make them better adapted to nocturnal feeding than black-footed albatrosses [86]. Consistent with this idea, we found that Laysan albatrosses landed more frequently at night than black-footed albatrosses during all breeding stages, but both species actively foraged at night. Overall, there was a bias towards foraging during daylight hours for both species, which is supported by other studies [124, 125]. During brooding, time in flight and daytime landing rates were similar between species, and although landing
<figure><img src="image_12.png" /><figcaption>**Fig. 8** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Fig. 8** Effects of covariates in final GAMMs for black-footed albatrosses during incubation (a) brooding (b) and chick-rearing (c). The contribution of each retained covariate to the linear predictor is plotted on the scale of the link function (y-axes); the plots can therefore be interpreted as population-level habitat preferences [80]. Dashed lines indicate approximate 95 % confidence intervals. GAMM: generalized additive mixed model; SST: sea surface temperature; SSHa: sea surface height anomaly; TZCF: Transition Zone Chlorophyll Front; EKE: eddy kinetic energy; and PP: primary productivity

rates differed at night, this was the case during all breeding stages. Therefore, we did not find evidence that Laysan and black-footed albatrosses changed their behavior during brooding as a means of reducing competitive interactions when spatial overlap was greatest [126].

Hawaiian albatrosses generally spend more time in flight [74-91 %; current study, [124], particularly during the brooding period, than other albatross species [44-69 %; [127, 128]. This may be adaptive for foraging in an oligotrophic environment, where prey resources are likely unpredictable [123]. Tropical seabirds that feed in lower productivity environments generally have a greater flight proficiency to allow travel between dispersed prey patches [38]. Available research on energetics of Laysan albatrosses does not suggest lower energy expenditure rates compared to other albatross species during the incubation period [129], however, a comparison to energy expenditure rates during brooding could provide insight as to whether increased time in flight reduces energetic costs when foraging in a low-productivity environment.

## Habitat preference

We combined satellite tracking data with geolocation data when modeling habitat preference, in order to obtain coverage of the entire breeding period. However, this could lead to potential biases given the different accuracies of each tag type (satellite tags are accurate to <10 km [53]; geolocation tags are accurate to ~200 km [54]). To address this potential bias, we investigated the effects of tag type on habitat models by comparing intercept-only GAMMs for each species and reproductive stage to models including tag type as a candidate variable. Because AIC values were not improved by the addition of this term, differences in tag accuracy should not have affected results of habitat preference modeling. We also did not find an effect of gender on habitat preference, which agrees with previous studies of these species [28, 29, 130].

SST accounted for much of the variability in final habitat models during incubation and chick-rearing, when both species preferentially selected cooler waters compared to availability. Similarly, when using a different technique to model habitat use of Hawaiian albatrosses during incubation, Kappes et al. [32] found that SST was the most important environmental variable related to first passage time [131], a measure of area-restricted

search [132]. Overall, black-footed albatrosses preferred slightly warmer marine habitats than Laysan albatrosses; this likely relates to differences in diet and spatial separation of prey resources along the latitudinal SST gradient in the NPTZ [86, 90].

We hypothesized that each species would use the same environmental covariates to select habitats throughout the breeding season, but that reproductive constraints and seasonal cycles in marine habitats would influence the composition of utilized habitats. SST accounted for the greatest proportion of variability in black-footed albatross habitat models during all reproductive stages and for Laysan albatross habitat models during incubation and chick-rearing; however, SST was not retained in the final habitat model for Laysan albatrosses during brooding. PP was the only environmental covariate retained in the brooding habitat model for Laysan albatrosses, but the majority of variability in habitat preference was explained by $d_c$. Both species demonstrated less distinct habitat preferences during brooding; preference increased only slightly with PP for Laysan albatrosses, and black-footed albatrosses had a smaller peak in preference in relation to SST compared to incubation and chick-rearing. This reflects the narrower range of environmental conditions available for selection during brooding, and helps explain why model goodness of fit was lower during brooding than during other breeding stages (Table 3).

Given the constraint of central place foraging during all stages of breeding, we included $d_c$ as a candidate covariate when modeling habitat preference. All final habitat models included $d_c$, suggesting that accessibility decreased at a rate greater than $1/d_c$ [80-82], and that the null usage model could be further refined. The inclusion of $d_c$ in final habitat models could also indicate a preference for resources closer to the colony that are not linked to the covariates examined. Albatrosses use gliding flight on sea surface winds for locomotion [98], so an important consideration may be the effect of wind on habitat accessibility [133-136]. Because Hawaiian albatrosses tend to avoid facing headwinds during long-distance movements [137], accounting for differences in wind direction might improve the null usage model. Overall accessibility of habitats may be affected by prevailing winds, but the fact that wind variables were not retained in final habitat models suggests that Laysan and black-footed albatrosses do not select foraging locations
based on wind conditions. Alternatively, wind variables examined may not sufficiently represent the complex manner in which albatrosses use surface winds.

Together, SST and $d_c$ accounted for most of the explained deviance in final habitat models for both species. For black-footed albatrosses during chick-rearing, $d_{TZCF}$ and EKE also accounted for considerable (>10%) proportions of explained deviance. Habitat preference declined after distances of ~2000 km from the TZCF, indicating that this response may be related to the northerly bias in tracking locations rather than a tight association with the TZCF itself. Although we expected that a higher degree of mesoscale activity would enhance foraging opportunities, black-footed albatrosses preferred habitats with lower EKE during chick-rearing. Their use of the pelagic northeast Pacific, a region of relatively low EKE [138], suggests that black-footed albatrosses do not rely on prey associated with high eddy activity during the chick-rearing period; leatherback turtles (*Dermochelys coriacea*) in the North Pacific also demonstrate a negative relationship between habitat use and EKE [139]. However, low accuracy of geolocation data may have limited our ability to detect responses to mesoscale activity.

The environmental covariate with the greatest influence on habitat preference of Laysan and black-footed albatrosses was SST. Sea floor depth, SSHa, PP, $d_{TZCF}$ and EKE were also retained as covariates in final habitat models, however, these terms represented a smaller proportion of explained deviance. This suggests that after accounting for SST preference, the location of large-scale fronts ($d_{TZCF}$), mesoscale activity (EKE and SSHa), and bathymetric domain do not have a strong influence on habitat selection of Hawaiian albatrosses. Other environmental covariates describing frontal characteristics (dySST) and wind-driven oceanic upwelling or downwelling (wind stress curl) were not retained in final models, despite characterizing physical features or processes that could drive prey availability for Hawaiian

albatrosses [88–91,93]. Studies of other marine top predators have also found the greatest behavioral responses to SST compared to other marine habitat variables (Simmons et al., 2007; Teo et al., 2007). Environmental covariates examined might also be improved by incorporating time lags, as the response of prey resources may trail behind increases in productivity.

## Breeding in an oligotrophic marine environment

Laysan and black-footed albatrosses differ from other albatross species in that they breed in an oligotrophic marine environment [140, 141]. They leave this environment during incubation and chick-rearing to forage in cooler, more northern waters, but are restricted to a low-productivity environment during brooding when energy requirements are greatest [21]. For wandering albatrosses, the brooding period coincides with seasonal increases in prey abundance in the Southern Ocean [15]. Little is known about temporal variability of prey consumed by Laysan and black-footed albatrosses, but it has been suggested that these species breed during winter to match peaks in prey abundance [142], or alternatively to avoid the period of maximum summer heat [113]. PP near the Tern Island colony is lower during brooding compared to the rest of the year (Fig. 9), so assuming PP is a good indicator of prey abundance, local prey resources may in fact be limited during the brooding period. Timing of breeding may instead be related to the relative proximity of preferred thermal habitats; the latitudinal position of the TZCF, and the cooler waters associated with it, are closest to the Tern Island colony during the brooding and early chick-rearing periods [50].

Body mass of adult Laysan albatrosses declines throughout the breeding season, from the start of the incubation period until the second half of the chick-rearing period [143]. This suggests a high cost of central place foraging from breeding sites in the Hawaiian Islands; other albatross species do not demonstrate equivalent declines in adult mass during breeding

<figure><img src="image_13.png" /><figcaption>**Fig. 9** Primary productivity (PP) near the Tern Island colony. Median monthly PP was calculated within the mean foraging range of Laysan and black-footed albatrosses during brooding (370 km); brooding periods are indicated by solid lines at the bottom of the figure</figcaption></figure>
[40, 144, 145]. The change to positive mass gain in mid-May may be related to seasonal increases in PP close to the colony (Fig. 9) and the closer proximity of preferred habitats of the NPTZ at this time [50]; the shorter commuting distance to preferred habitats could help adults restore body condition.

During chick-rearing, many Procellariform species alternate between short foraging trips that maximize energy delivery to young, with long trips that restore adult condition [146, 147]; this is reflected in higher food delivery rates after short trips compared to long trips [reviewed by [148]. Previous research has demonstrated that Hawaiian albatrosses mix short and long trips during chick-rearing [28, 29], though we did not observe strong bimodality in foraging trip duration during this stage (Fig. 3b). Contrary to expectation, Laysan albatrosses fledge chicks in better condition when they perform fewer short trips for every long trip during chick-rearing [149]. This suggests that shorter trips are less profitable, and that longer trips to preferred habitats of the NPTZ are required for successful chick-rearing. It is unclear, however, to what extent low food delivery rates affect chick survival. Along with food shortages, storm events and nest desertion are common causes of chick mortality [150, 151].

In addition to seasonal changes, the latitudinal position of the TZCF varies interannually, particularly in response to El Niño-Southern Oscillation events [50]. When the wintertime position of the TZCF is significantly north of its seasonal norm, Laysan and black-footed albatrosses experience dramatic breeding failures [152, 153]. Laysan albatrosses also demonstrate lower reproductive success and mass gains in years when they range farther and spend longer durations at sea [32]. This suggests that changes in the proximity of preferred habitats during breeding can lead to demographic effects in these species. Changes in SST in the North Pacific due to global climate change [154, 155] will likely have considerable impacts on Hawaiian albatrosses, given the importance of SST in foraging habitat selection. In addition, poleward shifts in westerly winds associated with climate change [156–158] may lead to changes in the positioning of the North Pacific Current [52] and access to preferred habitats. In the Southern Ocean, the poleward shift of the westerlies has positively affected breeding performance of wandering albatrosses, but a continuing shift may make the location of breeding colonies less optimal in the future [159]. A northward shift in the positioning of the TZCF due to climate change could have negative effects on Laysan and black-footed albatross populations, if preferred habitats become more distant from the colony during critical portions of the breeding season. Comparisons among breeding locations within the Hawaiian islands would provide insight into

these interactions [160], particularly at Midway Atoll and Laysan Island, where the greatest numbers of breeding pairs occur for both species [151].

## Conclusions

For central place foragers, the location of suitable breeding sites can have large effects on behavior and habitat use. Because Hawaiian albatrosses breed in an oligotrophic marine environment, they forage in lower-productivity waters than other albatross species during the energetically costly brooding period. They also spend a greater proportion of time in flight, a strategy that may be adaptive for locating unpredictable food resources. Additional research on the foraging energetics of these species could help resolve how Hawaiian albatrosses meet energetic demands during brooding, when spatial overlap between species is greatest and primary productivity near the colony is at an annual low. During incubation and chick-rearing, Laysan and black-footed albatrosses make long distance movements north of the colony, selecting cooler water temperatures and making use of the regionally productive NPTZ. A poleward shift in the TZCF due to climate change could negatively affect Laysan and black-footed albatross populations, given breeding failures in years when the TZCF is shifted significantly farther north. A comprehensive analysis of Hawaiian albatross foraging behavior and reproductive success in relation to the position of the TZCF will be necessary to predict population-level effects of climate change.

## Availability of supporting data

The tracking datasets supporting the results of this article are available by request from the BirdLife Tracking Ocean Wanderers database, [http://seabirdtracking.org].

## Additional file

**Additional file 1: Appendix A Laysan and black-footed albatross habitat preference modeling using an alternative null usage model.**

Habitat preference was modeled using a more restrictive null model of usage to reflect the northern bias of tracking locations; for each species-stage, minimum convex polygons of tracking locations were used to limit the southern extent of control locations, and maximum observed ranges were used to limit the northern extent of control locations. (DOCX 25 kb)

## Abbreviations

dE: Distance to the colony; dTZCF: Distance to the Transition Zone Chlorophyll Front; dySST: Latitudinal gradient in sea surface temperature; EKE: Eddy kinetic energy; GAMM: Generalized additive mixed model; NPTZ: North Pacific Transition Zone; PP: Primary productivity; SSHa: Sea surface height anomaly; SST: Sea surface temperature; TZCF: Transition Zone Chlorophyll Front; UD: Utilization distribution.

## Competing interests

The authors declare that they have no competing interests.
## Authors' contributions

All authors participated in the design of the study and contributed to the writing of the manuscript. DC and SS conceived of and coordinated the study. MK and YT collected tracking data in the field. MK, DF, DP, and SB worked on acquisition and interpretation of ocean data and development of habitat models; MK performed statistical analyses. All authors read and approved the final manuscript.

## Acknowledgements

The authors dedicate this manuscript to the memory of our colleague and friend, David Foley. We thank Angela Anders, Jill Awkerman, Daniel Barton, Sarah Chisholm, Melinda Conners, Kirsten Lindquist, Frank Mayer, Scott Seganti, and many U.S. Fish and Wildlife Service volunteers for invaluable assistance in the field. We also thank the Hawaiian Islands National Wildlife Refuge, U.S. Fish and Wildlife Service, Department of the Interior, for permission to conduct research on Tern Island and logistical support. Lindsay Young conducted the gender analysis. Alan Swithenbank and Steve Teo provided the SST-correction algorithm for light-based geolocation data. Jason Matthiopoulos and Ewan Wakefield provided helpful suggestions for defining accessibility and selecting control locations. SeaWiFS chlorophyll-a and photosynthetically available radiation were provided by NASA's Goddard Space Flight Center. QuikSCAT wind data were provided by NASA's Jet Propulsion Laboratory at the California Institute of Technology. Pathfinder SST data were obtained from NOAA's National Oceanographic Data Center. AMSR-E SST data were obtained from Remote Sensing Systems Inc. GOES SST data were obtained from NOAA Satellites and Information. Sea surface height anomaly data were obtained from AVISO/Centre National d'Études Spatiales, and ETOPO2 data from NOAA's National Geophysical Data Center. This research was part of the Tagging of Pacific Pelagics (TOPP) program, funded in part by the National Ocean Partnership Program (N00014-02-1-1012), the Office of Naval Research (N00014-00-1-0880 & N00014-03-1-0651), the Gordon and Betty Moore, David and Lucille Packard, and Alfred P. Sloan Foundations. M.A.K. was supported by a UCSC Regents Fellowship and Chancellor's Dissertation-Year Fellowship, a Friends of Long Marine Lab Research Award, and a STEPS Institute Graduate Research Award. All protocols employed in this study were approved by the Institutional Animal Care and Use Committees at UCSC. Henri Weimerskirch and Bruce Lyon provided thoughtful comments and suggestions which helped improve the quality of this manuscript. M.A.K. dedicates this work to the memory of her father, Stanislav Antolos.

## Author details

Received: 1 December 2014 Accepted: 9 September 2015

Published online: 29 September 2015

## References

1. Emen JM. The role of time and energy in food preference. Am Nat. 1966;100:611-7.

2. MacArthur RH, Pianka ER. On optimal use of a patchy environment. Am Nat. 1966;100:603-9.

3. Levins R. Evolution in changing environments: some theoretical explorations. Princeton University Press: Princeton; 1968.

4. Pyke GH. Optimal foraging theory: a critical review. Annu Rev Ecol Syst. 1984;15:523-75.

5. Svärdson G. Competition and habitat selection in birds. Oikos. 1949;1:157-74.

6. Charnov EL. Optimal foraging, the marginal value theorem. Theor Popul Biol. 1976;9:129-36.

7. Alerstam T, Högstedt G. Bird migration and reproduction in relation to habitats for survival and breeding. Ornis Scand. 1982;13:25-37.

8. Kamil AC, Roitblat HL. The ecology of foraging behavior: implications for animal learning and memory. Annu Rev Psychol. 1985;36:141-69.

9. Huey RB. Physiological consequences of habitat selection. Am Nat. 1991;137:S91-S115.

10. Rosenzweig ML. Habitat selection and population interactions: the search for mechanism. Am Nat. 1991;137:S5-S28.

11. Morris DW. Toward an ecological synthesis: a case for habitat selection. Oecologia. 2003;136:1-13.

12. Orians GH, Pearson NE. On the theory of central place foraging. In: Horn DJ, Stairs GR, Mitchell RD, editors. Analysis of ecological systems. Columbus, Ohio: Ohio State University Press; 1979. p. 155-77.

13. Schoener TW. Generality of the size-distance relation in models of optimal feeding. Am Nat. 1979;114:902-14.

14. Costa DP. The relationship between reproductive and foraging energetics and the evolution of the Pinnipedia. Symp Zool Soc Lond. 1993;66:293-314.

15. Salamolard M, Weimerskirch H. Relationship between foraging effort and energy requirement throughout the breeding season in the wandering albatross. Funct Ecol. 1993;7:643-52.

16. McCafferty DJ, Boyd IL, Walker TR, Taylor RL. Foraging responses of Antarctic fur seals to changes in the marine environment. Mar Ecol Prog Ser. 1998;166:285-99.

17. Hochscheid S, Godley BJ, Broderick AC, Wilson RP. Reptilian diving: highly variable dive patterns in the green turtle *Chelonia mydas*. Mar Ecol Prog Ser. 1999;185:101-12.

18. Shaffer SA, Costa DP, Weimerskirch H. Foraging effort in relation to the constraints of reproduction in free-ranging albatrosses. Funct Ecol. 2003;17:66-74.

19. Humphreys EM, Wanless S, Bryant DM. Stage-dependent foraging in breeding black-legged Kittiwakes *Rissa tridactyla*: distinguishing behavioural responses to intrinsic and extrinsic factors. J Avian Biol. 2006;37:436-46.

20. Boersma P, Rebstock G. Foraging distance affects reproductive success in Magellanic penguins. Mar Ecol Prog Ser. 2009;375:263-75.

21. Ricklefs RE. Some considerations on the reproductive energetics of pelagic seabirds. Stud Avian Biol. 1983;8:84-94.

22. Costa DP. Reproductive and foraging energetics of high latitude penguins, albatrosses and pinnipeds: implications for life history patterns. Am Zool. 1991;31:111-30.

23. Costa DP, Gales NJ. Energetics of a benthic diver: seasonal foraging ecology of the Australian sea lion, Neophoca cinerea. Ecol Monogr. 2003;73:27-43.

24. Beauplet G, Dubroca L, Guinet C, Cherel Y, Dabin W, Gagne C, et al. Foraging ecology of subantarctic fur seals *Arctocephalus tropicalis* breeding on Amsterdam Island: seasonal changes in relation to maternal characteristics and pup growth. Mar Ecol Prog Ser. 2004;273:211-25.

25. Shaffer SA. Annual energy budget and food requirements of breeding wandering albatrosses (*Diomedea exulans*). Polar Biol. 2004;27:253-6.

26. Weimerskirch H, Salamolard M, Sarrazin F, Jouventin P. Foraging strategy of wandering albatrosses through the breeding season: a study using satellite telemetry. Auk. 1993;110:325-42.

27. Guinet C, Koudil M, Bost C-A, Durbec JP, Georges JY, Mouchot MC, et al. Foraging behaviour of satellite-tracked king penguins in relation to sea-surface temperatures obtained by satellite telemetry at Crozet Archipelago, a study during three austral summers. Mar Ecol Prog Ser. 1997;150:11-20.

28. Fernández P, Anderson DJ, Sievert PR, Huyvaert KP. Foraging destinations of three low-latitude albatross (*Phoebastria*) species. J Zool. 2001;254:391-404.

29. Hyrenbach KD, Fernandez P, Anderson DJ. Oceanographic habitats of two sympatric North Pacific albatrosses during the breeding season. Mar Ecol Prog Ser. 2002;233:283-301.

30. Inchausti P, Guinet C, Koudil M, Durbec J-P, Barbraud C, Weimerskirch H, et al. Inter-annual variability in the breeding performance of seabirds in relation to oceanographic anomalies that affect the Crozet and the Kerguelen sectors of the Southern Ocean. J Avian Biol. 2003;34:170-6.

31. Pinaud D, Cherel Y, Weimerskirch H. Effect of environmental variability on habitat selection, diet, provisioning behaviour and chick growth in yellow-nosed albatrosses. Mar Ecol Prog Ser. 2005;298:295-304.

32. Kappes MA, Shaffer SA, Tremblay Y, Foley DG, Palacios DM, Robinson PW, et al. Hawaiian albatrosses track interannual variability of marine habitats in the North Pacific. Prog Oceanogr. 2010;86:246-60.

33. Whittow GC. Black-footed albatross (*Diomedea nigripes*). In: Poole A, Gill F, editors. The Birds of North America. Volume 65. Philadelphia: The Academy of Natural Sciences; Washington, D.C.: The American Ornithologists' Union; 1993.

34. Whittow GC. Laysan albatross (*Diomedea immutabilis*). In: Poole A, Gill F, editors. The Birds of North America. Volume 66. Philadelphia: The

<footnote n="1">¹Department of Ecology and Evolutionary Biology, University of California Santa Cruz, 100 Shaffer Road, Santa Cruz, California 95060, USA.</footnote>

<footnote n="2">²Environmental Research Division, Southwest Fisheries Science Center, NOAA Fisheries, 1352 Lighthouse Avenue, Pacific Grove, California 93950, USA.</footnote>

<footnote n="3">³Joint Institute for Marine and Atmospheric Research, University of Hawai'i at Manoa, 1000 Pope Road, Honolulu, Hawai'i 96822, USA. ⁴Present address: Department of Fisheries and Wildlife, Oregon State University, 104 Nash Hall, Corvallis, Oregon 97331, USA. ⁵Present address: Department of Biological Sciences, San Jose State University, One Washington Square, San Jose, California 95192, USA. ⁶Present address: Institut pour la Recherche et le Développement, UMR 212 IRD-IFREMER-UM2. Av Jean Monnet 34200, Sète, France.</footnote>
Academy of Natural Sciences; Washington, D.C.: The American Ornithologists' Union; 1993.

35. Tickell WLN. Albatrosses. New Haven: Yale University Press; 2000.

36. BirdLife International. Tracking ocean wanderers: the global distribution of albatrosses and petrels. Cambridge, UK: BirdLife International; Results from the Global Procellariiform Tracking Workshop, 1–5 September, 2003, Gordon's Bay, South Africa; 2004.

37. Ashmole NP. Seabird ecology and the marine environment. In: Farner DS, King JR, editors. Avian biology, vol. 1. London: Academic Press; 1971. p. 223–86.

38. Ballance LT, Pitman RL, Reilly SB. Seabird community structure along a productivity gradient: importance of competition and energetic constraint. Ecology. 1997;78:1502–18.

39. Awkerman JA, Fukuda A, Higuchi H, Anderson DJ. Foraging activity and submesoscale habitat use of waved albatrosses *Phoebastria irrorata* during chick-brooding period. Mar Ecol Prog Ser. 2005;291:289–300.

40. Weimerskirch H, Lys P. Seasonal changes in the provisioning behaviour and mass of male and female wandering albatrosses in relation to the growth of their chick. Polar Biol. 2000;23:733–44.

41. Ydenberg RC, Welham CVJ, Schmid-Hempel R, Schmid-Hempel P, Beauchamp G. Time and energy constraints and the relationships between currencies in foraging theory. Behav Ecol. 1994;5:28–34.

42. Schmelzer I. Seals and seascapes: covariation in Hawaiian monk seal subpopulations and the oceanic landscape of the Hawaiian Archipelago. J Biogeogr. 2000;27:901–914.

43. Rasmussen K, Palacios DM, Calambokidis J, Saborio MT, Dalla Rosa L, Secchi ER, et al. Southern Hemisphere humpback whales wintering off Central America: insights from water temperature into the longest mammalian migration. Biol Lett. 2007;3:302–5.

44. Weng KC, Foley DG, Ganong JE, Perle C, Shillinger GL, Block BA. Migration of an upper trophic level predator, the salmon shark *Lamna ditropis*, between distant ecoregions. Mar Ecol Prog Ser. 2008;372:253–64.

45. Roden GL. Subarctic-subtropical transition zone of the North Pacific: large-scale aspects and mesoscale structure. In: Wetherall JA, editor. Biology, oceanography, and fisheries of the North Pacific transition zone and subarctic frontal zone. NOAA Technical Report 105. 1991. p. 1–38.

46. Polovina JJ, Howell E, Kobayashi DR, Seki MP. The transition zone chlorophyll front, a dynamic global feature defining migration and forage habitat for marine resources. Prog Oceanogr. 2001;49:469–83.

47. Laurs RM, Fiedler PC, Montgomery DR. Albacore tuna catch distributions relative to environmental features observed from satellites. Deep-Sea Res I. 1984;31:1085–99.

48. Olson DB, Hitchcock GL, Mariano AJ, Ashjian CJ, Peng G, Nero RW, et al. Life on the edge: marine life and fronts. Oceanography. 1994;7:52–9.

49. Polovina JJ, Kobayashi DR, Parker DM, Seki MP, Balazs GH. Turtles on the edge: movement of loggerhead turtles (*Caretta caretta*) along oceanic fronts, spanning longline fishing grounds in the central North Pacific, 1997–1998. Fish Oceanogr. 2000;9:71–82.

50. Bograd SJ, Foley DG, Schwing FB, Wilson C, Laurs RM, Polovina JJ, et al. On the seasonal and interannual migrations of the transition zone chlorophyll front. Geophys Res Lett. 2004;31:L17204.

51. Polovina JJ, Howell E, Abecassis M. Ocean's least productive waters are expanding. Geophys Res Lett. 2008;35:L03618.

52. Sydeman WJ, Thompson SA, Field JC, Peterson WT, Tanasichuk RW, Freeland HJ, et al. Does positioning of the North Pacific Current affect downstream ecosystem productivity? Geophys Res Lett. 2011;38:L12606.

53. Costa DP, Robinson PW, Arnould JPY, Harrison A-L, Simmons SE, Hassrick JL, et al. Accuracy of ARGOS locations of pinnipeds at-sea estimated using Fastloc GPS. PLoS One. 2010;5:e8677.

54. Shaffer SA, Tremblay Y, Awkerman JA, Henry RW, Teo SLH, Anderson DJ, et al. Comparison of light- and SST-based geolocation with satellite telemetry in free-ranging albatrosses. Mar Biol. 2005;147:833–43.

55. Fridolfsson A-K, Leggren H. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 1999;30:116–21.

56. Wilson RP, Weimerskirch H, Lys P. A device for measuring seabird activity at sea. J Avian Biol. 1995;26:172–5.

57. Phillips RA, Xavier JC, Croxall JP. Effects of satellite transmitters on albatrosses and petrels. Auk. 2003;120:1082–90.

58. McConnell BJ, Chambers C, Fedak MA. Foraging ecology of southern elephant seals in relation to the bathymetry and productivity of the Southern Ocean. Antarct Sci. 1992;4:393–8.

59. Suryan RM, Sato F, Balogh GR, Hyrenbach KD, Sievert PR, Ozaki K. Foraging destinations and marine habitat use of short-tailed albatrosses: A multi-scale approach using first-passage time analysis. Deep-Sea Res II. 2006;53:370–86.

60. Spear LB, Ainley DG. Flight speed of seabirds in relation to wind speed and direction. Ibis. 1997;139:234–51.

61. Tremblay Y, Shaffer SA, Fowler SL, Kuhn CE, McDonald BJ, Weise MJ, et al. Interpolation of animal tracking data in a fluid environment. J Exp Biol. 2006;209:128–40.

62. Wilson RP, Ducamp JJ, Rees WG, Culik BM, Niekamp K. Estimation of location: global coverage using light intensity. In: Priede IG, Swift SM, editors. Wildlife telemetry: remote monitoring and tracking of animals. New York: Ellis Horwood, Inc.; 1992. p. 131–4.

63. Hill RD. Theory of geolocation by light levels. In: Le Boeuf BJ, Laws RM, editors. Elephant seals: population ecology, behavior, and physiology. Berkeley: University of California Press; 1994. p. 227–36.

64. Teo SLH, Boustany A, Blackwell S, Walli A, Weng KC, Block BA. Validation of geolocation estimates based on light level and sea surface temperature from electronic tags. Mar Ecol Prog Ser. 2004;283:81–98.

65. Rodionov SN. A sequential algorithm for testing climate regime shifts. Geophys Res Lett. 2004;31:L09204.

66. Pinheiro J, Bates D, DebRoy S, Sarkar D. R Core team. nlme: Linear and nonlinear mixed effects models. 2009.

67. Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biometr J. 2008;50:346–63.

68. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2010.

69. Berens P. CircStat: A MATLAB toolbox for circular statistics. J Stat Softw. 2009;31:1–21.

70. Worton BJ. Kernel methods for estimating the utilization distribution in home-range studies. Ecology. 1989;70:164–8.

71. Wood AG, Naef-Daenzer B, Prince PA, Croxall JP. Quantifying habitat use in satellite-tracked pelagic seabirds: application of kernel estimation to albatross locations. J Avian Biol. 2000;31:278–86.

72. Calenge C. The package "adehabitat" for the R software: A tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9.

73. Silverman BW. Density estimation for statistics and data analysis. London: Chapman & Hall; 1986.

74. Breed GA, Bowen WD, McMillan JJ, Leonard ML. Sexual segregation of seasonal foraging habitats in a non-migratory marine mammal. Proc R Soc Lond B. 2006;273:2319–26.

75. Weimerskirch H, Guionnet T, Martin J, Shaffer SA, Costa DP. Fast and fuel efficient? Optimal use of wind by flying albatrosses. Proc R Soc Lond B. 2000;267:1869–74.

76. Shaffer SA, Costa DP, Weimerskirch H. Behavioural factors affecting foraging effort of breeding wandering albatrosses. J Anim Ecol. 2001;70:864–74.

77. Lewin-Koh NJ, Bivand R. maptools: Tools for reading and handling spatial objects. 2010.

78. Richards JA, Jia X. Remote sensing digital image analysis. New York: Springer; 1999.

79. Qinke Y, Jupp D, Rui L, Wei L. Re-scaling lower resolution slope by histogram matching. In: Zhou Q, Lees B, Tang G, editors. Advances in digital terrain analysis. Berlin: Springer; 2008. p. 193–210.

80. Aarts G, MacKenzie M, McConnell B, Fedak M, Matthiopoulos J. Estimating space-use and habitat preference from wildlife telemetry data. Ecography. 2008;31:140–60.

81. Matthiopoulos J. The use of space by animals as a function of accessibility and preference. Ecol Model. 2003;159:239–68.

82. Wakefield ED, Phillips RA, Trathan PN, Arata J, Gales R, Huin N, et al. Habitat preference, accessibility and competition limit the global distribution of breeding black-browed albatrosses. Ecol Monogr. 2011;81:141–67.

83. Rosenberg DK, McKelvey KS. Estimation of habitat selection for central-place foraging animals. J Wildl Manage. 1999;63:1028–38.

84. Wood SN. Generalized additive models: an introduction with R. Boca Raton, Florida: Chapman & Hall; 2006.

85. Pinheiro JC, Bates DM. Mixed-effects models in S and S-Plus. New York: Springer Verlag; 2000.

86. Harrison CS, Hida TS, Seki MP. Hawaiian seabird feeding ecology. Wildlife Monogr. 1983;85:1–71.

87. Gould P, Ostrom P, Walker W. Trophic relationships of albatrosses associated with squid and large-mesh drift-net fisheries in the North Pacific Ocean. Can J Zool. 1997;75:549–62.
88. Gong Y, Kim S, An DH. Abundance of neon flying squid in relation to oceanographic conditions in the North Pacific. Int N Pac Fish Comm Bull. 1993;53:191–204.

89. Yatsu A, Shimada H, Murata M. Distributions of epipelagic fishes, squids, marine mammals, seabirds and sea turtles in the Central North Pacific. Int N Pac Fish Comm Bull. 1993;53:111–46.

90. Pearcy WG, Fisher JP, Anma G, Meguro T. Species associations of epipelagic nekton of the North Pacific Ocean, 1978–1993. Fish Oceanogr. 1996;5:1–20.

91. Ichii T, Mahapatra K, Sakai M, Inagake D, Okada Y. Differing body size between the autumn and the winter-spring cohorts of neon flying squid (*Ommastrephes bartramii*) related to the oceanographic regime in the North Pacific: a hypothesis. Fish Oceanogr. 2004;13:295–309.

92. Breder CM, Rosen DE. Modes of reproduction in fishes. New York: The American Museum of Natural History; 1966.

93. Bakun A. Fronts and eddies as key structures in the habitat of marine fish larvae: opportunity, adaptive response and competitive advantage. Sci Mar. 2006;70S2:105–22.

94. Johannes RE. Words of the Lagoon. Berkeley: University of California Press; 1981.

95. Oschlies A, Garçon V. Eddy-induced enhancement of primary production in a model of the north Atlantic Ocean. Nature. 1998;394:266–9.

96. Stammer D. On eddy characteristics, eddy transports, and mean flow properties. J Phys Oceanogr. 1998;28:727–39.

97. Waugh DW, Abraham ER, Bowen MM. Spatial variations of stirring in the surface ocean: a case study of the Tasman Sea. J Phys Oceanogr. 2006;36:526–42.

98. Pennycuick CJ. The flight of petrels and albatrosses (Procellariiformes), observed in South Georgia and its vicinity. Philos Trans R Soc Lond B. 1982;300:75–106.

99. Wolanski E, Hamner WM. Topographically controlled fronts in the ocean and their biological influence. Science. 1988;241:177–81.

100. Springer AM, McRoy CP, Flint MV. The Bering Sea Green Belt: shelf-edge processes and ecosystem production. Fish Oceanogr. 1996;5:205–23.

101. Powell BS, Arango HG, Moore AM, Di Lorenzo E, Milliff RF, Foley D. 4DVAR data assimilation in the Intra-Americas Sea with the Regional Ocean Modeling System (ROMS). Ocean Model. 2008;23:130–45.

102. Ullman DS, Cornillon PC, Shan Z. On the characteristics of subtropical fronts in the North Atlantic. J Geophys Res. 2007;112:C01010.

103. Behrenfeld MJ, Falkowski PG. Photosynthetic rates derived from satellite-based chlorophyll concentration. Limnol Oceanogr. 1997;42:1–20.

104. Sokal RR, Rohlf FJ. Biometry: the principles and practice of statistics in biological research. New York: W.H. Freeman and Co.; 1995.

105. Quinn GP, Keough MJ. Experimental design and data analysis for biologists. Cambridge: Cambridge University Press; 2002.

106. Wood SN. gamm4: Generalized additive mixed models using mgcv and lme4. 2010.

107. Crawley MJ. Statistical computing: an introduction to data analysis using S-Plus. West Sussex, England: John Wiley & Sons; 2002.

108. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R. New York: Springer; 2009.

109. Olivier F, Wotherspoon SJ. GIS-based application of resource selection functions to the prediction of snow petrel distribution and abundance in East Antarctica: Comparing models at multiple scales. Ecol Model. 2005;189:105–29.

110. Barclay RMR. The effect of reproductive condition on the foraging behavior of female hoary bats, Lasiurus cinereus. Behav Ecol Sociobiol. 1989;24:31–7.

111. Clarke J, Emmerson LM, Otahal P. Environmental conditions and life history constraints determine foraging range in breeding Adélie penguins. Mar Ecol Prog Ser. 2006;310:247–61.

112. Larimer SC, Fritzsche P, Song Z, Johnston J, Neumann K, Gattermann R, et al. Foraging behavior of golden hamsters (*Mesocricetus auratus*) in the wild. J Ethol. 2011;29:275–83.

113. Rice DW, Kenyon KW. Breeding cycles and behavior of Laysan and black-footed albatrosses. Auk. 1962;79:517–67.

114. Legaard KR, Thomas AC. Spatial patterns in seasonal and interannual variability of chlorophyll and sea surface temperature in the California Current. J Geophys Res. 2006;111:C06032.

115. Henson SA, Thomas AC. Interannual variability in timing of bloom initiation in the California Current System. J Geophys Res. 2007;112:C08007.

116. Fischer KN, Suryan RM, Roby DD, Balogh GR. Post-breeding season distribution of black-footed and Laysan albatrosses satellite-tagged in

Alaska: inter-specific differences in spatial overlap with North Pacific fisheries. Biol Conserv. 2009;142:751–60.

117. Glass GE, Slade NA. The effect of *Sigmodon hispidus* on spatial and temporal activity of *Microtus ochrogaster*: evidence for competition. Ecology. 1980;61:358–70.

118. Bocher P, Cherel Y, Hobson KA. Complete trophic segregation between South Georgian and common diving petrels during breeding at Iles Kerguelen. Mar Ecol Prog Ser. 2000;208:249–64.

119. González-Solís J, Croxall JP, Afanasyev V. Offshore spatial segregation in giant petrels *Macronectes spp.* differences between species, sexes and seasons. Aquat Conserv. 2008;17:522–36.

120. Navarro J, Forero MG, González-Solís J, Igual JM, Bécares J, Hobson KA. Foraging segregation between two closely related shearwaters breeding in sympatry. Biol Lett. 2009;5:545–8.

121. Xavier JC, Louzao M, Thorpe SE, Ward P, Hill C, Roberts D, et al. Seasonal changes in the diet and feeding behaviour of a top predator indicate a flexible response to deteriorating oceanographic conditions. Mar Biol. 2013;160:1597–606.

122. Pitman RL, Walker WA, Everett WT, Gallo-Reynoso JP. Population status, foods and foraging of Laysan albatrosses *Phoebastria immutabilis* nesting on Guadalupe Island, Mexico. Mar Ornithol. 2004;32:159–65.

123. Weimerskirch H, Gault A, Cherel Y. Prey distribution and patchiness: factors in foraging success and efficiency of wandering albatrosses. Ecology. 2005;86:2611–22.

124. Fernández P, Anderson DJ. Nocturnal and diurnal foraging activity of Hawaiian albatrosses detected with a new immersion monitor. Condor. 2000;102:577–84.

125. Hyrenbach KD, Dotson RC. Assessing the susceptibility of female black-footed albatross (*Phoebastria nigripes*) to longline fisheries during their post-breeding dispersal: an integrated approach. Biol Conserv. 2003;112:391–404.

126. Phillips RA, Silk JRD, Phalan B, Catry P, Croxall JP. Seasonal sexual segregation in two *Thalassarche albatross* species: competitive exclusion, reproductive role specialization or foraging niche divergence? Proc R Soc Lond B. 2004;271:1283–91.

127. Weimerskirch H, Guionnet T. Comparative activity pattern during foraging of four albatross species. Ibis. 2002;144:40–50.

128. Phalan B, Phillips RA, Silk JRD, Afanasyev V, Fukuda A, Fox J, et al. Foraging behaviour of four albatross species by night and day. Mar Ecol Prog Ser. 2007;340:271–86.

129. Pettit TN, Nagy KA, Ellis HI, Whittow GC. Incubation energetics of the Laysan albatross. Oecologia. 1988;74:546–50.

130. Hyrenbach KD, Keiper C, Allen SG, Ainley DG, Anderson DJ. Use of marine sanctuaries by far-ranging predators: commuting flights to the California Current System by breeding Hawaiian albatrosses. Fish Oceanogr. 2006;15:95–103.

131. Fauchald P, Tveraa T. Using first-passage time in the analysis of area-restricted search and habitat selection. Ecology. 2003;84:282–8.

132. Kareiva P, Odell G. Swarms of predators exhibit "preytaxis" if individual predators use area-restricted search. Am Nat. 1987;130:233–70.

133. Felicísimo AM, Muñoz J, González-Solis J. Ocean surface winds drive dynamics of transoceanic aerial movements. PLoS One. 2008;3:e2928.

134. Wakefield ED, Phillips RA, Matthiopoulos J, Fukuda A, Higuchi H, Marshall GJ, et al. Wind field and sex constrain the flight speeds of central-place foraging albatrosses. Ecol Monogr. 2009;79:663–79.

135. Adams J, Flora S. Correlating seabird movements with ocean winds: linking satellite telemetry with ocean scatterometry. Mar Biol. 2010;157:915–29.

136. Raymond B, Shaffer SA, Sokolov S, Woehler EJ, Costa DP, Einoder L, et al. Shearwater foraging in the Southern Ocean: the roles of prey availability and winds. PLoS One. 2010;5:e10960.

137. Suryan RM, Anderson DJ, Shaffer SA, Roby DD, Tremblay Y, Costa DP, et al. Wind, waves, and wing loading: morphological specialization may limit range expansion of endangered albatrosses. PLoS One. 2008;3:e4016.

138. Pascual A, Faugère Y, Larnicol G, Le Traon P-Y. Improved description of the ocean mesoscale variability by combining four satellite altimeters. Geophys Res Lett. 2006;33:L02611.

139. Benson SR, Eguchi T, Foley DG, Forney KA, Bailey H, Hitipeuw C, et al. Large-scale movements and high-use areas of western Pacific leatherback turtles, *Dermochelys coriacea*. Ecosphere. 2011;2:art84.

140. Longhurst AR. Ecological geography of the sea. San Diego, California: Academic Press; 1998.
141. Hardman-Mountford NJ, Hirata T, Richardson KA, Aiken J. An objective methodology for the classification of ecological pattern into biomes and provinces for the pelagic ocean. Remote Sens Environ. 2008;112:3341-52.

142. Harrison CS. Seabirds of Hawaii. Ithaca, New York: Cornell University Press; 1990.

143. Fisher HL. Body weights in Laysan albatrosses *Diomedea immutabilis*. Ibis. 1967;109:373-82.

144. Prince PA, Ricketts C, Thomas G. Weight loss in incubating albatrosses and its implications for their energy and food requirements. Condor. 1981;83:238-42.

145. Weimerskirch H. Regulation of foraging trips and incubation routine in male and female wandering albatrosses. Oecologia. 1995;102:37-43.

146. Chaurand T, Weimerskirch H. The regular alternation of short and long foraging trips in the blue petrel *Halobaena caerulea*: a previously undescribed strategy of food provisioning in a pelagic seabird. J Anim Ecol. 1994;63:275-82.

147. Weimerskirch H, Chastel O, Ackermann L, Chaurand T, Cuenot-Chaillet F, Hindermeyer X, et al. Alternate long and short foraging trips in pelagic seabird parents. Anim Behav. 1994;47:472-6.

148. Baduini CL, Hyrenbach KD. Biogeography of Procellariform foraging strategies: does ocean productivity influence provisioning? Mar Ornithol. 2003;31:101-12.

149. Sprague RS. Glucocorticoid physiology and behavior during life history transitions in Laysan albatross (*Phoebastria immutabilis*). PhD Thesis, University of Montana; 2009.

150. Fisher HL. Mortality and survival in the Laysan albatross, *Diomedea immutabilis*. Pac Sci. 1975;29:279-300.

151. Naughton MB, Romano MD, Zimmerman TS. A Conservation Action Plan for black-footed albatross (*Phoebastria nigripes*) and Laysan albatross (*P. immutabilis*). 2007.

152. Dearborn DC, Anders AD, Flint EN. Trends in reproductive success of Hawaiian seabirds: is guild membership a good criterion for choosing indicator species? Biol Conserv. 2001;101:97-103.

153. Seki MP, Flint EN, Howell E, Ichii T, Polovina JJ, Yatsu A. Transition Zone. In: Alexander V, Perry IR, editors. PICES report on marine ecosystems of the North Pacific. North Saanich. British Columbia: North Pacific Marine Science Organization; 2004. p. 191-8.

154. Levitus S, Antonov J, Boyer T. Warming of the world ocean, 1955-2003. Geophys Res Lett. 2005;32:L02604.

155. Climate Change IPCC. Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press; 2007. p. 2007.

156. Yin JH. A consistent poleward shift of the storm tracks in simulations of 21st century climate. Geophys Res Lett. 2005;32:L18701.

157. Bengtsson L, Hodges KI, Roeckner E. Storm tracks and climate change. J Clim. 2006;19:3518-43.

158. Archer CL, Caldeira K. Historical trends in the jet streams. Geophys Res Lett. 2008;35:L08803.

159. Weimerskirch H, Louzao M, de Grissac S, Delord K. Changes in wind pattern alter albatross distribution and life-history traits. Science. 2012;335:211-4.

160. Young LC, Vanderlip C, Duffy DC, Afanasyev V, Shaffer SA. Bringing home the trash: do colony-based differences in foraging distribution lead to increased plastic ingestion in Laysan albatrosses? PLoS One. 2009;4:e7623.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central