Introduction

Describing seasonal and migratory movements in many large marine vertebrates is challenging, largely due to the complexities of tracking individuals in water for sufficient time to obtain meaningful insight into migratory behaviour (Hammerschlag et al. 2011). Advances in satellite tracking technologies and attachment techniques now allow for observations of movements and space use over extended timescales (Hazen et al. 2012a; Hussey et al. 2015). These advance our ability to describe life-history events (Block et al. 2011; Hussey et al. 2015) and reveal an extraordinary array of behaviours from ocean basin migrations (Bonfil et al. 2005), to remarkable individual dive profiles into depths greater than 2500 m (Schorr et al. 2014).

Knowledge of vertical movements can provide insight into how animals use their environment (Musyl et al. 2011). Some parameters such as light, pressure, and temperature vary with depth, which can partition species ranges due to morphological and physiological constraints (Hussey et al. 2015). Animal movements between habitats may affect trophic dynamics such as nutrient transport between depth zones, which can represent a significant contribution to ecosystem nutrient cycling (Frid et al. 2009; Schmitz et al. 2010; Matich et al. 2011; Thorrold et al. 2014). For heavily exploited species, and species of conservation concern, coupling knowledge on vertical with horizontal movement behaviour is essential to fully understand the potential for overlap with specific fishing gear types and activity. With sharks in particular, reducing susceptibility to fisheries as bycatch is a major priority (Barker and Schluessel 2005; Hetherington et al. 2015; Oliver et al. 2015; Dulvy et al. 2017).

Historically, basking sharks were heavily exploited by directed fisheries in the north-east Atlantic, resulting in marked declines in local populations (Kunzlik 1988). This led to the IUCN listing the north-east Atlantic stocks as Endangered, with global populations assessed as Vulnerable (Fowler 2005). Incidental catch can still occur—with sharks caught annually in the bottom gill-net fishery in the Celtic Sea (Berrow 1994), and in trawl and set net fisheries in New Zealand (Francis and Duffy 2002). Basking sharks are protected under national (such as the Wildlife and Countryside Act (1981)) and international treaties, including being listed in Appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), and Appendices I and II in the Convention of Migratory Species (CMS) requiring international trade to be monitored.

Previous studies in the north-east Atlantic on vertical space use of basking sharks have shown that the majority of time is spent occupying shallow surface waters during the summer (Sims et al. 2003; Stéphan et al. 2011) moving into deeper, mesopelagic waters as winter approaches (Sims et al. 2003; Gore et al. 2008; Stéphan et al. 2011; Doherty et al. 2017a). Many of these studies, however, have been limited in tracking duration to see fine-scale vertical space use throughout winter. These patterns have also been described in the north-west Atlantic (Skomal et al. 2009; Braun et al. 2018). Sims et al. (2003) suggested that vertical movements were consistent with those associated with foraging, with animals likely feeding year round on zooplankton. These behaviours, however, differed with habitat type: when in deep stratified waters of the continental shelf-edge, basking sharks exhibited normal diel vertical migration (DVM; dusk ascent-dawn descent), but when in shallow, inner-shelf waters basking sharks conducted reverse diel vertical migration (rDVM; dusk descent-dawn ascent; Sims et al. 2005; Stéphan et al. 2011; Queiroz et al. 2017). Basking sharks were also shown to switch behaviours to a tidal rhythm when encountering boundaries between thermally stratified and mixed waters (Shepard et al. 2006).

Analysing fine-scale depth data can help identify patterns in vertical space use and gain insight into the underlying factors that influence habitat selection and behaviour (Vaudo et al. 2014; Afonso and Hazin 2015; Queiroz et al. 2017), particularly for species of conservation concern, such as basking sharks. Using data from archival satellite tags deployed on basking sharks, our aims were to (1) quantify basking shark depth and temperature distributions and vertical space use patterns, (2) investigate whether basking shark vertical space use varies throughout the year, and (3) investigate potential factors that might influence vertical behaviour.

Methods

Tag attachment and specification

Satellite tags were attached to basking sharks from a boat, using a titanium M-style dart (Wildlife Computers, Redmond, California, USA) inserted at the base of the first dorsal fin with a modified pole spear. Tags were attached via a tether consisting of heat-shrink covered stainless steel flexible cable, a swivel and monofilament line attached to the tag. Sharks were sexed when possible using a pole-mounted camera, and total body length estimated against the length of the survey vessel (10 m).

