Balkan Chamois ( Rupicapra rupicapra balcanica ) Avoids Roads, Settlements, and Hunting Grounds: An Ecological Overview from Timﬁ Mountain, Greece

: Balkan chamois ( Rupicapra rupicapra balcanica ) is a protected species with an Inadequate-Bad (U2) conservation status in Greece. Our study explores its seasonal range use pattern, demography and habitat selection in a site of the Natura 2000 network, Timﬁ Mountain. To this aim, we examined 1168 observations obtained from six seasonal surveys (2002: four seasons, 2014 and 2017: autumn) and performed an ecological-niche factor analysis (ENFA), using 16 environmental and human-disturbance variables. The species had an annual range of 6491 ha (25% of the study area), followed the typical range-use pattern, and presented the minimum core area during the rutting season (autumn). Timﬁ Mt hosted 469 individuals in 2017 (the largest population in Greece), increasing by 3.55 times since 2002. The species selected higher altitudes during summer and autumn, pinewoods over broad-leaved woods as winter grounds, and it avoided south-facing slopes. Our results supported the anthropogenic risk avoidance hypothesis; the species always selected remote areas away from roads, human settlements, and hunting grounds. In Greece, 40% of its distribution area falls within hunting ban areas (16.5% of the country). A national conservation policy is needed towards maintaining and increasing roadless areas and hunting-ban areas within Balkan chamois range nationwide.


Introduction
Wildlife shares the planet with an exploding human population [1] that transforms the land in such a way that further rapid wildlife habitat loss is expected, leading to even higher species extinction risks [2]. Although humans can use technological means to colonize and adapt to every corner of the Earth, complex evolutionary and ecological processes interact to determine the distribution patterns of wildlife species. Habitat selection is a fundamental spatiotemporal ecological process: animals select habitats that maximize their fitness, in a continuous trade-off between selecting the best foraging and breeding habitats that would increase their survival probability and reproduction success, while minimizing environmental stresses and anthropogenic risks [3]. The challenge, therefore, lies in pinpointing the underlying factors that restrict animal distribution, assessing in particular the role of human-induced factors as compared to natural ones in the habitat selection process, under the scope to minimize negative anthropogenic impacts through effective conservation management actions and policy decisions [4]. This challenge has to be met in a species-by-species approach at the local scale, where conservation action usually applies, rendering the conservation endeavour highly demanding. which the chamois distribution area coincides with hunting ban areas nationwide. The ultimate scope is to provide a science-based conservation and policy guideline as a paradigm fostering conservation action for the Balkan chamois in Greece and in the Balkans.

Study Area
The study area is located in Timfi Mt (39.877347 Ν, 20.784790 Ε) at the north-western part of Pindus mountain range in Greece, covering 25,896 ha ( Figure 1). The climate is intermediate between the Mediterranean and the Continental type: total annual precipitation of 1560 mm, average monthly temperature from 2.5 °C (January) to 24.4 °C (July), with a minimum in February (−8.6 °C) and a maximum in August (36 °C) (Data from Papingo Meteorological Station, 920 m: 2010-2013). It is a mountainous area (50% above 1000 m) within an elevation range between 412 m and 2497 m; the geological bedrock consists of limestone and flysch. The dominant habitat types in the study area are mountainous grasslands, pine forests (endemic black and endemic Mesogean pines), broadleaved forests (beech, mixed fir-beech, oak and mixed thermophilous forests), and screes (Table S1).  Table S1.
The study area coincides with the Natura 2000 site "Oros Timfi" (GR2130009), covering parts of two more overlapping Natura 2000 sites: "Vikos-Aoos National Park" (GR2130001:87%) and "Kentriko Tmima Zagoriou" (GR2130004:8%), managed by the Management Authority of Northern Pindus National Park. It is a sparsely populated area, including six settlements. The main human activities include livestock breeding, hunting and tourism. Transhuman livestock breeding takes place from June until October above 1300 m, but permanent breeding is confined around the settlements. Tourism is an essential component of the local economy including mountaineering, climbing and hiking, which take place mainly in the Vikos-Aoos National Park. Chamois presence was documented in the study area since the Late Upper Pleistocene, playing a crucial role as a food resource to prehistoric human communities [20]. The population in Timfi Mt is the largest population  Table S1.
The study area coincides with the Natura 2000 site "Oros Timfi" (GR2130009), covering parts of two more overlapping Natura 2000 sites: "Vikos-Aoos National Park" (GR2130001:87%) and "Kentriko Tmima Zagoriou" (GR2130004:8%), managed by the Management Authority of Northern Pindus National Park. It is a sparsely populated area, including six settlements. The main human activities include livestock breeding, hunting and tourism. Transhuman livestock breeding takes place from June until October above 1300 m, but permanent breeding is confined around the settlements. Tourism is an essential component of the local economy including mountaineering, climbing and hiking, which take place mainly in the Vikos-Aoos National Park. Chamois presence was documented in the study area since the Late Upper Pleistocene, playing a crucial role as a food resource to prehistoric human communities [20]. The population in Timfi Mt is the largest population in the country, being one of the 12 sub-populations that form the northern population block of the Pindus mountain range [6]. Pindus is the backbone of Greece, including three chamois population blocks that share similar genetic Diversity 2020, 12, 124 4 of 18 characteristics [21]. Chamois natural predators in the area are currently grey wolf, bear, red fox and golden eagle.

