Skip to main content

Advertisement

Effects of reproduction and environmental factors on body temperature and activity patterns of wolverines

Article metrics

Abstract

Background

Mammals in the far north are exposed to extreme seasonal changes in environmental conditions, such as temperature and photoperiod, which have notable effects on animal physiology and behaviour. The wolverine (Gulo gulo) is a carnivore with a circumpolar distribution and well-adapted to extreme environmental conditions. Still, ecophysiological studies on free-ranging wolverines are lacking. In this study, we used abdominally implanted body temperature loggers in combination with GPS collars with acceleration sensors on 14 free-ranging wolverines in northern Sweden to study daily and seasonal variation in body temperature and activity patterns. We used generalized additive mixed modelling to investigate body temperature patterns over time and Lomb-Scargle periodogram analysis to analyse circadian rhythms.

Results

We found that wolverines have an average core body temperature of 38.5 ± 0.2 °C with a daily variation of up to 6 °C. Body temperature patterns varied between reproductive states. Pregnant females showed a distinct decrease in body temperature during gestation. Wolverines were active both in day and night, but displayed distinct activity peaks during crepuscular hours. However, body temperature and activity patterns changed seasonally, with a gradual change from a unimodal pattern in winter with concentrated activity during the short period of day light to a bimodal pattern in autumn with activity peaks around dusk and dawn. Wolverines were less likely to display 24-h rhythms in winter, when hours of day light are limited.

Conclusions

The combination of different biologging techniques gave novel insight into the ecophysiology, activity patterns and reproductive biology of free-ranging wolverines, adding important knowledge to our understanding of animals adapted to cold environments at northern latitudes.

Background

