Thermal performance responses in free-ranging elasmobranchs depend on habitat use and body size

Predicting responses to climate change using a joint species, spatially dependent physiologically guided abundance model

Predicting the effects of warming temperatures on the abundance and distribution of organisms under future climate scenarios often requires extrapolating species-environment correlations to climatic conditions not currently experienced by a species, which can result in unrealistic predictions. For poikilotherms, incorporating species' thermal physiology to inform extrapolations under novel thermal conditions can result in more realistic predictions. Furthermore, models that incorporate species and spatial dependencies may improve predictions by capturing correlations present in ecological data that are not accounted for by predictor variables. Here, we present a joint species, spatially dependent physiologically guided abundance (jsPGA) model for predicting multispecies responses to climate warming. The jsPGA model uses a basis function approach to capture both species and spatial dependencies. We apply the jsPGA model to predict the response of eight fish species to projected climate warming in thousands of lakes in Minnesota, USA. By the end of the century, the cold-adapted species was predicted to have high probabilities of extirpation across its current range-with 10% of lakes currently inhabited by this species having an extirpation probability >0.90. The remaining species had varying levels of predicted changes in abundance, reflecting differences in their thermal physiology. Though the model did not identify many strong species dependencies, the variation in estimated spatial dependence across species suggested that accounting for both dependencies was important for predicting the abundance of these fishes. The jsPGA model provides a new tool for predicting changes in the abundance, distribution, and extirpation probability of poikilotherms under novel thermal conditions.

Upper thermal tolerance and population implications for the Magdalena River stingray Potamotrygon magdalenae

Knowledge of thermal tolerance limits provides important clues to the capacity of a species to withstand acute and chronic thermal changes. Climate models predict the increase and intensification of events such as heat waves, therefore understanding the upper thermal limits that a species can tolerate has become of utmost importance. We measured the upper thermal tolerance of the endemic Magdalena river stingray Potamotrygon magdalenae acclimated to experimental conditions, and then used critical thermal methodology to find the temperature at which an organism reaches a critical endpoint where locomotory activity becomes disorganized and the animal loses its ability to escape from conditions that will promptly lead to its death. We also describe the behavioral response of individuals to acute thermal stress and infer the possible consequences of temperature increases in the habitats of P. magdalenae populations. There were no significant differences between sexes in temperature tolerance or behavior. The critical thermal maximum (39°C) was 5.9°C above the maximum recorded temperature for the study area. Although P. magdalenae was tolerant to high temperature and currently is not living at its upper thermal limit, its survival in Guarinocito Pond will be threatened if temperatures continue to increase, considering the warming scenarios predicted for tropical regions due to climate change, even including short-term climate phenomena such as El Niño.

Intraspecific variation in muscle growth of two distinct populations of Port Jackson sharks under projected end-of-century temperatures

Although pervasive, the effects of climate change vary regionally, possibly resulting in differential behavioral, physiological, and/or phenotypic responses among populations within broadly distributed species. Juvenile Port Jackson sharks (Heterodontus portusjacksoni) from eastern and southern Australia were reared at their current (17.6 °C Adelaide, South Australia [SA]; 20.6 °C Jervis Bay, New South Wales [NSW]) or projected end-of-century (EOC) temperatures (20.6 °C Adelaide, SA; 23.6 °C Jervis Bay, NSW) and assessed for morphological features of skeletal muscle tissue. Nearly all skeletal muscle properties including cellularity, fiber size, myonuclear domain, and satellite cell density did not differ between locations and thermal regimes. However, capillary density was significantly influenced by thermal treatment, where Adelaide sharks raised at current temperatures had a lower capillarity than Jervis Bay sharks raised at ambient or projected EOC temperatures. This may indicate higher metabolic costs at elevated temperatures. However, our results suggest that regardless of the population, juvenile Port Jackson sharks may have limited acclimatory potential to alter muscle metabolic features under a temperature increase, which may make this species vulnerable to climate change.

