The ability to control hydration state is essential for terrestrial species, especially amphibians, which are highly susceptible to dehydration. Here, we examined how temperature (17°C versus 22°C) influenced behavioural hydroregulation in spotted salamanders (Ambystoma maculatum) using a laboratory humidity gradient. Salamanders defended a constant vapour pressure deficit (VPD) between temperatures by targeting higher relative humidity at 22°C than at 17°C, possibly to compensate for increased evaporative demand at warmer temperatures. Individuals selecting higher VPDs experienced greater evaporative water loss (EWL), with larger salamanders losing more water than smaller ones after accounting for temperature. Together, these results highlight a trade-off among body size, humidity preference and desiccation tolerance. Salamanders also rehydrated faster at 22°C than 17°C, highlighting temperature-dependent water uptake rates. Our finding that salamanders regulated a constant driving force of evaporation between temperatures suggests they are able to detect rates of EWL. Local evaporative cooling of the skin is a plausible mechanism: if moist-skinned ectotherms show local evaporative cooling on exposed surfaces but not on those in contact with the substrate, then the temperature gradient between dorsal and ventral skin could act as an effective cue to detect the drive for EWL. Ultimately, our study underscores the complexity of amphibian hydroregulation and emphasises the role of behaviour in maintaining hydration state.

The transition from aquatic to terrestrial environments was one of the most significant milestones in vertebrate evolution (Laurin, 2010). On land, water availability fluctuates both spatially and temporally, meaning that organisms face constant challenges related to the maintenance of water balance (Takei, 2015). Therefore, the ability to control hydration state is essential for terrestrial species. Effective hydroregulation may be achieved through a combination of non-mutually exclusive physiological and behavioural strategies. Physiological processes include the enhancement of skin resistance to water loss (Tracy et al., 2008), the presence of waterproofing lipids in the integument (Sadowski-Fugitt et al., 2012) or on the integument (Lillywhite, 2006), or the modulation of kidney and bladder function to increase water retention (Schmidt-Nielsen and Schmidt-Nielsen, 1952; Shoemaker and Bickler, 1979). Hydroregulatory behaviours may comprise the adoption of water-conserving postures (Pough et al., 1983), aggregation (Broly et al., 2014) or burrowing (Bulova, 2002). Often, however, behavioural hydroregulation is oversimplified to the selection of humid microhabitats (but see Dezetter et al., 2023). As a result, we still have limited information about how animals perceive and process changes in habitat humidity, and how behaviour can be used to mitigate potential disruptions to the maintenance of water balance (Hillyard and Willumsen, 2011; Hillyard et al., 1998).

Among vertebrate ectotherms, amphibians in particular face elevated hydroregulatory challenges. Amphibians actually have a relatively thin skin that is supported by blood vessels in the epidermis, which together facilitate cutaneous gas exchange (Tattersall, 2007). However, the ability to exchange gases through the skin also makes amphibians highly susceptible to evaporative water loss (EWL), especially at high temperatures (Lertzman-Lepofsky et al., 2020). The extent to which temperature impacts EWL depends on the vapour pressure deficit (VPD), which measures the drying power of air, and is determined by the difference between the saturated water vapour pressure (WVPsat) and the actual water vapour pressure (WVPact) (Anderson, 1936). Since warm temperatures typically produce high VPD, amphibians facing these conditions should suffer high evaporative demand and pay high hydroregulatory costs (Riddell et al., 2019, 2024). Because of this, many have argued that thermal constraints over hydroregulation explained the recurrent pattern of low body temperatures (Tb) in amphibians, on the basis that sustaining low Tb would minimise EWL (Brattstrom, 1979). This view has been recently expanded into a framework that posits that amphibians constantly balance thermoregulation and hydroregulation to sustain performance (Greenberg and Palen, 2021; Titon and Gomes, 2017). However, while this framework has both thermosensation and hydrosensation components, our current knowledge of the former outweighs the latter (Ortega et al., 2023; Tracy et al., 2008), particularly in terms of the behaviours that contribute to the maintenance of hydration state.

One possible explanation for this knowledge gap is the difficulty in disentangling temperature and humidity effects, since most laboratory studies test amphibians in combined humidity and thermal gradients (hereafter ‘hydrothermal gradients’) (Bundy and Tracy, 1977; Galindo et al., 2018). While one may argue that the use of hydrothermal gradients provides a realistic scenario, any response measured in such an apparatus is by definition confounded by temperature and humidity (Mitchell and Bergmann, 2016). Because one cannot properly gauge whether the observed behaviours are innately hydroregulatory, the extensive use of hydrothermal gradients resulted in the simplistic yet prevalent view that, for amphibians, ‘wetter is better’ (see Brattstrom, 1979). However, if one is interested in assessing humidity preference and the behaviours that underpin hydroregulation, one should test amphibians in a temperature-controlled humidity gradient (Heatwole, 1962; Shelford, 1913). For example, Spotila (1972) demonstrated interspecific differences in relative humidity (RH) preference (RH selection range: 60.4–90% at 15°C) across 12 closely related species of salamanders. Humidity preference was associated with microhabitat use, suggesting that salamanders were able to detect differences in air moisture content and respond in a manner that maximised hydration state (Spotila, 1972). Despite their potential, temperature-controlled humidity gradients have seldom been used in amphibian research since the 1970s (but see Chapman and Bidwell, 2023). As such, our understanding of behavioural hydroregulation remains limited.

In this study, we used a laboratory humidity gradient to assess how temperature (17°C versus 22°C) affected behavioural hydroregulation in the spotted salamander (Ambystoma maculatum). We focused on two main behavioural metrics: VPD selection and RH selection. By measuring VPD selection, one is able to obtain behavioural information based on a temperature-independent humidity gradient and gauge behavioural hydroregulation sensu stricto (Dezetter et al., 2023; Wolcott and Wolcott, 2001). In contrast, RH selection provides information on how changes in temperature impact the behaviours that contribute to the maintenance of hydration state (Galindo et al., 2018; Spotila, 1972). Based on this, we predicted that A. maculatum would alter RH selection to maintain a constant VPD between temperatures (Spotila, 1972). To further understand how behaviour contributes to the maintenance of hydration state, we tested for a functional trade-off among EWL, VPD selection, activity level and temperature. We predicted that highly active salamanders that also selected high VPD would be the ones with the highest EWL; these effects should be greater at 22°C than at 17°C (Heatwole and Newby, 1972; Spight, 1968). Finally, we also tested whether rehydration rates (ReRs) increased as a function of EWL and test temperature (Claussen, 1969; Spight, 1967). With this integrative approach, we aim to provide a broader understanding of the role of behaviour in the maintenance of hydration state in amphibians.

