ORIGINAL RESEARCH article
Sec. Aquatic Physiology
A Baseline Model For Estimating the Risk of Gas Embolism in Sea Turtles During Routine Dives
- 1Department of Research, Fundación Oceanogràfic de la Comunidad Valenciana, Valencia, Spain
- 2National Oceanic and Atmospheric Administration, National Marine Fisheries Service, Office of Protected Resources, University of Florida (duty station), Washington, DC, United States
- 3Virginia Aquarium and Marine Science Center, Virginia Beach, VA, United States
- 4Instituto de Biología de Organismos Marinos (IBIOMAR-CCT CONICET-CENPAT), Puerto Madryn, Argentina
- 5The Leatherback Trust, Fort Wayne, IN, United States
- 6Northeast Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration, Woods Hole, MA, United States
- 7Department of Clinical Sciences and Center for Marine Sciences and Technology, North Carolina State University, Raleigh, NC, United States
- 8Coonamessett Farm Foundation, East Falmouth, MA, United States
- 9Department of Comparative, Diagnostic, and Population Medicine, College of Veterinary Medicine, University of Florida, Gainesville, FL, United States
- 10Global Diving Research, Inc., Ottawa, ON, Canada
Sea turtles, like other air-breathing diving vertebrates, commonly experience significant gas embolism (GE) when incidentally caught at depth in fishing gear and brought to the surface. To better understand why sea turtles develop GE, we built a mathematical model to estimate partial pressures of N2 (PN2), O2 (PO2), and CO2 (PCO2) in the major body-compartments of diving loggerheads (Caretta caretta), leatherbacks (Dermochelys coriacea), and green turtles (Chelonia mydas). This model was adapted from a published model for estimating gas dynamics in marine mammals and penguins. To parameterize the sea turtle model, we used values gleaned from previously published literature and 22 necropsies. Next, we applied this model to data collected from free-roaming individuals of the three study species. Finally, we varied body-condition and cardiac output within the model to see how these factors affected the risk of GE. Our model suggests that cardiac output likely plays a significant role in the modulation of GE, especially in the deeper diving leatherback turtles. This baseline model also indicates that even during routine diving behavior, sea turtles are at high risk of GE. This likely means that turtles have additional behavioral, anatomical, and/or physiologic adaptions that serve to reduce the probability of GE but were not incorporated in this model. Identifying these adaptations and incorporating them into future iterations of this model will further reveal the factors driving GE in sea turtles.
Until the early 2000s, it was assumed that most air-breathing marine vertebrates were not susceptible to gas embolism (GE) or its associated diseases (Jepson et al., 2003). However, several studies have now provided direct evidence of GE in cetaceans, pinnipeds, and sea turtles (Van Bonn et al., 2011, 2013; Dennison et al., 2012a,b; Crespo-Picazo et al., 2020). Under natural conditions, most air-breathing marine vertebrates likely manage the risk of GE through a combination of behavioral, anatomical, and physiological adaptations (García-Párraga et al., 2018b; Fahlman et al., 2020, 2021). Yet there appears to be significant increase in the potential for GE if animals make abnormally rapid ascents after diving (Fahlman et al., 2021). This could be caused by animals being caught incidentally in fishing gear (Fahlman et al., 2017) or if they are disturbed by acoustic sources such as mid-frequency sonar (Tal et al., 2015). Severe cases of GE can even lead to immediate or delayed mortality after several days (Parga et al., 2020). Consequently, GE likely poses an underestimated threat to marine animals that are subject to high levels of fisheries bycatch such as sea turtles (Wallace et al., 2010, 2013).
The formation of gas emboli is driven by changes in gas dynamics associated with the variable hydrostatic pressure that an animal experiences while diving. During the descent, hydrostatic pressure increases in proportion to depth. For air-breathing vertebrates, the increased pulmonary gas tensions cause the gaseous N2, O2, and CO2 in the lungs to dissolve into solution within the pulmonary circulatory system (Fahlman et al., 2006, 2021). The gases are then transported via the vascular system throughout the various tissues of the body (Fahlman et al., 2006, 2009). The amount of gas that dissolves into the tissues during a dive is, therefore, influenced by the level of gas exchange in the lungs, cardiac output (heart rate and stroke volume), blood flow distribution, and tissue perfusion (Fahlman et al., 2006, 2009, 2018). Initially, the rate at which gas diffuses across the lungs increases with depth as higher hydrostatic pressures increase the quantity of gas forced into solution in the blood (Scholander, 1940; Berkson, 1967). However, if hydrostatic pressures continue to increase, the lung alveoli (or ediculi and faveoli in reptiles) compress and eventually collapse (Bostrom et al., 2008; Fahlman et al., 2009; Mcdonald and Ponganis, 2012). When all the alveoli have collapsed, gas exchange in the lungs ceases. Many diving animals also further reduce the rate at which gases dissolve in the their tissues through changes in circulation and perfusion (Fahlman et al., 2006). Specifically, some animals reduce cardiac output and blood flow to the lungs while simultaneously redistributing blood flow to O2-demanding tissues (Scholander, 1940).
On ascent, diminishing hydrostatic pressures decrease gas tensions and so the gases come out of solution. If this process occurs faster than the animal can transfer the excess blood gas to the lungs, intravascular gas bubbles will form (Fahlman, 2017; Fahlman et al., 2021). While minor GE can be tolerated, serious cases can damage neural and more vascularized tissues. In turn, this can lead to impaired motor skills, loss of consciousness, paralysis, and even death (van Hulst et al., 2003; Dennison et al., 2012b; García-Párraga et al., 2014).
Several recent studies have confirmed that GE in sea turtles is a wide-spread phenomenon associated with fisheries bycatch (García-Párraga et al., 2014; Crespo-Picazo et al., 2020; Parga et al., 2020) and can also occur in other interactions with submerged gear such as Hopper dredges (Harms et al., 2020). It has also been shown that the depth of capture is a key component in determining the likelihood of GE. Loggerhead turtles (Caretta caretta) caught by Mediterranean trawlers and gillnets at depths exceeding 65 m were 50% more likely to suffer from fatal GE than those caught at shallower depths (Fahlman et al., 2017). To understand why turtles are susceptible to fisheries-induced GE, we must first understand tissue gas dynamics during normal diving behavior. For example, it is possible that sea turtles live with elevated tissue N2 levels, comparable to that of a saturation diver (Brubakk et al., 2011). If so, even minor changes to routine diving patterns or disruption to those anatomical or physiological adaptations governing cardiac output and blood flow could lead to GE.
Our goal in this study was to develop a baseline model to estimate the risk of GE formation in sea turtles. To build this model, we built upon a published model for estimating gas tissue N2, O2, and CO2 tensions of marine mammals and penguins (Fahlman et al., 2006, 2007, 2009, 2018, 2021) and adapted it for use with loggerhead turtles, green turtles (Chelonia mydas), and leatherback turtles (Dermochelys coriacea). While there are many anatomical and physiological differences between marine mammals/penguins and sea turtles, building off these earlier models provided us with an initial framework for modeling gas dynamics that had been previously validated for use with air-breathing vertebrates. We parameterized the model based on estimates of key tissue compartment volumes in sea turtles derived from necropsies and computed tomography. We also gleaned other parameters on sea turtle metabolic rate and cardiac output data from the published literature. Next, we used this model to estimate partial pressures of N2 (PN2), O2 (PO2), and CO2 (PCO2) in various body compartments during routine dive patterns for each turtle species. Finally, we used this model to investigate whether (1) the percentage of body-fat (as a proxy for body-condition) and (2) cardiac output influenced the risk of GE. We chose to investigate body-fat as it changes the dynamics of the body N2 stores and has been considered as a potential risk variable in GE in humans (Lam and Yau, 1989; Schellart et al., 2012) and cetaceans (Fahlman et al., 2021). In addition, the volume of body fat is strongly influenced by reproduction in sea turtles (Davenport et al., 2011). We chose to investigate cardiac output as modeling efforts for marine mammals have shown this to be a key component influencing the risk of GE (Fahlman et al., 2006, 2009, 2021; Hooker et al., 2009).
Materials and Methods
We analyzed time-depth recorder (TDR) datasets from three loggerhead, six leatherback, and four green turtles (Table 1). The TDRs deployed on loggerhead turtles were Satellite Relay Data Loggers (SRDL–Sea Mammal Research Unit, United Kingdom), while MK10 tags (Wildlife Computers, United States) were used for both leatherback and green turtles. The loggerhead TDRs were attached to the top of the carapace using epoxy (for details, see Patel et al., 2018). For leatherback and green turtles, the TDRs were tethered to the posterior position of the carapace (for details, see Robinson et al., 2016). The loggerhead turtles were sub-adults sampled at foraging areas in the Mid-Atlantic Bight, while the leatherback and green turtles were inter-nesting females sampled on the east coast of South Africa and the northwest coast of Costa Rica, respectively (Table 1). The TDRs recorded depth to within 1 m accuracy at time-intervals of either 4 or 10 s. We standardized the sampling interval for all TDRs by interpolating the data to 1 s intervals.
Table 1. Time-depth data used in gas dynamics model. Data were analyzed from three loggerhead, six leatherback, and four green turtles.
TDR datasets ranged in duration 187–476 days for loggerhead turtles, 10–34 days for leatherback turtles, and 11–54 days for green turtles. As loggerheads are known to modify their diving behavior in response to seasonal variation in water temperatures (Iverson et al., 2019) and conduct extended “hibernation” dives during winter (Hochscheid et al., 2005), we only analyzed loggerhead data from the first 46 to 139 days (Table 1). This ensured that all data were collected during the summer months (May to November) when tag-derived sea surface temperature values were between 19 and 26°C. In contrast, tracking duration was short enough for leatherback and green turtles to assume that there were no significant changes in diving behavior due to seasonal factors within the available datasets. Moreover, all the leatherback and green turtles’ data were collected from tropical or sub-tropical habitats, where the sea surface temperature values only varied between 24 and 29°C. The loggerhead TDR data have been previously published in Patel et al. (2018). The leatherback and green turtle data have been published in Blanco et al. (2013), Robinson et al. (2017), and Clyde-Brockway et al. (2019).
To estimate gas dynamics in diving sea turtles, we adapted a model that was developed for estimating blood and gas tissue N2, O2, and CO2 tensions in marine mammals and penguins (Fahlman et al., 2006, 2007, 2009, 2018, 2021). In the model, gas exchange first occurs between the respiratory system and arterial blood, and then between the arterial blood and four compartments: (1) the brain; (2) fat and bone; (3) the central circulatory system, which included the heart, kidney, liver, and digestive tract but not blood; and (4) muscle, which included muscle, skin, connective tissue, and all other tissues and organs that were not included in the other compartments. Gas exchange next occurs between the four compartments and the mixed venous blood, and finally between the mixed venous blood and the respiratory system.
Gas exchange was driven by partial pressure/tension gradients as detailed in Fahlman et al. (2006), and a physiologic pulmonary shunt was allowed at increasing pressures due to alveolar compression (Bostrom et al., 2008; Fahlman et al., 2009). While reptilian lungs have faveoli or ediculi instead of the alveoli in mammals, the similarities in their structures mean that they likely respond similarly under pressure differentials. For ease of use, we hereafter refer to these differing features under the singular term faveoli. It should also be noted that sea turtles have the capacity to perform a right-to-left intracardiac shunt. As the role of this cardiac shunt in regulating N2 is unknown in aquatic reptiles (Burggren et al., 2020), we therefore adjusted the parameters for the modeled pulmonary shunt, so that it mirrored the arterial PN2 levels that previously reported for a sea turtles during forced submergence (Berkson, 1967).
The model assumed that breathing began the instant the animal reached the surface (i.e., a depth of 0 m). However, due to a combination of measurement error, the use of tethered tags on both leatherback and green turtles, and the original sampling frequency of the TDRs, not all breathing events were associated with a depth measurement of 0. While many studies address this issue by defining a dive as the period when a turtle descends beyond a pre-determined depth (e.g., 2 m), we chose to not use this dive definition as previous models have shown the risk of GE increases rapidly as a diver approaches the surface (Fahlman et al., 2018). Instead, we subtracted 2 m from all depth values. While this ensured that all breathing events were associated with a measurement of 0 depth, we acknowledge that this also slightly reduced the maximum depth of the dives.
When breathing, the arterial blood gas tensions of N2, O2, and CO2 were set to 0.741 atmospheres absolute (ATA), 0.164 ATA, and 0.033 ATA, respectively, as have been previously reported for loggerhead turtles (Lutcavage and Lutz, 1991) with water vapor being 0.062 ATA (Lutcavage and Lutz, 1991; Fahlman et al., 2008). We also assumed that arterial blood gas tensions while breathing were equal to the partial pressures in the faveoli, and that all CO2 that exchanged for O2 remained in gas phase in the lungs and did not dissolve into the lung parenchyma (Fahlman et al., 2006).
Lung Compression and Pulmonary Shunt
We used the model published by Bostrom et al. (2008) to estimate faveolar volume at depth and compared the estimated faveolar collapse depth from this model with empirical data previously published in green turtles forcibly submerged inside hyperbaric chambers (Berkson, 1967). This model assumed that the total lung capacity (TLC) includes both the maximum faveolar volume (VfAmax) and the volume of dead space in the trachea and bronchi (VD; Fahlman et al., 2009, 2018). We obtained measurements of TLC for each species from the relevant literature (Table 2; Berkson, 1967; Table 3; Lapennas and Lutz, 1982; Lutz and Bentley, 1985; Lutcavage et al., 1990; Lutcavage and Lutz, 1991; Hochscheid et al., 2007). However, as we were only able to find measurements of VD for green turtles (Berkson, 1967; Gatz et al., 1987), we assumed that VD was equivalent (7% of TLC) in both loggerhead and leatherback turtles. We assumed that gas exchange only occurred in the faveoli, and that gas exchange stopped immediately upon faveolar collapse (i.e., when the faveolar volume was zero, VfA = 0).
Table 2. Estimated model variables for mass-specific total lung capacity (TLC, ml • kg−1), O2 stores (ml • kg−1) of the lung (l-O2), blood (b-O2), and muscle (m-O2), hemoglobin concentration ([Hb], g • dl−1), packed cell volume (PCV, %), myoglobin concentration ([Mb], g Mb • kg−1 muscle).
Table 3. Body compartment size as a percentage of body mass in loggerhead, green, and leatherback turtles combined.
To account for the relationship between pulmonary shunt and faveolar collapse, we used Eq. 4 in Fahlman et al. (2009). We also estimated the faveolar volume (Eq. 4 in Bostrom et al., 2008), using the parameters as previously defined for human alveoli (a = 1.04, b = 0.20, and c = 1.21). Finally, we estimated dead space in the upper airways (Eq. 5 in Bostrom et al., 2008) and applied the parameters as previously defined for a compliant human trachea (Kp = −6.44, n = 0.74). We used the same parameters for all three turtle species.
We generated estimates of the relative size of the four compartments in the model using computed tomography of five loggerhead turtles ranging in body mass from 3.2 to 48.7 kg (Appendix 1). To support these results, we also performed mass dissections of two loggerheads, one leatherback, and nine green turtles ranging in body mass from 1 to 355.9 kg and encompassing both sexes. Dissections were conducted by the Fundación Oceanogràfic de la Comunidad Valenciana which is registered as a research unit under the official ID number: ES460250001024 under the collaboration agreement with the Conselleria d’Agricultura, Medi Ambient, Canvi Climàtic i Desenvolupament Rural of the Valencian Regional Goverment, the National Oceanic and Atmospheric Administration (Florida Fish and Wildlife Conservation Commission Marine Turtle Permit 20–081), and North Carolina State University (North Carolina Endangered Species Permit 20ST42). Organs were removed intact and as much blood as possible drained from the tissues prior to weighing. To weigh the heart, the great vessels were severed at the pericardium. Due to the difficulties separating fat and muscle, we did not attempt to estimate the fat/bone and muscle compartments during dissections.
To estimate blood volumes, we used a value midway between those measured in leatherback and green turtles (Thorson, 1968; Lutcavage et al., 1992; Table 3). Blood was separated into arterial and venous, respectively, 33 and 67% of the total blood volume (Fahlman et al., 2006). After the arterial blood gases were exchanged with the tissue, the venous side from all compartments joined immediately to form mixed venous blood, as previously detailed (Fahlman et al., 2006, 2009).
As we observed higher variability between individuals than between species and proportions in compartment sizes among species were similar, we combined the results from all the species to characterize compartment sizes for a “generic” sea turtle (Table 3). This turtle was composed of 7.00% blood, 0.06% brain, 30.00% fat/bone, 9.00% central circulation, and 53.94% muscle.
Cardiac Output and Blood Flow Distribution
We defined cardiac output () at the surface (i.e., depth = 0) based on measurements from green turtles (body mass = 1.2 kg; West et al., 1992). To account for allometric differences between species in body mass, we used the following equation:
where 8.4 was the (l • min−1) and sis the mass-specific (Davis and Kanatous, 1999). At the surface, the model transferred 60% of to the central circulation compartment, 34% to the muscle compartment, 4% to the brain compartment, and 2% to the fat and bone compartment. As no empirical data on this exist for sea turtles, we chose these values based on comparable data from marine mammal studies (Fahlman et al., 2018). While diving, the distribution of blood flow was tested iteratively to optimize O2 utilization and thus maximize the aerobic dive duration (Fahlman et al., 2006, 2009). This resulted in 92.8, 5.0, 2.0, and 0.2% of s being directed to the central circulation, muscle, brain, and fat and bone, respectively, while diving.
When at the surface, we used a metabolic rate of 4.23 ml O2 • min−1 • kg−1 for each compartment. This is comparable to those reported in the literature for exercising leatherback, loggerhead, and green turtles on land under air temperatures between 26 and 30°C (Jackson and Prange, 1979; Butler et al., 1984; Lutcavage et al., 1987, 1990, 1992; West et al., 1992). We used values for exercising MR instead of resting MR, as cardiac output is often higher at the surface to recover from previous dives (Southwood et al., 1999). When diving, we assumed that the metabolic rate was 5% of that at the surface, i.e., 0.21 ml O2 • min−1 • kg−1 based on data from a leatherback turtle conducting a 34-min dive (Southwood et al., 1999). This value also lies within the range of metabolic rates reported by Jackson and Prange (1979), Butler et al. (1984), Lutcavage et al. (1987, 1990, 1992), and West et al. (1992), and it also ensured that the modeled turtles did not run out of O2 during long dives.
Hemoglobin and Myoglobin Characteristics
We used published data on blood packed cell volume (PCV, %), hemoglobin ([Hb]), and muscle myoglobin ([Mb], Table 2) concentrations for each species (Lapennas and Lutz, 1982; Lutz and Bentley, 1985; Lutcavage et al., 1990). A P50 of the O2 dissociation curve for each species was gleaned from published literature (Table 3). We used a Hill coefficient of 2.7 for all species (Berkson, 1967; Wood et al., 1984; Lutcavage et al., 1990).
Gas Dynamics Computations
We used O2, CO2, and N2 gas balance equations from Fahlman et al. (2006) to estimate gas transport in the lungs and tissues over short-time intervals and small blood volumes (package volume). For these equations, the partial pressure/tension of each gas was used as the driving force for gas exchange, and each compartment was modeled as a single homogenous tissue with a one-to-one perfusion rate and tissue gas solubility (Fahlman et al., 2006).
We ran several model variations on each of the turtle datasets to test the effects of varying: (1) body-condition and (2a) the surface sand (2b) diving s on venous PN2
(1) To assess how changes in body-fat affect venous PN2, we varied the relative volume of the fat/bone compartment in the control model (30%) to either 20 or 40% to represent extremes of nutritional condition. When changing the relative size of the fat/bone compartment, the contribution of the other compartments was adjusted accordingly to ensure that they maintained the same overall mass (Table 3)
(2a and b) To assess the effect of varying the s, we adjusted the surface sfrom the control model (sof 7 ml • min−1 • kg−1) to either 2.5 ml • min−1 • kg−1 (Surface-Low) or 10 ml • min−1 • kg−1 (Surface-High). We also examined the effect of varying the sduring diving from the control model (5.0% of the control s) to either 10.0% (dive-low) or 3.3% (dive-high).
Finally, to evaluate how changes in body-fat or cardiac output affected the risk of GE, we calculated the end dive PN2 saturation of the mixed venous blood during the last 5 s of each dive. We then measured the percentage change between the mean (over all dives) and the maximum end dive PN2 for the different models using the formula:
Model Validation and Analysis
While there was no published data available on arterial PN2 levels in free-diving sea turtles, Berkson (1967) measured arterial PN2 in green turtles during forced submergence. Thus, to assess the accuracy of our model for estimated arterial PN2, we compared the modeled results to those observed in Berkson (1967) during comparable dives.
As these are modeling based on the physiological assumptions from a limited number of individuals, we opted to only perform limited statistical analyses on the results. As such, these results should be considered largely descriptive and are simply to provide a framework for developing testable hypotheses in future studies.
Dive summaries for each turtle are reported in Table 4. We divided these summaries between shallow (0–30 m), medium (30–90 m), and deep (>90 m) dives. We selected these thresholds according to the values in Berkson (1967) as there should be limited impact of pressure on lung diffusion at 30 m or less and there should be almost no lung diffusion deeper than 90 m. The table includes the average (±SD) dive duration, the maximum depth reached during each dive, and the average depth for each dive (Table 4).
When comparing data on arterial PN2 during forced submergence from Berkson (1967) to the model output, the modeled data were never more than 1.5 ATA from the measured data (Figure 1). Furthermore, the modeled data strongly reflected the observed changes in arterial PN2 at different pressures.
Figure 1. Measured (red circles) and estimated (broken black line) arterial N2 tension (PN2) in a Pacific green turtle (Chelonia mydas agassizii) during a step-wise forced dive to 19.4 ATA (Berkson, 1967). The model used a surface sof 3.42 ml • min−1 • kg−1 and a diving s0.17 ml • min−1 • kg−1 and a diving lung volume 50% of total lung capacity.
The model output showed that the accumulation of N2 into the central circulation occurred faster than it did in the fat (Figure 2). Because of this, PN2 within the central circulation compartment tended to reflect the ambient pressure, while fat PN2 was more influenced by the dive profiles over a series of hours to days (Figure 3). As such, central circulation only briefly contributes to GE risk during the ascent as N2 was rapidly removed and the supersaturation decreased. The fat compartment, on the other hand, slowly varied with time and depended on the previous diving pattern. Furthermore, unlike the central circulation, supersaturation occured close to the surface and during the entire surface interval. Nevertheless, the overall fat PN2 saturation seldom exceeded 1 ATA and therefore does not contribute appreciably to GE risk.
Figure 2. Model result for control model, showing ambient pressure (Pamb, ATA), and tissue gas tension (PN2) for central circulation (CC), fat, and mixed venous blood (MV) in a single loggerhead (Turtle ID 14293), leatherback (Turtle ID 10A1059a), and green sea turtles (Turtle ID 37800, Table 2). Data are shown for 100, 10, and 55 days for loggerhead, leatherback, and green turtles, respectively. The Pamb was used to generate the model output for the various tissues and the difference in tissue PN2 between CC (fast tissue) and fat (slow tissue), with the PN2 in MV is the composite of all tissues and their relative blood flow.
Figure 3. Model result for control model, showing depth (m), and tissue gas tension (PN2) in central circulation (CC), fat, and mixed venous blood (MV) in loggerhead (Turtle ID 14293), leatherback (Turtle ID 10A1059a), and green sea turtles (Turtle ID 37800, Table 2). The data presented are a subset of those in Figure 2 to show specific dive patterns and blood and tissue PN2 levels. The theoretical supersaturation [(mixed venous PN2 – ambient pressure) • ambient pressure−1, in ATA] for each of the tissues and the mixed venous blood are shown in color. The color bar shows the level of supersaturation, where negative values (blue) indicate N2 uptake and positive values indicates the whole body, or mixed venous PN2, exceeds the ambient pressure, i.e., risk of Gas Emboli formation (GE) and pathology (GEP). Thus, increasing value, i.e., yellow to red, indicates increasing risk of GEP. The results show increasing risk of GEP during ascent and how ascent rate or how sudden stops during ascent reduces risk, e.g., ascent from second dive in loggerhead or leatherback turtles.
Using the control model and combining the results from all three species, mean end-dive PN2 increased with depth for central circulation (χ2 = 32, 2 df, p < 0.01), muscle (χ2 = 22, 2 df, p < 0.01), brain (χ2 = 20, 2 df, p < 0.01), and mixed venous (χ2 = 34, 2 df, p < 0.01) but not for the fat compartment (χ2 = 2.4, 2 df, p > 0.2). When comparing dives within the shallowest dive bin (<30 m) to the deepest dive bin (>90 m), mean end-dive PN2 increased by 245% for central circulation, 11% for muscle, 189% for brain, and 221% for the mixed venous compartment. When the three species were considered separately, leatherback turtles had the lowest mean end-dive PN2, with green turtles having intermediate values, and loggerhead turtles having the highest (χ2 = 8.49, 2 df, p = 0.014; End-dive PN2 = 2.73 (±0.31)−0.56 (±0.28) × leatherback + 0.54 (±0.41) × loggerhead−0.63 (±0.30) medium−1.71 (±0.29) shallow).
With data from all three species combined, maximum end-dive PN2 increased with depth for central circulation (χ2 = 11.5, 2 df, p < 0.01), brain (χ2 = 8.8, 2 df, p = 0.012), and mixed venous blood (χ2 = 12.1, 2 df, p < 0.01), but decreased for fat (χ2 = 7.9, 2 df, p = 0.019), and did not change for muscle (χ2 = 5.1, 2 df, p = 0.079). Comparing dives in the shallow depth range to those in the deep depth bin, maximum end-dive PN2 increased by 80% for central circulation, by 68% for brain, and by 81% for mixed venous blood. For brain, both depth range and species helped to explain the variation (χ2 = 7.39, 2 df, p = 0.025). Overall, mean end-dive PN2 for brain increased with depth and leatherback turtles had the lowest maximum end-dive PN2, while green turtles had intermediate values, and loggerhead turtles had the highest (change in end-dive brain PN2 (%) = 4.23 (±0.65) − 0.50 (±0.66) × leatherback + 1.95 (±0.92) × loggerhead – 1.03 (±0.54) medium −1.86 (±0.53)). The maximum end-dive PN2 for mixed venous blood was 9.40 ATA, 6.17 ATA, and 5.03 ATA for leatherback, green, and loggerhead turtles, respectively, representing supersaturation ratios [(mixed venous PN2 – ambient pressure) • ambient pressure−1, or M-ratios] of 12.7, 8.3, and 6.8.
Model Variation – Body-Fat
When body fat was increased from 30 to 40% of body mass or decreased from 30 to 20% of body mass, there was no effect observed on mean or maximum end-dive mixed venous PN2 for either shallow or intermediate depths when all three species were combined (Figure 5). However, for the deepest depth bin, increasing or decreasing body-fat reduced the mean end-dive mixed venous PN2 (χ2 > 10.0, 2 df, p < 0.01; Figures 4A,B) and the maximum end-dive mixed venous PN2 (Figures 5A,B). Thus, a model that included depth, species, and body-fat was warranted; % change in maximal end-dive mixed venous PN2 = −13.3 (±1.2) − 12.4 (±2.6) leatherback – 5.2 (±3.8) logger + 12.5 (±3.0) medium depth + 11.0 (±2.9) shallow depth + 8.6 (±2.3) emaciated (χ2 > 14.4, 1 df, p < 0.01).
Figure 4. Data showing percent changes in mean end-dive PN2 for variation in gas dynamics model as compared with the control model [(new model – control model) • control model–1 • 100].
Figure 5. Data showing percent changes in maximal end-dive PN2 for variation in gas dynamics model as compared with the control model [(new model – control model) • control model–1 • 100].
Model Variation – Cardiac Output at Surface
When swas either increased or decreased, there were no consistent differences in mean mixed venous end-dive PN2 among species at the shallow depth bin (p > 0.3 for all, Figures 4C,D). For both medium and deep depths, when swas increased to 10 ml • min−1 • kg−1 (Surface-High) this increased mean end-dive PN2 by 10.6% (Figure 4D), while a decrease in sto 2.5 ml • min−1 • kg−1 (Surface-Low; Figure 4C) caused a 4.5% decrease in mean end-dive PN2 [χ2 > 26.5, 1 df, p < 0.01, change in end-dive PN2 (%) = 10.6 (±1.9) – 15.1 (±2.7) × (Surface-Low)].
For maximum end-dive mixed venous PN2, the changes were similar, with 43% increase in surface s(Surface-High) causing an increase in maximum end-dive PN2 by 13.4% (Figure 5D), and 64% decrease (Surface-Low) in scausing a 5.4% decrease [Figure 5C, χ2 > 31.0, 1 df, p < 0.01, change in end-dive PN2 (%) = 13.4(±2.1) – 18.7 (±3.0) × (Surface-Low)].
Model Variation – Cardiac Output at Depth
When altering diving s, there was no statistically significant effect of dive depth on mean end-dive mixed venous PN2 when all three species were combined (p > 0.5 for all, Figures 4E,F). Considering all species and depths together, an increase in diving s from 5 to 10% lead to 20.6% increase in mean end-dive mixed venous PN2 (Figure 5E). A decrease in diving sfrom 5 to 3.3% (Figure 4F) lead to a decrease of 13.2% in the mean end-dive mixed venous PN2 [χ2 > 52.8, 1 df, p < 0.01, change in end-dive PN2 (%) = 20.6(±2.7) – 33.9 (±3.8) × (1/30)].
When varying diving s, there was uniform effect on the maximal end-dive PN2, however, there were differences between species. Specifically, the maximum end-dive mixed venous PN2 was higher for leatherback turtles than either loggerhead or green turtles at all depths. In addition, the lower diving ssignificantly reduced the maximal end-dive PN2 (Figures 5E,F); [χ2 > 6.21, 1 df, p = 0.013, change in maximal end-dive PN2 (%) = 32.2 (±3.1) – 44.9 (±3.8) × (1/30) – 9.6 (±3.8) loggerhead/green].
GE associated with fisheries bycatch and other submerged gear interactions could be a significant threat to sea turtles worldwide (García-Párraga et al., 2014; Fahlman et al., 2017; Crespo-Picazo et al., 2020; Harms et al., 2020; Parga et al., 2020). To better understand how sea turtles manage the risk of GE during natural diving behavior, we adapted a theoretical model to estimate blood and tissue N2 tensions in diving marine mammals to fit loggerhead, leatherback, and green sea turtles. We also applied the new model to assess how body-fat and cardiac output influence the potential for GE during normal diving behavior in sea turtles. Demonstrating the accuracy of the model, it provided a reasonable estimate for blood and tissue PN2 from previously conducted forced diving experiments (Berkson, 1967). By varying the model, we observed that higher fat reserves can reduce end-dive PN2 by up to 25%. We also observed that an increase in surface s, as likely occurs as the turtle reaches the surface to help to increase O2 uptake and CO2 removal (Okuyama et al., 2020), led to an increase in end-dive PN2. Finally, the model showed that leatherback turtles have lower mean end-dive blood and tissue PN2, but higher maximum end-dive values, than loggerhead or green turtles during routine diving behavior. While higher maximum end-dives values may be expected considering that leatherback turtles dive considerably deeper than the other two species (Robinson and Paladino, 2015), the fact that they had lower mean end-dive blood and tissue PN2 also suggests that their increased body-size helps to reduce the risk of GE.
Model Assumptions and Limitations
Our goal was to construct the initial framework for a model to understand gas dynamics in diving sea turtles. As no such models previously existed, a logical first step was to adapt a previously constructed and validated model for another air-breathing vertebrate. While we chose to adapt a model originally designed for marine mammals (Fahlman et al., 2006, 2009, 2021), there are several key anatomical and physiological differences between marine mammals and sea turtles that have not yet been incorporated into the new model. Firstly, and of considerable importance, is that sea turtles, unlike marine mammals, are at least partially ectotherms. Specifically, even though several studies have shown that sea turtles can maintain body-temperature elevated above ambient conditions (Paladino et al., 1990; Sato et al., 1995), their body temperatures still fluctuate over wider ranges than marine mammals. While we attempted to control for this by only analyzing diving data for sea turtles within a restricted range of sea surface temperatures (19–29°C), the effect of temperature on both the metabolic rate and the gas tensions/solubilities cannot be ignored (Lutz et al., 1989).
Another important difference between marine mammals and sea turtles is that sea turtles have muscular sphincters within their pulmonary arteries and a partially compartmentalized ventricle that allow for central intracardiac shunting (García-Párraga et al., 2018a; Burggren et al., 2020). These features could allow a complete shunt and cessation of gas exchange as the animal begins the dive. In fact, this could be the mechanism by which sea turtles avoid GE (García-Párraga et al., 2018b; Fahlman et al., 2021). In this model, we assumed that the pulmonary shunt develops due to the passive compression of the terminal air spaces and the structural properties of the conducting airways and faveoli (Bostrom et al., 2008; Fahlman et al., 2009). While this is not necessarily the case in reptiles, the exact factors influencing the cardiac shunt in sea turtles are not well-established (Burggren et al., 2020). Without having a solid understanding of the functioning of the cardiac shunt, we adjusted the model so that the arterial PN2 reflected those of forced-diving sea turtles in Berkson (1967). This meant that the physiological shunt and the potential anatomical shunt were expressed together as a pulmonary physiological shunt. In addition, we acknowledge that forced-diving experiments in Berkson (1967) would likely elevate blood and tissue N2 levels beyond normal values, and our model likely thus provides a conservative estimate of blood and tissue N2 levels. However, if we consider exercise as one of the main driving forces in the cardiac output adjustment (Okuyama et al., 2020), physically restrained animals could minimize N2 uptake reflecting lower N2 values compared to free swimming or entrapped individuals with restricted movement.
Other significant assumptions with the model resulted from differences between the datasets available for each species. For example, we utilized diving data from subadult loggerhead in foraging habitats, while the data for leatherback and green turtles were from adult females in inter-nesting habitats. Indeed, it is known that foraging and inter-nesting turtles exhibit differing dive behaviors (Robinson and Paladino, 2015). Other factors including life-stage and temperature are also known to influence diving behavior but are once again not standardized between our turtle diving datasets.
As is the case for all theoretical physiological models, the complexity of the entire system cannot be incorporated holistically into a mathematical model. While such arguments are often used to disparage the value of such exercises, we stress that we aimed to build the foundation upon which more detailed assessments of how factors, such as behavior, life-stage, or life-stage, can influence the risk of GE in sea turtles. As such, the goal of this exercise was to try and to identify which factors and assumptions need to be prioritized in future models to gain a more accurate understanding of gas dynamics in sea turtles.
General Patterns in PN2 and Inter-Specific Variation
In the basic model, the mean end-dive PN2 differed between species and maximum depth of the previous dive. For shallow and deep dives, the end-dive PN2 values were, respectively, between 38–111% and 193–342% higher as compared with the surface equilibrium PN2. Interestingly, maximum PN2 increased with dive depth in all turtle species. This suggests that the risk of GE increases with depth as observed in previous studies (Fahlman et al., 2017).
Effect of Body-Fat on Gas Emboli Formation
While we observed no effect of body-fat on end-dive mixed venous PN2 at depths less than 90 m, significant differences were observed beyond 90 m. Interestingly, end-dive mixed venous PN2 was highest at intermediate body-fat levels, and it dropped if body-fat was either increased to 40% or dropped to 20%. Moreover, these differences were up to 25% lower in leatherback turtles relative to loggerhead or green turtles. These findings suggest that the higher levels of fat in leatherback turtles could perhaps reduce the effects of GE and play a role in the ability of leatherbacks to dive several times deeper than any other sea turtle species (Robinson and Paladino, 2015). That said it is unclear why lower body-fat values also result in reduced end dive PN2.
Effect of Cardiac Output on Gas Emboli Formation
Past studies have shown that variation in blood flow, through changes in cardiac output, and the level of gas exchange are the physiological variables that have the greatest effect on N2 uptake and removal (Fahlman et al., 2006, 2009; Hooker et al., 2009). Our model provided evidence that this may also be the case for sea turtles as a reduction of the diving s to 3.3% caused a mean 13.2% decrease in end-dive PN2, while increasing the diving s to 10% caused a mean 20.6% increase in end-dive PN2. Similar patterns were also observed for maximum end-dive PN2. Nevertheless, these patterns differed by species and leatherback turtles exhibited lower mean end-dive PN2 at both higher and lower diving s.
Implications for Gas Embolism Formation
Our model suggests that during “natural” diving behavior sea turtles experience blood and tissue N2 levels that would cause decompression sickness in land mammals of similar size. Indeed, the maximal end-dive PN2 values between species ranged from 5.03 ATA to 9.40 ATA, which is considerably higher than the level of end-dive PN2 that would result in severe DCS in 50% of similarly sized humans (Fahlman et al., 2020). As it is highly unlikely that sea turtles are perpetually suffering from GE during routine diving behavior, we propose that turtles, much like marine mammals, must have additional behavioral, anatomical, and/or physiological mechanisms to reduce N2 uptake that are not currently considered in this model (García-Párraga et al., 2018a; Burggren et al., 2020; Fahlman et al., 2021).
One mechanism that sea turtles may employ to lower the risk of GE is the Selective Gas Exchange hypothesis (García-Párraga et al., 2018b; Fahlman et al., 2021). This hypothesis, which is supported by theoretical, anatomical, and physiological studies (Fahlman et al., 2009, 2018, 2020; Olson et al., 2010; Hodanbosi et al., 2016; García-Párraga et al., 2018a), suggests that if ventilation and perfusion in the lung can be managed, it would allow exchange of O2 and CO2, with little or no exchange of N2 (García-Párraga et al., 2018b; Fahlman et al., 2021). In support of this, it has been shown that the pulmonary arterial sphincters contract and relax when exposed to parasympathetic and sympathetic neurotransmitters, respectively (García-Párraga et al., 2018a). If sea turtles exhibit Selective Gas Exchange, this would help to significantly reduce our modeled blood and tissue N2 levels and help to explain how turtles avoid GE during natural dives. It would also explain how the sympathetic stress response could result in elevated pulmonary blood flow to the still ventilated lung, yielding excessive N2 uptake and GE (Fahlman et al., 2021).
It is also possible that turtles couple physiological and behavioral mechanisms to further reduce the risk of GE. Indeed, sea turtles increase their heart rate, and likely cardiac output, during the ascent phase (Okuyama et al., 2020). In addition, several breath-hold diving vertebrates, including sea turtles, the beluga (Delphinapterus leucas; Martin and Smith, 1992), beaked whales (Hyperoodon ampullatus; Hooker and Baird, 1999), macaroni penguin (Eudyptes chrysolophus; Sato et al., 2004), and king penguins (Aptenodytes patagonicus; Sato et al., 2004) have been shown to reduce the ascent rate while approaching the surface. It has been suggested that this reduction may be an adaptation to reduce end-dive PN2 before reaching the surface (Fahlman et al., 2006; Fossette et al., 2010). When the reduction in ascent rate is coupled with increases in heart rate during the ascent, it may reduce end-dive PN2 by as much as 45% in marine mammals (Fahlman et al., 2006). A similar reduction in the supersaturation was seen in all the turtles in this study (Figure 4). It also appears that the tendency of sea turtles to conduct a short “pause” in their accent around 10–20 m, which we observed in all three species, may also help to reduce the supersaturation as the turtle returned to the surface (Fossette et al., 2010).
These adaptations to reduce PN2 during normal dives could also explain why sea turtles are predisposed to GE when caught by fisheries. For example, rapid accent to the surface in a fishing net coupled with disruption of the normal dive profile and physiological response to underwater entrapment may alter cardiac output and shunting, leading to rapid GE. It also suggests that slower retrieval of the fishing nets or trawls may help to significantly reduce blood and tissue N2 levels in bycatch sea turtles. If operationally possible, this may provide a simple and cost-effective strategy to minimize GE in turtle by-catch and could have the added benefit of reducing fish by-catch as well (Grothues et al., 2017). However, a significant concern with longer ascent times is increased risk of drowning, which may negate any benefits.
Our model provides a framework for future studies on the analysis of GE in sea turtles. While we adapted the model to test the effects of body-fat and cardiac output on GE, this model could be adapted in other beneficial ways. Indeed, incorporating a simulated capture event, whereby the speed of the ascent to the surface or the total time spent at depth is varied. The model could also be adapted to assess how variation in the cardiac shunt plays a role in GE during both routine dives and during incidental capture. The strong effect of cardiac output on GE also highlights that more information is needed on cardiac dynamics in diving sea turtles.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The animal study was reviewed and approved by Purdue Animal Care and Use Committee, NOAA Animal Care and Use Committee.
NR and AF designed the experiment with input from BS. NR and AF wrote the manuscript with input from all authors. AF built and ran the model. NR, GB, CC-B, HH, and SP provided the data from the data loggers. BS, DG-P, and AC supervised necropsies and computed tomography. All authors contributed to the article and approved the submitted version.
The model development was funded by the National Oceanic and Atmospheric Administration (NOAA). Loggerhead data logger deployments were funded the scallop industry Sea Scallop Research Set Aside program administered by the Northeast Fisheries Science Center and by the U.S. Department of the Interior, Bureau of Ocean Energy Management through Interagency Agreements M14PG00005, M10PG00075, and M19PG00007 with the U.S. Department of the Commerce, NOAA, National Marine Fisheries Service (NMFS), Northeast Fisheries Science Center (NEFSC). Leatherback and green turtle data logger deployments were funded by the Leatherback Trust. Nelson Mandela Metropolitan University provided essential support during the deployment of the leatherback data loggers. The staff at the Fundación Oceanogràfic de la Comunidad Valenciana, especially Jose Luis Crespo Picazo And Vicente Marco Cabedo; NOAA; and North Carolina State University provided essential support conducting necropsies.
Conflict of Interest
AF was employed without salary by the company Global Diving Research Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
We would like to thank the GIBI 2-30 group at the Instituto de Investigación la Fe de Valencia, especially Luis Marti and Amadeo Ten, for their support on CT imaging post-processing for segmentation and quantification of the different sea turtle body compartments.
ATA, Atmospheres Absolute; DVfA, Faveolar Volume; GE, Gas Embolism; Mb, Myoglobin; Pamb, Ambient Partial Pressure; PN2, Partial Pressure of N2; PO2, Partial Pressure of O2; PCO2, Partial Pressure of CO2; , Cardiac Output; TDR, Time-Depth Recorder; TLC, Total Lung Capacity; VD, Volume of Dead Space in the Trachea and Bronchi; VfAmax, Maximum Faveolar Volume; s, Mass-Specific Cardiac Output.
Blanco, G. S., Morreale, S. J., Seminoff, J. A., Paladino, F. V., Piedra, R., and Spotila, J. R. (2013). Movements and diving behavior of internesting green turtles along Pacific Costa Rica. Integr Zool 8, 293–306. doi: 10.1111/j.1749-4877.2012.00298.x
Bostrom, B. L., Fahlman, A., and Jones, D. R. (2008). Tracheal compression delays alveolar collapse during deep diving in marine mammals. Respir. Physiol. Neurobiol. 161, 298–305. doi: 10.1016/j.resp.2008.03.003
Butler, P. J., Milsom, W. K., and Woakes, A. J. (1984). Respiratory, cardiovascular and metabolic adjustments during steady state swimming in the green turtle, Chelonia mydas. J. Comp. Physiol. B. 154, 167–174. doi: 10.1007/BF00684141
Clyde-Brockway, C., Robinson, N. J., Blanco, G. S., Morreale, S. J., Spotila, J. R., Santidrián, T., et al. (2019). Local variation in the Internesting behavior of green turtles in the Gulf of Papagayo, Costa Rica. Chelonian Conserv. Biol. 18, 217–230. doi: 10.2744/CCB-1352.1
Crespo-Picazo, J. L., Parga, M., Bernaldo De Quirós, Y., Monteiro, D., Marco-Cabedo, V., Llopis-Belenguer, C., et al. (2020). Novel insights Into gas embolism in sea turtles: first description in three new species. Front. Mar. Sci. 7:442. doi: 10.3389/fmars.2020.00442
Davenport, J., Plot, V., Georges, J.-Y., Doyle, T. K., and James, M. C. (2011). Pleated turtle escapes the box – shape changes in Dermochelys coriacea. J. Exp. Biol. 214, 3474–3479. doi: 10.1242/jeb.057182
Dennison, S., Fahlman, A., and Moore, M. J. (2012a). The use of diagnostic imaging for identifying abnormal gas accumulations in cetaceans and pinnipeds. Front. Physiol. 3:00181. doi: 10.3389/fphys.2012.00181
Fahlman, A., Cozzi, B., Manley, M., Jabas, S., Malik, M., Blawas, A., et al. (2020). Conditioned variation in heart rate during static breath-holds in the bottlenose dolphin (Tursiops truncatus). Front. Physiol. 11:604018. doi: 10.3389/fphys.2020.604018
Fahlman, A., Crespo Picazo, J.-L., Sterba-Boatwright, B., Stacy, B. A., and Garcia-Parraga, D. (2017). Defining risk variables causing gas embolism in loggerhead sea turtles (Caretta caretta) caught in trawls and gillnets. Sci. Rep. 7:2739. doi: 10.1038/s41598-017-02819-5
Fahlman, A., Hastie, G. D., Rosen, D. A. S., Naito, Y., and Trites, A. W. (2008). Buoyancy does not affect diving metabolism during shallow dives in Steller Sea lions Eumetopias jubatus. Aquat. Biol. 3, 147–154. doi: 10.3354/ab00074
Fahlman, A., Hooker, S. K., Olszowka, A., Bostrom, B. L., and Jones, D. R. (2009). Estimating the effect of lung collapse and pulmonary shunt on gas exchange during breath-hold diving: the Scholander and Kooyman legacy. Respir. Physiol. Neurobiol. 165, 28–39. doi: 10.1016/j.resp.2008.09.013
Fahlman, A., Jensen, F., Tyack, P. L., and Wells, R. (2018). Modeling tissue and blood gas kinetics in coastal and offshore common bottlenose dolphins, Tursiops truncatus. Front. Physiol. 9:00838. doi: 10.3389/fphys.2018.00838
Fahlman, A., Moore, M. J., and Wells, R. S. (2021). How do Marine mammals manage and usually avoid gas emboli formation and gas embolic pathology? Critical clues From studies of wild dolphins. Front. Mar. Sci. 8:25. doi: 10.3389/fmars.2021.598633
Fahlman, A., Olszowka, A., Bostrom, B., and Jones, D. R. (2006). Deep diving mammals: dive behavior and circulatory adjustments contribute to bends avoidance. Respir. Physiol. Neurobiol. 153, 66–77. doi: 10.1016/j.resp.2005.09.014
Fossette, S., Gleiss, A. C., Myers, A. E., Garner, S., Liebsch, N., Whitney, N. M., et al. (2010). Behaviour and buoyancy regulation in the deepest-diving reptile: the leatherback turtle. J. Exp. Biol. 213, 4074–4083. doi: 10.1242/jeb.048207
García-Párraga, D., Crespo-Picazo, J. L., Bernaldo De Quirós, Y., Cervera, V., Martí-Bonmati, L., Díaz-Delgado, J., et al. (2014). Decompression sickness (“the bends”) in sea turtles. Dis. Aquat. Org. 111, 191–205. doi: 10.3354/dao02790
García-Párraga, D., Lorenzo, T., Ortiz, J.-L., Ortega, J., Crespo-Picazo, J.-L., Cortijo, J., et al. (2018a). Deciphering function of the pulmonary arterial sphincters in loggerhead sea turtles (Caretta caretta). J. Exp. Biol. 221:jeb.179820. doi: 10.1242/jeb.179820
García-Párraga, D., Moore, M., and Fahlman, A. (2018b). Pulmonary ventilation– perfusion mismatch: a novel hypothesis for how diving vertebrates may avoid the bends. Proc. R. Soc. B 285:20180482. doi: 10.1098/rspb.2018.0482
Harms, C. A., Boylan, S. M., Stacy, B. A., Beasley, J. F., García-Párraga, D., and Godfrey, M. H. (2020). Gas embolism and massive blunt force trauma to sea turtles entrained in hopper dredges in north and South Carolina, USA. Dis. Aquat. Org. 142, 189–196. doi: 10.3354/dao03542
Hochscheid, S., Mcmahon, C. R., Bradshaw, C. J., Maffucci, F., Bentivegna, F., and Hays, G. C. (2007). Allometric scaling of lung volume and its consequences for marine turtle diving performance. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 148, 360–367. doi: 10.1016/j.cbpa.2007.05.010
Hodanbosi, M., Sterba-Boatwright, B., and Fahlman, A. (2016). Updating a gas dynamics model using estimates for California sea lions (Zalophus californianus). Respir. Physiol. Neurobiol. 234, 1–8. doi: 10.1016/j.resp.2016.08.006
Hooker, S. K., and Baird, R. W. (1999). Deep-diving behaviour of the northern bottlenose whale, Hyperoodon ampullatus (Cetacea: Ziphiidae). Proc. R. Soc. London B Biol. Sci. 266, 671–676. doi: 10.1098/rspb.1999.0688
Hooker, S. K., Baird, R. W., and Fahlman, A. (2009). Could beaked whales get the bends? Effect of diving behaviour and physiology on modelled gas exchange for three species: Ziphius cavirostris, Mesoplodon densirostris and Hyperoodon ampullatus. Respir. Physiol. Neurobiol. 167, 235–246. doi: 10.1016/j.resp.2009.04.023
Iverson, A. R., Fujisaki, I., Lamont, M. M., and Hart, K. M. (2019). Loggerhead Sea turtle (Caretta caretta) diving changes with productivity, behavioral mode, and sea surface temperature. PLoS One 14:e0220372. doi: 10.1371/journal.pone.0220372
Lutz, P. L., Bergy, A., and Bergy, M. (1989). Effects of temperature on gas exchange and acid-base balance in the sea turtle Caretta caretta at rest and during routine activity. J. Exp. Biol. 144, 155–169. doi: 10.1242/jeb.144.1.155
Okuyama, J., Shiozawa, M., and Shiode, D. (2020). Heart rate and cardiac response to exercise during voluntary dives in captive sea turtles (Cheloniidae). Biology Open 9:bio049247. doi: 10.1242/bio.049247
Olson, K. R., Whitfield, N. L., Bearden, S. E., St. Leger, J., Nilson, E., Gao, Y., et al. (2010). Hypoxic pulmonary vasodilation: a paradigm shift with a hydrogen sulfide mechanism. Am. J. Physiology-Integrative, Comp. Physiol. 298, R51–R60. doi: 10.1152/ajpregu.00576.2009
Parga, M. L., Crespo-Picazo, J. L., Monteiro, D., García-Párraga, D., Hernandez, J. A., Swimmer, Y., et al. (2020). On-board study of gas embolism in marine turtles caught in bottom trawl fisheries in the Atlantic Ocean. Sci. Rep. 10:5561. doi: 10.1038/s41598-020-62355-7
Patel, S. H., Barco, S. G., Crowe, L. M., Manning, J. P., Matzen, E., Smolowitz, R. J., et al. (2018). Loggerhead turtles are good ocean-observers in stratified mid-latitude regions. Estuar. Coast. Shelf Sci. 213, 128–136. doi: 10.1016/j.ecss.2018.08.019
Robinson, N. J., Morreale, S. J., Nel, R., and Paladino, F. V. (2017). Movements and diving behaviour of inter-nesting leatherback turtles in an oceanographically dynamic habitat in South Africa. Mar. Ecol. Prog. Ser. 571, 221–232. doi: 10.3354/meps12136
Robinson, N. J., and Paladino, F. V. (2015). “Diving behavior and physiology of the leatherback turtle,” in The Leatherback Turtle. eds. J. R. Spotila and P. S. Tomillo (Baltimore, USA: Johns Hopkins University Press), 21–31.
Sato, K., Charrassin, J.-B. T., Bost, C.-A., and Naito, Y. (2004). Why do macaroni penguins choose shallow body angles that result in longer descent and ascent durations? J. Exp. Biol. 207, 4057–4065. doi: 10.1242/jeb.01265
Sato, K., Sakamoto, W., Matsuzawa, Y., Tanaka, H., Minamikawa, S., and Naito, Y. (1995). Body temperature independence of solar radiation in free-ranging loggerhead turtles, Caretta caretta, during internesting periods. Mar. Biol. 123, 197–205. doi: 10.1007/BF00353611
Schellart, N. A. M., van Rees Vellinga, T. P., van Dijk, F. H., and Sterk, W. (2012). Doppler bubble grades after diving and relevance of body fat. Aviat. Space Environ. Med. 83, 951–957. doi: 10.3357/ASEM.3189.2012
Southwood, A. L., Andrews, R. D., Lutcavage, M. E., Paladino, F. V., West, N. H., George, R. H., et al. (1999). Heart rates and diving behavior of leatherback sea turtles in the eastern pacific ocean. J. Exp. Biol. 202, 1115–1125. doi: 10.1242/jeb.202.9.1115
Tal, D., Shachar-Bener, H., Hershkovitz, D., Arieli, Y., and Shupak, A. (2015). Evidence for the initiation of decompression sickness by exposure to intense underwater sound. J. Neurophysiol. 114, 1521–1529. doi: 10.1152/jn.00466.2015
Van Bonn, W., Dennison, S., Cook, P., and Fahlman, A. (2013). Gas bubble disease in the brain of a living California Sea lion (Zalophus californianus). Front. Aqua. Physiol. 4:5. doi: 10.3389/fphys.2013.00005
Van Bonn, W., Montie, E., Dennison, S., Pussini, N., Cook, P., Greig, D., et al. (2011). Evidence of injury caused by gas bubbles in a live marine mammal: barotrauma in a California Sea lion Zalophus californianus. Dis. Aquat. Org. 96, 89–96. doi: 10.3354/dao02376
Wallace, B. P., Kot, C. Y., Dimatteo, A. D., Lee, T., Crowder, L. B., and Lewison, R. L. (2013). Impacts of fisheries bycatch on marine turtle populations worldwide: toward conservation and research priorities. Ecosphere 4:art40. doi: 10.1371/journal.pone.0001041
Wallace, B. P., Lewison, R. L., Mcdonald, S. L., Mcdonald, R. K., Kot, C. Y., Kelez, S., et al. (2010). Global patterns of marine turtle bycatch. Conserv. Lett. 3, 131–142. doi: 10.1111/j.1755-263X.2010.00105.x
West, N. H., Butler, P. J., and Bevan, R. M. (1992). Pulmonary blood flow at rest and during swimming in the green turtle, Chelonia mydas. Physiol. Zool. 65, 287–310. doi: 10.1086/physzool.65.2.30158254
Relative compartment size (blood, brain, fat, central circulation, and muscle) from seven loggerhead, eight green, and one leatherback turtles. Compartment size in five of the seven loggerhead turtles was assessed using computed tomography, while it was measured in-situ via necropsy in the remaining animals. *Blood was not measured directly and was estimated from values in the published literature.
Keywords: Physiology, ecological modeling, fisheries, decompression sickness, conservation, dive behavior
Citation: Robinson NJ, García-Párraga D, Stacy BA, Costidis AM, Blanco GS, Clyde-Brockway CE, Haas HL, Harms CA, Patel SH, Stacy NI and Fahlman A (2021) A Baseline Model For Estimating the Risk of Gas Embolism in Sea Turtles During Routine Dives. Front. Physiol. 12:678555. doi: 10.3389/fphys.2021.678555
Edited by:Fabrice Joulia, Université de Toulon, France
Reviewed by:Ghulam Nabi, Hebei Normal University, China
Rostain Jean-Claude, Aix Marseille Université, France
Copyright © 2021 Robinson, García-Párraga, Stacy, Costidis, Blanco, Clyde-Brockway, Haas, Harms, Patel, Stacy and Fahlman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Andreas Fahlman, firstname.lastname@example.org