Chamois Surveys
We conducted 47 transects (mountaineering trails) of 205 km length in total, in conjunction with a set of vantage points along them (winter 2001 to spring 2002) (Figure 1), in order to cover all chamois habitats across the elevation gradient. Chamois surveys were evenly repeated by conducting 47 field days per season as follows: winter (cold period: 10 December-9 March), spring (parturition: 10 March-9 June), summer (warm period: 10 June-9 September), and autumn (rutting period: 10 September-9 December) [22]. The dataset included direct (animal sightings) and indirect observations (chamois tracks and droppings found in a belt of 1 m on either side of the transects; older droppings that might refer to a previous season were not considered). A group of animals was considered as one direct observation when the distance from a neighboring group was greater than the average inter-individual distance within the group [23].
To estimate chamois population size and trend, we counted animals using the pointage flash method [24] during the rutting period (three successive days at the end of October) for the years: 2002, 2014, 2017. The same transects were walked by seven groups of two researchers, to scan the study area uniformly. We noted the coordinates, date, time, group size and demographic composition for each direct observation, to be able to detect and exclude duplicate observations. We recorded the age and sex of each animal as follows: kids (<1 year), yearlings (between 1-2 years), adult females (>2 years) and adult males (>2 years) [25].

Environmental Variables
We overlaid a grid of 100 m × 100 m in the study area and calculated the average values, or the percent frequency of occurrence of the following 16 environmental variables per grid cell. We used a digital elevation model (DEM) of 20 m resolution, provided by Anavasi Editions Company. The DEM was reclassified to match the resolution of the grid, to calculate seven topographic variables (elevation, slope, curvature and aspect). We used the available habitat map of the Natura 2000 site of Timfi Mt provided by the Management Agency of Northern Pindus National Park (year 2008). We considered six broader habitat types (grasslands, pine and broadleaved forests, rocks, scrubs, agricultural land), and calculated their cover (%) per grid cell (Table S1). We finally calculated three human disturbance variables: the hunting activity (binary), the distance from settlements and the distance from the nearest road per grid cell. Distances were computed accounting for the curvature of the surface. Spatial data for hunting ban areas (wildlife refuges, including the strict nature reserves of Northern Pindus national park), and settlements were openly available [26]. The spatial file of the road network for the broader study area was provided by Anavasi Editions Company in 2001 and it was cross-checked for 2014 and 2017 using aerial photography. Calculations were performed using ArcGIS 10.7 [27].

Annual and Seasonal Ranges
Using the data from the annual survey 2001-2002, we determined chamois seasonal ranges by employing the 95% fixed kernel density estimator (probability of species occurrence in its range greater than 95%) [28] and the extension "Animal Movement" [29] of ArcView. To define the core area as the area that is used more heavily in every season, excluding any random effect, we plotted the area of animal range (%) with the probability of area use (%) [30], as indicated in the upper right corner of Figure S1: the straight line reflects a random use of space, while the curve that sags below the straight line reflects a clumped use of space. The point where the curve's tangent becomes parallel to the line of random use is the threshold to define the core area. After defining the seasonal ranges and respective Diversity 2020, 12, 124 5 of 18 core areas, we merged them to define the annual range and core area, respectively. We calculated seasonal ranges overlap, using ArcGIS 10.7 [27].