Salamander collection and husbandry

In May 2022, we collected 44 post-breeding, adult Ambystoma maculatum (Shaw 1802) (22 females, 22 males) with a drift fence installed around Bat Lake, Algonquin Provincial Park, ON, Canada (45.5857°N, 78.5185°W). We recorded the sex of each salamander based on cloacal morphology (Petranka, 1998) and identified each individual based on their unique dorsal spot pattern. To transfer the salamanders to Brock University, we allocated them to plastic containers with ventilated fitted lids (34 cm×19.6 cm×12 cm). These containers had sphagnum moss, pine needles and water, and were placed in a transport box at ∼4°C to avoid overheating and dehydration during transportation.

In the lab, we housed salamanders in pairs in ventilated tanks containing coconut husk fibre, sphagnum moss, PVC pipe refuges and a water dish. These tanks were kept in a facility with controlled temperature, RH and photoperiod. We adjusted temperature and photoperiod seasonally to simulate native habitat conditions while keeping RH constant at 70% (Giacometti and Tattersall, 2024). We fed salamanders twice a week with mealworms, nightcrawlers and isopods dusted in multivitamin and calcium powder. We weighed salamanders weekly to the nearest 0.01 g using an analytical scale (Mettler Toledo, model PB602-S) to monitor changes in body mass (Mb) as a proxy of well-being; all animals maintained Mb throughout their period in the lab. The current study was conducted between late June and August 2024 with salamanders that had been acclimatised to summer conditions (14°C and a 14 h:12 h light:dark cycle) for at least 4 weeks.

Laboratory humidity gradient

To assess behavioural hydroregulation, we used an annular humidity gradient (total length=120 cm, outer radius=45 cm, inner radius=30 cm) built by Brock University Technical Services (Fig. S1). The thermally conductive floor of the gradient had copper pipes connected to a water bath (Haake™, model DC10) kept at either 17°C or 22°C. These temperatures represent, respectively, the median and maximum upper selected temperatures by A. maculatum in the summer (Giacometti and Tattersall, 2024). The gradient chamber was divided into an inner and outer lane by a 10.5 cm tall partition attached to the gradient lid, allowing the simultaneous study of two individuals per experiment (visually isolated). To create a humidity gradient, we delivered air of different moisture content to compartments (dry, intermediate and humid) in the gradient chamber. We created these compartments by attaching sponges (4 cm×3 cm×9 cm) to the gradient wall every 30 cm. To ensure the sponges remained in place and allowed for salamander navigation during the experiments, we sewed the hook side of Velcro tape onto the sponges and used aquarium safe silicon sealant (Adhesive Guru, product AG310) to attach the loop side of Velcro tape to the gradient wall. The gradient lid had two venting holes every 90 deg relative to the horizontal axis, such that the venting holes were positioned above the centre of the dry compartment at 0 deg, the intermediate compartments at 90 deg, and the humid compartment at 180 deg (Fig. S1).

To generate dry air at 0 deg, we used tubing to direct compressed air through a column of desiccant and then into the gradient at 1900 ml min−1 (split to direct evenly to both lanes). We generated humid air using a water vapour bubbling system connected to the venting holes at 180 deg. This system consisted of an 1800 ml flask filled with dechlorinated water and kept 3°C warmer than the prevailing experimental temperature. We aerated the water in the flask with an aquarium air stone (15 cm diameter) attached to tubing that received compressed air at 2100 ml min−1. The pressure created inside the flask directed water vapour from the flask to the venting holes in the gradient lid (split evenly to both lanes). To generate intermediate moisture levels at 90 deg, we used a manifold to create a subsample line that combined air from the dry and humid lines and directed mixed air into the gradient chamber at 1000 ml min−1. Since dry air is heavier than humid air (Haynes, 2016), we drew air at 250 ml min−1 from the dry line and 750 ml min−1 from the humid line. We verified flow rates with flow meters (Aalborg Instruments & Controls Inc., model PMR1-010972) at the beginning and end of each experiment to ensure proper functioning of the humidity gradient.

With this setup, we created a stable, near-linear humidity gradient at 17°C (R2=0.85) and 22°C (R2=0.84) (Fig. S2). We determined this with hygro-thermometers (Inkbird, model IBS-TH1; accuracy±3%) placed every 22.5 deg and set to record temperature and RH every 60 s for 24 h in the absence of any animals in the gradient. Since RH is a temperature-dependent measurement and therefore not readily comparable between 17°C and 22°C, we used VPD as a measure of gradient humidity between temperatures (Table S1). Following Haynes (2016), where T is temperature (°C), and vapour pressure units are kPa, we calculated:
(1)

We conducted our experiments in total darkness inside a temperature-controlled room (maintained at either 17 or 22°C). We placed the humidity gradient on top of 15 cm of insulation foam board and padding so that inadvertent building vibrations would not disturb the salamanders during the experiments. We positioned an infrared illuminator (wavelength=850 nm; TVPSii, model TP-IRBP15) and a high-resolution infrared webcam (Agama, model V-1325R) 180 cm above the centre of the gradient to record the salamanders via a time-lapse image acquisition software (HandyAVI®) set to capture an image every 30 s.

Experimental design

We allowed the salamanders a total of 12 h inside the humidity gradient, from ∼21:00 h to ∼09:00 h. Each salamander was randomly tested twice (once under each experimental temperature), with a minimum interval of 10 days between experiments. We considered the initial 3 h as the habituation period and the subsequent 9 h as the experimental period. All individuals were fasted for at least 1 week prior to experimentation to ensure a post-absorptive state (Secor and Boehm, 2006). We always allowed a minimum interval of 10 h between experiments after disinfecting the humidity gradient with 70% ethanol at the end of an experiment.

