Habitat use of golden jackals (Canis aureus) in riverine areas of northern Bosnia and Herzegovina

Recently confirmed expansion of golden jackals (Canis aureus) into countries without any previous records poses questions, one of them focusing on the species potential and possibly varying habitat use. In this study, we investigated the presence and distribution of golden jackals in northern Bosnia and Herzegovina, where knowledge about golden jackal distribution and habitat use is still scarce. We used bioacoustic stimulation as a non-invasive tool to gather data on golden jackal presence. Habitat structures potentially selected by the species were assessed at 92 calling stations and used as input for binary logistic regression models. Our study area covered approximately 1150 km2, and bioacoustic stimulation within this area resulted in an estimated minimum relative group density of 3.5 territorial groups per 100 km2. We found territorial groups at distances between 15 and 38 km southwards from the river Sava but always within a maximum range of 3 km to perennial watercourses. Habitat analysis identified shrub vegetation and pastures as structures with a significant effect on the presence of resident golden jackals. Probability that golden jackals answered at calling stations increased with increasing surface area covered with pastures and shrubs. Distances between golden jackal territories and the nearest human settlement were relatively small. Our results indicate that structures like transitional woodland-shrubs and pastures, together with other potential influencing factors like local agricultural practices, low hunting pressure, diverse natural and anthropogenic food sources, could have benefited the settlement of golden jackals in the northern lowlands of Bosnia and Herzegovina.


Introduction
The golden jackal (Canis aureus Linnaeus, 1758) is a medium-sized carnivore which, considering its morphology and behaviour, is a typical member of the genus Canis (Clutton-Brock et al. 1976;Jhala 2004). Due to its tolerance to various types of habitats and its omnivorous diet, the golden jackal was able to colonize large areas in Europe (Trouwborst et al. 2015;Krofel et al. 2017). The distribution of this species in the Balkan region and its surroundings underwent significant changes in the twentieth century ; from a first expansion wave in the 1950s, with new records for example in Slovakia and Slovenia (Brelih 1955;Ferianc 1955) to a near extinction of the Pannonian population (Tóth et al. 2009), to massive resettlement in the Balkans and another wave of expansion to central and most recently northern Europe since the 1980s (Spassov 1989;Kryštufek et al. 1997;Tóth et al. 2009;Trouwborst et al. 2015;Krofel et al. 2017;Spassov and Acosta-Pankov 2019).
Knowledge about golden jackal distribution in Bosnia and Herzegovina (B&H) is very scarce. Since there was no reliable data for this country, Arnold et al. (2012) assumed that at least vagrant individuals should be present. In Trouwborst et al. (2015), B&H was listed as a country with reproducing golden jackal groups. Trbojević et al. (2018) confirmed resident golden jackals and observed an increase in hunting bag data between 2000 and 2016. Within their studied area in northern B&H, golden jackal density was estimated to be a minimum of 0.33 groups per 10 km 2 . For the two neighbouring countries Serbia and Croatia, Šálek et al. (2014) estimated 1.1 and 0.9 groups per 10 km 2 , respectively. In Bulgaria, densities as high as 2.83 groups per 10 km 2 were recorded (Acosta-Pankov et al. 2018).
There are first indications that golden jackal core distribution is located along the Sava river in northern B&H (Trbojević et al. 2018), but it still remains unclear, which variables drive jackal habitat selection. Therefore, the goals of our study were (1) to identify essential habitat structures for golden jackals in the riverine areas of northern B&H and (2) to analyse the distances of detected golden jackals to the nearest human settlement. We hypothesised that golden jackals can be confirmed near habitat structures that provide cover as well as protection from human persecution, e.g. shrubs or dense forests (Admasu et al. 2004;Giannatos et al. 2005;Šálek et al. 2014;Hatlauf et al. 2016). Secondly, we hypothesised that golden jackals can be detected in close distance, less than 1 kilometre, from the nearest human settlement (Macdonald 1979;Giannatos et al. 2005). Jackal occurrence was assessed with bioacoustic stimulation, a non-invasive and effective monitoring tool that has been used for surveys of several canids like African wild dogs (Lycaon pictus Temminck, 1820), grey wolves (Canis lupus Linnaeus, 1758) and golden jackals (Robbins and McCreery 2003;Giannatos et al. 2005;Comazzi et al. 2016;Hatlauf et al. 2017;Passilongo et al. 2017;Suter et al. 2017).