Provisioning ecotourism does not increase tiger shark site fidelity

A perennial criticism of provisioning ecotourism is that it alters the natural behavior and ecology of the target species by providing an artificial food source. Here we evaluate its impact on the long-term site fidelity patterns of tiger sharks in French Polynesia. We hypothesized that a significant impact of provisioning would lead to (1) increases in individual site fidelity over time, and (2) an increase in the number of resident individuals over time. Of 53 individuals photo-identified and monitored during > 500 dives over five years, 10 individuals accounted for > 75% of all sightings, whereas 35 sharks were sighted very infrequently. Even the most frequently observed tiger sharks exhibited overall low fidelity at the site and showed no increase in site fidelity over time. Furthermore, the number of tiger sharks sighted during each dive did not increase. The observed patterns of tiger shark sightings were best explained by natural movements, including general roaming within home ranges along the coastline and seasonal migrations. Despite the apparent lack of impact of provisioning ecotourism on tiger shark ecology in Tahitian waters, it would be prudent to implement a strict code of conduct during any future provisioning activities to maximize the safety of participants and animals involved.

Thermal physiology integrated species distribution model predicts profound habitat fragmentation for estuarine fish with ocean warming

Species distribution models predict a poleward migration for marine ectotherms with ocean warming. However, a key limitation in current species distribution models (SDM) is that they do not account for population-specific heterogeneity in physiological responses to temperature change resulting from local adaptations and acclimatization. To address this gap, we developed a novel, Physiology Integrated BioClimate Model (PIBCM) that combines habitat-specific metabolic thermal physiological tolerance of a species into a bioclimate envelope model. Using a downscaling approach, we also established a fine-resolution coastal sea-surface temperature data set for 2050-2080, that showed a high degree of location-specific variability in future thermal regimes. Combining predicted temperature data with the PIBCM model, we estimated habitat distribution for a highly eurythermal intertidal minnow, the Atlantic killifish (Fundulus heteroclitus), a species that likely presents a best-case-scenario for coastal vertebrates. We show that the killifish northern boundary shifts southwards, while distinct habitat fragmentation occurs in the southern sub-population (due to migration of adjacent fish populations to the nearest metabolically optimal thermal habitat). When compared to current SDMs (e.g., AquaMaps), our results emphasize the need for thermal physiology integrated range shift models and indicate that habitat fragmentation for coastal fishes may reshape nursery habitats for many commercially and ecologically important species.

The distribution of manta rays in the western North Atlantic Ocean off the eastern United States

In 2018, the giant manta ray was listed as threatened under the U.S. Endangered Species Act. We integrated decades of sightings and survey effort data from multiple sources in a comprehensive species distribution modeling (SDM) framework to evaluate the distribution of giant manta rays off the eastern United States, including the Gulf of Mexico. Manta rays were most commonly detected at productive nearshore and shelf-edge upwelling zones at surface thermal frontal boundaries within a temperature range of approximately 20–30 °C. SDMs predicted highest nearshore occurrence off northeastern Florida during April, with the distribution extending northward along the shelf-edge as temperatures warm, leading to higher occurrences north of Cape Hatteras, North Carolina from June to October, and then south of Savannah, Georgia from November to March as temperatures cool. In the Gulf of Mexico, the highest nearshore occurrence was predicted around the Mississippi River delta from April to June and again from October to November. SDM predictions will allow resource managers to more effectively protect manta rays from fisheries bycatch, boat strikes, oil and gas activities, contaminants and pollutants, and other threats.

Ocean warming alters the distributional range, migratory timing, and spatial protections of an apex predator, the tiger shark (Galeocerdo cuvier)