We always handled the salamanders using nitrile gloves. Before introducing salamanders into the gradient, we placed them individually into a container filled with 30 ml dechlorinated water for 15 min so that water could be absorbed through the skin and animals could be introduced in a high hydration state. Then, we blotted the salamanders with paper towel, weighed them, and placed them into the gradient. We always introduced the salamanders into the gradient through the intermediate compartment, determining at random whether animals would be initially facing the dry or humid end, and whether they would be tested in the inner or outer gradient lanes. To minimise handling stress over hydroregulatory behaviours, we did not handle or manipulate the salamanders during the experiments. After finishing an experiment, we removed the salamanders from the gradient, blotted them with paper towel and weighed them. We considered the difference in Mb before and after an experiment divided by the length of our experiment as a proxy for evaporative water loss (EWL; ml h−1) (Ortega-Chinchilla et al., 2023). We determined ReR (ml h−1) by placing the salamanders individually in a container with 30 ml of dechlorinated water for 30 min. Then, we removed the salamanders from the containers, blotted them with paper towel, and reweighed them to determine ReR based on Mb gained per 30 min (Sherman and Stadlen, 1986). After this, we placed the salamanders back in their housing tanks.

Data processing

We obtained a total of 1440 images for each individual over the course of each 12 h experiment. We imported image sequences into Fiji (Schindelin et al., 2012) and for each experiment we recorded the identity (ID) and sex of the individual, as well as the test temperature (17°C or 22°C). We used the manual tracking plug-in in Fiji to track the position of the mid-body of the salamander within the gradient for each image in the sequence. As an output, we obtained a set of Cartesian (x,y) coordinates that had x,y=0,0 as the top left of the image. Since the coordinates in Fiji did not correspond to the centroid of the gradient, we had to re-centre (x,y=0,0 in the centre of the centroid) and convert Fiji coordinates into polar (θ) coordinates to convert salamander body position into selected gradient humidity; these calculations are defined in detail in Giacometti et al. (2021). We determined humidity selection in the gradient by calculating median selected VPD and median selected RH.

To quantify activity in the gradient, we measured the distance moved (dt) by each individual every 30 s based on the distance defined by the arc of a circle:
(2)
where r is the mean radius, and θ are polar coordinates between adjacent time points t and t–1 for all n time points. We then obtained total distance moved (Dt) by taking the sum of dt:
(3)
where n is the total number of images (n=1440).
Finally, we assessed the temperature dependency of ReR, using the Van't Hoff equation to calculate the temperature coefficient (Q10) for ReR between 17°C (T1) and 22°C (T2):
(4)

Data analyses