Study area
B&H is a mountainous country dominated by Dinaric Alps with lowland areas positioned along the major rivers. Moving north from the central parts of B&H, the landscape turns gradually into wide foothills and finally into lowlands. This area is characterized by a continental climate (Zurovec et al. 2015). Agricultural area makes 41.7% of the total surface area, and almost half of the arable land is left uncultivated. Less than 20% of the agricultural area, located mainly in the lowlands near the main river banks, is suited for intensive agriculture production (Volk 2010). Agriculture, in general, was massively influenced by the depopulation of rural areas induced by industrialisation (Zurovec et al. 2015). Small farms (up to 5 ha) represent about 80% of all agriculture farms. They are often split into more plots causing lower productivity and less efficiency. The most important crops are wheat and corn maize. Supported by the large proportion of grasslands in total (49%), livestock production is an essential segment of B&H agriculture (Volk 2010). The disposal of waste from livestock production is a commonly known problem in B&H, with almost 4000 t of animal tissue and 40 t of dead fish disposed yearly. Such waste is being buried near production sites, disposed with regular waste or in the rivers and surrounding area (Meško et al. 2011).
The vast majority of our study area was located in the flatlands of the northern part of B&H ( Fig. 1) along the valley of the river Sava with minor parts in hilly regions at a maximum distance of 29.35 km southwards to the river Sava. The whole study area covered approximately 1150 km 2 and stretched from the point where the river Sava enters B&H in the north-west up to the border with Serbia in the north-east.

Bioacoustic stimulation
Bioacoustic stimulation was performed in August and September 2016 on eight nights at 92 sampling points (calling stations). August is being considered one of the months with higher possible answer rates (Acosta-Pankov et al. 2018). Linear survey transects were created southwards of the river Sava using a geographical information system (QGIS 2012, version 2.16.0) and according to accessible roads. Since Croatia and Serbia, both neighbour countries of B&H, reported golden jackal presence (Kryštufek and Tvrtković 1990;Krofel 2008), our main and longest transect was positioned parallel to the border river Sava and contained 68 calling stations reaching from Donja Gradina in the north-west to the border with Serbia in the north-east. Furthermore, we positioned three short transects southwards from the main transect with the aim to inspect a potential dispersal of golden jackals towards areas with confirmed presence of grey wolves in central B&H (Trbojević 2016). The first one was located between the cities Derventa and Kotorsko and was composed of 12 calling stations. The second one encompassed seven calling stations and was located in the area of the Majevica mountain (918 m above sea level). The third transect contained five calling stations and was located between the cities Ormanica and Gracanica. These three transects had less calling stations because of difficult terrain, bad road connections and lack of time and personnel. We systematically set up a total of 92 calling stations along these transects at regular intervals of 3-4 km, following Giannatos et al. (2005), assuming the maximum human hearing distance is 2 km. Exceptions in the distance between single sampling points were made if cities or large villages interrupted the survey transect. In this case, calling stations were placed as close as possible, in front and after the interruption of the survey transect. In cases of landscape barriers like hills, the distance between calling stations was reduced to 1.5-2 km, since we considered these to function as sound barriers as well.
We used a GPS device (Garmin Ltd., 8200 Schaffhausen, Switzerland) to locate the sampling points in the field. Sometimes a small correction had to be made if the planned location was inappropriate (e.g. high-traffic roads, houses with dogs, minefields). The exact coordinates of the calling stations were saved to the GPS device. Bioacoustic monitoring was performed on the sampling points to count individual answers or group answers to the call.
Howling audio tracks (calls) were emitted with a 60-W Megaphone (MEG060, Vexus audio, 7602 Almelo, The Netherlands). We arrived at the first calling station 15 min before sunset and started the playback 5 min after sunset. We approached the location of every following calling station by car, turned the engine and lights off, exited the vehicle, and waited 5 min before playing the first track from the megaphone. The emitted tracks were recordings of a group of golden jackals (group call) and a single dominant male golden jackal (single call). The group call was played first, followed by a single call and lastly the group call again. A break of 5 min was kept between each of the three tracks to allow the audition of answer calls. At each calling station, we noted the following information: coordinates, date, time, type of answer (individual answer or group answer/s), direction of answer call/s, weather conditions, and additional information about potential sources of disturbance that could affect the audio monitoring (traffic, dogs barking, etc.). We assumed that groups that react to acoustic stimulation are territorial, as described by Giannatos et al. (2005).