Given climate change threats to ecosystems, it is critical to understand the responses of species to warming. This is especially important in the case of apex predators since they exhibit relatively high extinction risk, and changes to their distribution could impact predator-prey interactions that can initiate trophic cascades. Here we used a combined analysis of animal tracking, remotely sensed environmental data, habitat modeling, and capture data to evaluate the effects of climate variability and change on the distributional range and migratory phenology of an ectothermic apex predator, the tiger shark (Galeocerdo cuvier). Tiger sharks satellite tracked in the western North Atlantic between 2010 and 2019 revealed significant annual variability in the geographic extent and timing of their migrations to northern latitudes from ocean warming. Specifically, tiger shark migrations have extended farther poleward and arrival times to northern latitudes have occurred earlier in the year during periods with anomalously high sea-surface temperatures. A complementary analysis of nearly 40 years of tiger shark captures in the region revealed decadal-scale changes in the distribution and timing of shark captures in parallel with long-term ocean warming. Specifically, areas of highest catch densities have progressively increased poleward and catches have occurred earlier in the year off the North American shelf. During periods of anomalously high sea-surface temperatures, movements of tracked sharks shifted beyond spatial management zones that had been affording them protection from commercial fishing and bycatch. Taken together, these study results have implications for fisheries management, human-wildlife conflict, and ecosystem functioning.

Forecasting intraspecific changes in distribution of a wide-ranging marine predator under climate change

Globally, marine animal distributions are shifting in response to a changing climate. These shifts are usually considered at the species level, but individuals are likely to differ in how they respond to the changing conditions. Here, we investigate how movement behaviour and, therefore, redistribution, would differ by sex and maturation class in a wide-ranging marine predator. We tracked 115 tiger sharks (Galeocerdo cuvier) from 2002 to 2020 and forecast class-specific distributions through to 2030, including environmental factors and predicted occurrence of potential prey. Generalised Linear and Additive Models revealed that water temperature change, particularly at higher latitudes, was the factor most associated with shark movements. Females dispersed southwards during periods of warming temperatures, and while juvenile females preferred a narrow thermal range between 22 and 23 °C, adult female and juvenile male presence was correlated with either lower (< 22 °C) or higher (> 23 °C) temperatures. During La Niña, sharks moved towards higher latitudes and used shallower isobaths. Inclusion of predicted distribution of their putative prey significantly improved projections of suitable habitats for all shark classes, compared to simpler models using temperature alone. Tiger shark range off the east coast of Australia is predicted to extend ~ 3.5° south towards the east coast of Tasmania, particularly for juvenile males. Our framework highlights the importance of combining long-term movement data with multi-factor habitat projections to identify heterogeneity within species when predicting consequences of climate change. Recognising intraspecific variability will improve conservation and management strategies and help anticipate broader ecosystem consequences of species redistribution due to ocean warming.

Connecting post-release mortality to the physiological stress response of large coastal sharks in a commercial longline fishery

Bycatch mortality is a major factor contributing to shark population declines. Post-release mortality (PRM) is particularly difficult to quantify, limiting the accuracy of stock assessments. We paired blood-stress physiology with animal-borne accelerometers to quantify PRM rates of sharks caught in a commercial bottom longline fishery. Blood was sampled from the same individuals that were tagged, providing direct correlation between stress physiology and animal fate for sandbar (Carcharhinus plumbeus, N = 130), blacktip (C. limbatus, N = 105), tiger (Galeocerdo cuvier, N = 52), spinner (C. brevipinna, N = 14), and bull sharks (C. leucas, N = 14). PRM rates ranged from 2% and 3% PRM in tiger and sandbar sharks to 42% and 71% PRM in blacktip and spinner sharks, respectively. Decision trees based on blood values predicted mortality with >67% accuracy in blacktip and spinner sharks, and >99% accuracy in sandbar sharks. Ninety percent of PRM occurred within 5 h after release and 59% within 2 h. Blood physiology indicated that PRM was primarily associated with acidosis and increases in plasma potassium levels. Total fishing mortality reached 62% for blacktip and 89% for spinner sharks, which may be under-estimates given that some soak times were shortened to focus on PRM. Our findings suggest that no-take regulations may be beneficial for sandbar, tiger, and bull sharks, but less effective for more susceptible species such as blacktip and spinner sharks.

