DOI: 10.5061/dryad.p8cz8w9v0
title: “Breath holding” as a thermoregulation strategy in the scalloped hammerhead, a deep-diving tropical ectothermic shark

description: Fish moving between different thermal environments experience heat
 exchange via conduction through the body wall and convection from blood
 flow across the gills. Here we report a strategy of preventing convective
 heat loss at the gills during excursions into deep cold water by the
 tropical scalloped hammerhead shark (Sphryna lewini). Adult scalloped
 hammerhead sharks dive rapidly and repeatedly from warm (~26ºC) surface
 waters to depths exceeding 800 meters and temperatures as low as 5°C.
 Biologgers attached to adult sharks show that warm muscle temperatures
 were maintained throughout the deepest portion of each dive. Substantive
 cooling only occurred during the latter stages of the ascent phase, and
 once initiated, was rapid. Heat transfer coefficient modeling indicated
 that convective heat transfer was suspended, probably by suppressing gill
 function during deep dives. This previously unobserved strategy has broad
 similarities to marine mammal “breath hold” diving.
description: Note: (Refer to References section of publication for numbered
 citations) <strong>Measuring swimming
 performance, environment, and body temperature</strong>
 To measure scalloped hammerhead shark swimming performance,
 depth, ambient water and muscle temperatures we used an instrument package
 consisting of a tri-axial accelerometer tag and a depth and temperature
 archival tag housed in a syntactic foam float (2,000 m depth rating)
 equipped with a timed release mechanism and radio beacon transmitters to
 facilitate recovery.  The tri-axial accelerometer tag was either a
 TDR10-XB-340 (56 x 38 x 24mm 69g; Wildlife Computers., Redmond, WA) or a
 TDR10-Daily Diary-278 (74 x 57 x 36mm, 117g; Wildlife Computers., Redmond,
 WA). Tri-axial acceleration was sampled at either 16Hz or 32Hz.  Depth,
 water temperature and body temperature were sampled every 5 or 10 seconds
 using a MK9 archival tag with an 8cm thermistor stalk (Wildlife Computers,
 Redmond, WA). Each package also contained a SPOT5 or SPOT6
 satellite-linked transmitter (80 x 20 x 11mm, 30g; Wildlife Computers.,
 Redmond, WA) to indicate the package’s position when it floated to the
 surface following release from the tagged animal. Packages that contained
 a SPOT5 tag were also equipped with a VHF transmitter  (MM130B; 16mm
 diameter, 60mm length, 20g; ATS, USA).
 <strong>Shark capture and handling</strong>
 All sharks were caught using baited hooks on demersal longlines
 inside Kāneʻohe Bay (N 21.45°, W 157.80°), Oʻahu (Hawaiʻi, USA).  To
 ensure captured sharks were in good condition, lines were checked every 30
 minutes. Captured sharks were brought along the side of a 17ft skiff,
 secured with a rope around the caudal peduncle and provided with constant
 gill ventilation during measurement and instrumentation via a hose
 inserted into the mouth and connected to an in-water bilge pump. The tag
 package was attached to each shark via a fusible stainless steel cable tie
 (360 mm, 8 g; Little Leonardo Co., Tokyo, Japan) passed through two holes
 drilled through the base of the dorsal fin and secured around the
 syntactic foam float. The sensor stalk of the MK9 was inserted
 approximately 8 cm below the skin into the dorsal musculature while the
 package was being secured to the dorsal fin. Each package contained a
 timed-release mechanism (RT-4, 16mm diameter x 19mm length 10g; RT-5, 20mm
 diameter x 38mm length, 20g; Little Leonardo Co., Tokyo, Japan) set to
 release after 23 days. At release time, an integrated capsule severed the
 stainless steel band, releasing the package to float to the surface. We
 used a small vessel and handheld receivers equipped with directional
 antennas to home in on and recover the floating package locator beacons.
 Contact information was written on the packages in case they were found by
 members of the general public. Individual deployment durations ranged from
 23 to 29 days (Table S1). <strong>Swimming
 behavior data processing</strong> Archived
 accelerometer and depth data were downloaded from ten recovered tag
 packages. A total of 103 deep dives were recorded from 7 of the 10 sharks
 (Supplemental Table S1). Average maximum dive depth was 635.5m (± 91.8m;
 maximum 825m [HH10]). All 32Hz tri-axial acceleration data were resampled
 at 16Hz to facilitate analyses. Acceleration and depth data were analyzed
 using the ‘Ethographer’ package (<em>46)</em> in Igor Pro 8
 (WaveMetrics Inc., Portland, OR, USA).  A low pass filter of 0.3 Hz was
 used to estimate the static (gravitational) and dynamic (tail stroking)
 components of the acceleration signal for each axis. The static
 acceleration components from the x,y,z axes were used to calculate the
 roll and pitch angles of the shark, where x is the surge axis, y is the
 sway axis, and z is the heave axis (<em>47</em>):
 Eq. 2 Roll =
 arctan(y/x<sup>2</sup>+z<sup>2</sup>)<sup>1/2</sup>) (180/π) Eq. 3 Pitch = arctan(x/y<sup>2</sup>+z<sup>2</sup>)<sup>1/2</sup>) (180/π) The tag attachment angle was corrected to zero degrees centered for roll and to horizontal for pitch by quantifying tag pitch during upright swimming at a constant depth and subtracting this value from the archived data (<em>48</em>) (<em>49)</em>. The overall dynamic body acceleration (ODBA) was calculated by summing the absolute values of the dynamic acceleration from all three orthogonal axes (<em>23</em>) (<em>50</em>). We used continuous wavelet transformation to generate a spectrogram of the dynamic component of swaying acceleration, classified the dominant peak as tail beat cycles, and calculated tail beat frequency and amplitude of acceleration at every second. <strong>Thermal conductivity coefficient analysis</strong> Packages from three sharks (HH7, HH8, HH9) contained useful body temperature data for thermal coefficient analysis. HH11’s thermal probe was not inserted into the musculature and HH10 had unusable body temperature data, likely due to a malfunctioning thermal probe (see Table S1). HH1-3 and HH5 remained in shallow water and did not conduct deep dives following tagging. The package for HH4 was never recovered. To test the sensor response times of the ‘stalk’ thermistor (the one that measured body temperature) and the thermistor embedded in the body of the tag (the one that measured ambient temperature) we empirically tested the thermistor response times of the MK9 tags used for our study and compared our values with previously published values for MK9 tag thermistor response times in (<em>51</em>). Our empirical tests were similar to those of (<em>51</em>) and consisted of moving the MK9 tags back and forth between two buckets of water differing in temperature by 20 ºC. Our tests results concur with those of (<em>51</em>) and show a fast response in the ‘stalk’ thermistor that was inserted into the scalloped hammerhead musculature, and a minor lag in the thermistor embedded in the epoxy casing of the MK9. The lag response of the embedded MK9 thermistor yielded different depth-temperature profiles between the fast descents and the comparatively slower ascents of each dive. To account for this disparity, the depth-temperature profiles of the ascents were used to adjust the depth-temperature profile of the fast descents. To evaluate whether scalloped hammerhead sharks actively retain body heat or rely on simple thermal inertia during deep dives into cold water, predicted values of the whole-body heat transfer coefficient (<em>k</em>) and the rate of temperature change due to metabolic heat production (Τ<sub>o</sub>) were modeled to match the observed rates of body warming and cooling (sensu (<em>27</em>) (<em>28</em>) (Eq. 1)). The modeling used in (<em>27</em>) (<em>28</em>) (<em>30</em>) assumed the following conditions for <em>k</em>; (a) k  =  a constant value; (b) k = {<em>k<sub>1</sub></em>; <em>Ta</em>(t) &lt; <em>Tb</em> (t) <em>k<sub>2</sub></em>; <em>Ta</em>(t) ≥ <em>Tb</em>(t), where <em>k<sub>1</sub></em> and <em>k<sub>2</sub></em><sub> </sub>are two values for the whole-body heat transfer coefficient. Model (a) assumes the rate of heat exchange between the shark and the environment is constant and not altered by the shark. Model (b) assumes the rate of heat exchange is different when the shark is cooling (<em>k<sub>1</sub></em>) compared to warming (<em>k<sub>2</sub></em>). These conditions were not able to adequately model the body temperature of the sharks (all R2 fits &gt;0.60). A post-hoc analysis of swimming activity (ODBA, tailbeat frequency and amplitude) revealed high activity levels during each deep dive (Fig. S5). Body temperature was observed to increase throughout the deepest portion of each dive and then followed by a sharp cooling rate during ascent to the surface and subsequent decrease in swimming activity (Fig. 2, S5-8). These high activity rates will affect metabolic heat generation (Τ<sub>o</sub>) and possibly rates of blood flow through the core muscles or through the gills. These phenomena would ultimately affect the rate of change of body temperature (<em>dTb(t)/dt</em>) and whole-body thermal rate coefficient (<em>k</em>) (equation <em>3</em>) (<em>3</em>) (<em>31</em>). We therefore modeled predicted values of the whole-body heat transfer coefficient (<em>k</em>) and the rate of temperature change due to metabolic heat production (Τ<sub>o</sub>) based on different conditions associated with the periods of intense swimming activity during each dive. The following four conditions of <em>k</em> and Τ<sub>o </sub>were assumed: (1)  <em>k  </em>= <em> a constant value,</em> Τo = <em>a constant value</em>; (2)  <em>k </em>=<em> a constant value,</em> Τo = {Τ<sub>o1</sub>; <em>normal swimming</em> Τ<sub>o2</sub>; <em>deep diving</em>}; (3) <em>k </em>= <em>{k1; normal swimming k2</em>; deep diving}, Τ<sub>o</sub> = <em>a constant value</em>; (4)  <em>k </em>= <em>{k1; normal swimming k2; deep diving</em>}, Τo = {Τ<sub>o1</sub>; <em>normal swimming</em> Τ<sub>o2</sub>; <em>deep diving</em>}; where <em>k1<sub> </sub></em>and <em>k<sub>2</sub></em> are two values for the whole-body heat transfer coefficient. Model (1) assumes the rate of heat exchange between the shark and the environment (<em>k</em>) and rate of temperature change due to metabolic heat production (Τ<sub>o</sub>) is constant and not altered by the shark. Model (2) assumes rate of heat exchange between the shark and the environment (<em>k</em>) is constant while the rate of temperature change due to metabolic heat production (Τ<sub>o</sub>) alters between (Τ<sub>o1</sub>) during normal swimming and (Τ<sub>o2</sub>) during the high activity phases of a deep dive. Model (3) assumes (<em>k<sub>1</sub></em>) when the shark is swimming normally and (k2) during the high activity phases of a deep dive while the rate of temperature change due to metabolic heat production (Τo) is constant. Model (4) assumes (<em>k<sub>1</sub></em>) and (Τ<sub>o1</sub>) when the shark is swimming normally and (k2) and (Τ<sub>o2</sub>) during the high activity phases of a deep dive. The determination of switching (models 3 & 4) between (<em>k<sub>1</sub>)</em> (normal swimming) and (<em>k<sub>2</sub></em>) (deep diving) was based on established thresholds of vertical velocity between the descent phase and ascent phase of each dive. The same thresholds were used for the determination of switching (models 2 & 4) between (Τ<sub>o1</sub>) (normal swimming) and (Τ<sub>o2</sub>) (deep diving). The optimized parameters for each model were estimated using an iterative least-squares algorithm comparing the observed and predicted body temperature, and coefficient of determination (R-squared) computed to assess the fit across models. Models were run for each night that consisted of more than one deep dive for each shark. Start and stop times for each model duration were arbitrarily set to include all deep dives within a night duration and to include the durations where body-temperature reaches the ambient temperature after the conclusion of a series of deep dives. While both constant-k models yielded poor fits, both variable-k models performed well, with model (4) giving higher R2 values but also requiring an additional fit parameter (Τ<sub>o1</sub> and Τ<sub>o2</sub> vs only Τ<sub>o</sub>). Given the significantly enhanced tailbeat activity occurring during the rapid diving phase it is reasonable to assume higher metabolic heat generation during this period and that different Τo values should be used to reflect different metabolic heat generation during the low and high activity phases. <strong>Thermal conductivity of dead sharks</strong> To account for the rate of temperature change in the absence of convective heat transfer (blood flow across the gills), the body temperatures of two dead adult male scalloped hammerhead sharks with total lengths 240cm and 273cm were recorded by individually placing the sharks in an insulated tub and inducing abrupt alternations of water temperature from 23 ºC to 6 ºC and vice versa (<em>sensu (2) (30</em>)). Core muscle and ambient water temperatures were measured using the same TDR-Mk9 archival tag used in field experiments. Constant heat-transfer coefficients k were calculated using the heat budget model without metabolic heat production (Τ<sub>o</sub>) = 0 (Table S6). <strong>Allometric relationship of thermal coefficient rates</strong> To compare the influence of body size on the rate of heat exchange between the scalloped hammerheads from this study and other large fish species, the heat transfer coefficients from scalloped hammerheads HH7, HH8, and HH9 from model (4) were plotted with heating and cooling thermal coefficients against body mass published from other fish species collated in (<em>32</em>). For the scalloped hammerheads, normal swimming heat transfer coefficients (<em>k<sub>1</sub></em>) were incorporated into the warming regression and high activity deep dive phase (<em>k<sub>2</sub></em>) heat transfer coefficients from model (4) were incorporated into the cooling regression. Body mass estimates for each scalloped hammerhead were calculated using length-weight relationships from (<em>52</em>).
description: Thermal coefficient modeling of all scalloped hammerhead shark
 deep-diving was run using Matlab code "Script5.m" with Matlab
 version 2021b. (included in this dataset).  Thermal
 coefficient modeling of the dead scalloped hammerhead sharks was run using
 Matlab code "Script3.m" (included in this dataset) with Matlab
 version R2020a.  Telemtry was analyzed using IgroPro 8
 (Wavemetrics, Portland, OR).