We performed all analyses using R (version 4.4.1) in RStudio (version 2024-06-14; r-project.org) considering a significance level of 0.05. We built linear mixed-effects models using the ‘lmer’ function from the lme4 package (Bates et al., 2015), and we included ID as a random term in all models to account for repeated measures. We tested our hypotheses using data from the experimental period. We fit two individual models to test how humidity selection differed between temperatures. To this end, we considered either median VPD or RH as the response variables, experimental temperature as the predictor, and log-transformed Mb as a covariate. To assess how EWL differed between temperatures, we fitted a model that had log-transformed EWL as the response variable, experimental temperature, median VPD and Dt as the predictor variables and log-transformed Mb as a covariate. We also evaluated how ReR differed between temperatures with a model that had log-transformed ReR as the response variable, experimental temperature and log-transformed EWL as the predictors and log-transformed Mb as a covariate. For each model, we assessed residual autocorrelation with the ‘checkresiduals’ function from the forecast package (https://CRAN.R-project.org/package=forecast), and the ‘acf’ and ‘qqnorm’ functions from the R stats package. We evaluated model fit with the ‘check_model’ function from the performance package (Lüdecke et al., 2021) and computed marginal fixed effects through the ‘ggeffect’ function from the ggeffects package (Lüdecke, 2018). We created figures using the ggplot2 (https://CRAN.R-project.org/package=ggplot2), Thermimage (https://zenodo.org/records/11374916) and ggpubr (https://CRAN.R-project.org/package=ggpubr) packages.

During the experimental period, salamanders defended a constant VPD regardless of test temperature (Table 1; Fig. 1A). This effect was mediated by Mb, with smaller individuals selecting lower VPD than larger individuals (Table 2). Importantly, salamanders avoided high VPD at both test temperatures (Fig. S3). Salamanders selected different RH between temperatures (Table 1; Fig. 1B), with smaller individuals selecting higher RH than larger individuals (Table 2). Regardless of temperature and activity levels, salamanders with high EWL were those that selected higher VPD and were larger (Table 3; Fig. 2). Test temperature was the only factor impacting ReR (Q10=2.54±2.26), with salamanders rehydrating faster at 22°C than at 17°C (Table 4; Fig. 3). Salamanders moved similar total distances within the gradient regardless of temperature (Table 1).

Fig. 1.

Salamanders defend a constant vapour pressure deficit (VPD) through behavioural selection of different relative humidity (RH) levels between 17°C and 22°C. (A) VPD and (B) RH at 17 and 22°C. In both panels, large dots show mean marginal effects and bars indicate the corresponding 95% confidence intervals. Small dots show the mean predicted value for each individual salamander at a given temperature (17°C=blue, 22°C=red).

Fig. 1.

Salamanders defend a constant vapour pressure deficit (VPD) through behavioural selection of different relative humidity (RH) levels between 17°C and 22°C. (A) VPD and (B) RH at 17 and 22°C. In both panels, large dots show mean marginal effects and bars indicate the corresponding 95% confidence intervals. Small dots show the mean predicted value for each individual salamander at a given temperature (17°C=blue, 22°C=red).

Close modal
Fig. 2.

Relationship between log-transformed rates of evaporative water loss (logEWL) and log-transformed body mass (logMb) in Ambystoma maculatum tested in a humidity gradient at either 17°C or 22°C. Larger individuals had higher EWL and individuals that selected high VPDs were also the ones with higher EWL at both temperatures. In both panels, each dot represents an individual salamander. The solid lines and shaded areas indicate the predicted relationship between logEWL and logMb, and the 95% confidence interval, respectively. Vapour pressure deficit (VPD) is colour coded, with colder colours indicating wetter conditions and warmer colours indicating drier conditions.

Fig. 2.

Relationship between log-transformed rates of evaporative water loss (logEWL) and log-transformed body mass (logMb) in Ambystoma maculatum tested in a humidity gradient at either 17°C or 22°C. Larger individuals had higher EWL and individuals that selected high VPDs were also the ones with higher EWL at both temperatures. In both panels, each dot represents an individual salamander. The solid lines and shaded areas indicate the predicted relationship between logEWL and logMb, and the 95% confidence interval, respectively. Vapour pressure deficit (VPD) is colour coded, with colder colours indicating wetter conditions and warmer colours indicating drier conditions.

Close modal
Fig. 3.

Difference in log-transformed rehydration rates (logReR) between temperatures. Salamanders rehydrated faster at 22°C than at 17°C. Large dots show mean marginal effects and bars show the corresponding 95% confidence intervals. Small dots indicate the mean predicted value for each individual salamander between temperatures (17°C in blue and 22°C in red).

Fig. 3.

Difference in log-transformed rehydration rates (logReR) between temperatures. Salamanders rehydrated faster at 22°C than at 17°C. Large dots show mean marginal effects and bars show the corresponding 95% confidence intervals. Small dots indicate the mean predicted value for each individual salamander between temperatures (17°C in blue and 22°C in red).

Close modal
Table 1.

Behavioural, body size and physiological parameters of Ambystoma maculatum tested in a humidity gradient kept at either 17°C or 22°C

17°C (N=44)22°C (N=44)
Selected VPD (kPa) 0.56±0.16 (0.27, 0.95) 0.57±0.22 (0.12, 1.06) 
Selected RH (%) 71.20±8.31 (51.01, 85.88) 78.48±8.34 (60.09, 95.30) 
Mb (g) 12.00±3.91 (5.02, 23.70) 12.36±4.62 (4.66, 25.80) 
EWL (ml h−10.07±0.05 (0.01, 0.19) 0.06±0.06 (0.01, 0.30) 
ReR (ml h−1)  0.18±0.13 (0.02, 0.56) 0.31±0.25 (0.02, 1.14) 
Total distance moved (m) 40.00±20.80 (0.11, 90.62) 40.31±26.94 (0.34, 105.01) 
17°C (N=44)22°C (N=44)
Selected VPD (kPa) 0.56±0.16 (0.27, 0.95) 0.57±0.22 (0.12, 1.06) 
Selected RH (%) 71.20±8.31 (51.01, 85.88) 78.48±8.34 (60.09, 95.30) 
Mb (g) 12.00±3.91 (5.02, 23.70) 12.36±4.62 (4.66, 25.80) 
EWL (ml h−10.07±0.05 (0.01, 0.19) 0.06±0.06 (0.01, 0.30) 
ReR (ml h−1)  0.18±0.13 (0.02, 0.56) 0.31±0.25 (0.02, 1.14) 
Total distance moved (m) 40.00±20.80 (0.11, 90.62) 40.31±26.94 (0.34, 105.01) 

Values are shown as means±s.d. (min., max.)

VPD, vapour pressure deficit; RH, relative humidity; Mb, body mass; EWL, rate of evaporative water loss; ReR, rehydration rate.

Table 2.

Parameter estimates (β), 95% confidence intervals (95% CI) and P-values for the models testing how humidity selection differed between temperatures in A. maculatum (N=44)

Estimate95% CIP
VPD∼temperature+logMb+(1|ID)   
 Predictor Intercept 0.40 0.26–0.53 <0.001 
Temperature (22°C) 0.01 −0.06–0.07 0.883 
logMb 0.01 0.00–0.02 0.011 
 Random effect σ2 0.03 
τ00 ID 0.01 
ICC 0.22 
Marginal R2/conditional R2 0.10/0.30 
RH∼temperature+logMb+(1|ID) 
 Predictor Intercept 87.31 73.65–100.98 <0.001 
Temperature (22°C) 7.42 4.45–10.39 <0.001 
logMb −15.22 −27.96–−2.49 0.020 
 Random effect σ2 48.92 
τ00 ID 16.14 
ICC 0.25 
Marginal R2/conditional R2 0.22/0.42 
Estimate95% CIP
VPD∼temperature+logMb+(1|ID)   
 Predictor Intercept 0.40 0.26–0.53 <0.001 
Temperature (22°C) 0.01 −0.06–0.07 0.883 
logMb 0.01 0.00–0.02 0.011 
 Random effect σ2 0.03 
τ00 ID 0.01 
ICC 0.22 
Marginal R2/conditional R2 0.10/0.30 
RH∼temperature+logMb+(1|ID) 
 Predictor Intercept 87.31 73.65–100.98 <0.001 
Temperature (22°C) 7.42 4.45–10.39 <0.001 
logMb −15.22 −27.96–−2.49 0.020 
 Random effect σ2 48.92 
τ00 ID 16.14 
ICC 0.25 
Marginal R2/conditional R2 0.22/0.42 

We considered either VPD or RH as response variable, temperature and logMb as predictors and ID as a random term. Significant parameters are shown in bold.

σ2, residual variance; τ00 ID, individual variance; ICC, intraclass correlation coefficient.

Table 3.

Parameter estimates (β), 95% CIs and P-values for the models testing the trade-off between EWL and VPD selection between temperatures in A. maculatum (N=44)

logEWL∼VPD+temperature+logMb+ Dt+(1|ID)
Estimate95% CIP
Predictor Intercept −2.27 −2.84–−1.69 <0.001 
VPD 0.58 0.15–1.01 0.009 
Temperature (22°C) −0.08 −0.21–0.05 0.208 
logMb 0.64 0.10–1.18 0.021 
Total distance moved 0.00 −0.00–0.00 0.946 
Random effect σ2 0.09 
τ00 ID 0.02 
ICC 0.20 
Marginal R2/conditional R2 0.21/0.36 
logEWL∼VPD+temperature+logMb+ Dt+(1|ID)
Estimate95% CIP
Predictor Intercept −2.27 −2.84–−1.69 <0.001 
VPD 0.58 0.15–1.01 0.009 
Temperature (22°C) −0.08 −0.21–0.05 0.208 
logMb 0.64 0.10–1.18 0.021 
Total distance moved 0.00 −0.00–0.00 0.946 
Random effect σ2 0.09 
τ00 ID 0.02 
ICC 0.20 
Marginal R2/conditional R2 0.21/0.36 

We considered logEWL as the response variable, temperature, VPD, logMb and total distance moved (Dt) as the predictors, and ID as a random term. Significant parameters are shown in bold.

Table 4.

Parameter estimates (β), 95% CIs and P-values for the models testing how rehydration rates varied as a function of evaporative water loss and temperature in A. maculatum (N=44)

logReR∼logEWL+temperature+logMb+(1|ID)
Estimate95% CIP
Predictor Intercept −0.83 −1.57–−0.09 0.028 
logEWL 0.17 −0.06–0.39 0.147 
Temperature (22°C) 0.20 0.04–0.36 0.013 
logMb 0.18 −0.38–0.73 0.526 
Random effect σ2 0.13 
τ00 ID 0.00 
Observations 88 
Marginal R2/conditional R2 0.10/0.00 
logReR∼logEWL+temperature+logMb+(1|ID)
Estimate95% CIP
Predictor Intercept −0.83 −1.57–−0.09 0.028 
logEWL 0.17 −0.06–0.39 0.147 
Temperature (22°C) 0.20 0.04–0.36 0.013 
logMb 0.18 −0.38–0.73 0.526 
Random effect σ2 0.13 
τ00 ID 0.00 
Observations 88 
Marginal R2/conditional R2 0.10/0.00 

We considered logReR as the response variable, logEWL, temperature and logMb as the response variables, and ID as a random term. Significant parameters are shown in bold.

σ2=residual variance; τ00 ID=individual variance.

We used an integrative approach to evaluate the effect of temperature on behavioural hydroregulation in A. maculatum. Our findings suggest that salamanders behaviourally defended a constant VPD by altering RH selection between temperatures. By targeting higher RH at 22°C than at 17°C, salamanders potentially compensate for increased evaporative drive at warmer temperatures. We also found evidence for a trade-off among EWL, body size, and humidity selection in fully hydrated individuals. Specifically, larger salamanders selected higher VPDs and exhibited greater EWL than smaller ones, suggesting a size-dependent trade-off between humidity selection and evaporation. Temperature was the only factor explaining ReR, with salamanders rehydrating faster at 22°C than at 17°C. Together, our findings underscore the crucial role of behavioural strategies in the maintenance of hydration state in amphibians.

We found that salamanders defended a constant selected VPD between temperatures. Since VPD represents the driving force for EWL (Riddell et al., 2024), our data suggest that A. maculatum may be able to detect their rates of water loss and alter their behaviour patterns to control hydration state. This hypothesis, however, remains to be tested against other possible hypotheses that may include, for example, changes in body fluid composition or stretch receptors in the skin (Anderson et al., 2017). To maintain homeostasis of VPD, salamanders selected higher RH at 22°C than at 17°C. The maintenance of a constant VPD across temperatures indicates that amphibians may sense, process, and respond to moisture content in the air in a physiologically relevant manner that confers control over hydration state (Shelford, 1913). The prevailing paradigm in amphibian hydrothermal biology posits that individuals with high EWL must hydroregulate and will consequently face elevated thermoregulatory costs because of evaporative cooling (Tracy, 1976). Our work adds to this framework by demonstrating that A. maculatum not only showed evidence of humidity preference, but also employed behavioural adjustments to mitigate increased evaporative demand at warmer temperatures.

Our humidity preference values for A. maculatum were similar to those reported in the congener A. opacum (VPD=0.22 kPa or RH=74% at 13°C) (Marangio and Anderson, 1977). Humidity preference in Ambystoma species, however, appears to be lower than that reported for other amphibians (Bellis, 1962; Galindo et al., 2018; Shoemaker et al., 1992). The fossorial habit of most Ambystoma species coupled with their relatively large body size could explain this pattern. Fossorial amphibians typically bear behavioural and morphophysiological adaptations that allow for an increased independence from standing water in the environment compared with non-fossorial species (Shoemaker et al., 1992). For example, Ceratophrys stolzmanni uses soil humidity as a cue to adjust burrowing depth and maintain water balance (Székely et al., 2018) and A. maculatum can tolerate up to 41% of Mb loss through evaporation before losing function (Hall, 1922). Furthermore, our analyses showed that larger salamanders selected overall higher VPD than smaller ones at both test temperatures, suggesting that individuals with greater water content may be able to withstand greater VPDs (MacCracken and Stebbings, 2012). Small-bodied salamanders often coil to decrease the surface area available for EWL (e.g. Heatwole, 1960); however, we never observed A. maculatum assuming water-conserving postures despite the relatively wide range of body size in our study animals (4.66–25.80 g). Besides possessing a heightened desiccation tolerance (Giacometti and Tattersall, 2025) and being able to maintain hydration state behaviourally, A. maculatum may also benefit from the presence of cutaneous grooves in the venter and flanks which aid in transferring water to different parts of the body by capillarity (Lopez and Brodie, 1977; Toledo and Jared, 1993).

Our study also highlights how body size may influence behavioural hydroregulation in A. maculatum. Larger salamanders selected higher VPD and exhibited greater EWL than smaller salamanders, indicating a size-dependent trade-off between evaporation and humidity selection. Contrary to our prediction, however, the relationship between EWL and VPD was not mediated by temperature nor activity level. The allometric scaling between body size and EWL suggests that, relative to their size, smaller individuals should evaporate faster than larger individuals because small individuals have a high surface area-to-volume ratio (Jørgensen, 1997). That is, a salamander that weighs 5 g and loses 1 ml of water through evaporation would face greater hydroregulatory consequences than a salamander that weighs 15 g and loses the same volume of water. Therefore, smaller salamanders are predicted to experience greater hydric costs than larger salamanders, and behavioural adjustments should be crucial to mitigate EWL (Riddell et al., 2024). Our results supported this notion, as evidenced by the fact that smaller A. maculatum prioritised the selection of low VPDs that culminated in overall low EWL. By contrast, larger A. maculatum selected relatively high VPDs that resulted in increased EWL. In nature, size-specific habitat use is an important behavioural strategy to minimise dehydration in aboveground-dwelling salamanders (Heatwole, 1962; Spotila, 1972). However, the extent to which size-dependent habitat selection contributes to hydroregulation in fossorial amphibians still warrants further research.

Consistent with previous work (Claussen, 1969; Spight, 1967), we found that A. maculatum rehydrated faster at 22°C than at 17°C. This thermal dependency of rehydration highlights the dual role of temperature as both a driver of EWL (Giacometti and Tattersall, 2025) and a modulator of ReR. Additionally, we did not find a relationship between EWL and ReR in the current study. A possible reason for this lies in the fact that we allowed A. maculatum to hydroregulate at will prior to ReR assessments. Evidence from Rhinella marina indicates that toads may gain water at rates as fast as 30 ml h−1 when dehydrated to 30% of their initial Mb (Kosmala et al., 2020). However, water uptake is minimal when R. marina have unlimited access to water (Juarez et al., 2024). Our results back up this pattern, as rehydration in A. maculatum mostly consisted of salamanders regaining their initial Mb. The osmotic gradient between animal and water is one of the main drivers of ReR, with dehydrated animals with higher plasma osmolarity typically showing higher ReR (Anderson et al., 2017). Since A. maculatum lost a similar proportion of Mb between the two temperatures (6.69±3.57% at 17°C and 6.01±3.81% at 22°C), our reported values likely do not represent maximum ReRs. Despite the integumentary and vascular adaptations that favour water retention in salamanders, ReRs in this group are still lower than those of anurans (Marangio and Anderson, 1977; Spight, 1967; Toledo and Jared, 1993). Nonetheless, our calculated Q10 value for ReR suggests that at least some degree of temperature-dependent physiological regulation of water uptake exists in A. maculatum. The mechanisms behind such regulation are still unclear, but may involve hormonal control of skin permeability to enhance water uptake (Alvarado and Johnson, 1965), along with an increase in blood circulation at warmer temperatures (Jørgensen, 1997).

Conclusions

The ability to maintain hydration state under varying thermal conditions is crucial to amphibians (Ortega-Chinchilla et al., 2023). By integrating behavioural and physiological perspectives, our study demonstrates the critical role of temperature-compensating behaviours to the maintenance of hydration state in A. maculatum. Our finding that salamanders maintained a constant VPD between temperatures challenges the traditional use of RH measurements to assess hydroregulation and habitat suitability in amphibians (Galindo et al., 2018). Thus, to properly understand amphibian hydroregulation, future studies should incorporate information on how behaviour is used to manage water loss, rather than rely on simple physical measurements. Although our study provides insight into the behaviours involved in hydroregulation, we still have limited information about hydrosensation compared with thermal sensation. In this context, clarifying how amphibians use water-sensing units in the skin to detect water content in the air is invaluable (Hillyard, 1999; Spotila, 1972). Likewise, the observation that salamanders regulate a constant driving force of evaporation under changing thermal regimes suggests a mechanism for detecting rates of EWL. If moist-skinned ectotherms show local evaporative cooling on exposed surfaces but not on those in contact with the substrate (Andrade, 2016; Sinsch, 1984), then the temperature gradient between dorsal and ventral skin could act as a cue to effectively detect the drive for EWL. This hypothesis, however, remains to be tested and research is still necessary to determine the cues and sensors for EWL in amphibians. For example, future work may test the amount to which dorsal skin temperatures differ from ventral temperatures during any substantial exposure to drying conditions while measuring EWL concomitantly. Ultimately, our study underscores how behavioural adjustments are central to creating conditions conducive to the maintenance of a homeostatic hydration state in amphibians.

We thank Ontario Parks, the Ministry of Northern Development, Mines, Natural Resources, and Forestry, the Algonquin Wildlife Research Station, Kevin Kemmish, Njal Rollinson and Patrick Moldowan for facilitating access to Bat Lake and our study animals. We are grateful to Shawn Bukovac, Kristin Bray, Sarah Kehoe, Tom Eles and Natasha Hearn for assistance with animal care. We thank Stephen Renda, Mitch Sillaste and Art Reimer from Brock University's Machine Shop for building the apparatus used in our study. We would also like to thank Denis V. Andrade and the two anonymous reviewers whose suggestions greatly improved the quality of our manuscript.

Author contributions

Conceptualization: D.G., G.J.T.; Data curation: D.G.; Formal analysis: D.G.; Funding acquisition: D.G., G.J.T.; Investigation: D.G., G.J.T.; Methodology: D.G., G.J.T.; Project administration: D.G., G.J.T.; Resources: D.G., G.J.T.; Software: G.J.T.; Supervision: D.G., G.J.T.; Validation: D.G., G.J.T.; Visualization: D.G.; Writing – original draft: D.G., G.J.T.; Writing – review & editing: D.G., G.J.T.

Funding

Funding for this study was provided by a Doris White Memorial Graduate Scholarship and a Gervan Fearon Graduate Scholarship to D.G. and a Natural Sciences and Engineering Research Council of Canada Discovery Grant to G.J.T. (RGPIN-2020-05089). Open Access funding provided by Brock University. Deposited in PMC for immediate release.

Data and resource availability

Our data and code can be accessed from the Brock University Dataverse: https://doi.org/10.5683/SP3/IKHBTP.

Alvarado
,
R. H.
and
Johnson
,
S. R.
(
1965
).
The effects of arginine vasotocin and oxytocin on sodium and water balance in Ambystoma
.
Comp. Biochem. Physiol.
16
,
531
-
546
.
Anderson
,
D. B.
(
1936
).
Relative humidity or vapor pressure deficit
.
Ecology
17
,
277
-
282
.
Anderson
,
R. C. O.
,
Bovo
,
R. P.
,
Eismann
,
C. E.
,
Menegario
,
A. A.
and
Andrade
,
D. V.
(
2017
).
Not good, but not all bad: dehydration effects on body fluids, organ masses, and water flux through the skin of Rhinella schneideri (Amphibia, Bufonidae)
.
Physiol. Biochem. Zool.
90
,
313
-
320
.
Andrade
,
D. V.
(
2016
).
Temperature effects on the metabolism of amphibians and reptiles
. In
Amphibian and Reptile Adaptations to the Environment; Interplay Between Physiology and Behavior
(ed.
D. V.
de Andrade
,
C. R.
Bevier
and
J. E.
de Carvalho
), pp.
129
-
254
.
Boca Raton
:
CRC Press
.
Bates
,
D.
,
Mächler
,
M.
,
Bolker
,
B.
and
Walker
,
S.
(
2015
).
Fitting linear mixed-effects models using lme4
.
J. Stat. Softw.
67
,
1
-
48
.
Bellis
,
E. D.
(
1962
).
The influence of humidity on wood frog activity
.
Am. Midl. Nat.
68
,
139
-
148
.
Brattstrom
,
B. H.
(
1979
).
Amphibian temperature regulation studies in the field and laboratory
.
Am. Zool.
19
,
345
-
356
.
Broly
,
P.
,
Devigne
,
L.
,
Deneubourg
,
J. L.
and
Devigne
,
C.
(
2014
).
Effects of group size on aggregation against desiccation in woodlice (Isopoda: Oniscidea)
.
Physiol. Entomol.
39
,
165
-
171
.
Bulova
,
S. J.
(
2002
).
How temperature, humidity, and burrow selection affect evaporative water loss in desert tortoises
.
J. Therm. Biol.
27
,
175
-
189
.
Bundy
,
D.
and
Tracy
,
C. R.
(
1977
).
Behavioral response of American Toads (Bufo americanus) to stressful thermal and hydric environments
.
Herpetologica
33
,
455
-
458
.
Chapman
,
T. L.
and
Bidwell
,
J. R.
(
2023
).
Behavioral preference for microclimate conditions across elevation in Plethodon montanus
.
J. Exp. Zool. A Ecol. Integr. Physiol.
339
,
967
-
977
.
Claussen
,
D. L.
(
1969
).
Studies on water loss and rehydration in anurans
.
Physiol. Zool.
42
,
1
-
14
.
Dezetter
,
M.
,
Le Galliard
,
J.-F.
and
Lourdais
,
O.
(
2023
).
Behavioural hydroregulation protects against acute effects of drought in a dry-skinned ectotherm
.
Oecologia
201
,
355
-
367
.
Galindo
,
C. A.
,
Cruz
,
E. X.
and
Bernal
,
M. H.
(
2018
).
Evaluation of the combined temperature and relative humidity preferences of the Colombian terrestrial salamander Bolitoglossa ramosi (Amphibia: Plethodontidae)
.
Can. J. Zool.
96
,
1230
-
1235
.
Giacometti
,
D.
and
Tattersall
,
G. J.
(
2024
).
Seasonal variation of behavioural thermoregulation in a fossorial salamander (Ambystoma maculatum)
.
R. Soc. Open Sci.
11
,
240537
.
Giacometti
,
D.
and
Tattersall
,
G. J.
(
2025
).
Seasonal plasticity in the thermal sensitivity of metabolism but not water loss in a fossorial ectotherm
.
Oecologia
207
,
67
.
Giacometti
,
D.
,
Yagi
,
K. T.
,
Abney
,
C. R.
,
Jung
,
M. P.
and
Tattersall
,
G. J.
(
2021
).
Staying warm is not always the norm: behavioural differences in thermoregulation of two snake species
.
Can. J. Zool.
99
,
974
-
983
.
Greenberg
,
D. A.
and
Palen
,
W. J.
(
2021
).
Hydrothermal physiology and climate vulnerability in amphibians
.
Proc. R. Soc. B Biol. Sci.
288
,
20202273
.
Hall
,
F. G.
(
1922
).
The vital limit of exsiccation of certain animals
.
Biol. Bull.
42
,
31
-
51
.
Haynes
,
W. M.
ed. (
2016
).
CRC Handbook of Chemistry and Physics
.
Boca Raton
:
CRC Press
.
Heatwole
,
H.
(
1960
).
Burrowing ability and behavioral responses to desiccation of the salamander, Plethodon cinereus
.
Ecology
41
,
661
-
668
.
Heatwole
,
H.
(
1962
).
Environmental factors influencing local distribution and activity of the salamander, Plethodon cinereus
.
Ecology
43
,
460
-
472
.
Heatwole
,
H.
and
Newby
,
R. C.
(
1972
).
Interaction of internal rhythm and loss of body water in influencing activity levels of amphibians
.
Herpetologica
28
,
156
-
162
.
Hillyard
,
S. D.
(
1999
).
Behavioral, molecular and integrative mechanisms of amphibian osmoregulation
.
J. Exp. Zool.
283
,
662
-
674
.
Hillyard
,
S. D.
and
Willumsen
,
N. J.
(
2011
).
Chemosensory function of amphibian skin: integrating epithelial transport, capillary blood flow and behaviour
.
Acta Physiol.
202
,
533
-
548
.
Hillyard
,
S. D.
,
Von Hoff
,
K. S.
and
Propper
,
C.
(
1998
).
The water absorption response: a behavioral assay for physiological processes in terrestrial amphibians
.
Physiol. Zool.
71
,
127
-
138
.
Jørgensen
,
C. B.
(
1997
).
200 years of amphibian water economy: from Robert Townson to the present
.
Biol. Rev.
72
,
153
-
237
.
Juarez
,
B. H.
,
Quintanilla-Salinas
,
I.
,
Lacey
,
M. P.
and
O'Connell
,
L. A.
(
2024
).
Water availability and temperature as modifiers of evaporative water loss in tropical frogs
.
Integr. Comp. Biol.
64
,
354
-
365
.
Kosmala
,
G. K.
,
Brown
,
G. P.
,
Shine
,
R.
and
Christian
,
K.
(
2020
).
Skin resistance to water gain and loss has changed in cane toads (Rhinella marina) during their Australian invasion
.
Ecol. Evol.
10
,
13071
-
13079
.
Laurin
,
M.
(
2010
).
How Vertebrates Left the Water
.
Berkeley
:
University of California Press
.
Lertzman-Lepofsky
,
G. F.
,
Kissel
,
A. M.
,
Sinervo
,
B.
and
Palen
,
W. J.
(
2020
).
Water loss and temperature interact to compound amphibian vulnerability to climate change
.
Glob. Change Biol.
26
,
4868
-
4879
.
Lillywhite
,
H. B.
(
2006
).
Water relations of tetrapod integument
.
J. Exp. Biol.
209
,
202
-
226
.
Lopez
,
C. H.
and
Brodie
,
E. D.
(
1977
).
The function of costal grooves in salamanders (Amphibia, Urodela)
.
J. Herpetol.
11
,
372
-
374
.
Lüdecke
,
D.
(
2018
).
ggeffects: tidy data frames of marginal effects from regression models
.
J. Open Source Softw.
3
,
772
.
Lüdecke
,
D.
,
Ben-Shachar
,
M. S.
,
Patil
,
I.
,
Waggoner
,
P.
and
Makowski
,
D.
(
2021
).
performance: An R package for assessment, comparison and testing of statistical models
.
J. Open Source Softw.
6
,
3139
.
MacCracken
,
J. G.
and
Stebbings
,
J. L.
(
2012
).
Test of a body condition index with amphibians
.
J. Herpetol.
46
,
346
-
350
.
Marangio
,
M. S.
and
Anderson
,
J. D.
(
1977
).
Soil moisture preference and water relations of the marbled salamander, Ambystoma opacum (Amphibia, Urodela, Ambystomatidae)
.
J. Herpetol.
11
,
169
-
176
.
Mitchell
,
A.
and
Bergmann
,
P. J.
(
2016
).
Thermal and moisture habitat preferences do not maximize jumping performance in frogs
.
Funct. Ecol.
30
,
733
-
742
.
Ortega
,
Z.
,
Ganci
,
C. C.
and
Rivas
,
M. L.
(
2023
).
Thermoregulation and hydric balance in amphibians
. In
Evolutionary Ecology of Amphibians
(ed.
G.
Moreno-Rueda
and
M.
Comas
), pp.
103
-
119
.
Boca Raton
:
CRC Press
.
Ortega-Chinchilla
,
J. E.
,
Cabanzo-Olarte
,
L. C.
,
Vaz
,
R.
,
Lamadrid-Feris
,
F.
,
Bevier
,
C. R.
and
Navas
,
C. A.
(
2023
).
Behavioral models of hydrothermal regulation in anurans: a field study in the Atlantic Forest, Brazil
.
Biotropica
55
,
329
-
338
.
Petranka
,
J. W.
(
1998
).
Salamanders of the United States and Canada
.
Washington, DC
:
Smithsonian Institution Press
.
Pough
,
F. H.
,
Taigen
,
T. L.
,
Stewart
,
M. M.
and
Brussard
,
P. F.
(
1983
).
Behavioral modification of evaporative water loss by a Puerto Rican frog
.
Ecology
64
,
244
-
252
.
Riddell
,
E. A.
,
Roback
,
E. Y.
,
Wells
,
C. E.
,
Zamudio
,
K. R.
and
Sears
,
M. W.
(
2019
).
Thermal cues drive plasticity of desiccation resistance in montane salamanders with implications for climate change
.
Nat. Commun.
10
,
4091
.
Riddell
,
E. A.
,
Burger
,
I. J.
,
Muñoz
,
M. M.
,
Weaver
,
S. J.
and
Womack
,
M. C.
(
2024
).
Amphibians exhibit extremely high hydric costs of respiration
.
Integr. Comp. Biol.
64
,
366
-
376
.
Sadowski-Fugitt
,
L. M.
,
Tracy
,
C. R.
,
Christian
,
K. A.
and
Williams
,
J. B.
(
2012
).
Cocoon and epidermis of Australian Cyclorana frogs differ in composition of lipid classes that affect water loss
.
Physiol. Biochem. Zool.
85
,
40
-
50
.
Schindelin
,
J.
,
Arganda-Carreras
,
I.
,
Frise
,
E.
,
Kaynig
,
V.
,
Longair
,
M.
,
Pietzsch
,
T.
,
Preibisch
,
S.
,
Rueden
,
C.
,
Saalfeld
,
S.
,
Schmid
,
B.
et al.
(
2012
).
Fiji: an open-source platform for biological-image analysis
.
Nat. Methods
9
,
676
-
682
.
Schmidt-Nielsen
,
K.
and
Schmidt-Nielsen
,
B.
(
1952
).
Water metabolism of desert mammals
.
Physiol. Rev.
32
,
135
-
166
.
Secor
,
S. M.
and
Boehm
,
M.
(
2006
).
Specific dynamic action of Ambystomatid salamanders and the effects of meal size, meal type, and body temperature
.
Physiol. Biochem. Zool.
79
,
720
-
735
.
Shelford
,
V. E.
(
1913
).
The reactions of certain animals to gradients of evaporating power of air: a study in experimental ecology
.
Biol. Bull.
25
,
79
-
120
.
Sherman
,
E.
and
Stadlen
,
S. G.
(
1986
).
The effect of dehydration on rehydration and metabolic rate in a lunged and a lungless salamander
.
Comp. Biochem. Physiol. A
85
,
483
-
487
.
Shoemaker
,
V. H.
and
Bickler
,
P. E.
(
1979
).
Kidney and bladder function in a uricotelic treefrog (Phyllomedusa sauvagei)
.
J. Comp. Physiol. B
133
,
211
-
218
.
Shoemaker
,
V. H.
,
Hillman
,
S. S.
,
Hillyard
,
S. D.
,
Jackson
,
D. C.
,
McClanahan
,
L. L.
,
Withers
,
P.
and
Wygoda
,
M.
(
1992
).
Exchange of water, ions and respiratory gases in terrestrial amphibians
. In
Environmental Physiology of the Amphibians
(ed.
M. E.
Feder
and
W. W.
Burggren
), pp.
125
-
150
.
Chicago
:
The University of Chicago Press
.
Sinsch
,
U.
(
1984
).
Thermal influences on the habitat preference and the diurnal activity in three European Rana species
.
Oecologia
64
,
125
-
131
.
Spight
,
T. M.
(
1967
).
The water economy of salamanders: water uptake after dehydration
.
Comp. Biochem. Physiol.
20
,
767
-
771
.
Spight
,
T. M.
(
1968
).
The water economy of salamanders: evaporative water loss
.
Physiol. Zool.
41
,
195
-
203
.
Spotila
,
J. R.
(
1972
).
Role of temperature and water in the ecology of lungless salamanders
.
Ecol. Monogr.
42
,
95
-
125
.
Székely
,
D.
,
Cogălniceanu
,
D.
,
Székely
,
P.
and
Denoël
,
M.
(
2018
).
Dryness affects burrowing depth in a semi-fossorial amphibian
.
J. Arid Environ.
155
,
79
-
81
.
Takei
,
Y.
(
2015
).
From aquatic to terrestrial life: evolution of the mechanisms for water acquisition
.
Zool. Sci.
32
,
1
-
7
.
Tattersall
,
G.
(
2007
).
Skin breathing in amphibians
. In
Endothelial Biomedicine
(ed.
W. C.
Aird
), pp.
85
-
91
.
Massachusetts
:
Harvard University Press
.
Titon
,
B.
and
Gomes
,
F. R.
(
2017
).
Associations of water balance and thermal sensitivity of toads with macroclimatic characteristics of geographical distribution
.
Comp. Biochem. Physiol. A Mol. Integr. Physiol.
208
,
54
-
60
.
Toledo
,
R. C.
and
Jared
,
C.
(
1993
).
Cutaneous adaptations to water balance in amphibians
.
Comp. Biochem. Physiol. A
105
,
593
-
608
.
Tracy
,
C. R.
(
1976
).
A model of the dynamic exchanges of water and energy between a terrestrial amphibian and its environment
.
Ecol. Monogr.
46
,
293
-
326
.
Tracy
,
C. R.
,
Christian
,
K. A.
,
Betts
,
G.
and
Tracy
,
C. R.
(
2008
).
Body temperature and resistance to evaporative water loss in tropical Australian frogs
.
Comp. Biochem. Physiol. A Mol. Integr. Physiol.
150
,
102
-
108
.
Wickham
,
H.
(
2016
).
ggplot2: Elegant graphics for data analysis
.
New York
:
Springer-Verlag
.
Wolcott
,
T. G.
and
Wolcott
,
D. L.
(
2001
).
Role of behavior in meeting osmotic challenges
.
Am. Zool.
41
,
795
-
806
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information