Habitat analysis
The exact locations of all 92 calling stations were uploaded into the geographical information system (ESRI ArcMap 10.4.1). Then, we created buffers around these single points to analyse and quantify features of the landscape within the defined distance of a 2 km radius, with an area of about 12.5 km 2 that represents the possible home range size used by radio-tracked golden jackals (Giannatos 2004). Data on land cover types (see Table 1) were obtained from CORINE land cover 2012 (EEA 2018). We used the CORINE dataset as a classification of land cover types, but we did not verify these categories in the field, because landmine pollution did not allow detailed field observation within the whole study area. Finally, we quantified the spatial coverage of all land use categories (in km 2 ) within the investigated 92 buffer areas.

Data analysis
To analyse the effect of land cover on occurrence of golden jackals, we created a binary logistic regression using the statistical software R (version 3.3.1, R Core Team 2016). The binomial response variable answer indicated that either an individual golden jackal or a group was heard, or no answer was received. Before analysing data, the explanatory variables were tested for multicollinearity by calculating variance inflation factors (Zuur et al. 2010). Since correlation between land cover variables within the buffers around single calling stations (radius 2 km) was high (see Online Resource 1), and sample size within our study was limited (n = 92 calling stations), we had to remove variables before performing analyses. We excluded the following CLC classes from subsequent analyses: (a) beach, marsh, river and lakes, because most of our calling stations were located along the river Sava; (b) fruit, because only 4 out of the 92 calling stations were covered by the CLC class; and (c) indust and nonirigated, because they were highly correlated with discurban and complexcult, respectively (see Online Resource 1). We calculated a global model including the selected parameters and employed a multi-model inference approach (Burnham and Anderson 2002) to determine the most important variables. The covered surface area of land use categories (discurban, pastures, complexcult, agrinatur, forest, shrubs) and the distance to the nearest artificially surfaced areas (e.g. buildings, roads) were included as explanatory variables. Since calling stations were located in a riverine area, the mean distance to rivers was low (mean ± SE = 0.94 ± 0.63, range = 0.04-3.00 km) and the variable was not included in the model.
A set of models with all possible combinations of fixed effects was generated (R package MuMIn, Bartoń 2020) and graded according to their Akaike's information criterion with a correction for small sample sizes (AIC c ). A subset of models with ΔAIC c of less than 2 was obtained and model-averaged coefficients were extracted. In addition, relative importance (RI) ranging from 0 (not included in any model) to 1 (included Table 1 Overview and description of CORINE land cover (CLC) classes that are covering the study area. Variables in bold were finally included to analyse the habitat use of golden jackals (Heymann et al. 1994

Bioacoustic stimulation
Overall, we received golden jackal answers at 29 calling stations on a total area of 1150 km 2 , resulting in an estimated minimum relative group density of 3.5 territorial groups per 100 km 2 . Next to single howls (n = 7), we received answers from one up to four groups (mean ± SE = 1.8 ± 0.24 groups, n = 22) at single calling locations (see Table 2), resulting in an overall sampling success effort of 31.5%.

