Aquatic habitat use in a semi-aquatic mammal: the Eurasian beaver

Background: Semi-aquatic mammals exploit resources both on land and in water and may require both to meet their habitat requirements including foodand building resources, refuges, and for social interactions with conspecifics. Within this, the specific availability of both terrestrial and aquatic resources is expected to impact individual fitness. Beavers are highly dependent on water for movement and protection from predators. They are central place foragers and mostly forage on woody vegetation near water although aquatic vegetation may also be an important food resource. However, little is known about their use of aquatic habitats. We aimed to address this knowledge gap by dead-reckoning fine-scale movement tracks and classifying fine-scale diving events, which we then related to the spatial distribution of aquatic vegetation and habitat components within the territory. Results: Overall, there was a statistically clear decrease in probability that diving would occur at dawn and with increasing distance from territory borders. In addition, the distance from the lodge at which animals dived decreased through the night and during the spring/early summer. There was strong selection for diving habitats located closer to the riverbank, with stronger selection for these areas being observed in individuals with larger home ranges. We saw a higher selection for diving above clay sediment, and within 150 m from the lodge, presumably because mud and clay sediment tended to be located closer to the lodge than sand and rock sediment. Furthermore, we found a clear selection for diving in the presence of quillwort (Isoetes spp.), shoreweed (Littorella uniflora), and stonewort (Nitella spp.). Selection for these focal species was stronger among subordinate individuals. Individuals with lower body condition dived closer to the beaver lodge, and dives located further from the lodge were associated with high densities of aquatic vegetation. Conclusion: We provide new knowledge on the aquatic habitat use in a semi-aquatic mammal and show how energetic constraints may shape how beavers spatially use the aquatic environment, whereby short and shallow dives appear most beneficial. We show how aquatic habitats may have great importance for both foraging, building materials and safety, and discuss to how they may affect the fitness of individuals.