Population Demography
We assessed the annual chamois population density for the period 2001-2002 by dividing the estimated population size with the annual range. We considered only the autumn surveys (2002, 2014, and 2017) to compare the demographic characteristics of the population in terms of population structure and trend.

Habitat Selection
We performed an ecological-niche factor analysis (ENFA) [31] to explore the habitat selection by the species versus available environment [32] for the six survey datasets, namely the four seasons in 2002 and in autumn 2014 and 2017. In the case of direct observations, the absolute number of animals counted was assigned to each grid cell, whereas in the case of indirect observations, we assigned a value of one. ENFA assesses habitat selection by decomposing the ecological niche into two components within the ecological space: marginality and specialization. The marginality measures how different are the average conditions in the habitat used compared to those that are available in the environment. The specialization estimates habitat selection on a particular direction of the ecological space, i.e., the niche breadth of the habitat used. In the absence of high correlation among the candidate explanatory variables (Spearman's |r| ≤ 0.75), we performed ENFA considering all the 16 environmental variables, by contrasting the environmental conditions between the occupied grid cells with presence data and the available grid cells across the study area [33].
We first extracted one axis of marginality and then used the broken-stick method [34] to define the number of specialization axes on the basis of their eigenvalues, which express the amount of specialization they account for [31]. We then identified the degree to which candidate explanatory variables contributed to the marginality and specialization in terms of their respective contribution score to the corresponding axes, noting that high scores mean strong habitat selection and narrow niche. Finally, we tested the statistical significance of the above scores, using a Monte-Carlo randomization procedure (1000 iterations). All statistical analyses were conducted using R 3.6.1 [35] and the R-package adehabitatHS (v.0.3.14) [36].

Chamois Distribution and Hunting Ban Areas of Greece
As disturbance from hunting other game species is considered as high-impact threat for the Balkan chamois, we assessed the spatial overlap of the chamois distribution with the hunting ban area in Greece. According to the Greek legislation, hunting is allowed all over the country, except from the following cases: for reasons of public safety, hunting is banned at a distance less than 250 m of cities and villages, and 100 m from single houses, as well as at a distance less than 500 m from the eastern borderline, in archaeological sites and monuments, and in camping sites when operating. To estimate this hunting ban area approximately, we used Corine land-cover data [37], applied a buffer zone of 250 m around the polygons classified under the broad category of artificial surfaces, and calculated the buffer area of 500 m from the eastern borderline using open data. We also calculated the hunting ban areas for wildlife conservation, namely the area of wildlife refuges, core zones of the 10 old-standing national parks of the country, strict nature reserves and nature reserves, by merging the respective calculated areas, using open data from the European Environmental Agency. Using the same source of open data, we further estimated the merged area of the game breeding stations and the controlled hunting areas of Greece, where hunting is banned or controlled for game management reasons. We merged all the above areas to calculate the permanent hunting ban area, excluding overlapping areas. Finally, we screened all ministerial decisions temporarily banning hunting in specific areas for a given period (game protection, burnt areas etc.) and we considered their overall area for 2020 (37 decisions since 2016). It was not possible to estimate the further hunting ban areas, such as fenced private land or Diversity 2020, 12, 124 6 of 18 private agricultural land where hunting harms the harvest crop. We then considered the most recent distribution map of Balkan chamois, provided by the Greek Ministry of Environment and Energy, which is used in the action plan for the species that is currently under development, and we calculated the proportion of the distribution area that falls within the hunting ban area nationwide. Calculations were performed using ArcGIS 10.7, and percentages were calculated considering the area of Greece (132,029 km 2 ) and the distribution area of Balkan chamois according to the action plan (2706.67 km 2 ).

Results
Mountainous grasslands covered a great part of the study area (46%) and extended above the tree-line together with rocky habitats and screes (9%). The forest zone included pinewoods (15%) and broadleaved woods (14%), whilst scrubs (15%) and agricultural areas (<1%) lay in the lower altitudes ( Figure 1, Table S1). The overall road length was 174 km in the study area and has not changed since 2001. Hunting activity was excluded in 48% of the study area, namely in the core area of Vikos gorge and three wildlife refuges ( Figure 2). We collected 1168 observations overall during the six surveys. For the annual survey of 2001-2002, 36% were direct observations, and 47% and 17% were chamois droppings and tracks respectively (Table S2).
Diversity 2020, 12, x FOR PEER REVIEW 6 of 18 provided by the Greek Ministry of Environment and Energy, which is used in the action plan for the species that is currently under development, and we calculated the proportion of the distribution area that falls within the hunting ban area nationwide. Calculations were performed using ArcGIS 10.7, and percentages were calculated considering the area of Greece (132,029 km 2 ) and the distribution area of Balkan chamois according to the action plan (2706.67 km 2 ).

Results
Mountainous grasslands covered a great part of the study area (46%) and extended above the tree-line together with rocky habitats and screes (9%). The forest zone included pinewoods (15%) and broadleaved woods (14%), whilst scrubs (15%) and agricultural areas (<1%) lay in the lower altitudes ( Figure 1, Table S1). The overall road length was 174 km in the study area and has not changed since 2001. Hunting activity was excluded in 48% of the study area, namely in the core area of Vikos gorge and three wildlife refuges ( Figure 2). We collected 1168 observations overall during the six surveys. For the annual survey of 2001-2002, 36% were direct observations, and 47% and 17% were chamois droppings and tracks respectively (Table S2).

Figure 2.
Annual range and core area of Balkan chamois in the study area deriving from merging seasonal ranges and core areas, respectively; human settlements, roads and wildlife refuges (hunting ban areas) are indicated.

Annual and Seasonal Ranges
Seasonal ranges presented a minimum in winter and a maximum in summer, with the maximum aggregation within the core area observed in autumn (Table 1). We found an important overlap (above 49%) among seasonal ranges. Two main seasonal patterns emerged ( Figure S1). First, the summer-autumn pattern consisted of a single large range area and two core areas; the two ranges coincided to a great extent but less so their respective core areas. Second, the winter-spring pattern consisted of four range areas and one or two core areas; the core areas coincided to a great extent and less so their respective ranges ( Table 1). The annual range of the Balkan chamois in Timfi Mt covered Annual range and core area of Balkan chamois in the study area deriving from merging seasonal ranges and core areas, respectively; human settlements, roads and wildlife refuges (hunting ban areas) are indicated.

Annual and Seasonal Ranges
Seasonal ranges presented a minimum in winter and a maximum in summer, with the maximum aggregation within the core area observed in autumn (Table 1). We found an important overlap (above 49%) among seasonal ranges. Two main seasonal patterns emerged ( Figure S1). First, the summer-autumn pattern consisted of a single large range area and two core areas; the two ranges coincided to a great extent but less so their respective core areas. Second, the winter-spring pattern consisted of four range areas and one or two core areas; the core areas coincided to a great extent and less so their respective ranges ( Table 1). The annual range of the Balkan chamois in Timfi Mt covered an area of 6491 ha, accounting for 25% of the study area; the annual core area was about half of the annual range ( Figure 2, Table 1). Table 1. Seasonal ranges and core areas (ha) of the Balkan chamois population in Timfi Mt (2001)(2002) and respective proportions (Range: seasonal range/annual range, Core: core seasonal range/seasonal range), and overlap of the seasonal ranges and respective core areas (in parenthesis).

Season
Area (

Population Demography
In 2002, the population density of the species was estimated to be 2 individuals/100 ha. The population has shown an important increasing trend of the order of 9.5% per year, increasing from 132 individuals in 2002 to 469 individuals in 2017 (Table 2). Kids were about one fourth of the overall population, the sex ratio was always in favour of females, and the fecundity rate remained high (0.6 to 0.8).

Habitat Selection
The ENFA results indicated a pronounced habitat selection behaviour by chamois, as both marginality and specialization scores were statistically significant (Monte Carlo, p < 0.001). The environmental conditions in which chamois were found showed a clear deviation from the average conditions in the available grids of the study area, as indicated by the high marginality scores obtained (Table S3). The species was recorded in a rather narrow interval of environmental conditions, being sensitive to their change, as the specialization scores showed a much greater variance of the environmental conditions in the study area (50-126 times) than within the chamois ranges. The scatterplots showed that the position of the used habitat was clearly distinct from the available habitat across the seasons (Figure 3). Considering the environmental variables contributing most to the marginality axis we found that the chamois always selected areas away from the human settlements, roads, and areas where hunting was allowed, avoiding also south-facing slopes (Figure 3, Table 3 and Table S4). This pronounced avoidance pattern was persistent across seasons and years, although a few more factors contributed to the seasonal habitat selection. In winter and spring, chamois selected north-facing slopes in pine forests. In summer and autumn, it preferred higher elevations with a substantial cover of rocks and screes, favouring east-facing or north-facing slopes, respectively. On the basis of the specialization axis, chamois did not tolerate any cover of scrub or agricultural and anthropogenic land in its habitat (Figure 3, Table 3 and Table S4).  Table 3 for variable abbreviations and Table S4 for variable coefficients.

Chamois Distribution and Hunting Ban Area of Greece
The first rough estimation for the overall hunting ban area in Greece in 2019 was 21,842.37 km 2 accounting for 16.54% of Greece. The area where hunting was banned for reasons of wildlife conservation was 9.16% of Greece. About 40% of the Balkan chamois distribution area falls within the hunting ban area of Greece (Table 4).  Table 3 for variable abbreviations and Table S4 for variable coefficients.

Chamois Distribution and Hunting Ban Area of Greece
The first rough estimation for the overall hunting ban area in Greece in 2019 was 21,842.37 km 2 accounting for 16.54% of Greece. The area where hunting was banned for reasons of wildlife conservation was 9.16% of Greece. About 40% of the Balkan chamois distribution area falls within the hunting ban area of Greece (Table 4). Table 3. Average values (A) or percent frequency of occurrence (F) of the environmental variables in the occupied grid cells (100 m × 100 m) by chamois across the six surveys, and in the available grid cells of the study area. Shaded cells indicate the top five variables with the highest marginality coefficients (Table S4). Private land Private fenced land n/a n/a n/a Geospatial data not available Vineyards, cultivated land etc. in case that hunting harms crops. n/a n/a n/a Geospatial data not available

Range Pattern
Our results demonstrated that the Balkan chamois population in Timfi Mt was confined to a small part of the study area (25%), despite the large availability of typical chamois habitats there, such as mountainous grasslands, rocks and pine forests [5]. It exhibited two distinct seasonal spatial patterns (winter-spring and summer-autumn), presenting a minimum range in winter, a maximum range in summer and a minimum core area range in autumn. The species also seems to follow the "typical range-use pattern": the animals restrict to a small winter range and start to move from spring to summer grounds expanding gradually their range [38][39][40][41], while in autumn they aggregate again due to rutting. This pattern is well known for the chamois population in cooler climates, where winter is the main stress period for animal survival, due mainly to harsher climate and to restricted food resources [42][43][44][45][46]. We note, however, that in the southernmost limit of Balkan chamois distribution, Giona Mt in Greece, chamois follow exactly the inverse pattern (Mediterranean range use pattern) from the Timfi population; it exhibits a minimum range and core area in the hot and dry stress period of summer, possibly due to global warming that further restricts fresh palatable food resources [47].

Population Trend and Demography
According to our findings, Timfi Mt hosts a large population, substantially increasing during the last two decades, and accounting for about one third of the national population. The former population density (2002) was particularly low (2 individuals/100 ha), whilst the current population density is unknown. If we assume that the annual range remained stable, the population density would reach 7 individuals/100 ha in 2017, which is comparable to those in other European mountains, usually lying between 6 and 10 individuals/100 ha, and in some cases reaching up to 24 or 40 individuals/ 100 ha [48][49][50][51]. The current population density should be lower than 7 individuals/100ha, as the population increase is expected to result in gradual range expansion.
The population seems to be healthy in terms of its demographic characteristics. The fecundity rate found of the study population falls within the usual rate of other European populations (0.55 to 0.85) [52,53]. However, we found a sex ratio in favour of females. This could be explained by the theory of higher male mortality due to high competition for mating [54,55], but in other areas, similar survival rates of the two sexes in Rupicapra rupicapra are reported [42,56]. One possible explanation for our study area could be the selective poaching towards males, under the scope to maintain females for next year's harvest, or missing males that preferred to remain in lower forested habitats during the rutting season, and hence be less detectable during autumn surveys.

Seasonal Habitat Selection
Chamois selected higher altitudes in Timfi Mt during summer and autumn. They followed the typical circular seasonal migration pattern: moving to higher altitudes in spring as the snow recedes in search of sprouts and fresh palatable grass [22,57], occupying gradually the highest available habitats during summer and autumn, and moving slowly towards the lower altitudes in winter. We argue that it is possible that global warming will modify the above typical circular pattern in the future. For instance, we have recently noticed that some individuals remain at lower forested habitats during summer, as already reported in other cases [57][58][59][60]. In our case, the animals were observed drinking water in the rivers at the bottom of gorges, from 2017 onwards, either in the study area (Aoos gorge) or in the vicinal Tzoumerka Mt to the south (Aracthos gorge), potentially in search of water resources and cooler microhabitats. We have also noticed in the year 2017 a substantial delay in some cases in leaving the high-altitude autumn grounds of the study area towards winter grounds, potentially due to the milder weather conditions prevailing. Although the impact of global warming on chamois habitat selection and behaviour was not among the objectives of our study, this issue needs further investigation, in particular for the more vulnerable chamois populations of southern Europe. Furthermore, the species avoided south-facing slopes across all seasons, preferring eastern-facing slopes in summer and mostly northern-facing slopes in winter and spring. In other parts of Europe, animals tend to stay less often on slopes with a southerly aspect [61], or avoid them only in summer and prefer them in winter to cope with the harsh weather conditions [57,62], which does not seem to be the case in the milder winter conditions of Timfi Mt. We argue that the strong avoidance pattern of southern-facing slopes could be attributed to global warming, as the animals may search cooler microhabitats throughout seasons. Such responses are of local character and largely depend on prevailing environmental contexts. For instance, the aspect did not seem to play an important role for the most southern population of Balkan chamois in Giona Mt, where the population employed a more dynamic strategy to adapt to the arid stressful conditions of summer (Mediterranean range pattern) rather than using different aspects of the mountain [47]. However, we do not attribute the selection of norther-facing slopes to global warming, but rather to the availability of escape terrains mostly in the northern parts of the mountains, as these microhabitats can function as adequate escape terrains for the species to evade human disturbance [63]. In our study area, steep slopes lie mainly in the northern aspects of Timfi Mt and extend even at lower altitudes in the pine zone, and less in the broadleaved forest zone that is characterized by milder topography. The greater availability of such escape terrain microhabitats in the pinewood zone could also explain why chamois selected pinewoods over broadleaved woods as winter grounds, a preference that is known from the southwestern Pyrenees [60]. Therefore, we argue that the availability of escape terrains may play an important role in the habitat selection process, as an underlying driver explaining the selection of northern aspects and pinewoods in the study area, rendering this issue worthy of further investigation.

Anthropogenic Risk-Avoidance Pattern
Our results supported the anthropogenic risk-avoidance hypothesis, as the species clearly selected remote areas away from roads, human settlements, and hunting zones. The pattern was strong and permanent in time. When even a small new road penetrates into wilderness or previous remote areas, it can still trigger a cascade of disturbances and impairments of ecosystems, a phenomenon known as the "contagious effect" of roads [64]. The impact of roads on wildlife is negative, reducing habitat quality, whilst ungulates, including chamois, show a general avoidance pattern to road infrastructure in different ecosystems of the world [65][66][67][68]. However, the small forest roads in mountainous ecosystems cannot be considered detrimental per se for chamois. They cannot really act as barriers for animal dispersion and its subsequent free gene flow, and the collision-mortality problem is reported to be minimal [69], presumably concerning only paved roads. However, roads are surrogates of human disturbance, and more specifically of poaching and hunting activity. The animals in Timfi Mt avoid roads in order to fend off poachers and hunters. This is in line with the national report on the conservation status of chamois in Greece, highlighting roads, poaching and hunting as the main high-impact threats to the species nationwide [10]. Such interconnection of roads to hunting and poaching is widely recognized, because increased public access provided by roads increases the intensity of hunting and poaching of ungulates in several parts of the world [66,70,71]. Poaching is considered a common phenomenon in the country, enhanced by the poor law enforcement and poaching control by the Greek State [8,72]. For instance, the rangers of the Epirus hunting association denounced to the Court of Justice a poaching incidence (six chamois shot in the study area in February 2017), where barcoding techniques were also used to confirm the crime [73]. The litigation is still ongoing and if it ends up in a conviction of the defendants, this would be the first such case in Greece. Furthermore, our results support that legal hunting constitutes a pressure to the species, a fact that could be attributed to the disturbance from shooting and hunting dogs.

Chamois Distribution and Hunting Grounds
According to our results, the distribution area of Balkan chamois is quite well covered (40%) by hunting ban areas and particularly by wildlife refuges (29%). Given the negative impact of hunting disturbance on the chamois population, the designation of such refuges and reserves has been quite successful for chamois population conservation. Further research is recommended for the effectiveness of those hunting ban areas to benefit other game and non-game species, since the impact of hunting in terms of harvest or disturbance on other wildlife populations is still poorly studied [74,75]. At national scale, the hunting grounds of Greece were estimated to be about 83% of the country. This percentage is similar with other countries in Europe, where the hunting system is completely different, including hunting grounds leasing and the implementation of adequate local hunting management plans [76], which is not the case in Greece. Our findings include a degree of uncertainty. This is the first study that attempted to scan the current hunting legislation and estimate the hunting grounds in Greece, on the basis of which the overlapping distribution area of chamois was calculated. The legislation is scattered, and fragmented and spatial data are largely lacking for delineating the Greek hunting grounds accurately. The estimated hunting ban areas might be overestimated in our methodology, when considering the category of banning hunting for reasons of public safety, or underestimated, when considering the cases for which no data were available. Therefore, we suggest an update and improvement of the hunting legislation towards a single and explicit legal frame that would provide among others a clear delineation of the hunting grounds of Greece, to be available as open maps for hunters and the broad public.

Conservation Implications
We showed that road avoidance is a major determinant of the Balkan chamois space-use behaviour, which may apply to other populations in its distribution range in Greece and in the Balkans. Our findings, therefore, provide support for a national roadless policy, with special focus on the chamois distribution range. The policy should involve banning of road expansion in crucial chamois habitats and controlled access to existing mountainous roads, as road removal is a well-recognized strategy for restoring habitat for large mammals [77]. Such a roadless policy would be beneficial for chamois, largely diminishing the poaching pressure on the species, and fostering its fast population recovery at national scale. It would also benefit other protected large mammals that are negatively affected by roads in Greece [74], and would improve the environmental performance of the country, since the fragmentation pressure of Greece is high [16] and has been recently reported to be rapidly increasing [78]. Furthermore, our results highlighted the negative role of hunting disturbance on chamois space use in Timfi Mt but also showed that an important part of the national chamois population is protected within hunting ban areas. Therefore, we suggest the maintenance and expansion of the nature reserves and wildlife refuges all over the country, towards increasing their cover in chamois habitats. Such a policy would greatly contribute to the species population recovery, and it might be beneficial for other mountainous wildlife populations.
The two above policies should be implemented by priority in the Northern Pindus mountain range. Our study proved the importance of Timfi Mt for Balkan chamois population recovery at a national scale, hosting about one third of the national population. Roadless areas and new wildlife refuges or nature reserves should be established by priority in the mountains of the Northern Pindus population block, to link the respective chamois sub-populations through ecological corridors [21]. Located in the heart of the Northern Pindus mountain, the chamois population of Timfi Mt has the potential to recolonize previously occupied areas and secure gene flow towards vicinal small and declining populations, if such an upward trend persists [21], rendering the positive impact of the above two measures high for the national population. The animals have gradually started to recolonize adjacent areas, such as Vikos gorge in the study area or Mitsikeli Mt, where chamois was extinct.
Finally, our key research findings provide scientific evidence for measures and policies that could be integrated into the forthcoming national action plan for the Balkan chamois. The plan is under Diversity 2020, 12, 124 14 of 18 development by the Greek State and is expected to undertake a suite of well-coordinated actions and policies, in order to improve the conservation status of the Balkan chamois in the country, as an annexed species of the Habitats' Directive.
Supplementary Materials: The following are available online at http://www.mdpi.com/1424-2818/12/4/124/s1: Table S1: Classification of the 25 habitat types of the study area into six broad habitat types. For each habitat type, we present its area (ha) and the respective cover (%) of the study area, as well as the habitat area (ha) that falls within chamois annual range and core area. The habitat typology follows the Habitats Directive (92/43/EC), noting priority habitat types of Annex I (*). Rivers and lakes considered in the broadleaved woods category, settlements considered under the agriculture category, Table S2 Table S4: Coefficients of the explanatory environmental variables in the 2D environmental space formed by the marginality axis (M) and the first specialization axis (S) of the ecological-niche factor analysis (ENFA) across the six surveys. The top five variables with the highest marginality scores and the variable with the highest specialization score are shaded. Abbreviations as in Table 3, Figure S1: Balkan chamois seasonal ranges (Fixed Kernel Density Estimator: 95% probability of occurrence) and seasonal core areas respectively in Timfi Mountain in 2002. The upper right corner presents the method of defining the probability of occurrence in the core areas (80% in winter, 85% in spring, 70% in summer, and 75% in autumn).