HOME RANGE, MOBILITY AND HIBERNATION OF BROWN BEARS ( URSUS ARCTOS , URSIDAE) IN AREAS WITH SUPPLEMENTARY FEEDING

Supplementary feeding, although a common practice, is seldom studied in terms of its effect on non-targeted species, such as the brown bear ( Ursus arctos ). We conducted a GPS-GSM telemetry study on nine individuals (out of about 100 supposedly inhabiting that area) with the aim to inspect how supplementary feeding stations affect home range size, mobility and hibernation. We formulated three hypothesis: 1) there is a correlation between the home range size and the density of feeding stations; 2) t he influence of the artificial feeding stations is not changing during hyperphagia and outside the hyperphagia; 3) hibernation is affected by the density of feeding stations, regardless of the areas with trees in masting age present in the home range and forest age heterogeneity. Our analyses showed that the total home range averaged at 148.9 km 2 (range: 24.6–605.1 km 2 ) with 190.1 km 2 for males and 76.9 km 2 for females. Five out of nine brown bears never visited feeding stations and the density of feeding stations did not explain the home range size variability, thus not receiving evidence to support our first hypothesis. No evidence was found to reject the second hypothesis, while the third hypothesis was rejected for now with the proviso that our sample is not big enough for robust conclusions. The comparison of the home range size, as well as the average displacement distance during and outside the hyperphagia season showed no significant difference, contrary to the outcome of other published studies. Despite the variability of the masting areas proportions and age class heterogeneity, the lack of significant difference between the core area and the total home range for these two factors implied that these resources were not concentrated in the core area of the species, but were rather uniformly distributed throughout the whole territory. Neither density of feeding stations nor masting area proportions were found to be solely responsible for the formation of the hyperphagia home range. Additionally, the similarity in the sizes of the home range size during and outside the hyperphagia season suggested a balance between the variables forming the home range during these periods. Most of the brown bears in our study hibernate in their core area, predominantly in its very centre (0.1 percentile), as earlier hibernation than reported for this region from other studies was recorded. Further studies with a larger sample size on the role of the forest age heterogeneity and mast production on the formation of the home range are needed.


Introduction
The brown bear (Ursus arctos Linnaeus, 1758) is a species of the highest importance for ecosystems all over Europe. Due to its biology of a carnivore, it faces a lot of problems sharing the European landscapes with humans (Chapron et al., 2014). Human activities like forestry, game hunting, livestock farming and even tourism are having an impact on the brown bear distribution and behaviour. One of the human activities with a significant impact on the bears is the supplementary feeding (LCIE, 2018). Worldwide, this practice is widely used for herbivores , even in national parks (van Beeck Calkoen et al., 2020), where management interventions should be kept to a minimum. In Bulgaria, the intentional artificial feeding, targeting bears was banned with the National Action Plan (MOEW, 2008). This step proved not to be effective as the bears continued using the numerous ungulate feeding stations, in both state and private hunting grounds.
Three of the most studied topics in the brown bear-related spatial research are home range size, mobility and habitat use, with numerous publications worldwide (e.g., Seaman & Powell, 1990;McLoughlin et al., 2000;Dahle & Swenson, 2003;Ćirović et al., 2015;Gavrilov et al., 2015;Seryodkin et al., 2017). Yet, according to our knowledge, none of them were trying to relate the home range size to the effect of artificial feeding. Such a study was done on the red deer (Jerina, 2012).
Few publications deal with artificial feeding effect on bear movement (Selva et al., 2017), winter activity (Bojarska et al., 2019), diet (Kavčič et al., 2015), habituation/boldness to humans (Kojola & Heikkinen, 2012;Steyaert et al., 2016) and livestock depredation (Kavčič et al., 2013). The supplementary feeding is known to affect the territorial and social behaviour, mobility and activity of bears Selva et al., 2017), without proof that it is solving human dimension problems such as livestock depredation (Kavčič et al., 2013). Moreover, the provision of additional food might have some unforeseen negative effects on non-target species (Selva et al., 2014;Fležar et al., 2019).
The relationship between the home range size and the mast production is also poorly studied. So far, it was found that the mast production affected the home range use (Kozakai et al., 2011), movement behaviour (Koike et al., 2012) and daily activity time budgets (Kozakai et al., 2013) of the Japanese black bear (Ursus thibetanus Cuvier, 1823). This effect was also noted for specific changes in the Scandinavian brown bear behaviour, such as seasonal activity patterns, human settlement visits and habitat use (Hertel et al., 2019). For the black bear, some attempts were made to predict whether hard mast production would be related to bear visitation of baiting stations (Clark et al., 2005).
For the analysis in this study, we used the effects of two main variables -supplementary feeding stations presence and areas with acorn producing tree species to set up the following hypotheses: 1) There is a correlation between the home range size and density of the feeding stations, regardless of the habitat quality. Thus, the density of feeding stations plays an important role in shaping the home range.
2) The influence of the artificial feeding stations is not changing during hyperphagia and outside the hyperphagia, regardless of the habitat quality in sense of mast area distribution and size and forest age heterogeneity.
3) Hibernation, defined as start and duration, is affected by the density of feeding stations, regardless of distribution of areas with trees in masting age. This hypothesis was based on the biological assumption that after the natural food sources depletion, the bears usually go into hibernation. The presence of more artificial feeding stations will lead to a change in the natural hibernation cycle (earlier or later hibernation from the «typical» recorded for the country, or even no hibernation at all).