Knowledge of an animal’s body temperature (hereafter: Tb) profile provides information on thermoregulation, physiology and behaviour, as well as insight into animals` response to changing environmental conditions [1]. Despite large thermal fluctuations, endothermic species are able to regulate Tb within a species-specific narrow range and hence maintain relatively stable core Tb (homeothermy). In mammals, species-specific normal Tb varies from 30 to 40 °C [2] and often displays daily oscillations of 1–4 °C [3]. However, the greatest daily amplitude reported for a free-living large mammal that does not use torpor (Aardvark [Orycteropus afer]) was 8.6 °C, during reduced prey availability [4]. Changing environmental conditions not only impact seasonal and daily patterns of Tb in mammals [5,6,7,8] but also the circadian rhythmicity of activity [9]. Circadian rhythms are daily oscillations, which are at least partly endogenously controlled and recur in approximately 24-h intervals [10]. For instance, several species in the far north lack circadian rhythmicity in activity in summer, when the sun never sets [11, 12], indicating the importance of the light-dark cycle as a Zeitgeber to entrain circadian rhythms. Still, Arctic ground squirrels (Spermophilus parryii) maintain circadian rhythms under constant light conditions [13]. Studies of animals in laboratory conditions suggest that daily rhythms in Tb and activity are highly correlated [14,15,16], but this relationship is largely unexplored in free-ranging animals (but see [17, 18]). Studies in laboratory conditions have also found a profound influence of pregnancy on the physiology of mammals and Tb in particular [19,20,21]. This is expected to be present under natural conditions as well. For instance, pregnant female lions (Panthera leo) exhibit gestational hypothermia and decrease their Tb by 1.3 °C during gestation [22].

Species that are adapted to cold climates, such as the wolverine (Gulo gulo; [23,24,25]), may be particularly vulnerable to rising ambient temperatures associated with global warming [25,26,27]. The need for a framework incorporating ecology and physiology of cold-adapted mammals to identify vulnerability to winter climate change has been emphasized [28], as well as a need for long-term measurements of physiology of mammals in cold environments [29]. Current information on wolverine thermal physiology is limited to a single study of 4 captive individuals, where Tb was measured 4 times a day for 82 days during constant daylight and darkness [30]. Thus, information on year-round variation of wolverine Tb in their natural environment and light conditions is lacking. While wolverine reproductive biology has been studied in harvested [31], captive [32] and radio-collared animals [33,34,35,36], information on physiological characteristics during reproduction is lacking. Female wolverines display embryonic diapause [31, 37] after mating (April – August) [37] with implantation occurring from December to February. Literature suggests a gestation period of 30–50 days [32] and peak parturition from February to mid–March [36], which is earlier than for other northern non-hibernating carnivores [25].

In this study, we used intraperitoneally implanted Tb loggers to record core Tb in wolverines in northern Sweden (Fig. 1). In addition to Tb loggers, wolverines were equipped with global positioning system (GPS) collars with acceleration sensors to record locomotor activity. The combination of biotelemetry and biologging helps to link individual behaviour with physiology and energy status, while recognizing individual variation. This provides key information for conservation efforts of endangered species that are difficult to study with conventional methods [38] and are affected by human activities [39], such as the wolverine is in Sweden [40]. Wikelski and Cooke [41] proposed that conservation strategies would benefit from a better understanding of the physiological responses of organisms to a changing environment.

Fig. 1
figure1

The location of the study area in northern Sweden where wolverines where monitored from 2011 to 2014

The objectives of this study were to describe seasonal changes in wolverine Tb and to explore how reproductive state and environmental conditions influence Tb, as well as activity. We hypothesize that wolverines will lose circadian rhythmicity in Tb and activity in arctic winter and summer. Our aim was to increase the understanding of physiological and behavioural characteristics of a species that is considered highly adapted to cold climate (Fig. 2c).

Fig. 2
figure2

Daily averaged (a) body temperature (Tb, [°C]) and (b) activity of each individual wolverine (grey line) and all individuals together (pregnant individuals: black dots, non-pregnant individuals: hollow circles) and (c) daily averaged ambient temperature from 2011 to 2014 [°C] (black lines) and day length [hours] (orange line)

Methods

Study area

The study area encompasses approximately 7000 km2 and is located in and around the Laponia UNESCO World Heritage site (Kvikkjokk 67°00′N, 17°40′E, Fig. 1). The area is characterized by deep valleys at about 300 m above sea level (m.a.s.l.) and high mountainous plateaus of bare rock and glaciers with peaks up to 2000 m a.s.l. Tundra represents 57% of the study area and forest 36% [42]. Valleys are dominated by mountain birch (Betula pubescens), Scots pine (Pinus sylvestris), and Norway spruce (Picea abies), while mountain birch forms the tree line at 600–700 m a.s.l [43]. The climate is continental with distinct seasons and the ground is usually snow-covered from November to late May, but with large altitudinal variation. Monthly mean ambient temperatures range from − 12 °C in January to 12 °C in July (www.smhi.se; weather station Mierkenis; 66°68′ N, 16°6′E) and minimum day length (1.5 h) can be detected in December, whereas maximum day length of 24 h can be seen when the sun never sets from end of May until mid-July.

Data collection

Adult wolverines were captured by darting from helicopter in May, June and October 2011–2014. Juveniles were located by radiotracking the mother and subsequently captured on the ground [34, 44]. For individuals not captured as juveniles (i.e. of known age, n = 5), approximate age was estimated based on reproductive history [35] and/or genetic relatedness to individuals with known age [45]. Animals were handled, implanted with Tb loggers (Tb loggers: DST Centi, Star Oddi, Gardabaer, Iceland, 46 × 15 mm, 19 g) and VHF abdominal implants (IMP/150/L [21 g] and IMP/400/L [~ 95 g], Telonics, Mesa, Arizona, USA), and adults were fitted with GPS collars (Vectronics-aerospace, Berlin, Germany) according to established protocols [44, 46]. Adult wolverines were darted from helicopter with 4 mg medetomidine (Zalopine® or Domitor®, Orion Pharma Animal Health, Turku, Finland) + 100 mg ketamine (Narketan®, Chassot, Dublin, Irleand). Juveniles were immobilized with 0.1 mg/kg medetomidine + 5 mg/kg ketamine (i.m.). For further information on surgical procedures see Additional file 4. Tb loggers were individually calibrated by the manufacturer over the range of 5 °C to 45 °C and recorded Tb with a guaranteed accuracy of ±0.1 °C 1 year post calibration. Loggers were programmed to record Tb every 5, 10 or 15 min. All GPS collars were equipped with activity sensors measuring acceleration continuously 6–8 times per second in 2 orthogonal directions (X and Y axis). Acceleration values were averaged for each direction for a recording interval of 5 min, ranging from 0 to 255 for each axis and the sum was used to present overall activity, resulting in values from 0 to 510 [47]. Animals were recaptured in order to retrieve Tb loggers and GPS collars, and to downloaded data. In total, Tb data were obtained from 14 wolverines (13 females and 1 male; 4 months to 10 years old), after 2–24 months of data collection. For 10 of these wolverines, activity data was obtained for the same time period as Tb (9 females and 1 male). Data from the week of capture and fever events (daily mean Tb > = 40 °C) were excluded from all analyses. Three of the animals were found dead (2–16 months after capture) and data from the day of death were excluded. We identified 4 occasions in 4 individuals when when Tb dropped to unrealistic low Tb values (< 34 °C), followed by immediate return to baseline values, indicating Tb logger failure. These observations were identified and replaced by the average of the 3 previous and subsequent observations. Data for these four individuals were therefore included in analysis of daily mean Tb but days, with replaced data were excluded from analyses with daily variation as response variable. In total, we used 1,027,345 Tb measurements from 14 wolverines in 32 individual years and 922,673 activity measurements from 10 wolverines in 18 individual years.

Reproductive state

Female reproductive state was defined as either pregnant or non-pregnant based on daily mean Tb. Females were defined as pregnant if they showed a distinct pattern of continuously decreasing Tb in the beginning of the year (Additional file 1). This distinct pattern of decreasing Tb during gestation was evident in at least one breeding season for all 8 females older than 2 years, and in 10 out of 13 (77%) individual breeding seasons. Denning was documented subsequent to these 10 pregnancies in 4 cases and cubs were observed in June in 3 cases (J. Persson, unpublished data). Pregnancy was not detected in any of the 2 years-old females, although females can reproduce at this age [35]. The date of implantation was determined as the day when the trend of decreasing Tb began, using changepoint analysis for daily mean Tb in the R package ‘changepoint’ [48]. This method detects multiple change points in a time series using a pruned exact linear time (PELT) algorithm, which identifies the maximum number of segments a time series can be split into [49], and has been applied for changes in Tb [50]. Date of parturition was determined as the day when mean Tb showed a sudden increase (Additional file 1), resulting in values comparable to normal Tb, reached within a few days. Additionally, activity and location data from GPS collars were used to assess date of parturition. It was determined as the day in which daily mean activity showed a sudden decrease (Additional file 1) followed by low values over several days, and when location data showed that a female remained in a very restricted area for several consecutive days. For non-pregnant females, corresponding gestation period was determined as the period between the average date of implantation and parturition of pregnant females.

Statistical analysis

To investigate the annual variation in wolverine Tb and test for differences between reproductive states throughout the year we fitted Generalized additive mixed models (GAMM) with the function “bam” from the R package “mgcv” [51, 52]. GAMMs are particularly suitable for handling non-linear relationships between response and predictor and handling autocorrelation [53]. We used daily mean Tb as response variable, added a fitted smooth term for day of the year (DOY, 1–366), an ordered factor for reproductive state (pregnant vs. non-pregnant) and an ordered-factor-smooth-interaction term for day of the year multiplied by reproductive state. We included a random intercept and slope for each wolverine ID and added an autoregressive model (AR1) structure to account for detected residual temporal autocorrelation. We applied the “gam.check” function of the “mgcv” package to choose adequate basis dimension of k [51]. The difference of daily mean Tb between reproductive states was considered significant on days when the 95% confidence interval of the modelled difference did not overlap zero.

To investigate annual heterothermy we calculated the daily variation in Tb, defined as the individual range in wolverine Tb (individual daily maximum Tb – minimum Tb). We fitted GAMMs on the individual daily variation in Tb as response variable with the same combination of smooth terms, random structure and AR1 structure as in the model for annual variation in daily mean Tb.

To investigate the daily relationship between Tb, activity, photoperiod and ambient temperature (Ta) graphically, we averaged these parameters hourly for every individual, but excluded pregnant females during the time of gestation. We used a gam smoother with cubic spline function in the R package ggplot [54] on these hourly values averaged over each season to interpret visually how daily patterns in Tb and activity co-vary with Ta and photoperiod. Seasons were defined according to the Swedish Meteorological and Hydrological Institute (SMHI) based on ambient temperature climate indicator as the following: Winter (December – February), Spring (March – May), Summer (June – August), Autumn (September– November) [55]. Hourly Ta data were obtained from the SMHI weather station in Mierkenis, Sweden. Sunrise, sunset, dawn, dusk and day length for the study area were calculated using the R packages "rgeos" and "maptools" [56, 57] and defined as followed: sunrise and sunset are the points at which the top edge of the sun reaches the horizon (i.e. when the top of the sun < 1° below the horizon), dusk and dawn were defined as the onset and end of civil twilight (i.e. when the sun is 6° below the horizon in the morning and evening, respectively).

To investigate periodicity of circadian rhythms in Tb and activity we used Lomb-Scargle periodogram analysis in the R package "lomb" [58]. This method can detect periodicity in irregular spaced time series with missing values [58]. Analyses were performed on the raw data on a 15-day rolling window, i.e. we defined the periods by the DOY, selected 7 days on either side, performed the analysis and moved one DOY ahead. We tested for presence and significance of rhythms between 2 and 30 h, selected the highest significant peak and rounded it to the full hour. Significant periodicity was recognized when peaks in the Lomb-Scargle periodogram exceeded the 95% confidence limit. We determined the two most frequent period lengths (24 or 12 h), displayed by wolverines over the year and converted them into a binomial variable with 1 representing an approximate 24 or 12 h rhythm that day and 0 representing lack of a 24 or 12 h rhythm. To investigate if the probability of wolverines displaying a 24 or 12 h rhythm in Tb and activity changes over the course of a year, we fitted separate binomial GAMMs on the presence of the rhythm in question (binomial variable) as response variable. We added a smooth term for DOY and a random intercept and slope for each wolverine ID. We checked for autocorrelation in each model and decided that it was in an acceptable range and including an autocorrelation structure not necessary (Additional file 3).

Means are presented with associated standard deviation (mean ± SD) in the following text, if not otherwise stated.

Results

Annual variation in body temperature

On an annual basis, daily mean Tb for non-pregnant females (n = 6) was 38.5 ± 0.2 °C, and the range of individual daily variation was 0.9–6.0 °C (mean: 2.5 ± 0.7 °C). A steady decrease in daily mean Tb, indicating gestational hypothermia (Fig. 2a), revealed 10 pregnancies for 8 females > 2 years old and that all adult females were pregnant at least once. These 8 females were considered pregnant during 10 individual breeding seasons in further statistical analysis. Pregnancy was not revealed in either of the 2 – year old females (n = 2).

Non-pregnant wolverines showed a non-linear trend of increasing and decreasing daily mean Tb over the year. Lowest predicted values were found on day (DOY) 42 with 38.3 °C (95% CI: 38.2–38.4 °C), followed by an increase with predicted peak Tb on DOY 196 with 38.6 °C (95% CI: 38.5–38.7 °C) and a subsequent decrease towards the end of the year (Fig. 3). Daily mean Tb of pregnant females (n = 8) was significantly lower than for non-pregnant females from DOY 12–81 (i.e. modelled maximum difference of − 0.5 °C [95% CI: − 0.4 – − 0.6 °C] on DOY 41) but did not differ for the rest of the year (Fig. 3). At the day before parturition, pregnant females had a daily mean Tb of 37.5 ± 0.2 °C, which presented a decrease by 0.8 ± 0.3 °C when compared to Tb in December.

Fig. 3
figure3

Daily mean body temperature (Tb) [°C] for eight pregnant (black hollow circles) and six non-pregnant wolverines (red hollow circles) over the year (day of the year 1 = 1st January). Fitted values from GAMM model are presented for pregnant (black lines) and non-pregnant individuals (red lines). On the bottom, the black line shows the average modelled difference in daily mean Tb [°C] with 95% CI (dashed lines) between pregnant and non-pregnant wolverines, with a red line at 0 °C difference in the background as a reference level. The grey ribbon in the background shows the time span in which Tb of pregnant females is significantly different from Tb of non-pregnant wolverines

Daily variation in body temperature and activity

Similar to daily mean Tb, daily variation in Tb of pregnant and non-pregnant wolverines differed significantly from each other during the time of gestation, but showed similar patterns throughout the rest of the year. Pregnant females had significantly lower daily variation than non-pregnant individuals from DOY 7–70 with a modelled maximum difference of − 0.3 °C (95% CI: − 0.1 – − 0.5 °C) at DOY 40 (Fig. 4). Daily variation in Tb of non-pregnant females increased in the beginning of the year with a peak of 2.9 °C (95% CI: 2.8–3.0 °C) on DOY 85, followed by a steep decrease during summer with lowest value of 2.1 °C (95% CI: 2.0–2.2 °C) at DOY 251 and a second increase towards the end of the year (Fig. 4).

Fig. 4
figure4

Daily variation in body temperature (Tb) [°C, individual daily mamimum Tb – minimum Tb] for eight pregnant (black hollow circles) and six non-pregnant wolverines (red hollow circles) over the year (day of the year 1 = 1st January). Fitted values from GAMM model are presented for pregnant (black lines) and non-pregnant individuals (red lines). On the bottom, the black line shows the average modelled difference in daily variation in Tb [°C] with 95% CI (dashed lines) between pregnant and non-pregnant wolverines, with a red line at 0 °C difference in the background as a reference level. The grey ribbon in the background shows the time span in which daily variation in Tb of pregnant females is significantly different from Tb of non-pregnant wolverines

Daily patterns of both Tb and activity varied seasonally. In winter (December – February) wolverines showed an unimodal pattern in Tb. Activity peaked around the same time as Tb but drops faster than Tb (Fig. 5, Winter). During this period there was, on average, only 2–4 h of daylight and little variation in Ta during the day and wolverines concentrated their activity between dawn and midday. In spring (March – May) wolverines were least active during the time of daylight with lowest values in the afternoon but were active on a constant level during night. Tb varied little but followed the same trend (Fig. 5, Spring). In summer (June – August) Tb and activity followed the same trend as in spring. Highest values could be detected just after sunrise, followed by a gradual decrease with lowest values in the late afternoon, which coincided with highest Ta, followed by gradual increase towards the end of the day (Fig. 5, Summer). In autumn (September – November) wolverines showed a bimodal pattern in activity with highest values just after sunrise, followed by lowest values in late afternoon and another smaller peak just after sunset. Tb follows the same pattern with attenuated amplitude (Fig. 5, Autumn).

Fig. 5
figure5

Seasonal patterns in body temperature (Tb, [°C]) and activity patterns found in wolverines 2011–2014 in Sarek study area. Hourly and seasonal mean Tb (orange line), activity (grey line) and ambient temperature (Ta, [°C], blue line) are presented with 95% CI in shaded grey, applying a cyclic cubic regression spline smoother for time of the day. Grey lines in the background represent individual hourly and seasonal mean values in Tb and activity. Changes in photoperiod are indicated at the bottom of each plot (black bar = darkness, grey bar = dawn/dusk, white bar = sunlight)

Circadian rhythms in Tb and activity

In total, 24 h circadian rhythms were present in 54 and 50.6% of the time for Tb and activity, respectively. Significant 12 h rhythms in Tb were present in 23.5% and in activity for 23.8% of the time. This means that the majority of rhythms displayed by wolverines over the year consisted of 24 and 12 h rhythms (Tb: 77.5%, activity: 74.4%). The probability that wolverines displayed either of the two rhythms changed annually (Fig. 6). The highest probability that wolverines display a 24 h circadian rhythm in Tb was detected in April on DOY 111 with 61% [CI: 53–67%] (Fig. 6a). The highest probability for a 24 h circadian rhythm in activity was detected in July on DOY 194 with 67% [CI: 54–77%] (Fig. 6b). In January wolverines displayed the lowest probability of 24 h rhythm in both Tb (39% [CI: 32–47%]) and activity (36% [CI: 28–45%]). The probability of 12 h rhythms in both, Tb and activity followed an inverse trend of the 24 h rhythm and was generally lower (Fig. 6).

Fig. 6
figure6

Prediction plots of probability of 24 h (black line) or 12 h (grey line) rhythm in body temperature (a) and activity (b) over the course of a year for wolverines, monitored from 2011 to 2014 in Sarek study area. Dashed lines represent 95% confidence interval, orange lines represent associated day length and black dots represent binary results of Lomb-Scargle Periodgram analysis (100 = rhythm present, 0 = rhythm absent). Tb: n = 14 individuals, activity: n = 10 individuals

Discussion

We showed that wolverine Tb varies over the year and is influenced by reproductive state. Reproducing females had significantly lower daily mean Tb and daily variation in Tb than non-reproducing females from early January to early March. Seasonally varying Tb and activity patterns suggest that generally wolverines are active both day and night, with peaks during crepuscular hours. The probability of 24 and 12 h rhythms in Tb and activity changes annually and is lowest in winter time.

Mean wolverine Tb was 38.5 ± 0.2 °C, which is comparable to previous observations of Tb in captive wolverines [30]. Wolverines showed a mean daily variation in Tb of 2.5 ± 0.7 °C which can be considered normal daily oscillations for mammals (1–4 °C; [3]). Maximum daily variations in Tb exceeded 4 °C occasionally in 60% of the individual months (Additional file 2), which indicates a high level of heterothermy. The ability to cope with cold temperatures by being locally heterothermic and applying peripheral cooling [59,60,61] is widely accepted but evidence for heterothermy of core Tb in circumpolar species is lacking. Increased heterothermy in wolverines could be caused by regular events of reduced Tb, which in turn may be driven by an inadequate energy intake, as expected in cold environments [29]. Animals, which face inadequate energy intake may benefit from a reduced metabolic demand, accompanied by decreased Tb. We observed the highest daily variation and lowest mean Tb from December – April (Additional file 2). A decrease from 38.5 °C to 35.8 °C Tb in March (Additional file 2) at an average Ta of − 10 °C represents a 5.5% energy saved (Additional file 4). This period is characterized by the lowest ambient temperatures and relatively low and unpredictable food availability for wolverines, as migrating semi-domestic reindeer, the main prey, return, from winter grazing lands in late April and early May [62]. The daily variation in Tb is suggested to provide an index of any compromise experienced by a free-ranging large mammal, reflecting response to challenging environmental conditions or effects of endogenous factors [29]. From November – April most our study area is snow covered and wolverines are exposed to increased risk of mortality, due to poaching facilitated by use of snowmobiles [40, 63]. During this period, wolverines select stronger for steep and rugged habitats that hinders snowmobile use [63, 64], indicating an influence of human disturbance on behaviour and habitat use. Furthermore, human disturbance can have direct effects on the physiology of individuals by a rise in Tb, heart rate and/or metabolic rate [65,66,67]. The potential effect of human disturbance on both behaviour and physiology, suggests that snow season may represent the time of the year wolverines experience the biggest compromise, expressed by heterothermy.

Moreover, we found that the pattern of daily mean Tb of female wolverines differed between reproductive states, as reproductive females had significantly lower Tb from day 12–81 than non-reproducing females. Decreased Tb during gestation has been described in both domestic [20, 21] and wild mammals [5, 22, 68, 69]. Decreasing Tb during gestation in pregnant female brown bears is presumably caused by changes in progesterone levels, with maximum progesterone levels during implantation, followed by a decrease during gestation. Maintaining stable Tb during gestation and reducing maternal Tb may minimize hyperthermia events, which promotes fetal development [22]. Fetal Tb has been found to be dependent on and approximately 0.6 °C higher than maternal Tb [21]. Cell division stops at Tb higher than 40 °C [70], which makes a fetus potentially vulnerable to hyperthermia events [21].

Wolverines modified Tb and activity patterns over the year relative to photoperiod and ambient temperature. In general, activity was lowest in afternoon, similar to a previous study in the same area [71]. Furthermore, activity patterns changed with season, suggesting that wolverines may be active both day and night but with distinct peaks during crepuscular hours. Although changes in Tb were slightly delayed compared to changes in activity, both patterns were highly synchronized with peaks around the same time. The probability that wolverines exhibit a circadian rhythm in Tb and activity was lowest during polar nights and increased as the days became longer. This is contrary to our predictions that wolverines in northern Sweden will lose circadian rhythms in both polar night and day. Activity patterns in lynx (Lynx lynx) in northern Scandinavia were not influenced by daylight duration but rather by prey activity patterns [72] and circadian rhythms in sheep were disrupted due to environmental and dietary changes [73]. In late winter in northern ecosystems ungulate mortality increases due to starvation, avalanches or weakened animals [74]. Accordingly, wolverines are seasonally flexible in feeding strategy, shifting from mainly predation in spring and summer to scavenging in winter, as a response to increased carrion supply [75]. Changes in food availability feeding strategy may result in different activity patterns and attenuation of circadian rhythmicity in winter. Hence, photoperiod may not be the only influential factor, entraining circadian rhythms of wolverines in northern Sweden. Furthermore, circadian rhythmicity may also be influenced by a high individual component, which in turn may be influenced by e.g. demographic patterns and human disturbance.

Conclusion

The use of combined biologging techniques shed new light into the physiology, activity patterns and reproductive biology of wolverines. This adds new information on a species, adapted to cold climates, which could not be obtained with other methods. Knowledge on ecophysiological characteristics of species are needed to address further questions of interest, such as potential vulnerability to effects of climate change and human encroachment.

Availability of data and materials

Data used in this study are archived in the Dryad Digital Repository and will be accessible from one year following publication. Pending this time, data are available upon request from the wolverine project leader JP (jens.persson@slu.se).

Abbreviations

CI:

Confidence interval

DOY:

Day of the year

GAMM:

Generalised additive mixed effect model

GPS:

Global positioning system

i.m.:

Intramuscular

m.a.s.l:

Meters above sea level

SD:

Standard deviation

Ta :

Ambient temperature

Tb :

Body temperature

References

  1. 1.

    McCafferty DJ, Gallon S, Nord A. Challenges of measuring body temperatures of free-ranging birds and mammals. Animal Biotelemetry. 2015;3(1):33.

  2. 2.

    Clarke A, Rothery P. Scaling of body temperature in mammals and birds. Funct Ecol. 2008;22(1):58–67.

  3. 3.

    Refinetti R. The circadian rhythm of body temperature. Front Biosci. 2010;15:564–94.

  4. 4.

    Rey B, Fuller A, Mitchell D, Meyer LCR, Hetem RS. Drought-induced starvation of aardvarks in the Kalahari: an indirect effect of climate change. Biol Lett. 2017;13(7):5.

  5. 5.

    Mustonen AM, Asikainen J, Kauhala K, Paakkonen T, Nieminen P. Seasonal rhythms of body temperature in the free-ranging raccoon dog (Nyctereutes procyonoides) with special emphasis on winter sleep. Chronobiol Int. 2007;24(6):1095–107.

  6. 6.

    Hetem RS, Strauss WM, Fick LG, Maloney SK, Meyer LCR, Shobrak M, et al. Variation in the daily rhythm of body temperature of free-living Arabian oryx (Oryx leucoryx): does water limitation drive heterothermy? Journal of Comparative Physiology B. 2010;180(7):1111–9.

  7. 7.

    Signer C, Ruf T, Arnold W. Hypometabolism and basking: the strategies of Alpine ibex to endure harsh over-wintering conditions. Funct Ecol. 2011;25(3):537–47.

  8. 8.

    Turbill C, Ruf T, Mang T, Arnold W. Regulation of heart rate and rumen temperature in red deer: effects of season and food intake. J Exp Biol. 2011;214(6):963–70.

  9. 9.

    Williams CT, Barnes BM, Buck CL. Persistence, Entrainment, and Function of Circadian Rhythms in Polar Vertebrates. Physiology. 2015;30(2):86–96.

  10. 10.

    Halberg F. Chronobiology. Annu Rev. Physiol. 1969;31(1):675–726.

  11. 11.

    Reierth E, Stokkan KA. Activity rhythm in High Arctic Svalbard ptarmigan (Lagopus mutus hyperboreus). Can J Zool-Rev Can Zool. 1998;76(11):2031–9.

  12. 12.

    van Oort BEH, Tyler NJC, Gerkema MP, Folkow L, Blix AS, Stokkan KA. Circadian organization in reindeer. Nature. 2005;438(7071):1095–6.

  13. 13.

    Williams CT, Barnes BM, Buck CL. Daily body temperature rhythms persist under the midnight sun but are absent during hibernation in free-living arctic ground squirrels. Biol Lett. 2012;8(1):31–4.

  14. 14.

    Refinetti R. Phase relationship of the body temperature and locomotor activity rhythms in free-running and entrained rats. Biol Rhythm Res. 1997;28:19–24.

  15. 15.

    Weinert D, Waterhouse J. Diurnally changing effects of locomotor activity on body temperature in laboratory mice. Physiol Behav. 1998;63(5):837–43.

  16. 16.

    Refinetti R. Relationship between the daily rhythms of locomotor activity and body temperature in eight mammalian species. Am J Physiol-Regul Integr Comp Physiol. 1999;277(5):R1493–R500.

  17. 17.

    Hayes SR. Daily activity and body temperature of southern woodchuck (Marmota monax monax) in Northwestern Arkansas. Journal of Mammalogy. 1976;57(2):291–9.

  18. 18.

    Long RA, Martin TJ, Barnes BM. Body temperature and activity patterns in free-living arctic ground squirrels. Journal of Mammalogy. 2005;86(2):314–22.

  19. 19.

    Eliason HL, Fewell JE. Thermoregulatory control during pregnancy and lactation in rats. J Appl Physiol. 1997;83(3):837–44.

  20. 20.

    Fewell JE. Body temperature regulation in rats near-term of pregnancy. Can J Physiol Pharmacol. 1995;73(3):364–8.

  21. 21.

    Laburn HP, Mitchell D, Goelst K. Fetal and maternal body temperatures measured by radiotelemetry in near-term sheep during thermal stress. J Appl Physiol. 1992;72(3):894–900.

  22. 22.

    Trethowan PD, Hart T, Loveridge AJ, Haw A, Fuller A, Macdonald DW. Improved homeothermy and hypothermia in African lions during gestation. Biol Lett. 2016;12(11):4.

  23. 23.

    Hart J. Seasonal changes in insulation of the fur. Can J Zool. 1956;34(1):53–7.

  24. 24.

    Copeland JP, McKelvey KS, Aubry KB, Landa A, Persson J, Inman RM, et al. The bioclimatic envelope of the wolverine (Gulo gulo): do climatic constraints limit its geographic distribution? Can J Zool-Rev Can Zool. 2010;88(3):233–46.

  25. 25.

    Inman RM, Magoun AJ, Persson J, Mattisson J. The wolverine’s niche: linking reproductive chronology, caching, competition, and climate. Journal of Mammalogy. 2012;93(3):634–44.

  26. 26.

    Gilg O, Kovacs KM, Aars J, Fort J, Gauthier G, Gremillet D, et al. Climate change and the ecology and evolution of Arctic vertebrates. Ann N Y Acad Sci. 2012;1249(The Year in Ecology and Conservation Biology:166–90.

  27. 27.

    McKelvey KS, Copeland JP, Schwartz MK, Littell JS, Aubry KB, Squires JR, et al. Climate change predicted to shift wolverine distributions, connectivity, and dispersal corridors. Ecol Appl. 2011;21(8):2882–97.

  28. 28.

    Williams CM, Henry HAL, Sinclair BJ. Cold truths: how winter drives responses of terrestrial organisms to climate change. Biol Rev. 2015;90(1):214–35.

  29. 29.

    Hetem RS, Maloney SK, Fuller A, Mitchell D. Heterothermy in large mammals: inevitable or implemented? Biol Rev. 2016;91(1):187–205.

  30. 30.

    Folk GE, Folk MA, Craighead D. Comparison of body temperatures of least weasels and wolverines. Comparative Biochemistry and Physiology a-Physiology. 1977;58(3):229–34.

  31. 31.

    Wright P, Rausch R. Reproduction in the wolverine (Gulo gulo). Journal of Mammalogy. 1955;36(3):346–55.

  32. 32.

    Mead RA, Bowles M, Starypan G, Jones M. Evidence for pseudopregnancy and induced ovulation in captive wolverines (Gulo gulo). Zoo Biol. 1993;12(4):353–8.

  33. 33.

    Persson J. Female wolverine (Gulo gulo) reproduction: reproductive costs and winter food availability. Can J Zool. 2005;83(11):1453–9.

  34. 34.

    Persson J, Landa A, Andersen R, Segerstrom P. Reproductive characteristics of female wolverines (Gulo gulo) in Scandinavia. Journal of Mammalogy. 2006;87(1):75–9.

  35. 35.

    Rauset GR, Low M, Persson J. Reproductive patterns result from age-related sensitivity to resources and reproductive costs in a mammalian carnivore. Ecology. 2015;96(12):3153–64.

  36. 36.

    Aronsson M. ‘O neighbour, where art thou?’ Spatial and social dynamics in wolverine and lynx, from individual space use to population distribution [Ph.D thesis]. Uppsala: Swedish University of Agricultural Sciences; 2017.

  37. 37.

    Rausch RA, Pearson A. Notes on the wolverine in Alaska and the Yukon Territory. The Journal of Wildlife Management. 1972:249–68.

  38. 38.

    Cooke SJ. Biotelemetry and biologging in endangered species research and animal conservation: relevance to regional, national, and IUCN Red List threat assessments. Endang Species Res. 2008;4(1–2):165–185.

  39. 39.

    Krebs J, Lofroth E, Copeland J, Banci V, Cooley D, Golden H, et al. Synthesis of survival rates and causes of mortality in North American wolverines. J Wildl Manage. 2004;68(3):493–502.

  40. 40.

    Persson J, Ericsson G, Segerstrom P. Human caused mortality in the endangered Scandinavian wolverine population. Biol Conserv. 2009;142(2):325–31.

  41. 41.

    Wikelski M, Cooke SJ. Conservation physiology. Trends Ecol Evol. 2006;21(1):38–46.

  42. 42.

    Mattisson J, Andren H, Persson J, Segerstrom P. Influence of intraguild interactions on resource use by wolverines and Eurasian lynx. Journal of Mammalogy. 2011;92(6):1321–30.

  43. 43.

    Grundsten C. The Laponian area: a Swedish world heritage site. Fälths Tryckeri AB. Stockholm: Swedish Environmental Protection Agency; 1997.

  44. 44.

    Arnemo JM, Evans AL. Fahlman Å. Biomedical protocols for free-ranging brown bears, gray wolves, wolverines and lynx. Evenstad. Norway: Hedmark University College; 2012.

  45. 45.

    Hedmark E, Persson J, Segerstrom P, Landa A, Ellegren H. Paternity and mating system in wolverines Gulo gulo. Wildlife Biol. 2007;13:13–30.

  46. 46.

    Fahlman A, Arnemo JM, Persson J, Segerstrom P, Nyman G. Capture and medetomidine-ketamine anesthesia of free-ranging wolverines (Gulo gulo). J Wildl Dis. 2008;44(1):133–42.

  47. 47.

    Gervasi V, Brunberg S, Swenson JE. An individual-based method to measure animal activity levels: A test on brown bears. Wildl Soc Bull. 2006;34(5):1314–9.

  48. 48.

    Killick R, Eckley IA. changepoint: An R Package for Changepoint Analysis. J Stat Softw. 2014;58(3):1–19.

  49. 49.

    Killick R, Fearnhead P, Eckley IA. Optimal Detection of Changepoints With a Linear Computational Cost. J Am Stat Assoc. 2012;107(500):1590–8.

  50. 50.

    Evans AL, Singh NJ, Fuchs B, Blanc S, Friebe A, Laske TG, et al. Physiological reactions to capture in hibernating brown bears. Conserv Physiol. 2016;4:11.

  51. 51.

    Wood S. Mixed GAM computation vehicle with automatic smoothness estimation. R package vers 1.8–22. 2017.

  52. 52.

    van Rij J, Wieling M, Baayen RH, van Rijn HJRpv. itsadug: Interpreting time series and autocorrelated data using gamms. 2016;2:2.

  53. 53.

    Wood SN. Generalized additive models: an introduction with R. New York, USA: Chapman and Hall/CRC; 2017.

  54. 54.

    Wickham H, Chang W, Wickham MHJCsmRfhcr-powpggpGS. Package ‘ggplot2’. 2013.

  55. 55.

    Kjellstrom E, Boberg F, Castro M, Christensen JH, Nikulin G, Sanchez E. Daily and monthly temperature and precipitation statistics as performance indicators for regional climate models. Climate Research. 2010;44(2–3):135–50.

  56. 56.

    Bivand R, Lewin-Koh N. maptools: Tools for Reading and Handling Spatial Objects. R package version 0.9–2. 2017.

  57. 57.

    Bivand R, Rundel C. rgeos: Interface to Geometry Engine - Open Source (‘GEOS’). R package version 0.3–26. 2017.

  58. 58.

    Ruf T. The Lomb-Scargle periodogram in biological rhythm research: Analysis of incomplete and unequally spaced time-series. Biol Rhythm Res. 1999;30(2):178–201.

  59. 59.

    Arnold W, Ruf T, Reimoser S, Tataruch F, Onderscheka K, Schober F. Nocturnal hypometabolism as an overwintering strategy of red deer (Cervus elaphus). Am J Physiol-Regul Integr Comp Physiol. 2004;286(1):R174–R81.

  60. 60.

    Irving LJP. Heterothermous operation of warm blooded animals. Physiologist. 1959;2(2):18–32.

  61. 61.

    Scholander PF. Evolution of climatic adaption in homeotherms. Evolution. 1955;9(1):15–26.

  62. 62.

    Danell AC, Andren H, Segerstrom P, Franzen R. Space use by Eurasian lynx in relation to reindeer migration. Can J Zool-Rev Can Zool. 2006;84(4):546–55.

  63. 63.

    Rauset GR, Andren H, Swenson JE, Samelius G, Segerstrom P, Zedrosser A, et al. National Parks in Northern Sweden as Refuges for Illegal Killing of Large Carnivores. Conserv Lett. 2016;9(5):334–41.

  64. 64.

    Rauset GR, Mattisson J, Andren H, Chapron G, Persson J. When species’ ranges meet: assessing differences in habitat selection between sympatric large carnivores. Oecologia. 2013;172(3):701–11.

  65. 65.

    Regel J, Putz K. Effect of human disturbance on body temperature and energy expenditure in penguins. Polar Biol. 1997;18(4):246–53.

  66. 66.

    Macarthur RA, Johnston RH, Geist V. Factors influencing heart rate in free-ranging bighorn sheep- physiological approach to the study of wildlife harassment. Can J Zool-Rev Can Zool. 1979;57(10):2010–21.

  67. 67.

    Watts PD, Ferguson KL, Draper BA. Energetic output of subadult polar bears (Ursus maritimus)- resting, disturbance and locomotion. Comparative Biochemistry and Physiology a-Physiology. 1991;98(2):191–3.

  68. 68.

    Williams CT, Sheriff MJ, Schmutz JA, Kohl F, Toien O, Buck CL, et al. Data logging of body temperatures provides precise information on phenology of reproductive events in a free-living arctic hibernator. J Comp Physiol B-Biochem Syst Environ Physiol. 2011;181(8):1101–9.

  69. 69.

    Friebe A, Evans AL, Arnemo JM, Blanc S, Brunberg S, Fleissner G, et al. Factors Affecting Date of Implantation, Parturition, and Den Entry Estimated from Activity and Body Temperature in Free-Ranging Brown Bears. PLoS One. 2014;9(7):10.

  70. 70.

    López-Sáez J, Gimenez-Martin G. González-Fernández AJZfZuMA. Duration of the cell division cycle and its dependence on temperature. 1966;75(3):591–600.

  71. 71.

    Mattisson J, Andren H, Persson J, Segerstrom P. Effects of Species Behavior on Global Positioning System Collar Fix Rates. J Wildl Manage. 2010;74(3):557–63.

  72. 72.

    Heurich M, Hilger A, Küchenhoff H, Andrén H, Bufka L, Krofel M, et al. Activity patterns of Eurasian lynx are modulated by light regime and individual traits over a wide latitudinal range. PLoS One. 2014;9(12):e114143.

  73. 73.

    Wyse CA, Zhang X, McLaughlin M, Biello SM, Hough D, Bellingham M, et al. Circadian rhythms of melatonin and behaviour in juvenile sheep in field conditions: Effects of photoperiod, environment and weaning. Physiol Behav. 2018;194:362–70.

  74. 74.

    Wilmers CC, Stahler DR, Crabtree RL, Smith DW, Getz WM. Resource dispersion and consumer dominance: scavenging at wolf- and hunter-killed carcasses in Greater Yellowstone, USA. Ecol Lett. 2003;6(11):996–1003.

  75. 75.

    Mattisson J, Rauset GR, Odden J, Andren H, Linnell JDC, Persson J. Predation or scavenging? Prey body condition influences decision-making in a facultative predator, the wolverine. Ecosphere. 2016;7(8):14.

Download references

Acknowledgements

Body temperature loggers were provided by the French Space Agency. Captures were conducted by Peter and Einar Segerström.

Funding

The study was funded by the Norwegian Directorate for Nature Management, the Swedish Environmental Protection Agency, the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (Formas), Marie-Claire Cronstedt foundation and World Wildlife Fund for Nature, Sweden.

Author information

This study was conducted within an ongoing project led by JP. JMA, JP and ALE conceived this study and participated in the captures. JMA and ALE deployed and removed the devices. JP and MA compiled activity data and confirmed parturition based on GPS-data. BF compiled the body temperature data and contributed to the design of the modelling approach. AT analysed the data and drafted the manuscript. All authors contributed to writing of the manuscript and approved of the final version.

Correspondence to Alexandra Thiel.

Ethics declarations

Ethics approval

Captures, handling, implantation of body temperature loggers, very high frequency (VHF) implants and fitting of GPS collars were approved by the Swedish Ethical Committee on Animal Research for northern Sweden, Umeå (protocol number A13–09) and the Swedish Environmental Protection Agency.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Graph to illustrate pregnancy determination based on daily mean body temperature and activity. (PNG 23 kb)

Additional file 2:

Table of monthly averaged body temperature of non-pregnant and pregnant wolverines. (DOCX 13 kb)

Additional file 3:

Autocorrelation plots of GAMM models. (DOCX 592 kb)

Additional file 4:

Details about surgery and energy calculations. (DOCX 16 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Thiel, A., Evans, A.L., Fuchs, B. et al. Effects of reproduction and environmental factors on body temperature and activity patterns of wolverines. Front Zool 16, 21 (2019) doi:10.1186/s12983-019-0319-8

Download citation

Keywords

  • Biologging
  • Circadian rhythm
  • Ecophysiology
  • Gestation
  • Heterothermy
  • Northern ecosystem
  • Seasonality

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.