Habitat characteristics
Composition of land use types within the buffers around the calling sites was profoundly diverse (see Table 3). Areas within the 92 buffers were comprised of 2 to 9 different land use categories (mean ± SE = 5.7 ± 0.16). The vast majority of the studied area within buffers was covered by complexcult (24%), forest (22%), and nonirigated (20%). Pastures and shrubs covered only smaller proportions of the study area but were present in almost half of the buffers surrounding the calling stations (pastures: 6.6%, n = 46 calling stations, shrubs: 2.5%, n = 38 calling stations). Indust, fruit, marshes, beach and lake covered only small proportions (each below 1% of the whole study area). Land cover types occurring in most of the buffers surrounding the 92 calling stations were forest (n = 86), agrinatur (n = 85) and complexcult (n = 84). Distance from all the calling stations to the nearest urban settlements ranged between 0.0 and 20.6 km (mean ± SE = 3.39 ± 0.389 km). While the median distance to the nearest settlements was 1.68 km for negative calling stations, it was 3.1 km for the positive ones. The elevation at the chosen calling stations ranged between 77 m and 611 m above sea level. The subset of models explaining occurrence of golden jackals in the study area with ΔAIC c < 2 encompassed 3 models (see Table 4). The final averaged model included the variables pastures, shrubs, agrinatur, and forest. Discurban and complexcult were not useful predictors of occurrence of golden jackals. Probability that golden jackals answered at a calling station increased with increasing surface area covered with pastures (Estimate ± SE = 0.37 ± 0.15, RI: 1.00, 95% CI = 0.07, 0.66, see Fig. 2) and shrubs (Estimate ± SE = 1.02 ± 0.49, RI = 1.00, 95% CI = 0.05, 1.99, see Fig. 3). The relative importance of agrinatur (Estimate ± SE = 0.19 ± 0.16, RI = 0.33, 95% CI = −0.13, 0.50) and forest (Estimate ± SE = − 0.07 ± 0.11, RI = 0.20, 95% CI = −0.28, 0.14) was low and confidence intervals included zero, indicating little evidence that these predictors affect occurrence of golden jackals. The AUC values of the 3 models in the subset ranging between 0.74 and 0.76 indicated a fair discrimination between calling stations with golden jackal presence and absence.