Temporal niche partitioning as a novel mechanism promoting co-existence of sympatric predators in marine systems

Niche partitioning of time, space or resources is considered the key to allowing the coexistence of competitor species, and particularly guilds of predators. However, the extent to which these processes occur in marine systems is poorly understood due to the difficulty in studying fine-scale movements and activity patterns in mobile underwater species. Here, we used acceleration data-loggers to investigate temporal partitioning in a guild of marine predators. Six species of co-occurring large coastal sharks demonstrated distinct diel patterns of activity, providing evidence of strong temporal partitioning of foraging times. This is the first instance of diel temporal niche partitioning described in a marine predator guild, and is probably driven by a combination of physiological constraints in diel timing of activity (e.g. sensory adaptations) and interference competition (hierarchical predation within the guild), which may force less dominant predators to suboptimal foraging times to avoid agonistic interactions. Temporal partitioning is often thought to be rare compared to other partitioning mechanisms, but the occurrence of temporal partitioning here and similar characteristics in many other marine ecosystems (multiple predators simultaneously present in the same space with dietary overlap) introduces the question of whether this is a common mechanism of resource division in marine systems.

Effects of climate-change-driven gradual and acute temperature changes on shark and ray species

Climate change is altering distributions and abundances of marine species through both gradual and acute changes in temperature and productivity. Due to their high mobility and metabolic rates, elasmobranchs (sharks and rays) are likely to redistribute across latitudes and depths as they thermoregulate, but little is known about their responses to these climatic changes, which could vary widely across this diverse group of species. Here, we assessed how species with differing mobility and ecology responded to gradual changes in daily sea surface temperature (SST) and acute temperature anomalies, caused by the El Niño-Southern Oscillation (ENSO), at Cocos Island, Costa Rica, the site of multiple marine heatwaves. We used generalized linear mixed models to analyse 34,342 records of relative abundance or frequency of occurrence for seven shark and ray species collected in 27 years (1993-2019) by a dive company. We compared effect sizes for SST and the Oceanic Niño Index across the different species, which vary widely in body size and mobility. Large, mobile species responded strongly but inconsistently to temperature. For scalloped hammerhead sharks Sphyrna lewini, a 1℃ rise in SST reduced counts by over 14%, and dropped the occurrence of their large schools by almost one-fifth (19.4%). Mobula ray occurrence also declined substantially with a few degrees rise in SST, whereas tiger shark Galeocerdo cuvier occurrence sharply increased. These species also had divergent responses to the ENSO: S. lewini and G. cuvier were sighted with greater frequency during La Niña events, and their abundance dropped considerably during El Niño events-over a twofold decline between a strong La Niña and strong El Niño for S. lewini. In contrast, Mobula rays showed little response to ENSO. The smaller and sedentary Triaenodon obesus exhibited the weakest response of all species to both SST and the ENSO, reflecting its lower metabolic rates and mobility. Climate change will continue to impact elasmobranchs, even for smaller and more localized species, with the potential to impact the effectiveness of marine protected areas (MPAs). Our results compel further work on the diversity of elasmobranch responses to environmental change.

Large-scale space use of large juvenile and adult smalltooth sawfish Pristis pectinata: implications for management