Study area
The study was conducted in the central part of the Stara Planina Mnt., on the grounds of the Central Balkan National Park and several State Forestry and State Hunting Enterprises (Fig. 1). The State Hunting Enterprises are units not only for the management of forests and timber production like the Forestry Enterprises, but also include additional activities for intensive game management that include intensive supplementary feeding.
This study area was chosen as it is a part of the Central Balkan segment of the East Balkan population (Linnell et al., 2008;Kaczensky et al., 2013) with the highest registered density of the species in the country (1-2 bears/10 km 2 ) (MOEW, 2008). The estimated number of the population in the Central Balkan segment was reported to be around 100 bears, so we consider our sample representative for this relatively isolated part of the population.
The Central Balkan National Park is covering a larger part of our study area (716.7 km 2 ), with the highest altitude of 2376 m a.s.l (peak Botev). It was created to protect the largest non-urbanised and compact common beech (Fagus sylvatica L.) forest complex in the country, with an average age of the trees ~ 130 years. The rest of the forest complexes include the xero-thermophilic oak (Quercus sp.) woods, hornbeam (Carpinus sp.) and scattered patches of coniferous formations. More than 60% of the study area is covered predominantly with Fagus sylvatica. The rest are meadows, pastures, and agricultural land. The area is scarcely populated with humans, (density: 58.7 people/km 2 , range: 14.0-117.5).
In the study area, 106 supplementary feeding stations for ungulates (Cervus elaphus Linnaeus, 1758, Sus scrofa Linnaeus, 1758 and Capreolus capreolus Linnaeus, 1758), were registered (Fig.  1). These stations are located mainly in the hunting grounds of the State Hunting and Forestry Enterprises (Electronic Supplement 1). The feeding stations placement is based on the Regulation of Hunting and Game Conservation act (MAFF, 2001), according to which for Cervus elaphus and Capreolus capreolus there should be one feeding station per 10 individuals and additionally one feeding station per 3 km 2 hunting ground for Sus Fig. 1. Map of the study area and total home ranges of the GPS-collared individuals as calculated by the 90 th (0.9) percentile.
scrofa. Only 28 of the 106 feeding stations are situ-106 feeding stations are situ-feeding stations are situated in the Central Balkan National Park providing food mostly in the winter or on irregular intervals during the other seasons. The food provided on all the feeding stations is mainly corn (Zea mays L.). Brown bears are attracted by this food of high nutritional value, as Gunchev (1989) has confirmed a high percentage of corn intake by the species in the study area. Additionally, in the study area, brown bears were frequently registered by camera traps at supplemental feeding stations. On the hunting grounds, the food is provided more intensively, once per every 1-2 days in winter and spring and at least once per week in the other seasons (Electronic Supplement 2). During our study, 73.6% of these feeding stations were providing food regularly and 26.4% seasonally. Because of the variability of feeding stations management (differences in the frequency of feeding was found even within the same stations) for the purpose of this study, we assumed that all the feeding stations were attractants for the brown bears per se. Most of the stations had automatic feeders (not attractive for brown bears, because of the small amount of food provided periodically). However, almost all of these stations were supplemented with an additional amount of food, provided manually in various quantities, thus increasing their attractiveness to the animals. All of the feeding stations were visited for food or simply checked by brown bears, regardless of the intensity and type of feeding. None of the stations were providing enough food to sate one individual brown bear, especially with the competition of Sus scrofa and Capreolus capreolus.