animal movement is challenging as behavioural ecologists have to juggle with multiple interacting factors [20,21].
A major movement delimitator results from competition, when animals may exclude potential competitors from an area by being territorial [22,23]. As a result, some individuals end up defending habitats that do not have favourable combinations of essential resources, and thus have to trade-off resources within their ultimate habitat acquisition [24][25][26]. Consequently, some species benefit by being generalists rather than specialists, because the consequences of losing a resource are less severe since animals can exploit a greater variety of resources [27][28][29][30]. For example, semi-aquatic mammals have adapted to exploit resources on land and in water, and individuals of this group may express considerable plasticity to meet their habitat requirements according to the available food resources, shelter, and social interactions [31][32][33].
Habitat selection happens at various spatial scales [34] and is described as the use of resources (habitat) in a manner that is disproportionate to their availability [35]. Critically, resource-use may not necessarily be directly proportional to resource availability, but it may also be modulated by other ecological factors such as competition and predation [25,[36][37][38]. Although it has always been a challenge to quantify which habitats animals have available to them, and how much they use them, today, this information is important to inform space-value discussions so that confounding ecological variables affecting habitat value for the animal can be put into context [39][40][41].
In this study, we examine aquatic habitat use in Eurasian beavers (Castor fibre) to identify important characteristics of aquatic habitats within beaver territories and investigate potential differences in aquatic habitat use among individuals. To achieve our aim, we combine sophisticated animal-attached tags (GPS loggers and Daily Diary units), that allow determination of animal behaviour with locations, with a comprehensive assessment of aquatic habitat characteristics.
Technological developments have hugely enhanced what animal biotelemetry can do for us, elucidating, for example, fine-scale spatiotemporal location data on an increasing range of animals across various environments [40][41][42][43]. In particular, tri-axial accelerometers are increasingly being used to study wild animals [44,45], because they allow determination of an individual's behaviour [46]. They have been used to classify behaviour and activity level patterns in beavers, distinguishing seven behaviours with high precision, including swimming and diving [47,48]. These behaviours can be combined with fine-scale animal movement determined by dead-reckoning [49][50][51] to provide information on what animals do in the spaces they inhabit [52]. Altogether, fine-scale information on behaviour and movement, obtained from the use of acceleration and dead-reckoning, respectively, are being used increasingly to study wild animals that are hard to observe directly and/or without bias [44,49,[52][53][54][55][56].
Beavers are socially monogamous, monomorphic, nocturnal mammals that inhabit various freshwater bodies [38,57]. They live in family groups consisting of the dominant breeding pair, kits of the year, and older non-breeding offspring [33,58,59]. Beavers reach sexual maturity during their second winter [59], and give birth to one to five kits in mid-May at northern latitudes [57]. The kits emerge from the lodge in July when they start feeding on their own [33]. At around 2-3.5 years old, beavers tend to disperse from their natal territory to establish their own territory [15,60].
Beavers are central place foragers, mostly foraging near water and their lodges [61][62][63]. Their diet consists mainly of woody vegetation but varies seasonally, and comprised primarily bark from deciduous trees during winter to more nutritiously rich deciduous leaves, aquatic vegetation, and herbaceous plants in spring and summer [64][65][66][67][68][69][70][71][72]. In some areas, aquatic plants may seasonally account for up to 90% of the diet [65,66,73]. Aquatic vegetation may offer some nutritional benefits over terrestrial vegetation, including better digestibility, higher crude protein, and higher sodium and iron content [64,66,74,75]. Low concentrations of secondary compounds might also make aquatic vegetation more palatable [66], but this may vary with species [76]. Seasonally, rhizomes of aquatic plants can provide great nutritional value in winter and spring when plants store nutrients in the rhizomes in preparation for spring growth [77,78]. Diet variation may depend on nutrient content and digestibility of available forage as individual beavers attempt to maximize energy intake over time [66,76,[79][80][81]. Beaver foraging behaviour varies according to environmental factors that affect the distribution of food items [82], but ecological factors such as food plant density, human disturbance, presence of conspecifics, and predator activity may also affect their foraging choices and foraging locations [67,83,84]. No clear dietary differences have been found between sexes, ages, or social ranks in beavers [65,69,85], but several studies indicate that foraging behaviour may differ seasonally as territorial movements vary among individuals [86][87][88]. Furthermore, individuals may be affected by various ecological conditions during their lifetime, such as loss and acquisition of territories [15,89], that can affect their behavioural time-budgets and consequently their body condition, reproduction, and survival [57,90].
Water acts as a refuge for beavers [91] and is specifically used to minimize predation risk when foraging on land [92]. Even though beavers depend on the aquatic environment for movement and safety [83,86,87,93] and may even manipulate the environment to optimize aquatic conditions [82,94], aquatic behaviour and habitat use have not been well studied in these animals. Although beavers use aquatic habitats for foraging [64][65][66]71], protection from predators [92,95], and collecting resources for lodge-and dam-building [96,97], aquatic habitats may be particularly important when terrestrial vegetation is difficult to access, or is of low nutritional quality [38,64,66,75,98]. However, their spatial exploitation of the aquatic components of their territory is poorly understood. Research has found, however, that habitat use may differ between age groups and according to risk levels [87]. Diving behaviour has been studied using accelerometers, which has highlighted a preference for short (< 30 s) and shallow (up to 4 m, but most < 1 m) dives, which indicate some form of aquatic resource selection, although the link between diving and space use is vague [99]. Being only semi-aquatic [cf. 100], beavers may experience a higher cost exploiting aquatic resources [101,102] than a fully aquatic equivalent may do. This may explain why studies have found them to be diving for less than 3% of their nightly activity budget [99]. However, energy requirements have been reported to compare well with more fully aquatic mammals and birds [101].
Using fine-scaled dead-reckoned animal tracks to determine spatial and temporal locations of aquatic dives, we aim to examine important characteristics of aquatic habitat use by Eurasian beavers and investigate potential individual differences. Assuming that dives indicate aquatic habitat use in beavers, we hypothesize that habitat use vary temporally and spatially between individuals of various ages, sexes, social ranks, and by the composition of biotic and abiotic factors within their individual territories.

Study site
Our study site was located at the lower reaches of the river Sauar in Vestfold and Telemark County, southeastern Norway (Fig. 1). The river drains the lake Heddalsvatnet in the north and forms part of the catchment of the lake Norsjø in the south, stretching over approximately 13 km with a width of 45-250 m. The river sections are generally slow-flowing with stable water levels because of natural lakes and man-made impoundments along part of its length [103], although flooding events frequently occur. The river flows through small villages, farmlands and fields interspersed with riparian woodland that comprised mostly Norway spruce (Picea abies), Scots pine (Pinus sylvestris), birch (Betula spp.), grey alder (Alnus incana), aspen (Populus tremula) and mountain ash (Sorbus aucuparia) [62,103]. Fig. 1 a The location of the study site (red square) in Telemark and Vestfold County, Norway. b Overview of study river with random available sites (yellow circles) and identified beaver diving locations (red triangles) within beaver territories Beavers have inhabited the area since the 1920s when they recolonized the rivers [104]. The population is at carrying capacity, as territories of various sizes directly border each other [33]. Territory borders are identified based on scent mound concentrations, sight observations of known patrolling beavers backed up by GPS data. The main river contains ten distinguishable territories inhabited by approximately 32 individuals [33,57]. Predation pressure is low as wolves (Canis lupus) and bears (Ursus arctos) are functionally extinct in the area, and lynx (Lynx lynx) only occur at low densities [84,105].
The study river is part of a larger monitoring project where beavers in the area have been monitored through an extensive capture programme, the Norwegian Beaver Project (NBP), since 1997 [57]. The long-term monitoring project aims to capture all newcomers (kits and dispersers from outside the study site) annually, enabling identification of individuals at later encounters and family group sizes.

Capture and tracking protocol
Individuals were detected from a motorboat using searchlights and captured at night with large divingnets in shallow water or with land-nets [106]. Captured individuals were immobilized in cloth sacks, enabling easy handling without anaesthesia, and identified via microchips (PIT tag) and unique combinations of plastic and metal ear-tags. Beavers were weighed to the nearest 100 g. Body length was measured following the curvature of the spine from nose tip to the base of the tail. Tail length was measured from the base to the tip of the tail, and tail width was measured from edge to edge of the dorsal surface at the midpoint between tail base and tip. Measurements of body length and tail proportions in cm were used to calculate tail fat index ((tail length × tail width)/body length), representing the body condition of beavers [57,59].
Individuals were sexed based on the colour and viscosity of their anal gland secretion [107] and assigned a minimum age based on body mass at first capture [57,108]; minimum 2 years (subadult) when body mass was between 17 to 19.5 kg inclusive, and minimum 3 years (adult) when body mass was above 19.5 kg. Territorial dominance was in most cases attributed to adult territorial residents of each sex. Territorial dominance was verified by eventual dispersal of the alternative candidate, greatest body weight among same-sex group members or lactation in females (large nipples, i.e. > 0.5 cm). Individuals dispersing into a territory were posited to have achieved the dominant breeding position when the previous dominant same-sex individual had disappeared, or evidence outlined above was applicable. Unless proven otherwise, dominant individuals were assumed to maintain their social rank until they died or disappeared from the territory [109].
We captured and equipped nine beavers (five males, four females, Table 1) with GPS loggers (Gipsy-5, Techno Smart) and daily diary units (including accelerometer, magnetometer, thermometer, Wildbyte technologies [52]) in the spring and early summer of 2018 (one beaver) and 2019 (eight beavers). The data loggers were glued onto the fur on the lower back of the beavers, approximately 15 cm above the tail following the spine, and were removed again after 2 to 3 weeks if they had not fallen off by themselves [47]. To extend battery life, GPS loggers were programmed to take a fix position every 15 min between 7 p.m. and 7 a.m. to reduce numbers of unsuccessful GPS fix attempts from within the beaver lodges. Daily diary units logged continuously at 40 Hz. From the GPS positions, we identified territory borders to estimate territory size, expressed in bank length (km) and calculated 95% autocorrelated kernel density estimates (AKDE) to estimate the overall space use (home range) of each beaver [110]. Captured beavers were released near the capture site within their territory after approximately 40 min of handling time [106].

Identification of dives
Accelerometers in tandem with magnetometers can be used in dead-reckoning to accurately predict and reconstruct animals' fine-scale three-dimensional movement paths in space and time by sequentially integrating calculated travel vectors [49,51]. However, the estimated movement track accumulates error and therefore drifts over time, so it needs to be corrected through groundtruthing, e.g. correcting the track according to GPS fixes [49].
We calibrated the daily diary data in the software DDMT (Daily Diary Multiple Trace, Wildbyte Technologies). Using the acceleration and magnetism data, we dead-reckoned the movement track of each beaver in the software Framework4 [111]. The dead-reckoned movement tracks were hereafter corrected using the GPS positions as ground-truthing [49]. GPS positions were filtered to remove positions with horizontal dilution of precision (HDOP) values above five and with less than four available satellites to reduce the effects of imprecise GPS positions [112,113].
To identify diving locations, we divided the dead-reckoned movement tracks into ten second bursts and used the acceleration to assign behavioural activities to each burst based on the acceleration-based behavioural classification model by Graf et al. [47]. The classification model can clearly differentiate the acceleration between seven behaviours: swimming, diving, sleeping, feeding, standing, walking, and grooming. To furthermore filter out Table 1 Overview of tracked individuals with daily diaries and GPS loggers in a population of Eurasian beavers in southeastern Norway potentially falsely classified dives, we focused on diving sections where the previous 10-s bursts were additionally predicted to be swimming.

Assessment of aquatic habitats
We visited and assessed the aquatic habitat of all diving locations as well as random aquatic locations available within the territories of each beaver between June and October 2019 (Fig. 1b). We only sampled sites with a water depth of less than 10 m since vegetation growth here is more abundant due to light conditions [114], whilst also accounting for known diving depths in beavers [99].
For vegetation sampling at each location, sites were sampled using a 1 × 1 m quadrat with an aluminium frame. The quadrat was placed as close to the location as wind and currents allowed. The frame construction had a pyramidal shape, enabling a GoPro camera (GoPro Hero5) to be attached to the top, 0.8 m above the surface, which would keep the quadrat within the camera view. The quadrat was left at the bottom of the site after sediment settled, and pictures and films were recorded. When water depth allowed it, aquatic plants and physical characteristics were recorded in situ using an aqua scope. Aquatic plants were collected with a rake when identification required closer inspection. Plant species identification followed the database of Artsdatabanken [115]. Species abundance was quantified as coverage in percentage, rounding to the nearest 5% [116]. Plants with less than 5% cover were registered as 1% per species.
We categorized each site according to physical characteristics [water depth and sediment type (clay, mud, sand, and rock)], spatial characteristics (distance to riverbank, beaver lodge, and territory border), and characteristics of the aquatic vegetation. We characterized the aquatic vegetation by cover and species richness (number of species) to evaluate importance of quantity and diversity, respectively.

Statistical analysis
We used generalized linear mixed-effects models (GLMMs) with Poisson distribution, log link, and beaver ID as a random effect to investigate how the number of identified dives per night varied between individuals and components of the territories. We analysed the effects of sex (male, female), social rank (dominant, subordinate), age (years), body size (body mass, body length, and tail fat index), and home range size (95% AKDE), and the effects of territory size (bank length in km), mean water depth (m) and mean vegetation cover and species richness.
Using GLMMs with Bernoulli distribution, logit link, with beaver ID and tracking night as random effects, we also investigated how hourly diving probability (1 = a diving event was identified within a given hour, 0 = no diving events were identified within a given hour) varied through the night and between individuals. We analysed the effects of spatiotemporal variables (date, hour of the night, distance to riverbank, lodge, and territory border), individual differences (sex, social rank, age, body sizes, and home range size), and between components of the territories (territory size, mean water depth, mean vegetation cover, and mean species richness).
We used a detrended correspondence analysis (DCA) from the R package 'vegan' v. 2.5-7 [117] to assess the general distribution pattern of the aquatic vegetation in the study river in relation to physical (water depth and sediment type), spatial (distance to riverbank, beaver lodge, and territory border) and vegetation characteristics (vegetation cover and species richness) of the sites. A matrix including aggregated abundance per species was used for the DCA. However, species that only occurred in five sites or less were removed from the matrix. The correlation between the aquatic species compositions and physical, spatial, and vegetation characteristics of the sites were assessed by passively fitting them to the ordination (permutations = 999). From the ordination, comparing species with use, we could furthermore identify species that appeared to be of potential importance to the beavers when diving. For the subsequent resource selection functions, we included variables for the number of focal species present at a site and vegetation cover of these focal species.
We investigated the aquatic habitat selection within territories using GLMMs with Bernoulli distribution, logit link, and beaver ID as a random effect (1 = diving site, 0 = random available site within the territory) [34,35]. We analysed whether aquatic habitats located at varying water depths, sediment types, distances to riverbank, beaver lodge, and territory border, or with varying vegetation cover (overall and focal species), and species richness (overall and focal species) were used more than was generally available in the territories [35,118]. Additionally, we analysed the variations in diving selection among individuals in univariate models weighted by number of identified dives by fitting the resource selection function to each individual [119]. This enabled us to analyse how the individual diving selection coefficients varied between beavers of different age, sex, body size, and social rank that furthermore inhabit territories of different size and with varying amount of available water depth, vegetation cover (overall and focal species), and species richness (overall and focal species) [13,25,119]. As dives may have different purposes according to their spatial location, we furthermore investigated how log transformed distance to the lodge and riverbank varied temporally (date and hour of the night), by environmental characteristics (depth, sediment, vegetation cover, species richness), and among individuals of different age, sex, body size, social rank, territory size, and home range by using GLMMs with gaussian distribution and beaver ID as a random effect. To examine characteristics associated to the proximate area of the beaver lodge, we furthermore analysed how dives within 150 m of the beaver lodge varied according to the above variables using univariate GLMMs with Gaussian distribution and beaver ID as a random effect.
In all analyses, a list of candidate models was created using ecologically relevant combinations of fixed effects to account for variability in endogenous (such as sex, age, and social rank) and exogenous factors (such as territory size, vegetation composition) that may be important in describing the ecology of beavers (Fig. 2). Because of the sample size, individual effects (sex, age, social rank, territory size, and home range size) should be interpreted with care as they only imply possible ecological effects that should be investigated with more individuals in future studies. We included spatiotemporal interactions (between hour and distance to riverbank, lodge, and territory border, respectively) in the analysis for diving probability, but excluded interactions in all other analyses because of the limited sample size. Individual selection coefficients were similarly analysed in univariate models because of the limited sample size. The fixed effects used in all analyses were not correlated (Pearson r coefficients less than 0.5) and variance inflation factor values were less than 3 [120].
Model selection was based on Akaike's information criterion corrected for small sample size [121], and carried out using the R packages 'glmmTMB' v. 1.0.2.1 [122] and 'MuMIn' v. 1.43.17 [123]. The most parsimonious models within ΔAICc < 2 were chosen as the best models to describe the variation [121,124]. In each model, variables that included zero within their 95% confidence interval (CI) were considered uninformative and reported as unclear effects [124]. The best models were visually validated using the R package 'DHARMa' v. 0.4.1 [125] to plot standardized model residuals against the fitted values [120] and, when relevant, furthermore checked for zero-inflation. Top candidate models for all analyses can be found in the supplemental material (Additional file 1). All analyses were conducted in R 4.0.3 [126].

Results
Nine beavers were tracked with data loggers (Wildbyte technologies, Daily Diaries [52]) and GPS loggers (Techno Smart) affixed to the lower back for a total of 77 nights. Identified diving events lasted between 10 and 110 s, with the majority (80%) lasting 10 s or less. We identified on average (mean ± SD) 9.5 ± 3.1 dives per night for each beaver. We found no clear differences in the number of dives per night between males and females (10.9 ± 3.3 and 8.8 ± 3.0, respectively), among dominants and subordinates (12.9 ± 3.6 and 6.4 ± 2.5, respectively) or as a function of age, body size and tail fat index (Additional file 1). Furthermore, territory size, home range size, mean water depth in territory, and mean vegetation cover and species richness in territory did not have a clear effect on the number of identified dives per night (Additional file 1).

Diving probability
Hourly diving probability varied through the night, decreasing over the final hours of the night (Fig. 3) and with increasing distance from the territory borders (Table 2, Fig. 3). We found no differences in hourly diving probability between sexes, social ranks, age, or relative to body size and tail fat index. Furthermore, we found no clear effect of date, distance to riverbank, distance from the lodge, territory size, home range size, mean water depth in territory and mean vegetation diversity within the territory on nightly diving probability ( Table 2, Additional file 1).

Characterization of aquatic vegetation
The DCA ordination described up to 32.3% of the variation in the aquatic vegetation composition within the beaver territories with DCA1 and DCA2 describing 9.6 and 9.4% of the variation, respectively ( Table 3). The aquatic vegetation within the territories showed great variation and differed clearly with increasing water depth, vegetation cover, species richness, and sediment type (Fig. 4). Diving sites and random available sites within the beaver territories clearly differed in species composition. Diving sites were especially associated with varying Fig. 3 The predicted relationship ± 95% confidence interval between hourly diving probability, time of night, and distance from territory border among nine individuals in a Eurasian beaver population in southeastern Norway. Yellow boxes represent time of sunset and sunrise through the tracking period. Red boxes represent time of dusk and dawn through the tracking period    (Fig. 4). Furthermore, sites were predominantly more used for diving with decreasing water depths, increasing species richness, and decreasing vegetation cover. A list of all aquatic species can be found in Additional file 1.

Resource selection of aquatic habitats within the territory
We found a clear diving selection for aquatic locations closer to the riverbank but found no clear diving selection as a function of distance from lodge, distance from territory border, or water depth (Table 4, Fig. 5a, Additional file 1). Diving selection probability varied as a function of sediment type. Locations having either mud, sand, or rock sediment were less selected than locations with clay sediment ( Table 4, Fig. 5b). Furthermore, diving selection probability increased when several of either quillwort, shoreweed, and stonewort were present, but we found no effect of vegetation cover, overall species richness, or vegetation cover of the focal species on the diving selection probability (Table 4, Fig. 5c, Additional file 1). Individual diving selection coefficients varied among individuals and territories. Individuals exploiting larger home ranges (i.e. 95% AKDE) had a weaker selection for

Fig. 5
The predicted relationship ± 95% confidence interval between diving selection probability and a distance to riverbank, b sediment type and c number of focal vegetation species present among nine individuals in a Eurasian beaver population in southeastern Norway diving locations closer to the riverbank than individuals with smaller home ranges (Table 5, Fig. 6a). We found no clear context dependent effects on selection coefficients for the various sediment types (Additional file 1). Social rank had an effect on diving selection for number of focal species present which were stronger among subordinate individuals than dominant individuals (Table 5, Fig. 6b). We found no conditional effects of other variables on the selection coefficients (Additional file 1). We found that diving locations further from the riverbank had a deeper water depth and higher vegetation cover of focal species compared to diving locations closer to the riverbank ( Table 6, Fig. 7a, b). Diving locations located further from the beaver lodge had a higher vegetation cover than diving locations closer to the beaver lodge ( Table 6, Fig. 7c). Diving distance to the beaver lodge also increased with increasing tail fat index and decreased during the night and during the spring/early summer ( Table 6, Fig. 7d-f ). Focusing on dives within 150 m of the beaver lodge, we found a difference among sediment types with dives located on clay and mud sediment being closer to the lodge than dives located on sand and rock ( Table 6, Fig. 8).

Discussion
We provide new knowledge on aquatic habitat use in a semi-aquatic mammal, the Eurasian beaver, by examining finely resolved information on beaver movement and diving in relation to fine-scaled qualitative assessments of aquatic habitat characteristics within beaver territories. We observed clear spatiotemporal variations in hourly diving probability but found no differences among individuals or territories. Beavers selected for both spatial,  . 6 The predicted relationship ± 95% confidence interval between a selection coefficients for distance to riverbank and home range size (AKDE, autocorrelated kernel density estimate) and b selection coefficients for number of focal vegetation species present and social rank among nine individuals in a Eurasian beaver population in southeastern Norway physical and vegetation characteristics in their diving locations, highlighting the degree of choice they exercise for foraging behaviour although selection strength varied between individuals. Furthermore, spatial variations among dives indicate the energetic variability in aquatic habitat use. Often studies in freshwater-inhabiting semiaquatic mammals focus on the use of terrestrial habitat components, but we show how components of aquatic habitats similarly may be an important resource which can potentially have considerable fitness consequences for a semi-aquatic mammal like the beaver.

Diving patterns
The majority of our identified diving events were short, which matches previous findings for beavers [99,101] and other semi-aquatic mammals [127][128][129]. Beavers have previously been reported to spend less than 3% of their nightly activity budget on diving activities which corresponds well with their role as generalist herbivores that do not rely solely on aquatic foraging [99]. Similar diving patterns have been found among semi-aquatic generalist carnivores and may relate to semi-aquatic animals being less specialized for the aquatic environment [101,[129][130][131]. Beavers may preferably dive to bring large quantities of vegetation to the surface rather than consuming small amounts of vegetation underwater, saving energy by not diving repeatedly to exploit resources at depth, thereby also minimizing heat loss [132,133]. In fact, an extensive study on diving in beavers by Graf et al. [99] found that animals had high dynamic body acceleration (DBA) (a good proxy for movement-based energy expenditure [134]) for the descent, indicating work done against appreciable buoyancy, as has been noted for many birds with their air-filled plumage [133]. Curiously though, in stark contrast to birds, which use this buoyancy for passive ascents [135], beavers also had high DBA values during their return to the surface, which was suggested to be due to animals having to transport vegetation from the bottom to the surface for consumption [99]. Similar behaviour is commonly observed in water birds such as Eurasian coots (Fulica atra) when foraging on aquatic vegetation [136,137], but also in semi-aquatic carnivores such as American mink (Neovison vison) that occasionally consume aquatic food items at the surface or on the riverbank [128,138]. Critically, this behaviour is most advantageous when large amounts of food can be brought to the surface during one dive, which can then be consumed at leisure without the need for multiple, energetically onerous dives in repetitive feeding bouts such as those performed by carnivorous species [138,139]. We found a decreased hourly diving probability in the early morning which is contrary to the findings of Graf et al. [99], although they had appreciable variation. Other studies have found a peak in general activity (measured via overall body dynamic acceleration) in the middle of the beavers' principal activity period, suggesting increased activity in the middle of the night [48]. Diving generally has a high overall body dynamic acceleration compared to other behavioural activities and therefore a higher movement-based energy cost [47,134]. Diving patterns in beavers may be implicitly represented by this general activity pattern, peaking in the middle of the principal activity period, as aquatic habitat use may be costly for a semi-aquatic mammal [101,140]. Similar diving patterns peaking in the middle of the activity period have been found in other semi-aquatic mammals, such as American mink, that perform temporal niche shifts to avoid interspecific aggression from competitors [128,141].
Diving patterns may be structured according to the activity peaks of potential terrestrial predators [142]. When beavers bring aquatic resources onto the riverbank to be handled, their risk of predation increases [93,95]. Beavers are at a particular risk when on land because of their poor eyesight under low light conditions (i.e. they lack tapetum lucidum) [143] and dependence on Table 6 Effect size (β), standard error (SE), lower (LCI) and upper (UCI) 95% confidence interval of explanatory variables for the analyses of distance to riverbank and distance to beaver lodge among dives of nine individuals in a Eurasian beaver population in southeastern Norway Effects were modelled using a GLMMs with gaussian distribution. Beaver ID was included as random effect. Informative parameters are given in bold olfaction to detect potential risks [144]. Consequently, they may not detect potential predators, such as wolves (Canis lupus) and Eurasian lynx (Lynx lynx) that have advanced night vision [145,146], before being detected themselves. However, wolves have been observed to predate more at dawn, dusk, and during moonlit nights [142], which could increase the benefits of being most active (including diving) in the middle of the night when it is darkest. In North America, beavers make up a large proportion of wolf diet, especially during summer when wolves have been observed to ambush beavers at frequently used locations [95]. Natural predators are absent in our study area, but behavioural activities are known to be influenced by historical threats [84,147]. In addition, Fig. 7 The predicted relationship ± 95% confidence interval between distance to riverbank and beaver lodge and a water depth, b vegetation cover of focal species, c vegetation cover, d hour of the night, e Julian day, and f tail fat index among dives of nine individuals in a Eurasian beaver population in southeastern Norway. Points represent actual distances human activities both on land and on water, which are naturally reduced at night, may also influence the activity levels of beavers [147]. Consequently, diving for food resources at dusk and dawn may be perceived to be too risky for beavers. Diving is presumably not just shaped by risks, but also allows the animals to have access to important aquatic resources [4,81]. Hourly diving probability decreased with increasing distance from the territory borders which may indicate a possible depletion of aquatic resources near beaver lodges that are often located in central parts of the territory [cf. 148,149]. Although this may be expected to be related to territory size [33], we did not find this, possibly because our sample size was too small. Diving near the borders may also relate to territorial defence activities, which may help reduce aggressive territorial encounters when territorial intruders swim away unseen [89,150,151]. However, we were unable to quantify territorial behaviour.

Diving selection
We found a high selection for diving closer to the riverbank. Similar short and shallow dives have been found in other studies of semi-aquatic mammals which probably reflect the energetic constraints of bringing food resources to the riverbank for consumption and a preference for travelling along the riverbank [86,99,127,128]. Association with a riverbank can, in some senses, be treated as a central place [152], which makes clear the energetic costs of leaving the central place to acquire food. In that sense, the decreasing tendency to dive with increasing distance from the riverbank ties in with optimal foraging theory, which has animals maximizing reward and minimizing transit time and energy [153,154]. Distance to the riverbank may also reflect decreasing aquatic foraging options in terms of decreased macrophyte growth which is highly dependent on water depth and light penetration [114]. But short shallow dives may also be preferred as they are energetically cheaper for semi-aquatic animals like beavers that have high buoyancy [140]. In general, diving locations in semi-aquatic carnivores have been shown to follow distributions of aquatic food resources [127,155]. We noted that selection for diving locations near the riverbank was statistically stronger among individuals that exploited a larger home range: individuals that exploit a larger home range may have a greater need to stay closer to the riverbank due to higher territory patrolling efforts whereas individuals exploiting smaller home ranges may be able to forage further away from the riverbank [86]. Individuals using larger areas or inhabiting larger territories may also have reduced resource depletion so that, conversely, beavers restricted to smaller areas may be forced to exploit foraging areas further away from the riverbank [33,156]. Larger home ranges may also have a greater area of shallow water, which, depending on the time spent diving, can be energetically easier to exploit for a semi-aquatic animal [140]. Conversely, individuals with smaller home ranges may have to exploit all habitats to fulfil their energetic requirements [33]. We saw that water depth at the diving locations generally increased with increasing distance from the riverbank, which presumably represent higher energetic costs of diving away from the riverbank, but dives at longer distances from the riverbank also had higher amounts of quillwort, shoreweed, and stonewort present, possibly due to the Ashmole's halo effect [148,149] operating on animals preferentially associated with the shoreline. This indicates the interplay of depth and distance energetics in a semi-aquatic animal, which have to be balanced with the calorific value of the food plants and their location [132,133,140].
Beavers had a higher selection for diving at locations with clay sediment, which may be an important building material for lodges and dams [96], although dams are not present in our study site, and beavers may additionally make use of burrows dug into the riverbank whereby the clay sediment aids in enhancing the structural integrity of those burrows [157,158]. Lodges and dams are mostly repaired in the autumn [159,160], but may be repaired after flooding events too [161]. Mud is also widely used for beaver constructions [97], although mud substrates, Fig. 8 The predicted relationship ± 95% confidence interval between distance to beaver lodge and sediment type among dives within 150 m of the beaver lodge of nine individuals in a Eurasian beaver population in southeastern Norway. Points represent actual distances despite being highly abundant in the study area, seemed not to be selected. Beavers continuously need to apply fresh mud and other fine debris to seal their constructions as it is continuously washed away [96,97]. We found that dives on mud and clay sediment within 150 m of the beaver lodge were generally located closer to the beaver lodges than dives located on sand and rock sediment which, we believe, indicates the importance of fine sediment for building constructions. However, we cannot rule out that mud and clay sediment may contain advantageous foraging options, although our ordination did not show strong correlation between mud and clay sediment with aquatic species richness or vegetation cover.
We found a higher selection for diving locations with presence of species of quillwort, shoreweed, and stonewort, which may represent important food resources for the beaver. Other studies have similarly found a high preference for quillwort in the early summer [73], and algae like stonewort may be selected for its protein content and other nutrients [65]. Although beavers mainly forage on woody vegetation, aquatic vegetation seems to be seasonally important, with studies additionally reporting beavers foraging on, among others, water lilies (Nymphaea spp.), pondweed (Potamogeton spp.), water horsetail (Equisetum spp.), waterweed (Elodea spp.), and water lobelia (Lobelia dortmanna) [64-66, 68, 69, 71, 98]. Incorporating aquatic vegetation into a varied diet could be a strategy to minimize risk of nutrient deficiency [162], which may be seasonally beneficial to some individuals (e.g. females during lactation when energetic requirements increase) [163,164]. However, we did not find any selection differences between males and females which may be because of our limited number of individuals. Other studies on semi-aquatic mammals show how males and females in two species of shrew (Neomys fodiens and Sorex coronatus) use separate foraging habitats during the breeding season [165]. Males and females in our population have been shown to differ seasonally in aquatic foraging with peaks in the spring and late summer for females, whereas males only foraged on aquatic vegetation in the spring [166]. We found that selection for the focal plant species were statistically stronger among subordinate individuals, which may be linked to their higher energetic requirements resulting from their activities related to attempts to become dominant in a territory [57] (e.g. performing more extra territorial movements [15]). A higher use of aquatic vegetation may also be a risk-avoiding strategy, minimizing predation risk on land [92,95]. In other studies, adult beavers have been found to forage less on aquatic vegetation than subadult individuals (i.e. 2-year-old) [65] which may be less risk-willing as they potentially face higher fitness costs in terms of future reproductive success [167] and therefore adjust their foraging strategy accordingly [4,168]. Similarly, we found that beavers with lower tail fat index tended to dive closer to their lodge, which may be a consequence of territorial constraints as individuals with higher body condition may be better able to cope with the increased cost of patrolling and protecting territory borders [18,86]. Diving locations located further from the beaver lodge also had higher amounts of vegetation cover which indicate the energetic trade-offs between costs and benefits that a central place foraging individual experiences [4,169]. Diving distance to the lodge decreased through the night which may indicate a functional change in the purpose of the dive. As dives further from the lodge occurred at locations with more aquatic vegetation, these dives may be intended for foraging whereas dives later in the night may be used for building activities. We also saw that diving distance to the lodge decreased from early spring to summer, which may relate to increased parenting activities that require the beaver to stay closer to the lodge when the kits are born in mid-May [58]. Diving for food resources may be perceived as less risk-taking than going on land [87] and could be preferred by parenting individuals to ensure the growth and survival of their offspring [170].

Methodological limitations
Using dead-reckoning together with an accelerationbased behavioural classification model to identify the temporal and spatial distribution of clear diving events, we identified considerably fewer diving events than a previous study in beavers that identified typically 40 dives per night [99]. This may be related to seasonal variations as we only tracked beavers from April to June, whereas Graf et al. [99] also included observations from the autumn (September to October) where diving conditions may be more favourable, because of increased water temperature and life history patterns [15,57,99]. But dives may also be masked by our method. Despite a high classification accuracy of diving events from the model by Graf et al. [47], some events may not be identified because some behavioural activities, including behaviours not described by the model, can have similar patterns and mask each other [47]. Therefore, it is important to clearly define each behavioural activity in an acceleration-based classification model, but also to include enough variations of each activity from several individuals to improve the precision of the model. Different diving styles may be misclassified as other behavioural categories if the behavioural classification model is not trained on several variations of each behavioural category but only includes typical acceleration patterns for each behavioural category. For example, the ecological difference between a 'dive' and an 'almost dive' may be minimal when a beaver can access aquatic resources by just sticking its head underwater, but they may fall within two different behavioural categories because of variable acceleration patterns. The acceleration can also be affected by external environmental forces such as wave actions [134], which can be corrected by incorporating magnetism in the classification [50]. We also only gathered acceleration in 10 s bursts for the classification model, which may mask some of the diving events of shorter duration. A more detailed inspection of the fine-scale acceleration and body postures may improve the classification and provide more information on the actual behavioural activity [171,172]. Placing the behavioural activities into a larger context (i.e. what the animal did before and after an activity) would furthermore help understand the ecological significance of each activity. In addition, the ability of dead-reckoning procedures is spatially limited by the precision of the GPS positions that are used to ground-truth the deadreckoned movement tracks [49,51], which means that we will have introduced some spatial error in the locations of the diving events [112,113]. However, this potential error will be consistent along the tracking period making it less likely to bias our results. The high classification accuracy of the model together with the filtering of less likely diving locations (e.g. on land or not in combination with swimming) improve our confidence in our ability to classify relevant diving locations in our beaver population.

Conclusion
By coupling fine-scaled information on individual beavers' movement and diving with comprehensive qualitative assessments of aquatic habitat characteristics within beaver territories, we provided new knowledge on the aquatic habitat use by a freshwater semi-aquatic mammal. We showed how energetic constraints may shape beavers' spatial use of the aquatic environment, and how aquatic habitats may have great importance for both foraging, building materials and safety, even in absence of natural predators. However, future studies should investigate the importance of aquatic habitat use relative to terrestrial habitat use. Several groups of individuals experiencing various ecological conditions may benefit greatly from the use of aquatic resources, consequently affecting their body condition, reproduction, and survival [57,90], which should be investigated further in future studies including more individuals and populations.
Additional file 1: S1. The model selection results for the best candidate models investigating the number of dives per night among nine beavers in a Eurasian beaver population in south-eastern Norway. S2. The model selection results for the best candidate models investigating the diving probability among nine beavers in a Eurasian beaver population in south-eastern Norway. S3. List of aquatic species abundance in Saua river in south-eastern Norway. S4. The model selection results for the best candidate models investigating the diving habitat selection among nine beavers in a Eurasian beaver population in south-eastern Norway. S5. The model selection results for the best candidate models investigating the individual context-dependent diving selection for distance to riverbank among nine beavers in a Eurasian beaver population in south-eastern Norway. S6. The model selection results for the best candidate models investigating the individual context-dependent diving selection for clay sediment among nine beavers in a Eurasian beaver population in southeastern Norway. S7. The model selection results for the best candidate models investigating the individual context-dependent diving selection for mud sediment among nine beavers in a Eurasian beaver population in south-eastern Norway. S8. The model selection results for the best candidate models investigating the individual context-dependent diving selection for sand sediment among nine beavers in a Eurasian beaver population in south-eastern Norway. S9. The model selection results for the best candidate models investigating the individual context-dependent diving selection for rock sediment among nine beavers in a Eurasian beaver population in south-eastern Norway. S10. The model selection results for the best candidate models investigating the individual contextdependent diving selection for species richness of focal species among nine beavers in a Eurasian beaver population in south-eastern Norway. S11. The model selection results for the best candidate models investigating the diving distance to riverbank among nine beavers in a Eurasian beaver population in south-eastern Norway. S12. The model selection results for the best candidate models investigating the diving distance to beaver lodge among nine beavers in a Eurasian beaver population in south-eastern Norway. S13. The model selection results for the best candidate models investigating the diving distance to beaver lodge among dives within 150 m of the lodge of nine beavers in a Eurasian beaver population in south-eastern Norway.
with the relevant guidelines and regulations [174]. No short-term effects have been observed on the movement after tagging [48]. Body mass of dominant individuals have been observed to decrease with number of capture and handling events, but no statistically clear effects have been observed on survival or other body condition indices [57]. Number of capture and handling events was also observed to affect reproduction, but the population seemed habituated to repeated capture and handling in the long-term [57]. To minimize potential risks, the NBP prioritized capture and handling individuals that were necessary for the monitoring and experiments with clear objectives.

Consent for publication
Not applicable.