The smalltooth sawfish Pristis pectinata is an endangered species endemic to the Atlantic Ocean. The only known viable populations occur in the USA along both coasts of Florida and in the western Bahamas. Little is known about habitat use and movement ecology of large juvenile and adult smalltooth sawfish. Although Critical Habitat—a management designation in the USA—has been identified for small juveniles, it has yet to be identified for these life stages. Between May 2016 and April 2019, we used passive acoustic telemetry and 3 large data sharing networks of receivers to track movements of 43 large juvenile and adult smalltooth sawfish. During this study, 24 females and 19 males were implanted with transmitters with estimated 4 or 10 yr battery lives. These tagged individuals were detected off the southeastern USA on 461 receivers ranging from off the coast of Brunswick, Georgia, to the lower Florida Keys, and along the Gulf coast to Apalachee Bay, Florida. Seasonal migrations were undertaken by 58% (43% mature; 57% immature) of the tagged individuals, with the remainder being apparent residents of their tagging locations. Tagged sawfish from both size classes and of both sexes migrated, which indicates that neither sex nor length is a predictor of whether a sawfish will migrate or not. Although both coasts of Florida were used for migration, most individuals consistently used the same coast when they migrated. The areas surrounding Boca Grande, Cape Canaveral, and the lower Florida Keys were heavily visited sites that could be further evaluated as potential Critical Habitat for these life stages. Understanding the movement patterns of this Critically Endangered species is essential for creating policies to protect areas important for promoting growth of the population.

Wet season flood magnitude drives resilience to dry season drought of a euryhaline elasmobranch in a dry-land river

The increase in severity and occurrence of drought from environmental change poses a significant threat to freshwater ecosystems. However, many of the mechanisms by which periodic drought affects aquatic animals are poorly understood. Here we integrated physical, physiological, and behavioural measurements made in the field over a twelve-year period to provide a comprehensive understanding of the factors affecting the loss of body condition of fish in arid rivers, using the Critically Endangered freshwater sawfish (Pristis pristis) in the dryland Fitzroy River, Western Australia, as a model species. Sawfish lost condition throughout the long dry season in all years and had significantly poorer body condition throughout years characterized by low volumes of wet season flooding and little occurrence of overbank flooding. A mechanistic examination of factors leading to this loss of condition using measurements of body temperature, field energetics, and habitat use from telemetry techniques showed that the loss of condition throughout the season was likely due to substantial habitat compression and low productivity in drier years, while high rates of competition were more likely to drive this pattern in wetter years. This information can be used to forecast how climate change and water abstraction will affect aquatic fauna experiencing intermittent drought and can inform management decisions to help mitigate these threats.

Population variation in the thermal response to climate change reveals differing sensitivity in a benthic shark

Many species with broad distributions are exposed to different thermal regimes which often select for varied phenotypes. This intraspecific variation is often overlooked but may be critical in dictating the vulnerability of different populations to environmental change. We reared Port Jackson shark (Heterodontus portusjacksoni) eggs from two thermally discrete populations (i.e. Jervis Bay and Adelaide) under each location's present-day mean temperatures, predicted end-of-century temperatures and under reciprocal-cross conditions to establish intraspecific thermal sensitivity. Rearing temperatures strongly influenced ṀO_{2 Max} and critical thermal limits, regardless of population, indicative of acclimation processes. However, there were significant population-level effects, such that Jervis Bay sharks, regardless of rearing temperature, did not exhibit differences in ṀO_{2 Rest} , but under elevated temperatures exhibited reduced maximum swimming activity with step-wise increases in temperature. In contrast, Adelaide sharks reared under elevated temperatures doubled their ṀO_{2 Rest} , relative to their present-day temperature counterparts; however, maximum swimming activity was not influenced. With respect to reciprocal-cross comparisons, few differences were detected between Jervis Bay and Adelaide sharks reared under ambient Jervis Bay temperatures. Similarly, juveniles (from both populations) reared under Adelaide conditions had similar thermal limits and swimming activity (maximum volitional velocity and distance) to each other, indicative of conserved acclimation capacity. However, under Adelaide temperatures, the ṀO_{2 Rest} of Jervis Bay sharks was greater than that of Adelaide sharks. This indicates that the energetics of cooler water population (Adelaide) is likely more thermally sensitive than that of the warmer population (Jervis Bay). While unique to elasmobranchs, these data provide further support that by treating species as static, homogeneous populations, we ignore the impacts of thermal history and intraspecific variation on thermal sensitivity. With climate change, intraspecific variation will manifest as populations move, demographics change or extirpations occur, starting with the most sensitive populations.