Discussion
The bioacoustic method was successfully applied to analyse the occurrence of golden jackals in the study area in northern B&H (sampling effort success over 30%). Within the 92 calling stations, 29 were positive, 22 from which we received answers of groups (seven single answers). We detected a total of 40 presumed territorial groups in the survey area (the 7  single individual answers were not counted as a group). Minimum relative densities of territorial groups detected in this study are lower compared with those from northern Dalmatia (Krofel 2008), Peljesac peninsula (Krofel 2007), areas in Serbia and Romania (Šálek et al. 2014) or Bulgaria (Acosta-Pankov et al. 2018). Based on information from local hunters and residents, northern Bosnia seems to be the core area of golden jackal occurrence in this country with most of the groups located along the river Sava, similar to the denser jackal areas in the lowlands around the Danube in Serbia (Šálek et al. 2014). One possible reason for a diverse occurrence pattern in Bosnia could be the absence of wolves in the lowlands of northern B&H that might, on the other hand, be the limiting factor to jackal presence in other parts of the Balkan peninsula (Kryštufek and Tvrtković 1990;Krofel et al. 2017;Newsome et al. 2017). Still, we achieved to confirm six groups of golden jackals with territories southwards the river Sava, the furthest being 37 km towards central B&H. This dispersion could become a trend (southward) in the future due to habitat degradation by regular flooding events in the area along the Sava river (Boudaghpour et al. 2015).
Confirming our first hypothesis, the habitat analysis models show that areas covered with pastures and shrubs have a significant positive effect on the presence of golden jackals. Similar to previous habitat analyses from Serbia, that displayed an increased selection of shrubs as well, specifically when their coverage in the habitat was sparse (Šálek et al. 2014). The number of detected golden jackal groups decreased with an increasing proportion of arable land, whereas an increased presence of shrubs proved to influence golden jackal populations positively in Bulgaria (Šálek et al. 2014).
Our study shows that with only a slight increase of shrubs in the surveyed areas, the probability of golden jackal presence increased strongly (Fig. 2). Hence, if the area of shrubs surrounding the calling stations doubles from 2 km 2 (60% probability of presence) to 4 km 2 , the probability of golden jackal presence within this area rises to almost 100%. The lowlands in northern Bosnia are mainly used for agricultural purposes (Witmer and O'Loughlin 2009), already indicating potential habitat for golden jackals. After the war in B&H, which ended in 1995, almost 44% of the former arable land now lies idle (Volk 2010). Dense vegetation and single trees  . 2 The probability of golden jackal presence at calling stations was positively affected by pasture cover. For graphical illustration, simple logistic regression was performed (Estimate ± SE = 0.33 ± 0.15, p = 0.025) and the predicted probability was calculated afterwards Fig. 3 The probability of golden jackal presence at calling stations was positively affected by shrub cover. For graphical illustration, simple logistic regression was performed (Estimate ± SE = 0.87 ± 0.47, p = 0.062) and the predicted probability was calculated afterwards often overgrew these abandoned agricultural areas or have transformed into natural fallows (Witmer and O'Loughlin 2009;Volk 2010). Such habitats offer essential structures for golden jackals like thick cover for resting, hiding and building dens (Admasu et al. 2004;Giannatos et al. 2005). Furthermore, the species-rich natural pastures (Söderström et al. 2001;Hartel et al. 2013) as potential hunting grounds could be the reason why golden jackals seem to favour these areas. Anthropogenic food sources could explain the possible presence of golden jackals near human settlements, especially when shrub structures in the surrounding are present (Giannatos 2004). Our second hypothesis, regarding the possible detection of golden jackals within 1 km to human settlements, could be confirmed. Nevertheless, in more detail, our survey showed that while the median distance of negative calling stations was 1.68 km, it was 3.1 km for positive calling stations. This result might indicate a tendency to a stronger avoidance of human settlements compared to other regions, at least within studied area.
Our results reveal that golden jackals in northern B&H show preferences for habitat structures similar to those in surrounding countries. The shrub areas and pastures in the northern lowlands of B&H are a result of the agricultural policy and practices and seem to influence the presence of golden jackals positively by offering protection from humans (hunting pressure) and natural as well as anthropogenic food sources. Furthermore, the migration of the local population from rural areas, abandoned arable land, minefields, unsuitable hunting strategies, and improper waste management might represent further factors affecting the distribution of golden jackals in B&H directly or indirectly. We suggest the acoustic stimulation to be used country-wide to create a holistic picture (overview) of the golden jackal's distribution and relative density. Although golden jackals generally respond throughout the year, for future studies, we recommend the late summer as the best season for field surveys (next to a peak in answers from January to March), as it was suggested in the study of Acosta-Pankov et al. (2018) performed in Bulgaria. At this time of the year, wind conditions at night are favourable (low or no wind) and vegetation offers more cover for researchers, therefore less detection by golden jackals and hence a higher possibility of positive answers. Nevertheless, precautions should be made to avoid habituation to bioacoustic stimulation. For this reason, we suggest not to search for animals that could potentially move towards the calling station with a torch, as with this practice jackals might associate the presence of humans with the playback calls and avoid answering them in future surveys. Although costly, thermal imaging devices may serve as an alternative. Additionally, random visits of calling stations, only twice a year, are likely to avoid habituation.