Leatherback turtles in the Northwest Atlantic Ocean have a broad geographic range that extends from nesting beaches near the equator to seasonal foraging grounds as far north as Canada. The ability of leatherbacks to maintain core body temperature (Tb) higher than that of the surrounding water is thought to be a key element of their biology that permits them to exploit productive waters at high latitudes. We provide the first recordings of Tb from freely swimming leatherbacks at a northern foraging ground, and use these data to assess the importance of behavioral adjustments and metabolic sources of heat for maintenance of the thermal gradient (Tg). The mean Tb for individual leatherbacks ranged from 25.4±1.7 to 27.3±0.3°C, and Tg ranged from 10.7±2.4 to 12.1±1.7°C. Variation in mean Tb was best explained by the amount of time that turtles spent in the relatively warm surface waters. A diel trend in Tb was apparent, with daytime cooling suggestive of prey ingestion and night-time warming attributable to endogenous heat production. We estimate that metabolic rates necessary to support the observed Tg are ~3 times higher than resting metabolic rate, and that specific dynamic action is an important source of heat for foraging leatherbacks.
Leatherback turtles, Dermochelys coriacea (Vandelli 1761), forage on gelatinous zooplankton in temperate waters of the North Atlantic Ocean during the summer and autumn (James et al., 2006; Heaslip et al., 2012). The ability of leatherbacks to tolerate cold water temperatures experienced at high latitude is attributed to morphological features and physiological and behavioral adjustments that permit maintenance of a thermal gradient (Tg) between core body temperature (Tb) and ambient water temperature (Ta). Point measurements of Tb obtained from leatherbacks captured off Cape Breton, NS, Canada, indicate a Tg of 5.1–10.8°C is achievable in waters with an average sea surface temperature of 16.2°C (James and Mrosovsky, 2004).
Adult leatherbacks typically weigh 250–600 kg, and a small surface area-to-volume ratio combined with extensive fat deposits beneath the shell and throughout the head and neck region predisposes them to heat retention (Davenport et al., 2009a; Davenport et al., 2011). Vascular heat exchangers at the base of the flippers (Greer et al., 1973) and blood flow adjustments (Bostrom et al., 2010) permit tight control of heat loss at the extremities, and a vascular plexus lining the trachea minimizes respiratory heat loss (Davenport et al., 2009b). In order to maintain a steady Tb in cold water, the rate at which heat is lost to the environment must be matched by the rate of heat gain. The evidence to date suggests that metabolic heat production serves as the primary source of heat gain in leatherbacks (Paladino et al., 1990; Bostrom and Jones, 2007; Bostrom et al., 2010), and models indicate that the Tg observed in adult leatherbacks necessitates metabolic rates substantially higher than resting levels (Bostrom et al., 2010). Previous studies have demonstrated that juvenile leatherbacks (16–37 kg) increase swimming activity in response to a decrease in water temperature, resulting in a larger Tg in cool water compared with warm water (Bostrom et al., 2010). The relative importance of skeletal muscle heat production versus heat produced by other vital organs for thermoregulation in adult leatherbacks is less clear. Given high rates of prey intake (Heaslip et al., 2012), specific dynamic action (SDA, i.e. heat produced from ingestion, digestion, absorption and assimilation of a meal) (Secor, 2009) may be an important source of heat gain for adult leatherbacks at high latitude foraging grounds. Metabolic rate can increase by 1.6–2.1 times following meal ingestion in captive sea turtles (Jones and Seminoff, 2013).
Leatherbacks can use behavioral means to alter both sides of the heat gain=heat loss equation. Adjustments in swimming behavior allow leatherbacks to modulate rates of heat production (Bostrom and Jones, 2007; Bostrom et al., 2010), and adjustments in dive patterns allow leatherbacks to modulate rates of heat loss. For example, leatherbacks in tropical waters exhibit a significant negative correlation between dive depth and Tb (Southwood et al., 2005). Behavioral selection of cooler Ta at depth may act in conjunction with circulatory adjustments to permit effective heat transfer and prevent overheating while in tropical waters (Southwood et al., 2005; Wallace et al., 2005). Conversely, utilization of warmer surface waters while resident at high latitude foraging grounds may facilitate heat retention and maintenance of elevated Tb. Increased time at the surface would also allow for basking if environmental conditions were appropriate (James et al., 2005).
In this study, we document for the first time Tb of freely swimming adult leatherbacks during the resident foraging period and post-resident migration at high latitudes. Our primary goal was to investigate the behavioral and metabolic factors that contribute to the maintenance of elevated Tb under natural conditions. We deployed a combination of stomach temperature pills (STPs) and satellite-linked data recorders on leatherback turtles offshore Eastern Canada. The instrumentation provided data on dive patterns, Tb within the gastrointestinal tract, Ta, and location via the Argos Satellite System. Linear mixed models were used to explore the effects of dive behavior, diel status and residency status on Tb and the maximum change in Tb (ΔTb). We used our calculated values for mean Tg, a range of values for insulation thickness of the shell (L), and data and equations from the published literature to estimate rates
List of symbols and abbreviations
- surface area of shell
- curved carapace length
- specific heat capacity of prey
- thermal conductivity of shell
- insulation thickness
- mass of ingested prey
- meal energy
- metabolic rate
- platform transmitter terminal
- heat transferred across flippers
- heat transferred to ingested prey
- heat transferred across the shell
- heat loss
- specific dynamic action
- state space model
- stomach temperature pill
- ambient temperature
- body temperature
- thermal gradient
- change in body temperature
Satellite track durations ranged from 56 to 131 days, but for this study we focused only on data during the Tb monitoring period north of 42°N (range 12–60 days, Table 1). Turtles with Tb monitoring periods <45 days were assumed to have excreted the STP; cessation of Tb transmissions in all other turtles was assumed to reflect battery failure of the pill. The Tb transmissions for two turtles ceased during the resident period; Tb transmissions for five turtles extended into the post-resident period (Fig. 1). The mean rate of travel was 1.8 times higher during the post-resident period (2.7±0.7 km h−1) compared with the resident period (1.5±0.3 km h−1). Leatherbacks maintained a consistent dive pattern over the course of the entire monitoring period, with the same mean dive depth (20±6 m) during the resident and post-resident periods (Table 2). The mean percentage time at the surface was similar for the resident (41±8%) and post-resident (36±13%) periods (Table 2).
Mean Tb for individual turtles for the entire monitoring period ranged from 25.4±1.7 to 27.3±0.3°C. Mean Ta experienced by turtles over the course of the entire monitoring period ranged from 13.6±6.4 to 15.9±2.0°C, and mean Tg ranged from 10.7±2.4 to 12.1±1.7°C. Table 1 provides data on mean Tb, Ta and Tg for the resident and post-resident periods separately for each turtle. Fig. 2 provides an illustration of the daily variation in Tb for a female leatherback turtle (519 kg), and shows the stability of mean daily Tb in relation to Ta. The average mass of turtles was 467±76 kg (N=5). There was no significant correlation between mass and Tb (d.f.=3, r=0.809, P=0.097).
Results of the linear mixed models showed that the percentage time at the surface and diel status were the best predictors of mean Tb (Table 3). The 95% confidence intervals for the model estimates of percentage time at the surface (0.0326, 0.0902) and diel status (0.0014, 0.9888) excluded zero, which indicates a strong effect for both of these factors. The mean Tb increased with increased time at the surface, and mean night-time Tb tended to be higher than mean daytime Tb (Fig. 3). A strong diel effect was also observed for mean ΔTb. The 95% confidence interval for the model estimate of diel status (1.3365, 3.1808) excluded zero. Negative values for ΔTb (−1.2±1.1°C), indicative of cooling, occurred during the day and positive values for ΔTb (1.5±1.2°C), indicative of warming, occurred at night (Fig. 2B).
Fig. 4A illustrates estimated during the resident period of leatherbacks for which mass was measured. The mean decreases as insulation thickness (L) increases: 698±123 W at L=2 cm, 478±87 W at L=3.5 cm, 390±73 W at L=5 cm. The estimated mass-specific MR based on decreases accordingly as L increases (Fig. 4B): 1.498±0.122 W kg−1 at L=2 cm, 1.026±0.082 W kg−1 at L=3.5 cm and 0.837±0.066 W kg−1 at L=5 cm. As L increases, the proportion of attributable to heat loss across the body surface ( and ) decreases and the proportion due to increases (Fig. 4A). Mean is 186 W, or 0.397 W kg−1. The calculated rate of heat production due to SDA was 0.439 W kg−1.
Our results indicate that leatherbacks maintained Tg≥10°C in Canadian waters through a combination of both behavioral and physiological means. The amount of time that turtles spent in relatively warm (~15–17°C) surface waters had a strong influence on mean Tb (Table 3, Fig. 3). Individual leatherbacks spent 16–54% of time at the surface while north of 42°N. These results align with previous research at our study site, which documented surface times up to 41% for leatherbacks (James et al., 2005). Factors related to both foraging and thermoregulation are likely to contribute to the large amount of time spent at the surface. Leatherbacks have been observed to float at the surface to process large prey items (James and Mrosovsky, 2004; James et al., 2005), although video documentation from turtle-borne cameras suggests that the majority of prey can be manipulated and ingested in under 5 min (Heaslip et al., 2012). Modification of dive patterns to exploit warmer water at the surface and reduce rates of heat loss may be an important component of the leatherback's thermoregulatory strategy at high latitudes. Sightings of leatherbacks passively floating at the surface, flippers extended, during daylight hours at northern foraging grounds provide evidence that turtles may also bask to absorb solar radiation (James et al., 2005; Lindgren et al., 2014).
While basking may augment heat gain under suitable environmental conditions, the primary source of heat gain for leatherbacks derives from metabolic processes. Compelling evidence to support this assertion is found in the diel trends observed for Tb and ΔTb. It is important to note that we measured Tb from within the gastrointestinal tract; consequently, variation in Tb reflects cooling due to prey ingestion and warming due to SDA and transfer of heat via blood flow. We observed a trend for Tb of daytime cooling, suggestive of diurnal foraging, and night-time warming (Fig. 2). Model results did not show a strong correlation between dive parameters (depth and percentage time at the surface) and diel status, so alterations in dive behavior are unlikely to contribute to the night-time warming trend. We conclude that endogenous heat production warms gastrointestinal tract contents and accounts for the positive ΔTb observed during the night.
The non-surgical method used for deploying the STPs in this study limited our ability to know the precise location of the pills while monitoring Tb. It is possible that our measurements of Tb underestimate the true core temperature, as a result of decreases in Tb with prey ingestion. Alternatively, processes external to the gastrointestinal tract known to be involved in SDA [e.g. amino acid catabolism, anabolic processes and urea production by the liver (McCue, 2006)], as well as the localized effects of SDA within the stomach and small intestine, may result in Tb that is higher than temperatures in other parts of the core. The mean Tb for leatherbacks during the residency period in the current study (26.4±0.6°C, mass 391–589 kg) is 2.1°C higher than the mean cloacal temperature recorded for four leatherbacks captured at our study site and brought on-deck of a research vessel during a 2002–2003 field season (24.3±1.9°C, mass 315–380 kg) (James and Mrosovsky, 2004). Our calculations of Tg assume that Tb recorded in the gastrointestinal tract and averaged over a period of several days to weeks provides a good estimation of temperature throughout the body core.
Metabolic rates necessary to sustain observed Tg depend heavily on the thickness of the insulation layer underlying the shell (L) and rates of prey ingestion. Measurements of L obtained from stranded leatherbacks within our study area fall within the range of 2–5 cm (M.C.J., unpublished data). The mean mass-specific MR required to offset and maintain Tg for a foraging leatherback with L=3.5 cm is 1.026 W kg−1, 3.3 times higher than the average resting MR (0.312 W kg−1) (Wallace and Jones, 2008) but still considerably lower than maximal MR recorded for adult leatherbacks on nesting beaches (1.510 W kg−1) (Paladino et al., 1990). Insulation thickness is likely to increase over the course of the foraging season as leatherbacks store energy for migration, and this could result in a substantial decrease in and thermoregulatory costs (Fig. 4) (Davenport et al., 2011). Our estimates of MR based on field measurements of Tg are in agreement with outputs from previous models of leatherback thermoregulatory abilities based on laboratory studies (Bostrom and Jones, 2007; Bostrom et al., 2010).
Leatherbacks ingest large quantities of cold jellyfish during the foraging season (Heaslip et al., 2012), and a substantial amount of energy must be allocated to warm prey to core temperature (Davenport, 1998). The metabolic cost of warming food () is included in SDA, and can account for a large proportion of the total SDA response in cold environments (Secor, 2009). Based on our estimates for (0.397 W kg−1) and SDA (0.439 W kg−1), over 90% of the heat production attributed to SDA goes towards warming prey. Depending on the degree of insulation, SDA may account for ~30–50% of MR for leatherbacks during the foraging period (Fig. 4A). These estimates illustrate that visceral metabolic heat production plays an important role in maintaining stable Tb for leatherbacks foraging in cold water.
Cessation of area-restricted searching and an increase in the rate of travel occur as leatherbacks transition from resident foraging to post-resident migration. Surprisingly, model results indicated that residency status did not have a strong effect on mean Tb, and the pattern for ΔTb of daytime cooling and night-time warming persisted into the post-residency period, although it was slightly attenuated (Fig. 2). Furthermore, model results did not show a strong correlation between dive parameters (depth and percentage time at the surface) and residency status and there were no obvious shifts in dive behavior (Table 2). Taken together, these results provide evidence that leatherbacks continue to forage in the early stages of migration. The increase in the rate of travel with the onset of southward migration occurs without an appreciable change in dive patterns, which suggests that leatherbacks simply are moving in a more directed manner to cover greater distances in a given time. Under these circumstances, prey encounters may be opportunistic and prey ingestion rates are likely to be lower.
Measurements of Tb in freely swimming leatherbacks provide important insight into the factors that contribute to thermoregulation in this wide-ranging reptile. Our field data validate earlier models to assess leatherback thermoregulatory capacity (Paladino et al., 1990; Bostrom and Jones, 2007; Bostrom et al., 2010) and we provide additional support for endothermy in leatherbacks by illustrating the role of heat produced in visceral organs (SDA) for maintenance of Tg. As more data become available from remote monitoring studies conducted in both tropical and temperate waters, a clearer picture emerges of the integrative and versatile thermal strategy employed by leatherbacks. Mean Tb values recorded from leatherbacks foraging in Canadian waters (25.4–27.3°C) were remarkably similar to mean Tb recorded using the same techniques with leatherbacks during the nesting season in the Caribbean Sea (28.1–28.7°C) (Casey et al., 2010). The ability of leatherbacks to address thermal challenges and maintain high, stable Tb over a broad range of temperatures is the result of morphology and the orchestration of behavioral and physiological adjustments.
MATERIALS AND METHODS
Pre-calibrated stomach temperature pills (STP, 30 g in air, 63×24 mm, battery life 45–65 days, range 0–50°C, resolution 0.1°C; Wildlife Computers, Redmond, WA, USA) and platform transmitter terminals (PTT, 175 g in air, 93×52×26 mm, battery life 100–120 days; Wildlife Computers) were deployed on seven leatherbacks captured off Northern Cape Breton Island, NS, Canada (~47°N, 60°W), during August 2009 (Table 1). Turtles were live-captured at the surface, measured (curved carapace length, CCL), weighed, and equipped with identification tags (James et al., 2005). A PTT was attached using Tygon-coated stainless steel wires passed through drill holes (4 mm diameter) in the medial ridge (Casey et al., 2010). The turtle's mouth was opened using nylon straps and an STP was inserted ~40 cm into the esophagus using a lubricated, flexible tube (Casey et al., 2010). Turtles were released within 30 min of capture. All procedures were conducted in accordance with UNCW IACUC Protocol A0809-018 and Fisheries and Oceans Canada Licenses 323395 and 323398.
The ingested STP emitted a pulse-coded acoustic signal that corresponded to the coldest temperature detected by the instrument's thermistors (Tb). The Tb data were transmitted to and archived by the PTT at 1–2 min intervals. The effective range for communication between the STP and PTT was ≤2 m. Each PTT communicated with the Argos satellite system on a 24 h duty cycle for the first 40–44 days, and thereafter on a 48 h duty cycle. The PTT relayed archived time-series Tb data and binned data for Ta (range: −40 to +60°C, resolution 0.05°C) in 6 h blocks starting at 00:00 h (GMT). Data for Ta were transmitted as the total percentage of each 6 h block spent in a given temperature bin (time-at-temperature, TAT; see below). The PTTs also relayed dive behavior messages that included maximum depth (up to 1000 m, resolution 0.5 m), dive duration and post-dive surface time for dives with depths that exceeded 3 m. Turtles were considered to be at the surface at depths <3 m. Percentage time spent at the surface was calculated based on the dive duration and post-dive surface time for individual dives. The Argos satellite system provided Doppler-derived locations for the PTTs and an index of location accuracy from <150 to >1000 m.
Analysis and statistics
A state space model (SSM) interpolated latitude and longitude at 6 h intervals, using raw Argos locations (Jonsen et al., 2007). The SSM assigned a turtle's ‘state’ as either transiting or foraging based on rate of travel and turning frequency and angles. We interpreted a shift from resident to post-resident status when a turtle was no longer detected foraging north of 44°N. Only post-resident locations north of 42°N were included in our analyses to ensure that resident versus post-resident comparisons of Tb and Tg were made at similar Ta. Twenty rate of travel measurements, calculated as the distance between two consecutive SSM locations divided by the time lag between locations, were randomly sampled from each turtle's resident and post-resident period for descriptive statistics. A daily SSM location closest to 12:00 h (GMT) was used to determine the start/end time for daytime and night-time periods based on nautical twilight (http://aa.usno.navy.mil/data/docs/RS_OneYear.php). A random sample of 200 dive behavior messages from both the resident and post-resident periods were used to calculate descriptive statistics for dive parameters.
Analyses of time-series Tb data were constrained by the sporadic nature of STP transmissions. The Tb data were manually filtered for erroneous records, which occurred because of interrupted STP–PTT communication (OriginLab Corporation, Northampton, MA, USA). We then selected daytime and night-time periods that contained ≥30% of total possible Tb readings for diel analyses of mean Tb and maximum change in Tb (ΔTb). Random samples of 1500 Tb readings were drawn from the filtered data for the resident and post-resident periods to calculate mean Tb for each period and the entire monitoring period.
The manually filtered Tb data were also partitioned into 6 h blocks to correspond with TAT bin data for Ta. Only 6 h blocks with >40% of total possible Tb readings were analyzed. A single temperature value was assigned to each TAT bin for Ta: 0–5°C=5°C, 5–10°C=10°C, mean value for bins with intervals of 2°C, and >32°C=32°C. A proportional analysis was used to determine a representative Ta for each 6 h block based on the percentage time spent in each TAT bin. For example, if a turtle spent 60% of time in the 16–18°C bin and 40% of time in the 14–16°C bin, a Ta of 16.2°C [(0.6×17)+(0.4×15)] was assigned for that 6 h block. The Tg for each 6 h block was calculated using the mean Tb and Ta.
We used linear mixed models to investigate factors affecting Tb and ΔTb (lme4 package, R software, http://www.r-project.org). The mean percentage time at the surface, mean dive depth, diel status (day or night) and SSM state (resident or post-resident) were treated as fixed effects in the models, and turtle identity was treated as a random effect. Given the temporal limitations of the TAT data, we did not include Ta in the models. We assumed that Ta varies with depth. Models that exhibited subject variance ≥60% of the total variance were excluded from further analysis. We used second order AIC corrected for small sample size (AICc), AIC differences (Δi) and Akaike weights (wi) to evaluate models (Akaike, 1974). Confidence intervals calculated using Markov chain Monte Carlo simulations were used to assess the strength of fixed effects in models that exhibited the best fit. We did not include mass in the models, as we had just one measurement taken at the time of capture and could not predict mass change over the course of the foraging season (Davenport et al., 2011). The potential influence of mass on mean Tb for the entire monitoring period was investigated using Pearson's correlation.
To estimate metabolic heat production for leatherbacks, we assumed that (1) blood flow results in continual heat exchange throughout body core (Paladino et al., 1990) and (2) mean Tb reflects maintenance of an internal thermal steady state in relation to Ta. Under these conditions, heat production (i.e. metabolic rate, MR) is equivalent to total heat transfer () (Bostrom et al., 2010), i.e. the sum of heat transfer to ingested prey () and heat transfer to the external environment (i.e. heat loss) across the shell () and flippers (). Heat loss from the head and neck is assumed to be negligible as a result of extensive insulation (Davenport et al., 2009a). Heat loss due to excretion has not been well characterized for leatherbacks, but is assumed to be small relative to other sources of heat loss: (1) (Bostrom et al., 2010) where A is the surface area of the shell (m2), calculated as A=0.049mass0.69 (Bostrom et al., 2010); k is the thermal conductivity of the shell, 0.25 W m−1 K−1 (Bostrom et al., 2010); Tg is the mean thermal gradient (K); and L is insulation thickness (m): (2) as accounts for approximately 7% of total surface heat loss in cold water (Bostrom et al., 2010): (3) (Bostrom and Jones, 2007) where MP is prey mass consumed, equivalent to 73% of leatherback mass per day (Heaslip et al., 2012), or mass (8.45×10−6) for units of kg s−1; and CP is the specific heat capacity of prey (4186 J kg−1 K1) (Bostrom and Jones, 2007). We assume our calculated Tg provides a good approximation of the difference between core Tb and prey temperature, and warming of prey prior to ingestion is minimal (Heaslip et al., 2012).
Estimates for rate of heat production due to SDA were obtained by applying a general equation for SDA of reptiles: SDA=0.26ME–10.65, where ME is the meal energy (kJ) ingested per day based on published energy values for jellyfish prey (0.2 kJ g−1) (Doyle et al., 2007) and prey consumption equivalent to 73% of body mass per day (340,910 g of jellyfish for a 467 kg turtle) (Heaslip et al., 2012). The resulting value of SDA (kJ) for an average day's foraging effort was then converted to W kg−1, using our mean mass of 467 kg.
We are grateful to Molly Lutcavage for support of this research. We thank Bert Fricker, Blair Fricker, Josh Fricker, Kelly Fricker, and the Canadian Sea Turtle Network for logistical support, including assistance with tag deployments. We gratefully acknowledge Ian Jonsen for state-space modeling expertise. Leigh Anne Harden provided useful feedback on statistical analyses, and Christina Davis contributed to early analyses of dive behavior data.
J.P.C., M.C.J. and A.S.W. contributed to the conception and design of the experiment, interpretation of findings, and drafting and revision of the manuscript. J.P.C. and M.C.J. conducted field work, and A.S.W. conducted metabolic analyses.
The authors declare no competing financial interests.
This research was funded by the Large Pelagics Research Center (University of Massachusetts Amherst) [NA04NMF4550391 to A.S.W.], Environment Canada (M.C.J.), Fisheries and Oceans Canada (M.C.J.), and the Government of Canada Habitat Stewardship Program for Species at Risk (M.C.J.).
- © 2014. Published by The Company of Biologists Ltd