Thermal performance curves reveal shifts in optima, limits and breadth in early life

Understanding thermal performance at life stages that limit persistence is necessary to predict responses to climate change, especially for ectotherms whose fitness (survival and reproduction) depends on environmental temperature. Ectotherms often undergo stage-specific changes in size, complexity and duration that are predicted to modify thermal performance. Yet performance is mostly explored for adults, while performance at earlier stages that typically limit persistence remains poorly understood. Here, we experimentally isolate thermal performance curves at fertilization, embryo development and larval development stages in an aquatic ectotherm whose early planktonic stages (gametes, embryos and larvae) govern adult abundances and dynamics. Unlike previous studies based on short-term exposures, responses with unclear links to fitness or proxies in lieu of explicit curve descriptors (thermal optima, limits and breadth), we measured performance as successful completion of each stage after exposure throughout, and at temperatures that explicitly capture curve descriptors at all stages. Formal comparisons of descriptors using a combination of generalized linear mixed modelling and parametric bootstrapping reveal important differences among life stages. Thermal performance differs significantly from fertilization to embryo development (with thermal optimum declining by ∼2°C, thermal limits shifting inwards by ∼8-10°C and thermal breadth narrowing by ∼10°C), while performance declines independently of temperature thereafter. Our comparisons show that thermal performance at one life stage can misrepresent performance at others, and point to gains in complexity during embryogenesis, rather than subsequent gains in size or duration of exposure, as a key driver of thermal sensitivity in early life.

The power struggle: assessing interacting global change stressors via experimental studies on sharks

Ocean warming and acidification act concurrently on marine ectotherms with the potential for detrimental, synergistic effects; yet, effects of these stressors remain understudied in large predatory fishes, including sharks. We tested for behavioural and physiological responses of blacktip reef shark (Carcharhinus melanopterus) neonates to climate change relevant changes in temperature (28 and 31 °C) and carbon dioxide partial pressures (pCO₂; 650 and 1050 µatm) using a fully factorial design. Behavioural assays (lateralisation, activity level) were conducted upon 7–13 days of acclimation, and physiological assays (hypoxia tolerance, oxygen uptake rates, acid–base and haematological status) were conducted upon 14–17 days of acclimation. Haematocrit was higher in sharks acclimated to 31 °C than to 28 °C. Significant treatment effects were also detected for blood lactate and minimum oxygen uptake rate; although, these observations were not supported by adequate statistical power. Inter-individual variability was considerable for all measured traits, except for haematocrit. Moving forward, studies on similarly ‘hard-to-study’ species may account for large inter-individual variability by increasing replication, testing larger, yet ecologically relevant, differences in temperature and pCO₂, and reducing measurement error. Robust experimental studies on elasmobranchs are critical to meaningfully assess the threat of global change stressors in these data-deficient species.

Thermal tolerance and hypoxia tolerance are associated in blacktip reef shark (Carcharhinus melanopterus) neonates