Capture and handling of brown bears
As the brown bear is a strictly protected species in Bulgaria, it is under the regulation of Bulgarian Biodiversity Law, Annex 3. As such, the capture and handling of the brown bears were conducted with a permit from the Ministry of Agriculture, Food and Forestry, in which the handling procedure was prescribed and strictly followed. The capture sites and capture time were chosen based on the brown bears' presence, according to information from the Central Balkan National Park and SFE/SHE staff and camera trap data, which recorded the frequency of brown bear visits to the places later assigned as capture sites.
For the capture, we used standard Aldrich snares (Aldrich Snare Co., Washington). The trap placement followed the general procedures for using Aldrich snares (Jonkel, 1993). We used a brown bear capture alert system during the night, to minimise the brown bear stay in the trap to less than two hours.
During the study period (2007-2019) we captured 15 brown bears, nine of which were used in the current analysis ( Table 1). Two of the collared animals (CB5 and CB9) were second-year cubs, led by their mothers (present at the capture site) until finishing the study. This is why we consider both in the group of females since they represented their mother's territory. One of the brown bears was captured twice (CB2, Table 1). After the second capture, it was released in about 12 km straight line distance from the capture site. Another brown bear CB10 (also released from a poacher's snare) was released about 5 km straight line distance from the capture site within the boundary of the Central Balkan National Park. Most brown bears were captured in autumn when brown bears are in the hyperphagia period (n = 7). The capture site locations were recorded with a GPS device.
Upon capture, the animals were tranquilised remotely, weighed and measured. The status of the teeth was recorded for approximate assessment of the age.
Two different types of GPS-collars were used in this study: Vectronic (Vectronic GPS Pro Lite, Vectronic Aerospace GmbH, Germany) and Followit (Tellus GPS Medium Plus, Followit, Lindesberg AB, Sweden), equipped with a GSM module for data transfer.

Data collection and spatial analyses
In the current study, we used an individualbased approach to GPS telemetry data from the captured individuals to analyse the home range size, mobility and hibernation according to predefined variables.