Thirty-two archival satellite tags (Mini Pop-up Archival (MiniPAT, n = 10), and Pop-up Archival Transmitting with Fastloc™ GPS tags (PAT-F, n = 9; SPLASH-F, n = 13; Wildlife Computers, Washington, USA) were attached to basking sharks off the west coast of Scotland during July and August of 2012–2014. These tags were programmed to summarise depth data at 4 h intervals across twelve depth ranges of; 0–1 m, 1–5 m, 5–10 m, 10–25 m, 25–50 m, 50–75 m, 75–100 m, 100–250 m, 250–500 m, 500–750 m, 750–1000 m, and > 1000 m, along with 12 temperature ranges of 0 °C, 0–4 °C, 4–6 °C, 6–8 °C, 8–10 °C, 10–12 °C, 12–14 °C, 14–16 °C, 16–18 °C, 18–20 °C, 20–22 °C, and > 22 °C. Maximum reported daily depths of sharks were used to estimate vertical position within the water column during wide-ranging movements.

Location data processing

To observe movements into winter months, light geolocation data were obtained from archival tags with attachment durations into at least the January following tag attachment (> 165 days; n = 12; Doherty et al. 2017a). Hidden Markov Models (HMM; Neilson et al. 2014) were applied using light level, temperature, and depth data recorded by the tag to estimate the most likely track of individuals from high quality light measurements (Thygesen et al. 2009; by Collecte Localisation Satellites (CLS-Argos; www.argos-system.org)). For a detailed description of specific parameters of HMMs applied, see (Biais et al. 2017; Doherty et al. 2017a). In brief, to reconstruct the most probable track of basking sharks, HMMs were used to estimate the daily probability density (or Utilisation Distribution; UD) of the location of tracked animals (Thygesen et al. 2009). The HMM used a two-step process: first, at each sampling time a position prediction was made and, second, the predicted probability density was combined with latitude, longitude, SST, and depth recorded onboard the tag to produce the posterior distribution of the individual (Biais et al. 2017). These data were constrained by bathymetry (etopo2; https://www.ngdc.noaa.gov/mgg/global/etopo2.html), SST (GHRSST-OSTIA; https://www.ghrsst.org/), and known deployment and pop-off locations with daily location estimates were taken as the geographic mean of the grid locations weighted by their probability (Thygesen et al. 2009).

To examine variability in depth occupancy with respect to geographic location, polygon sampling grids [each cell was 50 km from grid cell centroid to perimeter; (cell area 8660 km2) encompassing mean error across all light-based geolocation estimates (98 km)] were created, bounded by the maximum limits of observed movement. These grids were spatially intersected with light geolocated tracking locations. The mean daily minimum and mean daily maximum depth within each grid cell were calculated for each shark and then averaged (mean) across all individuals calculated. Times of night and day were determined by applying sunrise and sunset times for mean latitude and longitude of occupancy within each month obtained using the suncalc package (Agafonkin and Thieurmel 2018) in R (R Core Team 2018).

A General Additive Model (GAM) was performed to investigate the relationship between the maximum daily depth occupied by basking sharks, latitude the depth was occupied at (fitting a cyclic cubic regression spline for each month), or month in which the depth occupancy occurred. Months included in the model were restricted to months outside of peak times of high occupancy of surface waters, retaining September to April. The response variable of daily maximum depth was loge transformed to fulfil assumptions of normality of data with individual shark identification number and year of deployment fitted as separate random effects. Models were fitted using the mgcv package (Wood 2017) and ranked by Akaike’s Information Criteria (AIC) using subset selection of the maximal model using the MuMIn package (Barton 2018). All terms within the maximal model were examined fixing month in the model to minimise bias of estimates from centred smooth terms. Top ranked models were defined as models ∆AIC ≤ 6 units of the best supported model (Harrison et al. 2018).

Sea surface temperature (SST) values were obtained from monthly averaged MODIS L3 nighttime remotely sensed rasters (°C; 4 km resolution, http://podaac.jpl.nasa.gov) prepared using Marine Geospatial Ecological Tools v0.8a64 (Roberts et al. 2010). Temporally specific locations were used to extract SST values and averaged for each month across years for the study period.

Data analyses were performed in R (R Core Team 2018) and all spatial analyses and maps created using Geospatial Modelling Environment (GME v 0.7.2.1; Beyer 2012) and ESRI ArcMap 10.1.

Results

General vertical space use

Thirty-two archival satellite tags (Table S1. PAT-F; n = 9, MiniPAT; n = 10, SPLASH-F; n = 13) were attached to basking sharks in 2012 (n = 9), 2013 (n = 14), and 2014 (n = 9). Satellite tags were attached to 12 females, 6 males and 14 individuals of unknown sex, measuring 4–5 m (n = 6), 5–6 m (n = 10), 6–7 m (n = 7), 7–8 m (n = 8), and 8–9 m (n = 1) estimated total length. These tags transmitted depth and temperature data over a cumulative 4489 days (mean 140 ± 97, range 10–292 days; Table S1). Sharks spent on average more time in shallower waters in the peak of summer [July; most occupied depth class = 0–1 m (21% of time), August; 25–50 m (25%)], particularly during the daytime [July; most occupied depth class in the day = 0–1 m (25% of time), compared to 10–25 m at night (23% of time), August; 0–1 m in the day (27%), compared to 25–50 m at night (28%)]. In the peak of the winter, sharks were generally in waters away from the surface [December; most occupied depth class = 50–75 m (28% of time), January; 75–100 m (19%)], but were shallower during the night-time (December; most occupied depth class in the day = 50–75 m (28% of time), compared to 25–50 m at night (28% of time), January; 100–250 m in the day (18%), compared to 50–75 m at night (24%); Tables S2 and S3; Figs. 1 and S1).

Fig. 1
figure 1

Monthly time-at-depth occupancy of basking sharks by time of day. Mean time-at-depth for basking sharks (n = 32) for each month of tag attachment. Bars represent mean percentage of time (summarised in 4-h intervals) occupying eleven depth ranges ± standard error (se) for day (yellow bars) and night (grey bars). Number of sharks contributing to each monthly plot shown. Dotted lines denote labelled ocean zones

Sharks exhibited evidence of spatially explicit depth occupancy in the north-east Atlantic. During the peak summer months (July–August), sharks occupied shallow waters off the west coast of Scotland and north Northern Ireland. These sharks moved further south along the west coast of Ireland and into the Irish and Celtic Seas during autumn (September–October) exhibiting maximum depth occupancy greater than 100 m. During the early winter (November–December), sharks began to move into waters off the continental shelf, west of the Bay of Biscay remaining in epipelagic waters. Late winter and early spring (January–March) showed the emergence of use of waters greater than 750 m, most often occurring west of the Bay of Biscay towards the Iberian Peninsula. During spring (April–May), sharks exhibited a contraction of spatial extent, within waters closer to the coasts of Ireland and the UK, showing increased occurrence of maximum depths in waters shallower than 250 m (Fig. 2). Occupancy of deeper waters appeared not to be governed by bathymetry as maximum depths were not at the seafloor and often remained within the epipelagic zone when in higher latitudes (Fig. 2).

Fig. 2
figure 2

Grid density enumeration of monthly minimum and maximum depths for satellite tracked basking sharks. Mean daily minimum and mean daily maximum depth of basking sharks on a hexagonal grid (cell edge size: 50 km; cell area: 8660 km2). Locations derived from modelled light geolocation positions from archival tags for basking sharks (n = 12). Country international two letter codes shown in regional map (UK United Kingdom, IE Ireland, FR France, ES Spain, PT Portugal, and MA Morocco). Broken black line denotes 200 m depth contour. Maps projected in Europe Albers Equal Area Conic coordinates system. Number of sharks contributing to each monthly plot shown

Individual variation in vertical space use

Analysis of vertical space use of basking sharks with extended tracking durations (> 165 days; n = 12) revealed occupancy of deep waters (> 500 m) during late winter–early spring (Fig. 3). Seven of the twelve sharks also used waters greater than 1000 m, with two sharks reaching depths of 1500 m (Fig. 3; Table S1).

Fig. 3
figure 3

Basking shark depth-use time-series data. Daily minimum (grey line) and daily maximum (black line) depths from modelled light geolocation positions from archival tags of 12 basking sharks. Latitude of greatest maximum daily depth is shown for sharks exhibiting vertical space use greater than 1000 m. Tag number, sex (F female, M male, and U unknown), and body length (m) shown, dotted line denotes 200 m depth

Five of the 12 archival tags with attachment durations greater than 165 days were physically recovered, allowing for high-resolution time-series (depth and temperature data recorded at 10–15 s intervals) profiles to be obtained (Table S1). Data indicated a seasonal shift in vertical space use from summer to autumn/winter where sharks moved into deeper waters, spent less time associated with surface waters, and experienced smaller fluctuations in water temperature (Fig. 4). These time-series data also revealed a conspicuous switch in vertical movements from a relatively uniform use of 50–250 m depths during the winter, to deep, rapid and repeated ‘bounce’ or ‘yo–yo’ dive behaviour from depth to surface waters in late winter–early spring. These movements were often carried out whilst making more directed movements off the continental shelf (Fig. 4).

Fig. 4
figure 4

Basking shark depth-use from high-resolution time-series data from physically recovered satellite tags. a PAT-F tag depth and temperature recorded at 10-s intervals, b SPLASH-F tag depth and temperature recorded at 15-s intervals, c PAT-F tag depth and temperature recorded at 10-s intervals, and d, e MiniPAT tag depths and temperatures recorded at 15-s intervals. Left panels: Maps showing movements from archival satellite tags for five basking sharks tracked for over 165 days, displaying modelled light geolocation positions (grey circles) with associated light geolocation error (light grey ellipses), and track end points (white stars). Broken black line denotes 200 m depth contour. Right panels: Depth (black lines) and temperature (grey lines) time-series data of the same individuals. Blue shaded sections in the left panels correspond to locations highlighted in blue circles in right panel, as regions of deep diving behaviour

Maximum daily depths occupied by basking sharks were shown to be influenced by the latitude in which the depth occupancy occurred and varied throughout the year. There was a trend of occupancy of deeper water occurring at lower latitudes (south of 50° N) with January to April showing evidence of movement into deeper waters (> 500 m; Table S2; Fig. 5).

Fig. 5
figure 5

Spatio-temporal depth-use predictions for basking sharks from a general additive model (GAM). Relationship between latitude and maximum daily depth occupied by basking sharks tracked for extended durations (> 165 days; n = 12) for months September to April. Solid black lines denote predictions from general additive model (GAM) with blue polygon representing standard error (se). Month and number of sharks contributing to the plot shown. Rug plots show the latitude at which data are present for each month

High-resolution time-series data from tags demonstrated basking sharks occupied waters of a mean temperature of 11.6 °C (± 0.24 °C; range 5.5–19.1 °C); however, a mean of 99% (± 2.8%; 88–100%) of time was spent in waters of between 8 and 16 °C, with water temperatures experienced by these sharks decreasing with the onset of winter (Fig. 4; Fig. S2). Sea surface temperatures in the same locations that sharks were tracked showed similar seasonal changes and a comparable temperature range to that experienced by basking sharks at depth (mean 11.8 °C ± 1.7; range 10–14 °C).

Discussion

The need to identify areas providing suitable conditions for large marine vertebrates to undertake key life-history events, such as foraging, mating, and parturition is critical to understanding the biology and ecology of species of conservation concern. Most marine organisms function within three-dimensional environments and a number of important factors including; light, pressure, temperature, and oxygen and salinity concentrations vary with depth (Hussey et al. 2015). These physical features of the ocean can act as barriers to species’ movement due to adaptations or limitations of physiology and morphology; however, large marine vertebrates can often occupy many distinct vertical habitats (Boustany et al. 2002; Thorrold et al. 2014).

The data demonstrate that basking sharks exhibit seasonal changes in vertical movement in the north-east Atlantic, with shallower vertical space use in summer compared to winter months, whereby sharks occupy epipelagic (0–200 m) to mesopelagic depths (200–1000 m), with some forays into the bathypelagic zone (1000–4000 m). Basking sharks are capable of finding dense prey patches (Sims et al. 2005, 2006), and activity at the surface during the summer months is likely in relation to these food sources, especially during the day in nearshore waters, with data from the continuous plankton recorder showing increased biomass of copepods in these regions (Richardson et al. 2006). Diel vertical migration (DVM) has been described for all three species of planktivorous sharks: megamouth shark [Megachasma pelagios (Nelson et al. 1997), whale sharks (Rhincodon typus; Graham et al. 2006, Wilson et al. 2006), and basking sharks (Sims et al. 2005; Stéphan et al. 2011)]. Reverse diel vertical migration (rDVM) can also occur, with movements switching towards the surface during the day. Both DVM and rDVM have been observed in whale sharks (Rowat et al. 2006) and basking sharks (Shepard et al. 2006; Stéphan et al. 2011) and may be determined by the heterogeneous environmental conditions of the water column they inhabit. Basking sharks in this study exhibited rDVM during the summer months, moving to relatively shallow waters during daylight hours when they spent the majority of time within the Sea of the Hebrides (Doherty et al. 2017b). During peak winter months, basking sharks switched to exhibiting DVM, likely due to the sharks moving from shallow, well-mixed areas with plankton blooms in summer, which breakdown at the onset of winter to deeper, more stratified, waters (Sims et al. 2005; Stéphan et al. 2011) sometimes at mesopelagic depths (Gore et al. 2008; Skomal et al. 2009; Doherty et al. 2017a). Copepods dominate zooplankton biomass in the North Atlantic (Gallienne et al. 2001), with basking shark abundance shown to positively correlate with increased abundance of copepods (Sims and Merrett 1997; Sims and Quayle 1998; Sims 1999). The annual life cycle of copepods often includes a period of active feeding and growth in surface waters during spring and summer and a period of diapause in deep water (Hirche 1998). Diapause is characterised by a reduction in development, feeding, and metabolic rates, and is associated with migration to deep water in summer and autumn and ascension towards surface waters in late winter (Hirche 1998). Seasonal cycles also show evidence of predator avoidance by copepods during summer months by reversing migration (rDVM) in relation to predator vertical position during the day (Irigoien et al. 2004). Previous studies have suggested the optimal emergence from diapause in the North Atlantic occurs when copepods overwinter at temperatures less than 5 °C (Saumweber and Durbin 2006). However, some species undergoing diapause tolerate a broad range of temperatures throughout the North Atlantic (−1 to 11 °C), with dormancy prevailing in higher latitudes (Krumhansl et al. 2018). It is possible that basking sharks have the ability to track the patterns in vertical movement exhibited by prey species by also adopting rDVM movement strategy in summer months when copepods are in surface waters in the daytime during the foraging stage of the annual cycle or avoiding predators. During the winter, it is unclear if basking sharks continue to actively feed, but may switch to DVM movements at greater depths to track winter prey daily cycles and/or vertical occupancy during diapause.

Basking sharks in the present study were recorded making repeated ‘yo–yo’ vertical movement between the surface and deeper waters in late winter–spring (February–April). This has been recorded in a wide range of other shark species including whale sharks (Brunnschweiler et al. 2009), white sharks [Carcharodon carcharias, (Klimley et al. 2002; Domeier and Nasby-Lucas 2008)], scalloped hammerhead sharks (Sphyrna lewini, (Jorgensen et al. 2009)), and tiger sharks (Galeocerdo cuvier (Nakamura et al. 2011)). These behaviours are generally attributed to feeding—when sharks begin to actively seek out foraging areas (horizontally or vertically); however, it is also possible that they are involved in thermoregulation or aid energy conservation (Holland et al. 1992; Klimley et al. 2002; Thums et al. 2013). This type of vertical movement may also serve as a means of orientation, as it has been suggested that sharks use cognitive maps of their home ranges at various spatial scales (Meyer et al. 2010; Papastamatiou et al. 2011).

The spatially explicit data in the present study show evidence of some individuals displaying directed movements off the continental shelf before exhibiting ‘yo–yo’ vertical behaviour between 6 and 13 days, before returning to continental shelf waters. These movements may be examples of foraging trips, failed foraging trips, or “sampling” of the water column for orientation and/or olfactory cues. Queiroz et al. (2017) showed that basking sharks exhibit “V”-shaped time-depth profiles (short durations at maximum depth prior to ascent) to a wide range of depths, further suggesting that these movements serve as an exploratory function. These authors also showed that when approaching an area of increased chlorophyll, maximum depth of V-shaped dives decreased as well as the variation in targeted depths. These repeated movements may represent the beginning of increased activity from a more torpid state during winter months if not actively feeding, or in anticipation of the onset of spring plankton blooms. Olfactory prey cues are considered to be more extensively distributed horizontally than vertically, due to current shear between water layers of different densities (Carey and Scharold 1990; Klimley et al. 2002). Therefore, repeated deep diving behaviour may facilitate sampling of multiple depth layers for chemical cues over short time periods, a behaviour documented for pelagic predatory sharks (Carey and Scharold 1990; Boustany et al. 2002) and is likely a common search strategy for both epipelagic planktivores and macropredators (Sims et al. 2003; Queiroz et al. 2017).

During the present study, we observed several sharks travelling to depths up to 1500 m, most often off the continental shelf to the west of the Bay of Biscay at the Celtic-Armorica shelf margin. This is comparable to the deepest depth recorded for a basking shark of 1504 m (Braun et al. 2018), with previous observations of basking sharks occupying waters at 1264 m (Gore et al. 2008), and between 750 and 1000 m (Francis and Duffy 2002; Sims et al. 2003; Skomal et al. 2009). Other planktivorous elasmobranchs have also been shown to perform vertical space use beyond 1000 m, such as whale sharks diving to 1928 m (Tyminski et al. 2015) and Chilean devil rays (Mobula tarapacana) diving to 1896 m, suggesting the possibility of foraging at depths providing an important ecological link between surface foragers and forage species occupying pelagic habitats below the euphotic zone (Thorrold et al. 2014).

Basking sharks in this study consistently occupied waters between 8 and 16 °C regardless of depth or time of year. This appears to be a relatively narrow thermal niche compared to wide ranging teleost fish (Atlantic cod (Gadus morhua) range − 1.5 to 20 °C (Righton et al. 2010); Atlantic bluefin tuna (Thunnus thynnus) range 3–31 °C (Block et al. 2001)) and sharks (salmon sharks (Lamna ditropis) range 2–24 °C (Weng 2005); shortfin mako sharks (Isurus oxyrinchus) range 5–30 °C (Vaudo et al. 2016); white sharks range 5–26 °C (Boustany et al. 2002); whale sharks range 4–31 °C (Graham et al. 2006)). This suggests that the vertical space use exhibited by basking sharks in the present study could help to maintain a consistent thermal envelope, which could be advantageous as many physiological rates and functions (e.g. metabolic rate and reproduction) are determined by body temperature (Angilletta et al. 2002; Whitney et al. 2016). We observed basking sharks moving into deeper waters in the winter; however, our results also suggest that conditions may be thermally suitable for basking sharks to occupy surface waters year-round in the north-east Atlantic with sea surface temperature values showing similar seasonal shifts to temperature-at-depth. We observed sharks occupying depths of 1500 m; however, we cannot be certain of the maximum depths these sharks can tolerate, which may be governed by temperature. Other variables such as oxygen levels may impact the depths at which these sharks can occupy (Nasby-Lucas et al. 2009; Prince et al. 2010; Comfort and Weng 2015). Species metabolic rates and, therefore, oxygen requirements may determine environments in which they can inhabit (Bernal et al. 2012). These results may provide evidence for basking sharks feeding on overwintering zooplankton in deeper waters during the winter, although we have no direct observations of this. Incorporating new knowledge of vertical space use into tag programming should be carried out in future research, allowing for depth class assignment to reflect that of what is utilised by sharks. Reducing the breadth of the deeper depth classes may allow for further investigation of timings and proportion of time at more extremes of vertical space use for these individuals.

Conclusion

Characterising vertical space use in basking sharks within areas of high relative importance can contribute to discussions regarding management options for the conservation of basking sharks. Basking sharks exhibit seasonality in their vertical space use, likely reacting to changes in environmental conditions in the areas they inhabit, with occupancy of surprisingly deep waters during the winter. Combining detailed understanding of space use in three dimensions (latitude, longitude, depth, and time) in response to seasonal biological and environmental conditions can elucidate the ecosystem role and environmental preferences of these sharks (Coffey et al. 2017). This spatio-temporally explicit depth information can strengthen our ability to estimate potential interactions with fisheries such as mid-water trawls in order to mitigate bycatch risk. Understanding patterns of habitat use and environmental preferences for species of conservation concern will best inform on implementation of sustainable conservation and management strategies, especially when considering projections of range shifts under future climate scenarios (Hazen et al. 2012b; Miller et al. 2015).