Thermal dependence of growth and metabolism can influence thermal preference and tolerance in marine ectotherms, including threatened and data-deficient species. Here, we quantified the thermal dependence of physiological performance in neonates of a tropical shark species (blacktip reef shark, Carcharhinus melanopterus) from shallow, nearshore habitats. We measured minimum and maximum oxygen uptake rates (ṀO2), calculated aerobic scope, excess post-exercise oxygen consumption and recovery from exercise, and measured critical thermal maxima (CTmax), thermal safety margins, hypoxia tolerance, specific growth rates, body condition and food conversion efficiencies at two ecologically relevant acclimation temperatures (28 and 31°C). Owing to high post-exercise mortality, a third acclimation temperature (33°C) was not investigated further. Acclimation temperature did not affect ṀO2 or growth, but CTmax and hypoxia tolerance were greatest at 31°C and positively associated. We also quantified in vitro temperature (25, 30 and 35°C) and pH effects on haemoglobin–oxygen (Hb–O2) affinity of wild-caught, non-acclimated sharks. As expected, Hb–O2 affinity decreased with increasing temperatures, but pH effects observed at 30°C were absent at 25 and 35°C. Finally, we logged body temperatures of free-ranging sharks and determined that C. melanopterus neonates avoided 31°C in situ. We conclude that C. melanopterus neonates demonstrate minimal thermal dependence of whole-organism physiological performance across a seasonal temperature range and may use behaviour to avoid unfavourable environmental temperatures. The association between thermal tolerance and hypoxia tolerance suggests a common mechanism warranting further investigation. Future research should explore the consequences of ocean warming, especially in nearshore, tropical species.

Divergent field metabolic rates highlight the challenges of increasing temperatures and energy limitation in aquatic ectotherms

Environments where extreme temperatures and low productivity occur introduce energetically challenging circumstances that may be exacerbated by climate change. Despite the strong link between metabolism and temperature in ectotherms, there is a paucity of data regarding how the metabolic ecology of species affects growth and fitness under such circumstances. Here, we integrated data describing field metabolic rates and body condition of two sympatric species of ectotherms with divergent lifestyles, the benthic freshwater (or largetooth) sawfish (Pristis pristis) and the epipelagic bull shark (Carcharhinus leucas) occurring in the Fitzroy River, Western Australia, to test the implications of their differing metabolic ecologies for vulnerability to rising temperatures. Over a temperature range of 18-34 °C, sawfish had lower field metabolic rates (63-187 mg O₂ kg^-0.86 h^-1) and lower temperature sensitivity of metabolic rates [activation energy (E_A) = 0.35 eV] than bull sharks (187-506 mg O₂ kg^-0.86 h^-1; E_A = 0.48 eV). Both species lost body mass throughout the dry season, although bull sharks significantly more (0.17% mass loss day^-1) than sawfish (0.07% mass loss day^-1). Subsequent bioenergetics modelling showed that under future climate change scenarios, both species would reach potentially lethal levels of mass loss during dry season periods before the end of the century. These results suggest that ectotherms with low metabolic rates may be better suited to extreme environmental conditions, and that even small increases in temperature due to climate change could have substantial impacts on the ability of ectotherms to grow and survive in harsh conditions, including high temperatures and energy-limiting circumstances.

Body temperature stability observed in the whale sharks, the world's largest fish

It is generally assumed that the body temperature of large animals is less likely to change due to their large body size, resulting in a high thermal inertia and a smaller surface area to volume ratio. The goal of this study was to investigate the stability of body temperature in large fish using data from field experiments. We measured the muscle temperatures of free-ranging whale sharks (Rhincodon typus), the largest extant fish globally, and investigated their ectothermic physiology and the stability of their body temperatures. The measured muscle temperature of the whale sharks changed substantially more slowly than the water temperature fluctuations associated with vertical movements, and the whole-body heat-transfer coefficients (HTC) of whale sharks estimated using heat-budget models were lower than those of any other fish species measured to date. The heat-budget models also showed that internal heat production does not contribute to changes in muscle temperature. A comparative analysis showed that the HTC at cooling in various fish species including both ectothermic and endothermic species ranging from 10-4 to 103 kg was proportional to body mass-0.63. This allometry was present regardless of whether the fish were ectothermic or endothermic, and was an extension of the relationship observed in previous studies on small fish. Thus, large fish have the advantage of body temperature stability while moving in environments with large temperature variations. Our results suggest that the large body size of whale sharks aids in preventing a decrease in body temperature during deep excursions to more than 1000 m depths without high metabolic costs of producing heat.