GPS-collars data
The collars were initially set to record locations at a varying time interval, ranging from 2 h to 12 h, depending on the season: every 2 h during spring, summer and autumn and every 12 h during winter (power save mode). During the breeding season, some of the collars were set to take more frequent fixes (10 min. per fix for one week). This called upon the need for sub-sampling to equal time frame, due to two reasons: 1). dealing with autocorrelation (De Solla et al., 1999), which is a pre-requisite for the kernel density estimation model, chosen for home range estimation in our study; 2) analysing and comparing the step length (Euclidean distance between subsequent locations, forming the animal displacement in space, hereafter called «displacement») across individuals throughout the whole study period. Twelve hours between each location fix (at 00:00 and 13:00 ) were chosen (Table 1) to compromise between the shortest step length possible, the coincidence of available GPS fixes from the collars, the continuousness of the data flow and the possible effect of correlated locations. Despite that in most of the published literature the daily (24 h) displacement distances were considered, many recent studies showed that the brown bear activity (including the alternation of the species rest and activity) did not support the idea of discrete 24-h samples (including for displacement) as biologically meaningful. Many behavioural studies recommend using data at a high as possible resolution. That is why in our study we used the highest possible resolution (12 h.). The same was used for building the utilisation (kernel) model of the home range. Note: * -re-capture (poacher's snare); ** -the cub considered in the female's group, due to representing its mother's area; *** -released at a different from the capture location.

Spatial analyses
The supplementary feeding station locations were mapped with a GPS device at the centroid of the feeding ground (with an average area of approximately 30 × 30 m). Since the acorns are the main part of brown bear's diet during hyperphagia in Bulgaria (Gunchev, 1989), the spatial analysis was focused also on hard mast producing species. As data on the mast production rates is deficient, the current study relied on modelling the forest dynamic using the age and origin of the trees to provide an open concept of the possible mast production. This method allowed us to have a general understanding and the closest to accurate information on possible mast production to use further in the analysis. The spatial distribution of areas with mast (nuts and acorns) capacity was extracted from the Database of Forest Inventory plans, provided by the Executive Forestry Agency/Ministry of Agriculture, Food and Forests. Three main variables (tree species, their age and their origin (from seeds or shoots)), were used to create an algorithm for trees in masting or potentially masting age class based on Delkov (1984). The following tree species, represented within the GPS-collared brown bears' home ranges, were considered: Fagus sylvatica and various oak (Quercus robur L., Q. frainetto Ten., Q. petraea Liebl., and Q. cerris L.) species. The extracted data were classified into two categories: masting/potentially masting areas and non-masting areas, as in the analyses we used only the proportions of masting/ potentially masting areas in % (also referred hereafter as masting area). The forest composition is important not only as a food source, but also for providing visibility cover and day resting places for the brown bears. Young forests (age classes from I st to III rd ) are preferred for those activities due to the higher density of trees, typical for forests in that age. Old forests might be used for avoidance of tourists, since territories with such high age class (VI th ) are usually found mainly in Protected Areas. To account for those brown bear activities, we extracted and analysed the forest age heterogeneity, expressed as proportions of forest age classes: 0-20 (I st ), 20-40 (II nd ), 40-60 (III rd ), 60-80 (IV th ), 80-100 (V th ) and over 100 (VI th ) (Hundeshagen, 1826). The non-forested areas are marked as Class 0. During the analyses, it has been taken into consideration, whether the trees are part of a coppice (low) or a high forest, thus adjusting the classification accordingly.
All spatial analyses were conducted with ArcGis Desktop 10.2.2 -ArcMap (ESRI). The resampling of the locations to 12-h intervals, the calculation of the home range and the individual displacement (step length) for each individual were done with ArcMET 10.2.2.v3 extension for ArcGIS Desktop (Wall, 2014).
The annual home ranges for each brown bear were calculated through fixed kernel density distributions (Worton, 1989) and «href» smoothing parameter. Kernel isopleths (i.e. probability contours) were calculated for each 10% increase in kernel density (hereafter called «percentiles» in decimal values). Percentiles were used for identification of denning site related to the home range. The 50% and 90% kernel isopleth range were further defined as a core area (defined as a portion of the total home range which had the most intensive (Samuel et al., 1985) or concentrated (Seaman & Powell, 1990) use and expressed as 50% kernel isopleth range, e.g. percentile 0.5) and a total home range area (90% kernel isopleth range, e.g. percentile 0.9). The locations outside the 90% kernel isopleth range were not considered to account for the numerous outliers. For the core and total areas, a comparison of the variables was made between home range size, density (and where appropriate, number) of feeding stations, average 12-h displacement distances (hereafter average displacement distances), proportion of the masting areas (in %) and forest age heterogeneity.
To account for the pre-denning conditions (e.g. food resources available), we calculated the home range size during hyperphagia (using fixed kernel density distributions on locations between 1 st September and 30 th November). The same was done for the period outside this season at 0.9 percentile to be able to test for differences in the home range size and variables during these two periods. The same hyperphagia period is used in the research of Pop et al. (2018) in Romania and by Štofik et al. (2013) in Slovakia, where the climate and respectively vegetation and tree species are similar. For both of these home ranges, we extracted the density of feeding stations, average displacement, and proportions of masting area and tested for differences between the hyperphagia and outside hyperphagia.
We extracted the start-end date of the hibernation and duration in days from the first, inter-mediate and, respectively, last location clustered within less than 100 m (expected GPS error under denning conditions) in the period when denning was expected (end of November -end of March). The denning sites were overlaid with the percentiles layers, identifying their locations within the home range percentiles. Additionally, the proportion of the masting areas and density of feeding stations during hyperphagia were compared to the hibernation data.

Statistical analyses
All statistical analyses and plots were carried out with Statgraphics Centurion 18. 1.11 (Statgraphics Technologies,Inc.). Due to deviation from normality, the non-parametric Mann-Whitney (Wilcoxon), W-test was used to test for differences in medians of the masting areas between the core area and total home range; between the home range sizes in hyperphagia and outside this season; between the displacement distance in the core area and total home range. The same was used for comparison of forest age heterogeneity proportions between core area and total home range. Spearman and Pearson correlations were used to test for relationships between several variables such as home range sizes, density of feeding stations, average displacement, and masting area size. The differences and correlations were tested with two-tailed probabilities at 0.05 significance, within the 95.0% confidence level.

Home range size and density of feeding stations
The shape and size of the individual home ranges (both in sense of core and total area) varied depending on the sex and age of the animals (Fig. 1, Table 2). The average size of the core areas was 37.0 km 2 (with the range being 6.5-124.2 km 2 ), differed between sexes as follows: 47.23 km 2 (median: 23.2 km 2 ) for males and 18.9 km 2 (median: 16.6 km 2 ) for females. The total home range averaged at 148.9 km 2 (range: 24.6-605.1 km 2 ) with 190.1 km 2 for males and 76.9 km 2 for females. The CB2 individual showed a larger home range after the recapture and translocation, despite it being tracked fewer days (330 days before recapture vs. 52 days after recapture). The movement of the subadult CB10 showed a clear dispersal pattern. This individual was tracked for 91 days within its first home range (CB10-HR1 in Table 2). After 11 days of oriented movement, it settled into its new home range (CB10-HR2 in Table 2), where it was tracked for another 244 days. The second home range had a 1.7 times larger core area, forming a 1.9 times larger total home range than the estimated size of the first home range.
All the home ranges of the studied individuals (with the exeption of the second home range of CB10) contained feeding stations. Yet, five out of nine brown bears (CB1, CB4, CB9, CB10 and CB12) never visited a feeding station during the study period. The core areas of six individuals (Table 2) were associated with a low density of feeding stations (one feeding station per 10 km 2 ), in three individuals (CB2-HR1, CB10 -both home ranges and CB12) none were observed. No correlation was found between the core area size and the density of feeding stations (Spearman correlation, P = 0.480). A clearly positive correlation was found between the number of feeding stations and the total home range size (Spearman's correlation ρ = 0.74, df = 10, P = 0.020). However, no correlation was found between the total home range size and the feeding stations density (Spearman correlation, P = 0.328).

Home range, masting areas and forest age heterogeneity
The quality of the habitat in the sense of masting area size showed a wider range in the core area ( Fig. 2A) than in the total home range. No statistically significant difference was found in the proportions of masting areas between the core area and total home range (Mann-Whitney test on medians, W = 62.0, P = 0.948). Around half or less of the core areas and total home ranges of most of the studied individuals was covered by a masting or potentially masting trees (Fig. 2B). For two individuals (CB5 and CB8), the core area has shown a concentration (higher density) of feeding stations in the core area, compared to the total home range.
The forest age class proportions for both core area and total home range are presented in Fig. 3A,B. The proportions of non-forested areas (Class 0) in the core area were more variable with a larger interquartile range compared to the total home range. The same is valid for the classes III (40-60 years) and V (80-100 years). There is no statistically significant difference between the median values of the age class proportions between the core area and total home range (Mann-Whitney test on medians, W = 25.5, P = 0.949).  Fig. 2. Proportions of the masting / non-masting areas (%) within the home range expressed as a decimal index. Designations: A -masting area size for all brown bears in the core and total home range area. (diamond -average value, horizontal line -median value, box -interquartile range); B -comparison of the masting area proportions and density of feeding stations in the core and total home range area per individual. (mast area% -CA -proportions of the core area (in %) covered with masting trees; mast area% -HR -proportions of the core area (in %) covered with masting trees; FSdens -CA (red dots)density of feeding stations in the core area; FSdens -HR (black triangles) -density of feeding stations in total home range; horizontal black line -equal proportion line (50%-50%).

Fig. 3.
Box-plot comparison between forest age classes' proportions distribution for core area and total home range. Designations: A -in core area; B -in total home range size; diamonds -average value, horizontal line -median value, boxinterquartile range, dots -outliers, dots with plus sign -extreme outliers.
We compared the home range size during hyperphagia (between September and November) and outside this season with the density of feeding stations present (Table 3). No statistically significant difference was found between the home range sizes during and outside hyperphagia (Mann-Whitney test on medians, W = 30.0, P = 0.871) and between the density of feeding stations during hyperphagia and outside it (Mann-Whitney test on medians, W = 27.0, P = 0.636). Both average and median values for feeding stations density were higher during hyperphagia than outside of it (average 0.135, median 0.140 per 1 km 2 vs. average 0.104, median 0.090 per 1 km 2 , respectively). No correlation was found between the home range size in hyperphagia and the density of feeding stations present in this range (Spearman correlation, P = 0.102), and between the home range size in hyperphagia and the masting area proportions (Spearman correlation, P = 0.488). A positive correlation was found outside the hyperphagia period between the home range size and the number of feeding stations (Pearson correlation ρ = 0.77, df = 7, P = 0.026). However, no correlation was found between the home range size and the density of feeding stations (Pearson correlation, P = 0.389).

Home range and mobility
All brown bears exhibited relatively similar 12-h displacement distances in the core area and total home range (Fig. 4A), despite the different density of feeding stations present (Fig. 4B). The average 12-h displacement distances varied between 292 m and 1410 m for the core area and between 297 m and 1519 m for the total area. The Mann-Whitney W-test on medians showed no statistically significant difference between the displacement distances in the core area and in total home range (Mann-Whitney test on medians, W = 67.0, P = 0.696). No significant correlation was found between the average displacement distance and the density of feeding stations in the core area (Spearman correlation, P = 0.995) and the average displacement distance and the density of feeding stations in the total home range (Spearman correlation, P = 0.625). 28.5% 0.08 0.08 Note: * -As a hyperphagia, the period between September and November is considered; ** -here are included only those individuals, with data both during and outside hyperphagia, so without CB4. Fig. 4. Mobility of the GPS-collared individuals in the core area and total home range, expressed as average displacement distances and feeding station density. Designations: A -box-plot comparison between average displacement distance for all brown bears in the core and total home range area (diamonds -average value, horizontal line -median value, box -interquartile range, dots = outliers); B -comparison between the average displacement distance and the density of feeding stations in the core and total home range area per individual (CA -average displacement in the core area; HR -average displacement in the total home range; FSdens -CA (red dots) -density of feeding stations in the core area; FSdens -HR (black triangles) -density of feeding stations in total home range).
The average displacement distances were significantly and positively correlated with the masting areas size within the core area (Spearman's correlation, ρ = 0.65, df = 10, P = 0.04). The same was found for the displacement in the total home range (Spearman correlation, P = 1) vs. masting area size (Spearman's correlation, ρ = 0.76, df = 10, P = 0.016). At the same time, the average displacement distances during and outside hyperphagia were similar (1455 m vs. 1475 m) as more variability was expressed in-between the individuals during hyperphagia, than outside the hyperphagia (Fig. 5). Yet, the compari-5). Yet, the compari-). Yet, the comparison of the medians did not show any statistically significant differences (Mann-Whitney test on medians, W = 51.0, P = 0.377). Two of the individuals (CB10 and CB13) were even more mobile during hyperphagia than outside this period.
CB13 expressed the highest mobility of all the brown bears (more than 1.2 km per 12 h). Most of the brown bears kept their 12-h displacement distances below 600 m. The two individuals with separate home ranges (CB2 and CB10) had a different mobility in the first and second home ranges. CB2 showed almost two times longer average displacement distances in the second home range (after translocation), while C10 had even a little lower average displacement distance after settling in his second home range.

Home range and hibernation
The hibernation duration varied between 28 days (CB10) and four months (CB4, 122 days). Almost all individuals have chosen their winter den sites in the core areas of their home ranges (Table  4). Four out of the nine brown bears showed preferences for the percentile 0.1, acting as a «nucleus» of the core area.
No significant correlation was found between the hibernation duration and masting area proportions (Pearson correlation, P = 0.502) and between the hibernation duration and density of feeding stations (Pearson correlation, P = 0.375) during hyperphagia. Yet, the brown bears which started their hibernation in November had slightly higher average masting area proportions (0.49) than those starting their hibernation in December (average masting area proportions: 0.39). The average density of feeding stations is slightly lower though (0.15 per 1 km 2 ) for the earlier hibernators, than for those starting hibernation in December (0.16 per 1 km 2 ). Most of the earlier hibernators have either a higher proportion of masting area during hyperphagia (CB12) or a proportion masting areas above 0.30 (30%) and a higher density of feeding stations (more than one feeding station per 10 km 2 ). The only individual with a late and the shortest hibernation (CB10, 28 days) had the lowest masting area proportion (0.25) in his hyperphagia home range and no feeding stations.

Discussion
First, we should acknowledge that more feeding stations might be present in our study area. Before and during the hunting seasons, the hunters set temporary unannounced feeding stations that were impossible to locate and map. Secondly, we acknowledge that our small sample size and the outcome of this study are relative only to this local area and cannot be representative for the whole country. Yet, these are, according to our knowledge, the first efforts of trying to identify the complex relations between the home range size and mobility on one side and feeding stations presence, masting area proportions and hibernation on the other.
In this study, we worked under the assumption of Powell (2000) that the core area size is not a pre-requisite for the total area size. The core area size for our studied individuals ranged between 20% and 31% of the total home range. The total area showed a high variability. This variability was not only due to the sex differences, but also due to the sub-adults or cubs recently separated from  their mothers (~ 2 years old). Such individuals (like CB1, CB2 and CB10) are more susceptible to segregation from larger and more dominant individuals. We also found that the males' home range size was more variable. Some literature sources stated that this might be due to segregation with an increased local population density (Nagy & Haroldson, 1990), as this effect was more prominent in males than females due to matrilineal assemblages (Dahle et al., 2006). Further on, some publications provide proof that despite the variability and generally larger home ranges in the less-productive areas, the home range size was more related to social interactions than food availability (Nagy & Haroldson, 1990;Dahle & Swenson, 2003;Dahle et al., 2006). Thus, the comparison of the home range size between populations and even between regions would be meaningless, if these sizes were not considered together with the sex and age structure of the local population and its density. These data, regretfully, were seldom available, reported or compiled with the home range size to be able to draw robust conclusions.
Due to supplementary feeding, corn is part of the brown bear's diet (more specifically in autumn) in many countries throughout the species range. It is usually obtained by the brown bears either at hunter's ungulate feeding stations or in agricultural fields (Landers et al., 1979;Rigg & Gorman, 2005). In Bulgaria, the usual approach is to provide corn as an appropriate supplement or substitute for the high energy acorns and nuts. Being concentrated and often abundant in resources, the artificial feeding stations were reported to increase the home range overlap, decrease home range size, increase group size, and reduce territoriality (Robb et al., 2008). This is also reported in other species, such as Cervus elaphus: the increase in density of feeding stations or the frequency of feeding, as well as the increase of species density, is directly leading to a decrease of the home range size (Jerina, 2012). A decrease of home range size with the increase of feeding stations density however, was not observed in our study (no correlation found both for core area and total home range). The lack of correlation between the average displacement distance and the density of feeding stations also supported the idea of low impact of the feeding stations presence on the home range formation. In some of the individuals an increase of the home range was observed possibly due to segregation by other more dominant conspecifics or/and the local management of the feeding stations (low and irregular food provided), complemented with the effect of a higher local density of the ungulates (competition). Usually the food is provided in too small quantities, not enough to sate a brown bear. Additionally, the ungulates often arrive earlier than the brown bears and consume most of the food provided, which force the brown bears to seek other sources (often, it means other feeding stations). Some brown bears were registered at several stations within duration of one night. Yet, this is not valid for all the studied individuals, as some never visited any feeding stations during the study period (provided, that visits were not omitted due to the time gaps between GPS fixes). This implies that the repeated observations of brown bears at feeding stations (recorded also with camera traps, see Electronic Supplement 1) were most probably due to habituated individuals to artificial feeding. Such habituation is often observed Steyaert et al., 2014), as brown bears are highly intelligent and adaptable and tend to utilise various food resources.
Brown bears that do not visit feeding stations might change their behaviour when the resources run low as a compensatory mechanism (Mattson et al., 1992;Oka et al., 2004;Schwartz et al., 2006). For example, a negative relationship between bait-station visitation and oak mast production was observed in a grizzly bear (Ursus arctos horribilis (Ord, 1815)) study, which set up a hypothesis that poor body condition «could increase bait station visitation by bears and vice versa» (Clark et al., 2005). So even if we did not observe such visits to the feeding stations during our study period, this did not mean that they will not start doing this during low resources. For example, CB9 never visited feeding stations, while he was with his mother. But a few months after dispersal (contact with the collar lost before that), this brown bear was found poached with a stomach full of corn, that in this area could be found only on a feeding station.
In the light of the finding that more than half of the studied individuals never visited feeding stations and no correlation was found between the density of feeding stations and the home range size, we looked at the other high energy source, acorns. The role of the mast was confirmed for the Ursus thibetanus which changed their habitat use in autumn (Kozakai et al., 2011). Other studies from North America and Japan had also supported the mast importance in autumn (Mattson et al., 1992;Schwartz et al., 2006) as bears rely on the mast during hyperphagia to increase their body mass in preparation for hibernation (Hertel et al., 2019). We hypothesised that the distribution of the masting areas and their proportions were important for a better understanding of the role of the feeding stations and the formation of the home range itself. Although the mast production is a seasonal event, we worked under the assumption that brown bears, like other animals, would memorise and integrate such sources of high nutritional value in their «cognitive map» (Powell, 2000) to maximise the chances for survival.
Forest management is suspected to affect both brown bear habitats and their behaviour through the created age class heterogeneity. Young forests (age class I) found in early succession areas, after a calamity, disturbance or timber harvesting, have a canopy in the process of closure, thus being an abundant source for berries, i.e. an essential part of the brown bear's diet. During the process of re-establishment of mature overstory (classes II and III), the tree density increases, making them suitable for day resting places. Mature forest (classes IV and V) provides the highest nutritional food, making them vital during the season of hyperphagia, while old forest (age class VI) provides opportunities for denning sites and low disturbance. All of the above requires additional analyses in the search of the unrevealed factors affecting the home range size and distribution (Goodrum et al., 1971;Stransky & Roese, 1984;Mitchell & Powell, 2003). Despite the variability of the masting areas proportions and age class heterogeneity, the lack of significant difference in these two factors between the core area and the total home range implied that these resources were not concentrated in the core area of the species, but were rather uniformly distributed throughout the whole area. This is a logical outcome, as the areas with mast potential are used only seasonally.
The brown bear mobility in hyperphagia can be explained by mast, as well as berry production. Despite the high value of the masting area size, it explained a very low percent of the displacement distance variability in our GPS-collared brown bears. During our study, the Fagus mast production was high and it was expected that it would result in shorter displacement. Instead, the mobility of the brown bears was not directly connected to the size of the masting areas. Even more, the lack of difference between the average displacement during and outside the hyperphagia season suggested that the brown bears were probably foraging on other food sources, such as artificial food.
Additionally, we observed a positive significant correlation between the home range size (core area and total home range) with the masting areas size, although the average displacement distances during and outside hyperphagia were similar. This was contradicting to other studies where during berry seasons a shift towards shorter displacement was noted in brown bears (Ordiz et al., 2012) or a negative relationship between movement distances and mast production was found in Ursus thibetanus (Koike et al., 2012). The specificity of the mast is the variability of its production and a synchronicity over large spatial areas (Silvertown, 1980), which is supposed to affect the overall bear movement. When low mast during hyperphagia occurs, Kozakai et al. (2011) found that U. thibetanus increased their movement activity (Koike et al., 2012), increased their home range size, and even made short-range excursions outside their home range in the search of food (Pelchat & Ruff, 1986;Koike et al., 2012). The latter was observed for individual CB13 in our study, who during the hyperphagia period spent several days in nearby areas, before returning to his established area.
The lack of correlation between home range size in hyperphagia with the density of feeding stations and with the masting area proportions and between average displacement distance and the density of feeding stations implied than none of these two factors (density of feeding stations and the masting area) were solely responsible for the formation of the home range in hyperphagia. Additionally, the similarity in the size of the home range during and outside the hyperphagia season suggested a balance between the variables forming it during these periods. This is a reason to reject the first hypothesis. The lack of statistical support for the difference of the role of the artificial feeding stations during hyperphagia and outside the hyperphagia allows us to confirm the second hypothesis.
The patterns in hibernation are still not entirely understood (Manchi & Swenson, 2005;Krofel et al., 2017). Most of the brown bears in our study hibernate in their core area, predominantly in its centre (0.1 percentile). As far as we know, there is no such study with which we can compare these findings.
After seven years of monitoring, Gunchev (1989) reported that brown bears in this region normally enter into hibernation in December. The level of artificial feeding and the masting potential in Gunchev's (1989) study were similar as today, due to the fact that the management has not changed. Hence, the data about the hibernation start might be considered valid today. Our results though point at an earlier hibernation, starting in November. No correlation was found between hibernation duration and masting area proportions/ density of feeding stations. But these results should be treated cautiously with the proviso that our sample is not big enough for robust conclusions. The insignificant, but yet present, difference during hyperphagia in average masting area proportions and the density of feeding stations between the «earlier (November)» and «normal (December)» imply that earlier hibernators probably rely more on mast production, while the normal hibernators might stay longer where more feeding stations are present. This is in compliance with Pigeon et al. (2016) who suggested that access to food supply can affect denning chronology and continuousness in terms of higher food availability resulting in later den entering. Research in North America and Southern Europe has even shown that weather conditions and food availability can result in bears remaining active during winter (Van Daele et al., 1990;Huber & Roth, 1997;Linnell et al., 2000). So, the outcomes from this part of our study give us the reason that currently we could reject the third hypothesis.

Conclusions
In conclusion, summarising the outcome and with the acknowledgement of the small sample size, we can state that the first and the third hypothesis were rejected, while the second hypothesis was confirmed. The reason for this is mainly because we did not found a solid proof that there was a dominant variable responsible for formation of the annual brown bear's home range and those during and outside hyperphagia in our study area. Further studies on forest age heterogeneity and mast production are needed with a larger sample size.

Supporting Information
Detailed characteristics of feeding stations (Electronic Supplement 1: Feeding stations details -GPS co-ordinates, forestry/national park unit, locality and frequency of feeding which is most commonly applied), and photos of brown bears on feeding stations (Electronic Supplement 2: Camera-trap photos of brown bears on feeding stations in the study area) may be found in the Supporting Information here.