Get 20M+ Full-Text Papers For Less Than $1.50/day. Subscribe now for You or Your Team.

Learn More →

Demography of a Eurasian lynx (Lynx lynx) population within a strictly protected area in Central Europe

Demography of a Eurasian lynx (Lynx lynx) population within a strictly protected area in Central... www.nature.com/scientificreports OPEN Demography of a Eurasian lynx (Lynx lynx) population within a strictly protected area in Central Europe 1* 2,3 3 4 4,5 Stefano Palmero , Elisa Belotti , Luděk Bufka , Martin Gahbauer , Christoph Heibl , 1,4,6 5 1,4,7 Joe Premier , Kirsten Weingarth‑Dachs & Marco Heurich Large carnivores promote crucial ecosystem processes but are increasingly threatened by human persecution and habitat destruction. Successful conservation of this guild requires information on long‑term population dynamics obtained through demographic surveys. We used camera traps to monitor Eurasian lynx between 2009 and 2018 in a strictly protected area in the Bohemian Forest Ecosystem, located in the core of the distribution of the Bohemian–Bavarian–Austrian lynx population. Thereby, we estimated sex‑specific demographic parameters using spatial capture– recapture (SCR) models. Over 48,677 trap nights, we detected 65 unique lynx individuals. Density increased from 0.69 to 1.33 and from 1.09 to 2.35 individuals/100 km for open and closed population SCR models, respectively, with corresponding positive population growth rates (mean = 1.06). Estimated yearly sex‑specific survival probabilities for the entire monitoring period were high (females 82%, males 90%) and per capita recruitment rate was low (females 12%, males 9%), indicating a low yearly population turnover. We ascertained an average number of recruits of 1.97 and a generation time of 2.64 years when considering resident reproducing females. We confirmed that reproduction in the study area took place successfully every year. Despite the overall increase in local lynx densities, the number of detected family groups remained constant throughout the study period. These results indicated that the strictly protected study area acts as a source for the multi‑use landscapes in its surroundings. In this first open population SCR study on lynx, we provide sex ‑specific demographic parameters that are fundamental information for lynx management in the study area as well as in similar contexts Europe‑wide. Large carnivores shape ecosystems through top-down control of herbivores and intraguild predation of meso- carnivores, which in turn can trigger trophic cascades . Their occurrence is therefore crucial for the functioning of ecosystems. However, large carnivore conservation is particularly challenging as carnivores’ food acquisition and large spatial requirements often lead them into conflict with human activities and their wide-ranging behaviour necessitate challenging transboundary management . Moreover, low densities and reproduction rates make these species vulnerable to the effects of human persecution . Consequently, a better understanding of the dynamics of large carnivore populations is imperative for successful conservation m anagement . Despite their importance, the temporal dynamics and environmental factors that drive populations are still poorly under- stood for many species. This is because only relatively few species are studied in the high detail needed for a thorough understanding. The necessity of sampling a large number of individuals over long periods makes it especially difficult to obtain a mechanistic understanding of why and how populations increase, decline, or go Faculty of Environment and Natural Resources, University of Freiburg, Tennenbacher Straße 4, 79106 Freiburg, Germany. Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Kamýcká 1176, 16521 Prague 6, Czech Republic. Department of Research and Nature Protection, Šumava National Park Administration, Sušická 399, 34192 Kašperské Hory, Czech Republic. Department of Visitor Management and National Park Monitoring, Bavarian Forest National Park, Freyunger Str. 2, 94481 Grafenau, Germany. Plant Biodiversity Research, Technische Universitӓt München, Emil-Ramann Straße 2, 85354 Freising-Weihenstephan, Germany. Leibniz Institute for Zoo and Wildlife Research (IZW), Alfred-Kowalke-Str. 17, 10315 Berlin, Germany. Inland Norway University of Applied Science Institute for Forest and Wildlife Management, Campus Evenstad, 2480 Koppang, Norway. email: st.palmero@gmail.com Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 1 Vol.:(0123456789) www.nature.com/scientificreports/ extinct. In most cases, even basic information, including demographic parameters such as abundance, density, survival and recruitment, are not available for populations of free-ranging animals, which makes conservation of threatened populations a hazardous game. Only a detailed time series of demographic parameters will allow a deeper understanding of population dynamics and reliable prediction of populations’ future developments, which is crucial for well-informed conservation management. Therefore, long-term data collection is required . Information on demographic parameters is challenging to collect for large carnivores because they usually exhibit low population densities, are primarily nocturnal and oe ft n live in areas rich in cover. The rising popu - larity of camera traps for wildlife monitoring in recent years has helped to overcome some of these d ifficulties . Moreover, such non-invasive devices allow simultaneous monitoring of different species and help to avoid stress- ful animal imm obilisation . Camera traps can provide high-quality pictures that enable the identification of naturally marked animals such as f elids . This has led to their extensive use in combination with capture–recap - ture (CR) methods for estimating demographic parameters of marked felids, e.g. tiger (Panthera tigris) , ocelot 10 11 12 (Leopardus pardalis) , jaguar (Panthera onca) and various lynx species (Lynx sp.) . Recently developed spatial capture–recapture (SCR) models oe ff r improvements over conventional non-spatial CR models for estimating demographic parameters because they also incorporate spatial information such as the location of individuals and traps, and habitat suitability . Closed population SCR models assume demographic closure, i.e. there is no emigration, immigration, mortality or reproduction, and are normally used to estimate abundance and density within one “session” . More recent open population SCR models can be applied across multiple sessions, which has the advantage of providing further parameters, such as survival, per capita recruitment and population growth rate , making them well suited for demographic analyses. In the last century, the Eurasian lynx (Lynx lynx), was eradicated across Central Europe but, following legal protection and population reintroductions, the species has since recolonised parts of its former ra nge . However, most of the reintroduced populations in Central Europe have remained isolated and small, mainly due to human disturbance and habitat f ragmentation . As a typical example, the Bohemian–Bavarian–Austrian population, which was reintroduced in the 1980s and has been considered stagnating in recent years, is today classified as "endangered". This is despite the availability of sufficient suitable habitat . From the perspective of a potential Central European metapopulation, the range of this population is located in a crucial area, but likely because of the low population size and spatial isolation, subadults cannot connect with neighbouring populations (e.g. 17 18 the Harz, Carpathian, and Alpine) , resulting in reduced genetic variability . Illegal killing is considered as the mortality cause most constraining the Bohemian–Bavarian–Austrian population and the protected areas in the region are therefore crucial factors for its persistence . The protection of source populations has been proposed as a strategy for recovering predator species, for example, the tiger (Panthera tigris) . However, whether pro- tected areas in Central Europe, including the study area, can host source populations with their limited size is poorly understood. In this study, we conducted a 10-year (2009–2018) demographic study in the core of the distribution of a Eurasian lynx (hereafter lynx) population within one of Central Europe’s largest strictly protected areas, the Bohemian Forest Ecosystem, using camera trapping and SCR methods. The results of this study help to improve our understanding of lynx demography in reintroduced populations. We used open population SCR models to estimate sex-specific demographic parameters such as abundance, density, survival probabilities, per capita recruitment rate and population growth rate. We expected higher density estimates for females since males have larger home r anges . We estimated density via both closed and open population SCR models to compare the reliability of the different methods for assessing population status and allow cross-comparison with existing studies that used conventional closed population methods. As lynx is a K-selected species, we expected generally high survival probabilities and a low per capita recruitment rate. Specifically, we expected lower male survival probabilities because they generally take higher risks, for example by getting closer to human activities to exploit higher prey densities in those a reas or patrolling their larger home ranges. Camera trapping data provided auxiliary information on reproductive parameters such as generation time and the average number of recruits, 22,23 which we predicted to be at around 2 years of age and 1.5 k ittens . Finally, we calculated the relative abundance index (RAI) to obtain information on the development of lynx’s main prey and a common mesopredator. Material and methods Study area. e B Th ohemian–Bavarian Forest is situated in Central Europe at the border between Austria, Czechia, and Germany, and includes two adjacent national parks: the Bavarian Forest National Park (BFNP) 2 2 (240  km ) in Germany and the Šumava National Park (SNP) (690  km ) in Czechia (Fig.  1). These protected areas are not divided by ecological barriers and therefore represent a continuous area in the core of the Bohe- mian–Bavarian–Austrian lynx population distribution. The BFNP is surrounded by the Bavarian Forest Natural 2 2 Park (3007  km ) and the SNP by the Šumava Protected Landscape Area (1000  km ), which together comprise the Bohemian Forest Ecosystem. Elevation ranges from 600 to 1456  m.a.s.l. and snow cover can persist for 5–8 months with greatest depths from January to March. e a Th rea is covered by a mixed mountainous forest composed mainly of Norway Spruce (Picea abies ), followed by European Beech (Fagus sylvatica) and Silver Fir (Abies alba) and hosts ungulate species such as roe deer (Capreolus capreolus), red deer (Cervus elaphus), wild boar (Sus scrofa) and moose (Alces alces) . In the BFNP roe 227 2 deer density ranges from 1.1 to 5 animals/km and red deer was estimated at 1.56 animals/km via coordinated counting at winter feeding stations. Densities of both species are higher in the S NP . Wildlife control within both national parks is conducted by trained staff and is limited to red deer and wild boar. Outside the national parks, roe deer, red deer and wild boar are hunted . e Th Bohemian-Bavarian-Austrian lynx population originates from 5 to 10 individuals illegally reintroduced in the 1970s and 17 individuals (11 males and 6 females) released officially in the 1980s . The total population Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 2 Vol:.(1234567890) www.nature.com/scientificreports/ CZ DE 2.7 x 2.7 Camera trapping grid Bavarian Forest Narional Park Šumava National Park Land cover Open area Forest cover Figure 1. Map of the study area with forest coverage including the Bavarian Forest National Park (BFNP) on the German side (DE), the Šumava National Park (SNP) on the Czech side (CZ) and the 2.7 × 2.7 km camera trapping grid in which sites were located in every second cell. The map was created using QGIS 3.4 https:// qgis. org/ it/ site/. size, which includes the wide surroundings of our study area (13,000  km ), was recently estimated as 97–143 29 2 (mean = 120) . Lynx density within the study area was estimated at 0.4–0.9 individuals/100  km in a preliminary study before the beginning of our study p eriod . Lynx is legally protected in all three countries but illegal kill- ing still occurs . The main prey species of lynx in the study area are roe deer and red deer with 80% and 17% of detected kills, r espectively . Camera trap monitoring. Data were collected between 2009 and 2018. e Th monitoring design was devel- oped specifically for lynx but also provided data on other species. The spatial organisation of camera traps, developed by Weingarth et  al. , consisted of a 2.7 × 2.7  km grid with camera traps situated in every second grid-cell and positioned along forest paths, roads and trails (Fig. 1). The design maximises detection probability and avoids gaps that might include female home ranges, i.e. a minimum of 122  km estimated locally via radio- telemetry . The camera trapping array occupied the entire BFNP and two-thirds of the SNP and measured 2 2 16 760  km , i.e. almost twice the mean male home range size (436  km ) measured through telemetry . This was in line with the camera trapping array size recommended to achieve adequate sample sizes for robust estimates . Since individual identification requires high-quality pictures of both animal flanks, most camera trapping sites included two opposing cameras (Cuddeback) with white flash. In the BFNP, camera traps functioned almost constantly for 10 “lynx years”. A lynx year is defined as the period from 01/05 to 30/04 the subsequent year, since kittens are typically born on M ay . In the SNP, most camera traps were only active from mid-September to the end of December and for technical reasons, camera trapping sites could not be active in year 2012. Camera trap- ping sites were sometimes slightly moved in the national parks at the beginning of each session and a few camera trapping sites were added in the study area over the entire monitoring period. We identified lynx individuals from images by comparing their unique coat patterns. When the individual identification posed difficulties, pictures were judged by at least one additional expert with long term experi- ence in this field. Poor quality images that precluded individual identification were discarded (~ 2% of captured images). Sex was determined by observing females with kittens or the genital area of the animal. Age could only be assessed for individuals first photographed as kittens and recaptured over the years. Additionally, we followed Weingarth et al. and assigned animals into two categories referring to their status: "juvenile" and "independent". The juvenile category included individuals < 1-year-old, i.e. kittens detected with their mother. Individuals in this category were excluded from the SCR analysis because of the high mortality of kittens . The independent category consisted of all individuals > 1-year-old (i.e. subadults and adults) and individuals of unknown age but with proof of independence. This category included floaters (non-resident or dispersing individuals), which were Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 3 Vol.:(0123456789) www.nature.com/scientificreports/ defined either as known juveniles of the previous year, in their second year of life, or as individuals of unknown origin that appeared for the first time in the study area whose age could not be determined. In both cases, these individuals held no territories and could either disperse from or settle in, the study area the following year. All individuals of unknown status were discarded from analyses. Spatial capture–recapture models. Demography of lynx was investigated using both open and closed population SCR modelling frameworks. Although the former method allows estimation of the full range of demographic parameters, we included the latter for backward compatibility with previous studies on lynx con- 35,36 ducted in Europe using closed population SCR models , as well as for comparison of both methods within one study. Following Weingarth et al. , we selected a time frame of 100 days from 15/09 to 24/12 of each year as a primary period to ensure demographic closure and a sufficient number of recaptures for robust estimates. The first primary periods were not placed in such an ideal period. Specifically, camera traps started functioning in November in both the two national parks in year 2009 and just in the SNP in year 2010 and 2011. In the fourth primary period (i.e. year 2012), as already mentioned, camera traps were not active in the SNP. We combined all primary periods according to a classical "robust design" , and thus performed open population SCR analysis. We defined one secondary period as one day and restricted the number of detections to at most one per site in any secondary period in line with a Bernoulli distribution, thereby reducing temporal autocorrelation. As such, our study included ten primary periods with 100 secondary periods each. The SCR method assumes the baseline detection probability g of any individual declines with the distance from its theoretical home range centre, the detection function scale σ, in the state space S, which should be at least 3σ, . We, therefore, created a rectangular state-space mask with a continuous bue ff r of 18,000  km around camera traps based on results of preliminary closed SCR analyses (i.e. σ ~ 3.5–4 km). We used a density-independent population growth model adapted for sex-specific demographic parameters . The population growth model was based on a spatial point process indicating the number and location of individuals at the initial population state (t = 1) and modelling abundance and distribution at the time t = 2 as a function of both survival probabilities and per capita recruitment rate. Using the R package “OpenPopSCR” , we ran ten pooled chains comprising 100,000 Markov chain Monte Carlo iterations each with an augmented observed population size M of 400 individuals, i.e. much greater than the overall number of independent individuals detected across all primary periods. We estimated combined and sex-specific yearly abundance and density by dividing abundance by the area of the state space. We used a spatially explicit movement model to estimate sex-specific yearly survival probabilities and per capita recruitment rate separately from emigration and immigration, respectively . For movement, we used a Markov activity centre relocation type with activity centres of individuals in primary period l + 1 centred around the activity centre in the first primary period l in the state-space according to a bivariate normal distribution. Estimates of per capita recruitment rate indicated the number of individuals per sex class added to total abundance each year. Realised sex-specific yearly population growth rate and sex ratio were estimated as derived parameters. The former was derived from the sex-specific abundances. The sex ratio indicated the probability of any individual being a female. Unknown sexes were considered as a latent covariate and estimated through the sex vector augmented to the length of M. The movement parameter indicated yearly activity centre relocation according to the Markov activity centre relocation type. All point estimates were obtained using posterior modes and interval estimates were calculated through 95% highest posterior density (HPD) intervals at the 2.5% and 97.5% quantiles of the posterior distribution . We used the Gelman-Rubin statistic to assess convergence of open population SCR parameters through the “gelman.diag” function from the R package “coda” . This calculates the potential scale reduction factor for each parameter with upper confidence limits (CI), whereby parameters with a 95% upper CI substantially above 1 are considered to lack c onvergence . Closed population SCR models in a maximum likelihood framework were fitted using the R package “secr” . 9 44 We defined one sampling occasion as 5  days and used detector type proximity according to a Bernoulli distri- bution to estimate the combined density of males and females. To allow method comparison, we used for each session a continuous buffer of 18,000  km as for open population SCR models and created rectangular state-space masks accordingly (Supplementary Figs. S1–S10). First, we fitted a model M which assumes baseline detection probability g and detection function scale σ to be equal for all individuals . Secondly, we included sex as a covariate for both parameters and fitted two different models for each session, i.e. one keeping baseline detection probability g constant and testing sex on the detection function scale σ and one testing sex on both parameters for each session. We compared Akaike’s Information Criterion corrected for small sample sizes (AICc) to deter- mine the best model for our data and kept those with ΔAICc < 2 for model-averaging . Uknown sexes were considered as a latent covariate by performing a hybrid mixture m odel . This model consists of a combination of latent (e.g. missing sexes) and known classes (e.g. male and female) and estimates the mixing proportion of the different classes by calculation of the parameter “pmix”, which corresponds to the sex ratio. Additionally, this parameter allows class-specific modelling of the detection parameters for investigating sex die ff rences, for example, in the detection probability. Statistical significance was evaluated using the 95% CI. Reproductive parameters and age distribution. Based on camera trapping data, we investigated reproductive parameters such as generation time and number of recruits. The first was calculated as the mean age of resident reproducing females (i.e. females with kittens) at their first documented reproduction. Age was assessed only for individuals first photographed as kittens and recaptured over time. We, therefore, only included resident reproducing females of known age in the calculation. In contrast, no resident reproducing females were discarded when calculating the average number of recruits. This was calculated as the overall mean number of kittens with the mother at the onset of winter (November–January) . Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 4 Vol:.(1234567890) www.nature.com/scientificreports/ Status Known mortality causes Independent (German side only) Ee ff ctive Female (with Total Camera trapping Unique lynx Traffic Session kittens) Male Unknown sex (floaters) Juvenile Unknown trapping sites nights events accidents Illegal killing 2009 6 (4) 3 0 9 (–) 8 3 55 1728 23 0 0 2010 6 (4) 6 1 13 (4) 8 0 62 3813 52 0 0 2011 8 (3) 7 1 16 (7) 5 4 66 4124 82 0 0 2012 7 (3) 5 1 13 (5) 5 0 31 2621 48 0 1 2013 8 (6) 7 1 16 (0) 8 0 66 5885 86 0 1 2014 9 (6) 8 0 18 (2) 9 0 65 6151 94 1 0 2015 11 (2) 10 0 22 (6) 5 0 64 6022 135 1 1 2016 12 (5) 11 3 27 (6) 11 0 65 5964 144 2 0 2017 11 (4) 14 3 29 (10) 6 0 69 5957 121 2 0 2018 11 (7) 13 1 25 (4) 15 0 69 6412 130 1 0 Table 1. Summary of the lynx information obtained every primary period, including: the number of independent females and resident reproducing females (with kittens), males, individuals of unknown sex, total number of independent individuals including floaters (– indicates not available), juveniles, individuals of unknown status, camera trapping sites, effective trapping nights, unique lynx events (at most one detection per site in any secondary period) of independent individuals and documented mortality cases with known causes. e fir Th st three sessions had reduced number of effective trap nights because the monitoring started later in November in both the national parks in year 2009 and just in the SNP in years 2010 and 2011. For the session 2012, due to technical reasons, lynx camera trapping only took place on the German side of the study area. e n Th umber of juveniles referred to the entire lynx year. Finally, we also included all the independent individuals of known age and sex detected during the primary periods for investigating the age distribution of the population. Relative abundance index. When species cannot be individually identified, camera trapping data can 47 23 be used to calculate abundance indices . The relative abundance index (RAI) , which expresses the number of captures per monitoring day, is widely used in wildlife research . We matched the session length for RAI to that of the SCR analysis and considered lynx’s main prey species, i.e. roe deer and red deer, and red fox (Vulpes vulpes). To avoid temporal autocorrelation, we defined capture events with a threshold of 10  min and bootstrapped the RAI estimates by randomly selecting one camera trap per site at a time over 5000 replicates to reduce camera and site capture biases. This aspect of the study was conducted in the BFNP only. Use of experimental animals. No animals were used or handled in the study. Results Camera trap monitoring. We detected 65 unique independent individuals (25 males, 28 females and 12 individuals of unknown sex) in a total of 48,677 trap nights (Table 1). For the reasons mentioned in the meth- ods, the standardised session length of 100 days was not achieved in the first three primary periods (i.e. year 2009, 2010 and 2011) and in the fourth primary period (i.e. year 2012) camera traps were not active in the SNP, potentially resulting in estimates with higher uncertainty. The overall number of trapping nights was partially reduced during each season by technical problems such as camera trap failure, theft or snowfall. Considering the primary periods, only ten pictures (between 0 and 2% depending on the primary period) did not allow individual identification because of poor quality, while the remaining pictures were identified with certainty. We observed an increase in numbers of juveniles and independent individuals including males, females and indi- viduals of unknown sex, while the number of resident reproducing females (i.e. females with kittens) remained quite stable. We could determine the number of floaters in most years, as well as their origin. Around 33% of the detected floaters were born in the study area while around 67% were born outside or had unknown origin. Compared to the other categories of independent lynx, the number of floaters was the value that fluctuated the most (from 0 in session 2013 to 10 in session 2017, see Table 1). Apart from three and four individuals in the first and third sessions, respectively, the status (i.e. juvenile or independent) was determined. The number of unique lynx events (at most one detection per site in any secondary period) including independent individuals increased by around six times over the monitoring period. For the German side, we were able to confirm human- related mortality causes of some detected individuals in the vicinity of the study area, which showed an increase in traffic accidents and occasional illegal killing. In particular, three adults, two subadults and two juveniles were killed in traffic accidents in the BFNP, while two adults and one juvenile of unknown sex were killed illegally outside the national park in Germany. Spatial capture–recapture models. e co Th mbined density estimates for each primary period using open population SCR models increased from 0.69 (95% HPD intervals 0.50–1.10) to 1.33 (95% HPD inter- Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 5 Vol.:(0123456789) www.nature.com/scientificreports/ Combined Male Female Closed SCR 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 Figure 2. Posterior modes and 95% highest posterior density (HPD) intervals of open population spatial capture–recapture (SCR) models referring to combined, male and female lynx densities. 95% HPD intervals Parameter Sex Estimate Lower Upper M 0.01 0.01 0.02 Detection probability F 0.01 0.01 0.01 M 4.07 3.80 4.38 Detection function scale F 3.96 3.66 4.32 M 0.90 0.79 0.96 Survival probabilities F 0.82 0.72 0.91 M 0.09 0.04 0.14 Per capita recruitment rate F 0.12 0.07 0.20 Movement C 4.75 4.03 5.65 Table 2. Posterior modes and 95% highest posterior density (HPD) intervals of open population spatial capture–recapture (SCR) model referring to male (M) and female (F) lynx yearly baseline detection probability g , detection function scale σ (km), survival probabilities and per capita recruitment rate. Per capita recruitment rate indicates the number of individuals of each sex added per year per total abundance. The yearly movement (km) was calculated for sexes combined (C). vals 1.05–1.79) individuals/100  km while, as expected, sex-specific density estimates were on average higher for females (Fig. 2). However, the overlapping 95% HPD intervals suggest there was no statistically significant difference between sexes. Combined abundance ranged from 38 (95% HPD intervals 26–60) to 75 (95% HPD intervals 58–100) individuals, that of males from 17 (95% HPD intervals 9–31) to 36 (95% HPD intervals 21–50) and that of females from 21 (95% HPD intervals 10–38) to 37 (95% HPD intervals 26–58). Concerning closed population SCR models, estimates of density ranged from 1.09 (SE 0.35) to 2.36 (SE 0.79) individuals/100  km (Table 1; Supplementary Table S2) and fluctuated strongly compared to open population SCR models (Fig.  2). Abundance ranged from 47.94 (SE 11.03) to 121.25 (SE 42.45) (Supplementary Table S2). Density estimates from each method fell just within the respective 95% CI or HPD of the other method, apart from the year 2012, indi- cating estimates were comparable and suggesting no statistically significant difference across methods. Further results of closed population SCR models are shown in the appendix (Supplementary Tables S2, S3). All Gelman–Rubin diagnostic statistics had a 95% upper CI < 1.1 indicating convergence was reached. The parameter posterior modes (Table 2) of the yearly baseline detection probability g were equal across sexes. The yearly per capita recruitment rate was higher for females. The yearly survival probabilities of males, i.e. 0.90 (95% HPD intervals 0.79–0.96), were higher than those of females, i.e. 0.82 (95% HPD intervals 0.72–0.91). Males showed a higher yearly detection function scale σ, due to larger home ranges. However, considering the overlap of 95% HPD intervals, none of the parameters suggested a statistically significant difference across sexes. The realised population growth rate was 1.06 (range 1.01–1.10), that of males was 1.07 (range 1.00–1.12) and that of females was 1.05 (range 1.00–1.14), suggesting there was no predominant sex in the population. The probability of any individual to be a female was quite stable over years ranging from 0.51 (95% HPD intervals 0.36–0.65) to 0.56 (95% HPD intervals 0.33–0.76), meaning the sex ratio was slightly skewed towards females, with, however, no apparent statistically significant die ff rence. Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 6 Vol:.(1234567890) Individuals/100km2 0.00.5 1.01.5 2.0 2.5 3.0 3.5 www.nature.com/scientificreports/ Events Session Roe deer Red deer Red fox Camera trapping sites Ee ff ctive trapping nights 2009 12 12 136 29 1238 2010 28 56 306 31 2786 2011 29 87 291 31 2914 2012 52 52 315 31 2621 2013 57 57 315 31 2865 2014 57 86 401 30 2865 2015 110 110 274 29 2745 2016 57 172 58 29 2871 2017 27 189 189 29 2698 2018 27 191 273 29 2732 Table 3. Summary of the available information about roe deer, red deer and red fox obtained every primary period, including number of events, camera trapping sites and effective trapping nights. The number of events refers to the mean number of events resulting from all iterations in which one camera trap per site was sampled at a time. roe deer red deer red fox 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 Figure 3. Estimates and standard deviations of the relative abundance index (RAI) of red fox and lynx’s prey species roe deer and red deer calculated for the BFNP only. Reproductive parameters and age distribution. Regarding reproductive parameters, we included 14 resident reproducing females of known age and calculated a mean generation time of 2.64 years. One individual was first detected with kittens when it was seven years old and thus considered an outlier. The average number of recruits was 1.97 (range 1–3) and was calculated from information on 22 detected females (Supplementary Table S4). Age and sex could be confirmed for a total of 24 independent individuals (15 females and 9 males) and were included in the age distribution pyramid (Supplementary Fig. S11). The oldest individual of known age was a 10-year-old male individual. The number of females at each age was equal to or higher than that of males except for ages 8–10. Relative abundance index. e n Th umber of camera trapping sites used to calculate the RAIs ranged from 29 to 31 and effective trapping nights was 26,335 (Table  3). The RAI of red fox ranged from 0.07 (SD 0.01) to 0.14 (SD 0.02), that of red deer from 0.01 (SD 0.00) to 0.07 (SD 0.01), that of roe deer from 0.01 (SD 0.00) to 0.04 (SD 0.01) (Fig. 3). Fox abundance increased from 2009 to 2014 and decreased in the remaining years, overall showing an oscillating trend. Red deer and roe deer slightly increased until 2015, aer w ft hich the former almost doubled while the latter decreased to earlier values. Discussion Our modelling approach revealed that the number of independent lynx detected in the Bohemian–Bavarian Forest has increased over the past decade, with concurrently increasing abundance, density and positive popu- lation growth rates. Survival probabilities were high and per capita recruitment rate was low indicating a low yearly population turnover. Reproductive parameters such as generation time and average litter size consider- ing all resident reproducing females indicated successful reproduction occurred every year, and the number of Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 7 Vol.:(0123456789) Relative abundance index (RAI) 0.00 0.05 0.10 0.15 0.20 www.nature.com/scientificreports/ reproducing females remained stable throughout the study period despite the increase in the total number of detected independents. All these findings indicate that the protected areas act as an important, stable source area for the Bohemian–Bavarian–Austrian lynx population. Comparison of spatial capture–recapture methods. To our knowledge, this is the first demographic study on lynx comprising annual camera trapping over a decade and using open population SCR models. These models can provide a wide range of demographic parameters useful for long-term lynx monitoring. In addition, estimates from open population SCR models were not prone to inter-annual fluctuations that ae ff cted closed population SCR models and were generally lower, due to lower abundances compared to those from the lat- ter (Table 1; Supplementary Table S2) divided by an equal buffer area. However, higher estimates with greater uncertainty, concerning the years 2009 and 2012, for which data were least complete, resulted from issues when accounting for varying efforts, i.e. the exact number of days the cameras were working. Allowing for varying effort has the benefit that detection parameters are related to a specific unit of effort thus being unbiased . In the open population SCR models, the activity of each camera trap was informed over the primary periods. Con- cerning closed population models, in the year 2009, it was not possible to include the period from September to November because no cameras were functioning, resulting in zero occasions. While in the year 2012, the activity of the cameras in the SNP could not be specified because the devices were not in the field in the yearly camera trapping array that year. Therefore, data gaps resulted in biased estimates in this case because some spatio-tem- poral requirements, namely adequate sample size and number of recaptures, were not fulfilled. Closed popula- tion SCR density estimates of primary periods for which it was possible to account for varying effort, including years 2010 and 2011 for which data were also not complete, were comparable to those from open population SCR models since their 95% CI and HPD intervals overlapped. One advantageous feature of open population SCR methods in a Bayesian framework is their ability to deal with incomplete detection (e.g. non-annual moni- toring). This is possible because they produce posterior distributions of the demographic parameters that incor - porate the uncertainty resulting from data gaps by using information derived from other primary periods . This highlights the robustness of this method, which we recommend for future studies and monitoring, especially when incomplete detections occur . Nonetheless, closed population SCR models still represent a reliable and conventional tool for abundance and density assessments . Spatial capture–recapture models. Our density estimates of both open (0.69–1.33 individuals/100  km ) and closed population SCR models (1.02–2.39 individuals/100  km ) were on the same order of magnitude as those resulting from other closed population SCR studies conducted on both reintroduced and autochthone populations across Europe. For example, the reintroduced populations of the French Jura and Vosges Mountains, 2 35 with the lowest density estimate reported in Europe (0.24–0.91 individuals/100  km ) , or the Swiss Alps (1.47 2 36 and 1.38 individuals/100  km ) , and the autochthone population of the Western Carpathians (0.26–1.85 indi- 2 50 viduals/100  km of suitable habitat) . The only outlier is a subpopulation in southwest Asia (Turkey) with a high lynx density of 4.20 individuals/100  km . However, lynx in Turkey live under different ecological conditions compared to those in Western and Central Europe, feeding mainly on lagomorphs instead of ungulates which leads to smaller home range sizes . In our study, the yearly baseline detection probability g was equal between sexes (Table 2). Although males are generally more active in order to patrol their larger home ranges, females with kittens hunt at a higher rate, which can result in augmented ac tivity , and might explain similarities across sexes for this parameter. The yearly detection function scale σ was not significantly higher for males compared to females (Table  2). This sug- 20 10 gests similar home range sizes, which stands in contrast to expectations for l ynx and most felids . However, the restricted and seasonal session period was chosen for demographic closure and not appropriate for annual home range estimation Furthermore, many detected individuals were not residents thus resulting in potentially biased estimates for the parameter in question. e Th combined yearly survival probabilities for independent individuals reached 85%. This high value is likely because our estimates come from strictly protected areas where lynx have a higher chance of survival. However, none of the detected individuals used the territory included in the protected areas exclusively, thus potentially exposing themselves to a higher risk. The survival estimated through open population SCR methods in our study site is among the highest reported for lynx. In the Western Carpathians, Dula et al. conducted an SCR survey on an autochthone lynx population across multiple seasons and found an overall apparent survival (i.e. including emigration) of 63% with high human-related mortality occurring in the area. The only other values of survival rates available for comparisons come from telemetry studies because almost all of the SCR surveys conducted on lynx in Europe based on camera trapping data used closed population models and thus did not provide information on this parameter. In Switzerland, Breitenmoser-Würsten et al. found an overall survival rate of 76% for adults and 53.3% for subadults. In Poland, the survival rate was only 63% , considering subadults and adults combined. In three Scandinavian study sites, Andrén et al. found lower survival rates for subadults (70, 77 and 71%) than adults (87, 91 and 84%), likely due to lower mortality related to vehicle collision and hunting in adulthood. Comparisons to SCR results are however difficult given the contrasting underlying meth- odologies. Against our expectations, the survival probabilities were higher for males, though not significantly. This was due to a higher overall number of apparent survival events for males, which consist of the number of consecutive detections over years including gaps during which the animal was alive but not detected. A similar open population SCR study, conducted on a low-density ocelot population in Belize, also showed no significant differences in sex-specific survival probabilities, despite they found probabilities of 0.86 for females and 0.78 for males, respectively . e a Th uthors suggested the statistical power was not enough to detect significant differences between sexes in these parameters, although they were able to determine sex for a large number of adult ocelots Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 8 Vol:.(1234567890) www.nature.com/scientificreports/ (n = 322). These are probably also the reasons why we could not prove statistically significant die ff rences in per capita recruitment rate between sexes (Table 2). Camera trapping does not allow assessment of the fates of all disappearing individuals and thus does not provide information on natural mortality. However, our auxiliary findings on dead animals show an increasing number of traffic accidents in recent years (Table  1), indicating that the prevention of lynx vehicle collisions represents an important management action for the population. We could not assess the actual impact of illegal killing since carcasses are seldom found, but considering the high survival of the lynx individuals monitored in the study area, illegal killing seems to play no important role within the study site. However, a high poaching rate in the surroundings of our study area was suggested by a modelling approach, which underlines the impor- tance of the protected areas for the survival of the p opulation . These results are consistent with other areas in 5 54 Europe e.g. S witzerland , and P oland were traffic accidents and poaching have been found as the main cause of lynx mortality. Concerning limitations of the open population SCR models, the movement parameter had poor mixing, which is in line with Gardner et al. , meaning Markov chain Monte Carlo iterations slowly converged to unbiased posterior distributions. This is the reason why the parameter was estimated for both sexes combined and it was necessary to run relatively long chains (100,000 iterations). We attempted to separate survival from emigration and recruitment from immigration using open popula- tion SCR models with a spatially explicit model for activity centre relocation . However, the discrimination ability of these models depends on how accurately the movement model can describe the actual activity centre relocation between primary periods . This would require a larger camera trapping array to detect dispersing movements of large carnivores. We, therefore, could not reliably separate emigration and immigration from survival and recruitment, respectively. However, it is reasonable to assume that true survival and recruitment are at least as high as the value we estimated when the underestimation due to the inclusion of emigration and immigration is considered. Reproductive parameters and age distribution. Generation time was calculated as the mean age of resident reproducing females at their first documented reproduction, potentially resulting in overestimated val- ues if previous litters went undetected by camera traps. This may be the case for individuals living close to the boundaries of the monitored area, such as the outlier we detected. In Scandinavia, generation time of lynx was investigated across different environments and ranged from two to more than three years with the highest values in the northern territories due to slower life cycles and a different feeding regime and thus body mass . Similar results were found in other areas of S candinavia . We were not able to fully assess the proportion of reproduc- tive females because sex was not determined for all detected individuals. However, our mean generation time of 2.64 years fall in the range of the Scandinavian studies. We found a rare case of a reproducing female first reproducing at one year of age (Supplementary Table S3). The first reported case of a 1-year-old female Eurasian lynx breeding in the wild was recently observed in the Bohemian–Bavarian–Austrian population range, outside the study area . The authors hypothesised this could be related to very high prey density or high turnover due to poaching. The average number of recruits of 1.97 (range 1–3) was comparable to that reported in Scandinavia. Gaillard et al. , investigated the number of recruits related to litter size and found higher values for multiparous females (i.e. females with multiple reproduction events) ranging approximately between 0.5 and 2 with a litter size between one and four kittens. However, we were not able to determine litter size through camera trapping. Regarding age distribution, the maximum age determined with certainty in the study area was 10 years, which was constrained by the duration of the study. In Switzerland, Breitenmoser-Würsten et al. found individuals 14–15 years of age through telemetry. The exact age can only be determined by camera trapping if an individual was first photographed as a kitten. As such, some of the individuals already classified as independent in the first monitoring session reached higher ages than those reported, even if their exact age remained unknown. Relative abundance index. According to the RAI results, roe deer abundance in the German part of the study site remained stable even though they were not hunted since 2012. Contrastingly, the red deer population strongly increased (Fig. 3). As lynx in the study area mainly feed on roe deer and are normally only able to kill red deer calves and subadults, they have a stronger effect on the abundance of roe deer than red deer . The stable RAI for roe deer might therefore indicate a limiting effect of lynx, a result consistent with observations from the 58 59 same and other a reas . Concerning red fox, intraguild predation of lynx on this species has already been reported in die ff rent areas 60–62 of Scandinavia . The fox RAI shows a decrease with increasing lynx numbers, but the trend is not clear. In the study area, the red fox was found as only 1% of lynx k ills meaning other factors such as food availability and diseases could be driving red fox dynamics. RAI results were not related statistically to SCR models as they could not be used as an individual covariate or as a trap covariate since this aspect of the study was limited to the BFNP only. Therefore the RAI results should be interpreted cautiously. Although RAI can be ae ff cted by biases attributed to changes in detection p robability , we believe this method well represented the abundance trends of the species in question since the study design was uniform over the years. Population development and prospects. For lynx, we found increasing abundance and density, as well as positive population growth rates, high survival probabilities and low per capita recruitment, resulting in a low yearly population turnover. This suggests that lynx is not strongly ae ff cted by human-related mortality within the study area. This contrasts with Duľa et al. , who found fluctuating density estimates in the Western Carpathians due to illegal killing and traffic collisions. The comparatively stable development in our study area was accompanied by a relatively constant number of family groups (Table 1), i.e. females accompanied by their Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 9 Vol.:(0123456789) www.nature.com/scientificreports/ kittens, which form the most stable structure of the population and are therefore of particular importance. On the other hand, the number of individuals belonging to the least stable part of the population, i.e. floaters likely dispersing through the study area, highly fluctuated throughout the study period (from a minimum of 0 in 2013 to a maximum of 10 in 2017). Furthermore, the increase (or decrease) of the total number of detected individuals from one monitoring season to the following mainly corresponded with an increase (or decrease, respectively) in the number of detected floaters (Table  1). This suggests that the number of floaters visiting the study area in the given year played a major role in determining local abundance lynx for that year. Among all floaters detected in this study, only a lower proportion of individuals originated inside the study area (~ 33%). This reveals potential changes in the Bohemian–Bavarian–Austrian lynx population. Based on coordinated transboundary camera trapping within an area of 13,000  km , the estimated entire population size showed a positive trend and confirmed lynx presence and regular reproduction in a large part of the entire population range, including its outskirts . Thus, fluctuations in the number of floaters may be related to what takes place outside of the protected study area. In particular, increases in the number of floaters reaching the study area from 1 year to the following may be related to a local increase in survival of juveniles and yearlings in part of the Bohemian–Bavarian–Austrian population range. This may be partially related to a local decrease in poaching, possibly as a result of law enforcement and long-term public relations campaigns aimed at improving accept- ance in the wake of high-profile poaching incidents. These findings suggest that the overall conditions of the Bohemian–Bavarian–Austrian population have at least partially improved in recent years and this do not only concern the protected areas lying within this population’s range. This is fundamental for the long term genetic viability of a potential future Central European metapopulation as the conservation of species with large spatial requirements cannot rely on protected areas a lone . Since data on poaching are almost absent, the causes lead- ing to this positive trend of the population remain speculation. However, other potential explanations such as immigration from other populations or variability in reproduction rates within our study area can be excluded since the Bohemian–Bavarian–Austrian is an isolated population . Conclusions High survival probabilities and regular reproduction thanks to local protection indicate that the protected areas act as a stable source of lynx for the wider distribution range of the Bohemian–Bavarian–Austrian population, which is in line with previous studies . Thus our results show the high value of the protected areas for the persistence of the Bohemian–Bavarian–Austrian As proposed for other large carnivores , we conclude that protected areas can be considered as strategic areas for lynx recovery across its range because full protection increases survival probabilities and ensure stable reproduction. However, we must consider that the management of the surrounding landscapes is fundamental for the long term survival of the population, as protected areas in Central Europe are too small to encompass enough animals to maintain healthy populations over the long term. Our study revealed crucial demographic parameters of a lynx population that can be used for improving con- servation and management plans. Finally, our results stress the importance of long-term systematic monitoring as a basis for the understanding of the population dynamics of large carnivore populations and recommend the use of open population SCR methods to achieve such aims. Data availability Data supporting the findings of this study are not included due to conservation concerns. Data may, however, be available from the authors upon reasonable request and with permission of both Administrations of both the BFNP and SNP. Received: 16 February 2021; Accepted: 22 September 2021 References 1. Ripple, W. J. et al. Status and ecological ee ff cts of the world’s largest carnivores. Science 343, 1241484–1241484 (2014). 2. Treves, A. & Karanth, K. U. Human–carnivore conflict and perspectives on carnivore management worldwide. Conserv. Biol. 17, 1491–1499 (2003). 3. Linnell, J. D. C. & Boitani, L. Building biological realism into wolf management policy: The development of the population approach in Europe. Hystrix Ital. J. Mammal. 23, 80–91 (2011). 4. Heurich, M. et al. Illegal hunting as a major driver of the source-sink dynamics of a reintroduced lynx population in Central Europe. Biol. Conserv. 224, 355–365 (2018). 5. Breitenmoser-Würsten, C., Vandel, J.-M., Zimmermann, F. & Breitenmoser, U. Demography of lynx Lynx lynx in the Jura Moun- tains. Wildl. Biol. 13, 381–392 (2007). 6. Clutton-Brock, T. & Sheldon, B. C. Individuals and populations: The role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends Ecol. Evol. 25, 562–573 (2010). 7. O’Connell, A., Nichols, J. D. & Karanth, K. U. Camera Traps in Animal Ecology: Methods and Analyses. (Springer Tokyo, 2011). 8. Noss, A. J. et al. A Camera trapping and radio telemetry study of lowland tapir (Tapirus terrestris) in Bolivian Dry Forests. Tapir Cons. 12, 9 (2003). 9. Karanth, K. U. & Nichols, J. D. Estimation of tiger densities in India using photographic captures and recaptures. Ecology 79, 11 (1998). 10. Satter, C. B., Augustine, B. C., Harmsen, B. J., Foster, R. J. & Kelly, M. J. Sex‐specific population dynamics of ocelots in Belize using open population spatial capture–recapture. Ecosphere 10, e02792 (2019). 11. Silver, S. C. et al. The use of camera traps for estimating jaguar Panthera onca abundance and density using capture/recapture analysis. Oryx 38, 148–154 (2004). 12. Zimmermann, F., Breitenmoser-Würsten, C., Molinari-Jobin, A. & Breitenmoser, U. Optimizing the size of the area surveyed for monitoring a Eurasian lynx (Lynx lynx) population in the Swiss Alps by means of photographic capture–recapture. Integr. Zool. 8, 232–243 (2013). 13. Royle, J. A., Chandler, R. B., Sollmann, R. & Gardner, B. Spatial Capture–Recapture. (Elsevier, 2014). Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 10 Vol:.(1234567890) www.nature.com/scientificreports/ 14. Chandler, R. B. & Clark, J. D. Spatially explicit integrated population models. Methods Ecol. Evol. 5, 1351–1360 (2014). 15. Kaczensky, P . et al. Status, management and distribution of large carnivores—Bear, lynx, wolf and wolverine in Europe (Europe- anCommission, 2013). 16. Magg, N. et al. Habitat availability is not limiting the distribution of the Bohemian–Bavarian lynx Lynx lynx population. Oryx 50, 742–752 (2016). 17. Müller, J. et al. Protected areas shape the spatial distribution of a European lynx population more than 20 years ae ft r reintroduc - tion. Biol. Conserv. 177, 210–217 (2014). 18. Bull, J. K. et al. The effect of reintroductions on the genetic variability in Eurasian lynx populations: The cases of Bohemian–Bavar - ian and Vosges–Palatinian populations. Conserv. Genet. 17, 1229–1234 (2016). 19. Walston, J. et al. Bringing the tiger back from the brink—The six percent solution. PLoS Biol. 8, e1000485 (2010). 20. Schmidt, K., Jędrzejewski, W. & Okarma, H. Spatial organization and social relations in the Eurasian lynx population in Bialowieza Primeval Forest, Poland. Acta Theriol. (Warsz.) 42, 289–312 (1997). 21. Bunnefeld, N., Linnell, J. D. C., Odden, J., van Duijn, M. A. J. & Andersen, R. Risk taking by Eurasian lynx (Lynx lynx) in a human- dominated landscape: Ee ff cts of sex and reproductive status. J. Zool. 270, 31–39 (2006). 22. Gaillard, J.-M., Nilsen, E. B., Odden, J., Andrén, H. & Linnell, J. D. C. One size fits all: Eurasian lynx females share a common optimal litter size. J. Anim. Ecol. 83, 107–115 (2014). 23. Nilsen, E. B., Linnell, J. D. C., Odden, J., Samelius, G. & Andrén, H. Patterns of variation in reproductive parameters in Eurasian lynx (Lynx lynx). Acta Theriol. (Warsz.) 57, 217–223 (2012). 24. O’Brien, T. G., Kinnaird, M. F. & Wibisono, H. T. Crouching tigers, hidden prey: Sumatran tiger and prey populations in a tropical forest landscape. Anim. Conserv. 6, 131–139 (2003). 25. Cailleret, M., Heurich, M. & Bugmann, H. Reduction in browsing intensity may not compensate climate change effects on tree species composition in the Bavarian Forest National Park. For. Ecol. Manag. 328, 179–192 (2014). 26. Heurich, M. et al. Country, cover or protection: What shapes the distribution of red deer and roe deer in the Bohemian Forest Ecosystem?. PLoS ONE 10, e0120960 (2015). 27. van Beeck Calkoen, S. T. S. et al. The blame game: Using eDNA to identify species-specific tree browsing by red deer (Cervus elaphus) and roe deer (Capreolus capreolus) in a temperate forest. For. Ecol. Manag. 451, 117483 (2019). 28. Wöl, fl M. et al. Distribution and status of lynx in the border region between Czech Republic, Germany and Austria. Acta e Th riol. 46, 181–194 (2001). 29. Mináriková, T. et al. Lynx monitoring report for Bohemian–Bavarian–Austrian lynx population for lynx year 2017 (INTERREG Central Europe, 2019). 30. Weingarth, K. et al. First estimation of Eurasian lynx (Lynx lynx) abundance and density using digital cameras and capture–recap- ture techniques in a German national park. Anim. Biodivers. Conserv. 35, 197–207 (2012). 31. Belotti, E. et al. Patterns of lynx predation at the interface between protected areas and multi-use landscapes in Central Europe. PLoS ONE 10, e0138139 (2015). 32. Tobler, M. W. & Powell, G. V. N. Estimating jaguar densities with camera traps: Problems with current designs and recommenda- tions for future studies. Biol. Conserv. 159, 109–118 (2013). 33. Zimmermann, F., Breitenmoser-Würsten, C. & Breitenmoser, U. Natal dispersal of Eurasian lynx ( Lynx lynx ) in Switzerland. J. Zool. 267, 381 (2005). 34. Andrén, H. et al. Survival rates and causes of mortality in Eurasian lynx (Lynx lynx) in multi-use landscapes. Biol. Conserv. 131, 23–32 (2006). 35. Gimenez, O. et al. Spatial density estimates of Eurasian lynx (Lynx lynx) in the French Jura and Vosges Mountains. Ecol. Evol. 9, 11707–11715 (2019). 36. Pesenti, E. & Zimmermann, F. Density estimations of the Eurasian lynx (Lynx lynx) in the Swiss Alps. J. Mammal. 94, 73–81 (2013). 37. Weingarth, K. et al. Hide and seek: Extended camera-trap session lengths and autumn provide best parameters for estimating lynx densities in mountainous areas. Biodivers. Conserv. 24, 2935–2952 (2015). 38. Pollock, K. H. A capture–recapture design robust to unequal probability of capture. J. Wildl. Manag. 46, 752 (1982). 39. Augustine, B. benaug/OpenPopSCR. (2019). https:// github. com/ benaug/ OpenP opSCR. 40. Ergon, T. & Gardner, B. Separating mortality and emigration: Modelling space use, dispersal and survival with robust-design spatial capture–recapture data. Methods Ecol. Evol. 5, 1327–1336 (2014). 41. Schaub, M. & Royle, J. A. Estimating true instead of apparent survival using spatial Cormack–Jolly–Seber models. Methods Ecol. Evol. 5, 1316–1326 (2014). 42. Plummer, M., Best, N., Cowles, K. & Vines, K. CODA: Convergence diagnosis and output analysis for MCMC. R News 6, 7–11 (2005). 43. Brooks, S. P. & Gelman, A. General methods for monitoring convergence of iterative simulations. J. Comput. Graph. Stat. 7, 434–455 (1998). 44. Efford, M. secr 4.1—Spatially explicit capture–recapture in R. (2019). https://cra n.micr os oft.co m/sn apsh ot/2019- 12- 24/ w eb/p acka ges/ secr/ vigne ttes/ secr- overv iew. pdf. 45. Burnham, K. P. & Overton, W. S. Robust estimation of population size when capture probabilities vary among animals. Ecology 60, 927–936 (1979). 46. Burnham, K. P., Anderson, D. R. & Burnham, K. P. Model Selection and Multimodel Inference: A Practical Information–e Th oretic Approach (Springer, 2002). 47. O’Brien, T. G. Abundance, density and relative abundance: A conceptual framework. In Camera Traps in Animal Ecology (eds O’Connell, A. F. et al.) 71–96 (Springer Japan, 2011). https:// doi. org/ 10. 1007/ 978-4- 431- 99495-4_6. 48. Rovero, F. & Zimmermann, F. Camera Trapping for Wildlife Research (Pelagic Publishing Ltd, 2016). 49. Augustine, B. C. et al. Sex-specific population dynamics and demography of capercaillie (Tetrao urogallus L.) in a patchy environ- ment. Popul. Ecol. 62, 80–90 (2020). 50. Duľa, M. et al. Multi-seasonal systematic camera-trapping reveals fluctuating densities and high turnover rates of Carpathian lynx on the western edge of its native range. Sci. Rep. 11, 9236 (2021). 51. Avgan, B., Zimmermann, F., Güntert, M., Arıkan, F. & Breitenmoser, U. The first density estimation of an isolated Eurasian lynx population in southwest Asia. Wildl. Biol. 20, 217–221 (2014). 52. Mengüllüoğlu, D., Ambarlı, H., Berger, A. & Hofer, H. Foraging ecology of Eurasian lynx populations in southwest Asia: Conserva- tion implications for a diet specialist. Ecol. Evol. 8, 9451–9463 (2018). 53. Heurich, M. et al. Activity patterns of Eurasian lynx are modulated by light regime and individual traits over a wide latitudinal range. PLoS ONE 9, e114143 (2014). 54. Jedrzejewski, W . et al. Population dynamics (1869–1994), demography, and home ranges of the lynx in Bialowieza Primeval Forest (Poland and Belarus). Ecography 19, 122–138 (1996). 55. Gardner, B., Sollmann, R., Kumar, N. S., Jathanna, D. & Karanth, K. U. State space and movement specification in open population spatial capture–recapture models. Ecol. Evol. 8, 10336–10344 (2018). 56. López-Bao, J. V. et al. Eurasian lynx t fi ness shows little variation across Scandinavian human-dominated landscapes. Sci. Rep. 9, 8903 (2019). 57. Engleder, T. et al. First breeding record of a 1-year-old female Eurasian lynx. Eur. J. Wildl. Res. 65, 17 (2019). Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 11 Vol.:(0123456789) www.nature.com/scientificreports/ 58. Heurich, M. et al. Selective predation of a stalking predator on ungulate prey. PLoS ONE 11, e0158449 (2016). 59. Andrén, H. & Liberg, O. Large impact of Eurasian lynx predation on roe deer population dynamics. PLoS ONE 10, e0120570 (2015). 60. Elmhagen, B. & Rushton, S. P. Trophic control of mesopredators in terrestrial ecosystems: Top-down or bottom-up?. Ecol. Lett. 10, 197–206 (2007). 61. Wikenros, C. et al. Fear or food—Abundance of red fox in relation to occurrence of lynx and wolf. Sci. Rep. 7, 9059 (2017). 62. Helldin, J. O., Liberg, O. & Glöersen, G. Lynx (Lynx lynx) killing red foxes (Vulpes vulpes) in boreal Sweden? Frequency and population ee ff cts. J. Zool. 270, 657–663 (2006). 63. Sollmann, R., Mohamed, A., Samejima, H. & Wilting, A. Risky business or simple solution—Relative abundance indices from camera-trapping. Biol. Conserv. 159, 405–412 (2013). 64. Linnell, J. D. C., Kaczensky, P., Wotschikowsky, U., Lescureux, N. & Boitani, L. Framing the relationship between people and nature in the context of European conservation: Relationship between people and nature. Conserv. Biol. 29, 978–985 (2015). Acknowledgements The study was funded by the German Academic Exchange Service (DAAD) and allowed by the Administra- tions of both the BFNP and SNP. Financial support was also provided by EU Program Interreg IV (Objective 3 Czech Republic—the Independent State of Bavaria, Grant no. 18) and by the World Wide Fund for Nature Section Germany (Grant no. 533500/53000/10120). We greatly thank Ben Augustine who kindly shared open population SCR models in the R package OpenPopSCR and supported the following analysis as well as helped interpret results. We sincerely thank all the colleagues of both the BFNP and the SNP who facilitated the study. In particular, we thank Nicolò Borgianni who helped with part of the analyses. Author contributions This submission represents original work carried out by the authors. All authors have contributed substantially to the manuscript and have approved the manuscript and agree to its submission. Conceptualisation: S.P., J.P., M.H., methodology: S.P., C.H., J.P., M.H., data collection: M.G., K.W., E.B., L.B., formal analysis and investigation: S.P., C.H., writing—original draft preparation: S.P., J.P., writing—review and editing: J.P., E.B., M.H., funding acquisition: S.P, E.B., L.B., M.H. Funding Open Access funding enabled and organized by Projekt DEAL. Competing interests The authors declare no competing interests. Additional information Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598- 021- 99337-2. Correspondence and requests for materials should be addressed to S.P. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. © The Author(s) 2021 Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 12 Vol:.(1234567890) http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

Demography of a Eurasian lynx (Lynx lynx) population within a strictly protected area in Central Europe

Loading next page...
 
/lp/springer-journals/demography-of-a-eurasian-lynx-lynx-lynx-population-within-a-strictly-l21cHPEV0G

References (72)

Publisher
Springer Journals
Copyright
Copyright © The Author(s) 2021
eISSN
2045-2322
DOI
10.1038/s41598-021-99337-2
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports OPEN Demography of a Eurasian lynx (Lynx lynx) population within a strictly protected area in Central Europe 1* 2,3 3 4 4,5 Stefano Palmero , Elisa Belotti , Luděk Bufka , Martin Gahbauer , Christoph Heibl , 1,4,6 5 1,4,7 Joe Premier , Kirsten Weingarth‑Dachs & Marco Heurich Large carnivores promote crucial ecosystem processes but are increasingly threatened by human persecution and habitat destruction. Successful conservation of this guild requires information on long‑term population dynamics obtained through demographic surveys. We used camera traps to monitor Eurasian lynx between 2009 and 2018 in a strictly protected area in the Bohemian Forest Ecosystem, located in the core of the distribution of the Bohemian–Bavarian–Austrian lynx population. Thereby, we estimated sex‑specific demographic parameters using spatial capture– recapture (SCR) models. Over 48,677 trap nights, we detected 65 unique lynx individuals. Density increased from 0.69 to 1.33 and from 1.09 to 2.35 individuals/100 km for open and closed population SCR models, respectively, with corresponding positive population growth rates (mean = 1.06). Estimated yearly sex‑specific survival probabilities for the entire monitoring period were high (females 82%, males 90%) and per capita recruitment rate was low (females 12%, males 9%), indicating a low yearly population turnover. We ascertained an average number of recruits of 1.97 and a generation time of 2.64 years when considering resident reproducing females. We confirmed that reproduction in the study area took place successfully every year. Despite the overall increase in local lynx densities, the number of detected family groups remained constant throughout the study period. These results indicated that the strictly protected study area acts as a source for the multi‑use landscapes in its surroundings. In this first open population SCR study on lynx, we provide sex ‑specific demographic parameters that are fundamental information for lynx management in the study area as well as in similar contexts Europe‑wide. Large carnivores shape ecosystems through top-down control of herbivores and intraguild predation of meso- carnivores, which in turn can trigger trophic cascades . Their occurrence is therefore crucial for the functioning of ecosystems. However, large carnivore conservation is particularly challenging as carnivores’ food acquisition and large spatial requirements often lead them into conflict with human activities and their wide-ranging behaviour necessitate challenging transboundary management . Moreover, low densities and reproduction rates make these species vulnerable to the effects of human persecution . Consequently, a better understanding of the dynamics of large carnivore populations is imperative for successful conservation m anagement . Despite their importance, the temporal dynamics and environmental factors that drive populations are still poorly under- stood for many species. This is because only relatively few species are studied in the high detail needed for a thorough understanding. The necessity of sampling a large number of individuals over long periods makes it especially difficult to obtain a mechanistic understanding of why and how populations increase, decline, or go Faculty of Environment and Natural Resources, University of Freiburg, Tennenbacher Straße 4, 79106 Freiburg, Germany. Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Kamýcká 1176, 16521 Prague 6, Czech Republic. Department of Research and Nature Protection, Šumava National Park Administration, Sušická 399, 34192 Kašperské Hory, Czech Republic. Department of Visitor Management and National Park Monitoring, Bavarian Forest National Park, Freyunger Str. 2, 94481 Grafenau, Germany. Plant Biodiversity Research, Technische Universitӓt München, Emil-Ramann Straße 2, 85354 Freising-Weihenstephan, Germany. Leibniz Institute for Zoo and Wildlife Research (IZW), Alfred-Kowalke-Str. 17, 10315 Berlin, Germany. Inland Norway University of Applied Science Institute for Forest and Wildlife Management, Campus Evenstad, 2480 Koppang, Norway. email: st.palmero@gmail.com Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 1 Vol.:(0123456789) www.nature.com/scientificreports/ extinct. In most cases, even basic information, including demographic parameters such as abundance, density, survival and recruitment, are not available for populations of free-ranging animals, which makes conservation of threatened populations a hazardous game. Only a detailed time series of demographic parameters will allow a deeper understanding of population dynamics and reliable prediction of populations’ future developments, which is crucial for well-informed conservation management. Therefore, long-term data collection is required . Information on demographic parameters is challenging to collect for large carnivores because they usually exhibit low population densities, are primarily nocturnal and oe ft n live in areas rich in cover. The rising popu - larity of camera traps for wildlife monitoring in recent years has helped to overcome some of these d ifficulties . Moreover, such non-invasive devices allow simultaneous monitoring of different species and help to avoid stress- ful animal imm obilisation . Camera traps can provide high-quality pictures that enable the identification of naturally marked animals such as f elids . This has led to their extensive use in combination with capture–recap - ture (CR) methods for estimating demographic parameters of marked felids, e.g. tiger (Panthera tigris) , ocelot 10 11 12 (Leopardus pardalis) , jaguar (Panthera onca) and various lynx species (Lynx sp.) . Recently developed spatial capture–recapture (SCR) models oe ff r improvements over conventional non-spatial CR models for estimating demographic parameters because they also incorporate spatial information such as the location of individuals and traps, and habitat suitability . Closed population SCR models assume demographic closure, i.e. there is no emigration, immigration, mortality or reproduction, and are normally used to estimate abundance and density within one “session” . More recent open population SCR models can be applied across multiple sessions, which has the advantage of providing further parameters, such as survival, per capita recruitment and population growth rate , making them well suited for demographic analyses. In the last century, the Eurasian lynx (Lynx lynx), was eradicated across Central Europe but, following legal protection and population reintroductions, the species has since recolonised parts of its former ra nge . However, most of the reintroduced populations in Central Europe have remained isolated and small, mainly due to human disturbance and habitat f ragmentation . As a typical example, the Bohemian–Bavarian–Austrian population, which was reintroduced in the 1980s and has been considered stagnating in recent years, is today classified as "endangered". This is despite the availability of sufficient suitable habitat . From the perspective of a potential Central European metapopulation, the range of this population is located in a crucial area, but likely because of the low population size and spatial isolation, subadults cannot connect with neighbouring populations (e.g. 17 18 the Harz, Carpathian, and Alpine) , resulting in reduced genetic variability . Illegal killing is considered as the mortality cause most constraining the Bohemian–Bavarian–Austrian population and the protected areas in the region are therefore crucial factors for its persistence . The protection of source populations has been proposed as a strategy for recovering predator species, for example, the tiger (Panthera tigris) . However, whether pro- tected areas in Central Europe, including the study area, can host source populations with their limited size is poorly understood. In this study, we conducted a 10-year (2009–2018) demographic study in the core of the distribution of a Eurasian lynx (hereafter lynx) population within one of Central Europe’s largest strictly protected areas, the Bohemian Forest Ecosystem, using camera trapping and SCR methods. The results of this study help to improve our understanding of lynx demography in reintroduced populations. We used open population SCR models to estimate sex-specific demographic parameters such as abundance, density, survival probabilities, per capita recruitment rate and population growth rate. We expected higher density estimates for females since males have larger home r anges . We estimated density via both closed and open population SCR models to compare the reliability of the different methods for assessing population status and allow cross-comparison with existing studies that used conventional closed population methods. As lynx is a K-selected species, we expected generally high survival probabilities and a low per capita recruitment rate. Specifically, we expected lower male survival probabilities because they generally take higher risks, for example by getting closer to human activities to exploit higher prey densities in those a reas or patrolling their larger home ranges. Camera trapping data provided auxiliary information on reproductive parameters such as generation time and the average number of recruits, 22,23 which we predicted to be at around 2 years of age and 1.5 k ittens . Finally, we calculated the relative abundance index (RAI) to obtain information on the development of lynx’s main prey and a common mesopredator. Material and methods Study area. e B Th ohemian–Bavarian Forest is situated in Central Europe at the border between Austria, Czechia, and Germany, and includes two adjacent national parks: the Bavarian Forest National Park (BFNP) 2 2 (240  km ) in Germany and the Šumava National Park (SNP) (690  km ) in Czechia (Fig.  1). These protected areas are not divided by ecological barriers and therefore represent a continuous area in the core of the Bohe- mian–Bavarian–Austrian lynx population distribution. The BFNP is surrounded by the Bavarian Forest Natural 2 2 Park (3007  km ) and the SNP by the Šumava Protected Landscape Area (1000  km ), which together comprise the Bohemian Forest Ecosystem. Elevation ranges from 600 to 1456  m.a.s.l. and snow cover can persist for 5–8 months with greatest depths from January to March. e a Th rea is covered by a mixed mountainous forest composed mainly of Norway Spruce (Picea abies ), followed by European Beech (Fagus sylvatica) and Silver Fir (Abies alba) and hosts ungulate species such as roe deer (Capreolus capreolus), red deer (Cervus elaphus), wild boar (Sus scrofa) and moose (Alces alces) . In the BFNP roe 227 2 deer density ranges from 1.1 to 5 animals/km and red deer was estimated at 1.56 animals/km via coordinated counting at winter feeding stations. Densities of both species are higher in the S NP . Wildlife control within both national parks is conducted by trained staff and is limited to red deer and wild boar. Outside the national parks, roe deer, red deer and wild boar are hunted . e Th Bohemian-Bavarian-Austrian lynx population originates from 5 to 10 individuals illegally reintroduced in the 1970s and 17 individuals (11 males and 6 females) released officially in the 1980s . The total population Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 2 Vol:.(1234567890) www.nature.com/scientificreports/ CZ DE 2.7 x 2.7 Camera trapping grid Bavarian Forest Narional Park Šumava National Park Land cover Open area Forest cover Figure 1. Map of the study area with forest coverage including the Bavarian Forest National Park (BFNP) on the German side (DE), the Šumava National Park (SNP) on the Czech side (CZ) and the 2.7 × 2.7 km camera trapping grid in which sites were located in every second cell. The map was created using QGIS 3.4 https:// qgis. org/ it/ site/. size, which includes the wide surroundings of our study area (13,000  km ), was recently estimated as 97–143 29 2 (mean = 120) . Lynx density within the study area was estimated at 0.4–0.9 individuals/100  km in a preliminary study before the beginning of our study p eriod . Lynx is legally protected in all three countries but illegal kill- ing still occurs . The main prey species of lynx in the study area are roe deer and red deer with 80% and 17% of detected kills, r espectively . Camera trap monitoring. Data were collected between 2009 and 2018. e Th monitoring design was devel- oped specifically for lynx but also provided data on other species. The spatial organisation of camera traps, developed by Weingarth et  al. , consisted of a 2.7 × 2.7  km grid with camera traps situated in every second grid-cell and positioned along forest paths, roads and trails (Fig. 1). The design maximises detection probability and avoids gaps that might include female home ranges, i.e. a minimum of 122  km estimated locally via radio- telemetry . The camera trapping array occupied the entire BFNP and two-thirds of the SNP and measured 2 2 16 760  km , i.e. almost twice the mean male home range size (436  km ) measured through telemetry . This was in line with the camera trapping array size recommended to achieve adequate sample sizes for robust estimates . Since individual identification requires high-quality pictures of both animal flanks, most camera trapping sites included two opposing cameras (Cuddeback) with white flash. In the BFNP, camera traps functioned almost constantly for 10 “lynx years”. A lynx year is defined as the period from 01/05 to 30/04 the subsequent year, since kittens are typically born on M ay . In the SNP, most camera traps were only active from mid-September to the end of December and for technical reasons, camera trapping sites could not be active in year 2012. Camera trap- ping sites were sometimes slightly moved in the national parks at the beginning of each session and a few camera trapping sites were added in the study area over the entire monitoring period. We identified lynx individuals from images by comparing their unique coat patterns. When the individual identification posed difficulties, pictures were judged by at least one additional expert with long term experi- ence in this field. Poor quality images that precluded individual identification were discarded (~ 2% of captured images). Sex was determined by observing females with kittens or the genital area of the animal. Age could only be assessed for individuals first photographed as kittens and recaptured over the years. Additionally, we followed Weingarth et al. and assigned animals into two categories referring to their status: "juvenile" and "independent". The juvenile category included individuals < 1-year-old, i.e. kittens detected with their mother. Individuals in this category were excluded from the SCR analysis because of the high mortality of kittens . The independent category consisted of all individuals > 1-year-old (i.e. subadults and adults) and individuals of unknown age but with proof of independence. This category included floaters (non-resident or dispersing individuals), which were Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 3 Vol.:(0123456789) www.nature.com/scientificreports/ defined either as known juveniles of the previous year, in their second year of life, or as individuals of unknown origin that appeared for the first time in the study area whose age could not be determined. In both cases, these individuals held no territories and could either disperse from or settle in, the study area the following year. All individuals of unknown status were discarded from analyses. Spatial capture–recapture models. Demography of lynx was investigated using both open and closed population SCR modelling frameworks. Although the former method allows estimation of the full range of demographic parameters, we included the latter for backward compatibility with previous studies on lynx con- 35,36 ducted in Europe using closed population SCR models , as well as for comparison of both methods within one study. Following Weingarth et al. , we selected a time frame of 100 days from 15/09 to 24/12 of each year as a primary period to ensure demographic closure and a sufficient number of recaptures for robust estimates. The first primary periods were not placed in such an ideal period. Specifically, camera traps started functioning in November in both the two national parks in year 2009 and just in the SNP in year 2010 and 2011. In the fourth primary period (i.e. year 2012), as already mentioned, camera traps were not active in the SNP. We combined all primary periods according to a classical "robust design" , and thus performed open population SCR analysis. We defined one secondary period as one day and restricted the number of detections to at most one per site in any secondary period in line with a Bernoulli distribution, thereby reducing temporal autocorrelation. As such, our study included ten primary periods with 100 secondary periods each. The SCR method assumes the baseline detection probability g of any individual declines with the distance from its theoretical home range centre, the detection function scale σ, in the state space S, which should be at least 3σ, . We, therefore, created a rectangular state-space mask with a continuous bue ff r of 18,000  km around camera traps based on results of preliminary closed SCR analyses (i.e. σ ~ 3.5–4 km). We used a density-independent population growth model adapted for sex-specific demographic parameters . The population growth model was based on a spatial point process indicating the number and location of individuals at the initial population state (t = 1) and modelling abundance and distribution at the time t = 2 as a function of both survival probabilities and per capita recruitment rate. Using the R package “OpenPopSCR” , we ran ten pooled chains comprising 100,000 Markov chain Monte Carlo iterations each with an augmented observed population size M of 400 individuals, i.e. much greater than the overall number of independent individuals detected across all primary periods. We estimated combined and sex-specific yearly abundance and density by dividing abundance by the area of the state space. We used a spatially explicit movement model to estimate sex-specific yearly survival probabilities and per capita recruitment rate separately from emigration and immigration, respectively . For movement, we used a Markov activity centre relocation type with activity centres of individuals in primary period l + 1 centred around the activity centre in the first primary period l in the state-space according to a bivariate normal distribution. Estimates of per capita recruitment rate indicated the number of individuals per sex class added to total abundance each year. Realised sex-specific yearly population growth rate and sex ratio were estimated as derived parameters. The former was derived from the sex-specific abundances. The sex ratio indicated the probability of any individual being a female. Unknown sexes were considered as a latent covariate and estimated through the sex vector augmented to the length of M. The movement parameter indicated yearly activity centre relocation according to the Markov activity centre relocation type. All point estimates were obtained using posterior modes and interval estimates were calculated through 95% highest posterior density (HPD) intervals at the 2.5% and 97.5% quantiles of the posterior distribution . We used the Gelman-Rubin statistic to assess convergence of open population SCR parameters through the “gelman.diag” function from the R package “coda” . This calculates the potential scale reduction factor for each parameter with upper confidence limits (CI), whereby parameters with a 95% upper CI substantially above 1 are considered to lack c onvergence . Closed population SCR models in a maximum likelihood framework were fitted using the R package “secr” . 9 44 We defined one sampling occasion as 5  days and used detector type proximity according to a Bernoulli distri- bution to estimate the combined density of males and females. To allow method comparison, we used for each session a continuous buffer of 18,000  km as for open population SCR models and created rectangular state-space masks accordingly (Supplementary Figs. S1–S10). First, we fitted a model M which assumes baseline detection probability g and detection function scale σ to be equal for all individuals . Secondly, we included sex as a covariate for both parameters and fitted two different models for each session, i.e. one keeping baseline detection probability g constant and testing sex on the detection function scale σ and one testing sex on both parameters for each session. We compared Akaike’s Information Criterion corrected for small sample sizes (AICc) to deter- mine the best model for our data and kept those with ΔAICc < 2 for model-averaging . Uknown sexes were considered as a latent covariate by performing a hybrid mixture m odel . This model consists of a combination of latent (e.g. missing sexes) and known classes (e.g. male and female) and estimates the mixing proportion of the different classes by calculation of the parameter “pmix”, which corresponds to the sex ratio. Additionally, this parameter allows class-specific modelling of the detection parameters for investigating sex die ff rences, for example, in the detection probability. Statistical significance was evaluated using the 95% CI. Reproductive parameters and age distribution. Based on camera trapping data, we investigated reproductive parameters such as generation time and number of recruits. The first was calculated as the mean age of resident reproducing females (i.e. females with kittens) at their first documented reproduction. Age was assessed only for individuals first photographed as kittens and recaptured over time. We, therefore, only included resident reproducing females of known age in the calculation. In contrast, no resident reproducing females were discarded when calculating the average number of recruits. This was calculated as the overall mean number of kittens with the mother at the onset of winter (November–January) . Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 4 Vol:.(1234567890) www.nature.com/scientificreports/ Status Known mortality causes Independent (German side only) Ee ff ctive Female (with Total Camera trapping Unique lynx Traffic Session kittens) Male Unknown sex (floaters) Juvenile Unknown trapping sites nights events accidents Illegal killing 2009 6 (4) 3 0 9 (–) 8 3 55 1728 23 0 0 2010 6 (4) 6 1 13 (4) 8 0 62 3813 52 0 0 2011 8 (3) 7 1 16 (7) 5 4 66 4124 82 0 0 2012 7 (3) 5 1 13 (5) 5 0 31 2621 48 0 1 2013 8 (6) 7 1 16 (0) 8 0 66 5885 86 0 1 2014 9 (6) 8 0 18 (2) 9 0 65 6151 94 1 0 2015 11 (2) 10 0 22 (6) 5 0 64 6022 135 1 1 2016 12 (5) 11 3 27 (6) 11 0 65 5964 144 2 0 2017 11 (4) 14 3 29 (10) 6 0 69 5957 121 2 0 2018 11 (7) 13 1 25 (4) 15 0 69 6412 130 1 0 Table 1. Summary of the lynx information obtained every primary period, including: the number of independent females and resident reproducing females (with kittens), males, individuals of unknown sex, total number of independent individuals including floaters (– indicates not available), juveniles, individuals of unknown status, camera trapping sites, effective trapping nights, unique lynx events (at most one detection per site in any secondary period) of independent individuals and documented mortality cases with known causes. e fir Th st three sessions had reduced number of effective trap nights because the monitoring started later in November in both the national parks in year 2009 and just in the SNP in years 2010 and 2011. For the session 2012, due to technical reasons, lynx camera trapping only took place on the German side of the study area. e n Th umber of juveniles referred to the entire lynx year. Finally, we also included all the independent individuals of known age and sex detected during the primary periods for investigating the age distribution of the population. Relative abundance index. When species cannot be individually identified, camera trapping data can 47 23 be used to calculate abundance indices . The relative abundance index (RAI) , which expresses the number of captures per monitoring day, is widely used in wildlife research . We matched the session length for RAI to that of the SCR analysis and considered lynx’s main prey species, i.e. roe deer and red deer, and red fox (Vulpes vulpes). To avoid temporal autocorrelation, we defined capture events with a threshold of 10  min and bootstrapped the RAI estimates by randomly selecting one camera trap per site at a time over 5000 replicates to reduce camera and site capture biases. This aspect of the study was conducted in the BFNP only. Use of experimental animals. No animals were used or handled in the study. Results Camera trap monitoring. We detected 65 unique independent individuals (25 males, 28 females and 12 individuals of unknown sex) in a total of 48,677 trap nights (Table 1). For the reasons mentioned in the meth- ods, the standardised session length of 100 days was not achieved in the first three primary periods (i.e. year 2009, 2010 and 2011) and in the fourth primary period (i.e. year 2012) camera traps were not active in the SNP, potentially resulting in estimates with higher uncertainty. The overall number of trapping nights was partially reduced during each season by technical problems such as camera trap failure, theft or snowfall. Considering the primary periods, only ten pictures (between 0 and 2% depending on the primary period) did not allow individual identification because of poor quality, while the remaining pictures were identified with certainty. We observed an increase in numbers of juveniles and independent individuals including males, females and indi- viduals of unknown sex, while the number of resident reproducing females (i.e. females with kittens) remained quite stable. We could determine the number of floaters in most years, as well as their origin. Around 33% of the detected floaters were born in the study area while around 67% were born outside or had unknown origin. Compared to the other categories of independent lynx, the number of floaters was the value that fluctuated the most (from 0 in session 2013 to 10 in session 2017, see Table 1). Apart from three and four individuals in the first and third sessions, respectively, the status (i.e. juvenile or independent) was determined. The number of unique lynx events (at most one detection per site in any secondary period) including independent individuals increased by around six times over the monitoring period. For the German side, we were able to confirm human- related mortality causes of some detected individuals in the vicinity of the study area, which showed an increase in traffic accidents and occasional illegal killing. In particular, three adults, two subadults and two juveniles were killed in traffic accidents in the BFNP, while two adults and one juvenile of unknown sex were killed illegally outside the national park in Germany. Spatial capture–recapture models. e co Th mbined density estimates for each primary period using open population SCR models increased from 0.69 (95% HPD intervals 0.50–1.10) to 1.33 (95% HPD inter- Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 5 Vol.:(0123456789) www.nature.com/scientificreports/ Combined Male Female Closed SCR 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 Figure 2. Posterior modes and 95% highest posterior density (HPD) intervals of open population spatial capture–recapture (SCR) models referring to combined, male and female lynx densities. 95% HPD intervals Parameter Sex Estimate Lower Upper M 0.01 0.01 0.02 Detection probability F 0.01 0.01 0.01 M 4.07 3.80 4.38 Detection function scale F 3.96 3.66 4.32 M 0.90 0.79 0.96 Survival probabilities F 0.82 0.72 0.91 M 0.09 0.04 0.14 Per capita recruitment rate F 0.12 0.07 0.20 Movement C 4.75 4.03 5.65 Table 2. Posterior modes and 95% highest posterior density (HPD) intervals of open population spatial capture–recapture (SCR) model referring to male (M) and female (F) lynx yearly baseline detection probability g , detection function scale σ (km), survival probabilities and per capita recruitment rate. Per capita recruitment rate indicates the number of individuals of each sex added per year per total abundance. The yearly movement (km) was calculated for sexes combined (C). vals 1.05–1.79) individuals/100  km while, as expected, sex-specific density estimates were on average higher for females (Fig. 2). However, the overlapping 95% HPD intervals suggest there was no statistically significant difference between sexes. Combined abundance ranged from 38 (95% HPD intervals 26–60) to 75 (95% HPD intervals 58–100) individuals, that of males from 17 (95% HPD intervals 9–31) to 36 (95% HPD intervals 21–50) and that of females from 21 (95% HPD intervals 10–38) to 37 (95% HPD intervals 26–58). Concerning closed population SCR models, estimates of density ranged from 1.09 (SE 0.35) to 2.36 (SE 0.79) individuals/100  km (Table 1; Supplementary Table S2) and fluctuated strongly compared to open population SCR models (Fig.  2). Abundance ranged from 47.94 (SE 11.03) to 121.25 (SE 42.45) (Supplementary Table S2). Density estimates from each method fell just within the respective 95% CI or HPD of the other method, apart from the year 2012, indi- cating estimates were comparable and suggesting no statistically significant difference across methods. Further results of closed population SCR models are shown in the appendix (Supplementary Tables S2, S3). All Gelman–Rubin diagnostic statistics had a 95% upper CI < 1.1 indicating convergence was reached. The parameter posterior modes (Table 2) of the yearly baseline detection probability g were equal across sexes. The yearly per capita recruitment rate was higher for females. The yearly survival probabilities of males, i.e. 0.90 (95% HPD intervals 0.79–0.96), were higher than those of females, i.e. 0.82 (95% HPD intervals 0.72–0.91). Males showed a higher yearly detection function scale σ, due to larger home ranges. However, considering the overlap of 95% HPD intervals, none of the parameters suggested a statistically significant difference across sexes. The realised population growth rate was 1.06 (range 1.01–1.10), that of males was 1.07 (range 1.00–1.12) and that of females was 1.05 (range 1.00–1.14), suggesting there was no predominant sex in the population. The probability of any individual to be a female was quite stable over years ranging from 0.51 (95% HPD intervals 0.36–0.65) to 0.56 (95% HPD intervals 0.33–0.76), meaning the sex ratio was slightly skewed towards females, with, however, no apparent statistically significant die ff rence. Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 6 Vol:.(1234567890) Individuals/100km2 0.00.5 1.01.5 2.0 2.5 3.0 3.5 www.nature.com/scientificreports/ Events Session Roe deer Red deer Red fox Camera trapping sites Ee ff ctive trapping nights 2009 12 12 136 29 1238 2010 28 56 306 31 2786 2011 29 87 291 31 2914 2012 52 52 315 31 2621 2013 57 57 315 31 2865 2014 57 86 401 30 2865 2015 110 110 274 29 2745 2016 57 172 58 29 2871 2017 27 189 189 29 2698 2018 27 191 273 29 2732 Table 3. Summary of the available information about roe deer, red deer and red fox obtained every primary period, including number of events, camera trapping sites and effective trapping nights. The number of events refers to the mean number of events resulting from all iterations in which one camera trap per site was sampled at a time. roe deer red deer red fox 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 Figure 3. Estimates and standard deviations of the relative abundance index (RAI) of red fox and lynx’s prey species roe deer and red deer calculated for the BFNP only. Reproductive parameters and age distribution. Regarding reproductive parameters, we included 14 resident reproducing females of known age and calculated a mean generation time of 2.64 years. One individual was first detected with kittens when it was seven years old and thus considered an outlier. The average number of recruits was 1.97 (range 1–3) and was calculated from information on 22 detected females (Supplementary Table S4). Age and sex could be confirmed for a total of 24 independent individuals (15 females and 9 males) and were included in the age distribution pyramid (Supplementary Fig. S11). The oldest individual of known age was a 10-year-old male individual. The number of females at each age was equal to or higher than that of males except for ages 8–10. Relative abundance index. e n Th umber of camera trapping sites used to calculate the RAIs ranged from 29 to 31 and effective trapping nights was 26,335 (Table  3). The RAI of red fox ranged from 0.07 (SD 0.01) to 0.14 (SD 0.02), that of red deer from 0.01 (SD 0.00) to 0.07 (SD 0.01), that of roe deer from 0.01 (SD 0.00) to 0.04 (SD 0.01) (Fig. 3). Fox abundance increased from 2009 to 2014 and decreased in the remaining years, overall showing an oscillating trend. Red deer and roe deer slightly increased until 2015, aer w ft hich the former almost doubled while the latter decreased to earlier values. Discussion Our modelling approach revealed that the number of independent lynx detected in the Bohemian–Bavarian Forest has increased over the past decade, with concurrently increasing abundance, density and positive popu- lation growth rates. Survival probabilities were high and per capita recruitment rate was low indicating a low yearly population turnover. Reproductive parameters such as generation time and average litter size consider- ing all resident reproducing females indicated successful reproduction occurred every year, and the number of Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 7 Vol.:(0123456789) Relative abundance index (RAI) 0.00 0.05 0.10 0.15 0.20 www.nature.com/scientificreports/ reproducing females remained stable throughout the study period despite the increase in the total number of detected independents. All these findings indicate that the protected areas act as an important, stable source area for the Bohemian–Bavarian–Austrian lynx population. Comparison of spatial capture–recapture methods. To our knowledge, this is the first demographic study on lynx comprising annual camera trapping over a decade and using open population SCR models. These models can provide a wide range of demographic parameters useful for long-term lynx monitoring. In addition, estimates from open population SCR models were not prone to inter-annual fluctuations that ae ff cted closed population SCR models and were generally lower, due to lower abundances compared to those from the lat- ter (Table 1; Supplementary Table S2) divided by an equal buffer area. However, higher estimates with greater uncertainty, concerning the years 2009 and 2012, for which data were least complete, resulted from issues when accounting for varying efforts, i.e. the exact number of days the cameras were working. Allowing for varying effort has the benefit that detection parameters are related to a specific unit of effort thus being unbiased . In the open population SCR models, the activity of each camera trap was informed over the primary periods. Con- cerning closed population models, in the year 2009, it was not possible to include the period from September to November because no cameras were functioning, resulting in zero occasions. While in the year 2012, the activity of the cameras in the SNP could not be specified because the devices were not in the field in the yearly camera trapping array that year. Therefore, data gaps resulted in biased estimates in this case because some spatio-tem- poral requirements, namely adequate sample size and number of recaptures, were not fulfilled. Closed popula- tion SCR density estimates of primary periods for which it was possible to account for varying effort, including years 2010 and 2011 for which data were also not complete, were comparable to those from open population SCR models since their 95% CI and HPD intervals overlapped. One advantageous feature of open population SCR methods in a Bayesian framework is their ability to deal with incomplete detection (e.g. non-annual moni- toring). This is possible because they produce posterior distributions of the demographic parameters that incor - porate the uncertainty resulting from data gaps by using information derived from other primary periods . This highlights the robustness of this method, which we recommend for future studies and monitoring, especially when incomplete detections occur . Nonetheless, closed population SCR models still represent a reliable and conventional tool for abundance and density assessments . Spatial capture–recapture models. Our density estimates of both open (0.69–1.33 individuals/100  km ) and closed population SCR models (1.02–2.39 individuals/100  km ) were on the same order of magnitude as those resulting from other closed population SCR studies conducted on both reintroduced and autochthone populations across Europe. For example, the reintroduced populations of the French Jura and Vosges Mountains, 2 35 with the lowest density estimate reported in Europe (0.24–0.91 individuals/100  km ) , or the Swiss Alps (1.47 2 36 and 1.38 individuals/100  km ) , and the autochthone population of the Western Carpathians (0.26–1.85 indi- 2 50 viduals/100  km of suitable habitat) . The only outlier is a subpopulation in southwest Asia (Turkey) with a high lynx density of 4.20 individuals/100  km . However, lynx in Turkey live under different ecological conditions compared to those in Western and Central Europe, feeding mainly on lagomorphs instead of ungulates which leads to smaller home range sizes . In our study, the yearly baseline detection probability g was equal between sexes (Table 2). Although males are generally more active in order to patrol their larger home ranges, females with kittens hunt at a higher rate, which can result in augmented ac tivity , and might explain similarities across sexes for this parameter. The yearly detection function scale σ was not significantly higher for males compared to females (Table  2). This sug- 20 10 gests similar home range sizes, which stands in contrast to expectations for l ynx and most felids . However, the restricted and seasonal session period was chosen for demographic closure and not appropriate for annual home range estimation Furthermore, many detected individuals were not residents thus resulting in potentially biased estimates for the parameter in question. e Th combined yearly survival probabilities for independent individuals reached 85%. This high value is likely because our estimates come from strictly protected areas where lynx have a higher chance of survival. However, none of the detected individuals used the territory included in the protected areas exclusively, thus potentially exposing themselves to a higher risk. The survival estimated through open population SCR methods in our study site is among the highest reported for lynx. In the Western Carpathians, Dula et al. conducted an SCR survey on an autochthone lynx population across multiple seasons and found an overall apparent survival (i.e. including emigration) of 63% with high human-related mortality occurring in the area. The only other values of survival rates available for comparisons come from telemetry studies because almost all of the SCR surveys conducted on lynx in Europe based on camera trapping data used closed population models and thus did not provide information on this parameter. In Switzerland, Breitenmoser-Würsten et al. found an overall survival rate of 76% for adults and 53.3% for subadults. In Poland, the survival rate was only 63% , considering subadults and adults combined. In three Scandinavian study sites, Andrén et al. found lower survival rates for subadults (70, 77 and 71%) than adults (87, 91 and 84%), likely due to lower mortality related to vehicle collision and hunting in adulthood. Comparisons to SCR results are however difficult given the contrasting underlying meth- odologies. Against our expectations, the survival probabilities were higher for males, though not significantly. This was due to a higher overall number of apparent survival events for males, which consist of the number of consecutive detections over years including gaps during which the animal was alive but not detected. A similar open population SCR study, conducted on a low-density ocelot population in Belize, also showed no significant differences in sex-specific survival probabilities, despite they found probabilities of 0.86 for females and 0.78 for males, respectively . e a Th uthors suggested the statistical power was not enough to detect significant differences between sexes in these parameters, although they were able to determine sex for a large number of adult ocelots Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 8 Vol:.(1234567890) www.nature.com/scientificreports/ (n = 322). These are probably also the reasons why we could not prove statistically significant die ff rences in per capita recruitment rate between sexes (Table 2). Camera trapping does not allow assessment of the fates of all disappearing individuals and thus does not provide information on natural mortality. However, our auxiliary findings on dead animals show an increasing number of traffic accidents in recent years (Table  1), indicating that the prevention of lynx vehicle collisions represents an important management action for the population. We could not assess the actual impact of illegal killing since carcasses are seldom found, but considering the high survival of the lynx individuals monitored in the study area, illegal killing seems to play no important role within the study site. However, a high poaching rate in the surroundings of our study area was suggested by a modelling approach, which underlines the impor- tance of the protected areas for the survival of the p opulation . These results are consistent with other areas in 5 54 Europe e.g. S witzerland , and P oland were traffic accidents and poaching have been found as the main cause of lynx mortality. Concerning limitations of the open population SCR models, the movement parameter had poor mixing, which is in line with Gardner et al. , meaning Markov chain Monte Carlo iterations slowly converged to unbiased posterior distributions. This is the reason why the parameter was estimated for both sexes combined and it was necessary to run relatively long chains (100,000 iterations). We attempted to separate survival from emigration and recruitment from immigration using open popula- tion SCR models with a spatially explicit model for activity centre relocation . However, the discrimination ability of these models depends on how accurately the movement model can describe the actual activity centre relocation between primary periods . This would require a larger camera trapping array to detect dispersing movements of large carnivores. We, therefore, could not reliably separate emigration and immigration from survival and recruitment, respectively. However, it is reasonable to assume that true survival and recruitment are at least as high as the value we estimated when the underestimation due to the inclusion of emigration and immigration is considered. Reproductive parameters and age distribution. Generation time was calculated as the mean age of resident reproducing females at their first documented reproduction, potentially resulting in overestimated val- ues if previous litters went undetected by camera traps. This may be the case for individuals living close to the boundaries of the monitored area, such as the outlier we detected. In Scandinavia, generation time of lynx was investigated across different environments and ranged from two to more than three years with the highest values in the northern territories due to slower life cycles and a different feeding regime and thus body mass . Similar results were found in other areas of S candinavia . We were not able to fully assess the proportion of reproduc- tive females because sex was not determined for all detected individuals. However, our mean generation time of 2.64 years fall in the range of the Scandinavian studies. We found a rare case of a reproducing female first reproducing at one year of age (Supplementary Table S3). The first reported case of a 1-year-old female Eurasian lynx breeding in the wild was recently observed in the Bohemian–Bavarian–Austrian population range, outside the study area . The authors hypothesised this could be related to very high prey density or high turnover due to poaching. The average number of recruits of 1.97 (range 1–3) was comparable to that reported in Scandinavia. Gaillard et al. , investigated the number of recruits related to litter size and found higher values for multiparous females (i.e. females with multiple reproduction events) ranging approximately between 0.5 and 2 with a litter size between one and four kittens. However, we were not able to determine litter size through camera trapping. Regarding age distribution, the maximum age determined with certainty in the study area was 10 years, which was constrained by the duration of the study. In Switzerland, Breitenmoser-Würsten et al. found individuals 14–15 years of age through telemetry. The exact age can only be determined by camera trapping if an individual was first photographed as a kitten. As such, some of the individuals already classified as independent in the first monitoring session reached higher ages than those reported, even if their exact age remained unknown. Relative abundance index. According to the RAI results, roe deer abundance in the German part of the study site remained stable even though they were not hunted since 2012. Contrastingly, the red deer population strongly increased (Fig. 3). As lynx in the study area mainly feed on roe deer and are normally only able to kill red deer calves and subadults, they have a stronger effect on the abundance of roe deer than red deer . The stable RAI for roe deer might therefore indicate a limiting effect of lynx, a result consistent with observations from the 58 59 same and other a reas . Concerning red fox, intraguild predation of lynx on this species has already been reported in die ff rent areas 60–62 of Scandinavia . The fox RAI shows a decrease with increasing lynx numbers, but the trend is not clear. In the study area, the red fox was found as only 1% of lynx k ills meaning other factors such as food availability and diseases could be driving red fox dynamics. RAI results were not related statistically to SCR models as they could not be used as an individual covariate or as a trap covariate since this aspect of the study was limited to the BFNP only. Therefore the RAI results should be interpreted cautiously. Although RAI can be ae ff cted by biases attributed to changes in detection p robability , we believe this method well represented the abundance trends of the species in question since the study design was uniform over the years. Population development and prospects. For lynx, we found increasing abundance and density, as well as positive population growth rates, high survival probabilities and low per capita recruitment, resulting in a low yearly population turnover. This suggests that lynx is not strongly ae ff cted by human-related mortality within the study area. This contrasts with Duľa et al. , who found fluctuating density estimates in the Western Carpathians due to illegal killing and traffic collisions. The comparatively stable development in our study area was accompanied by a relatively constant number of family groups (Table 1), i.e. females accompanied by their Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 9 Vol.:(0123456789) www.nature.com/scientificreports/ kittens, which form the most stable structure of the population and are therefore of particular importance. On the other hand, the number of individuals belonging to the least stable part of the population, i.e. floaters likely dispersing through the study area, highly fluctuated throughout the study period (from a minimum of 0 in 2013 to a maximum of 10 in 2017). Furthermore, the increase (or decrease) of the total number of detected individuals from one monitoring season to the following mainly corresponded with an increase (or decrease, respectively) in the number of detected floaters (Table  1). This suggests that the number of floaters visiting the study area in the given year played a major role in determining local abundance lynx for that year. Among all floaters detected in this study, only a lower proportion of individuals originated inside the study area (~ 33%). This reveals potential changes in the Bohemian–Bavarian–Austrian lynx population. Based on coordinated transboundary camera trapping within an area of 13,000  km , the estimated entire population size showed a positive trend and confirmed lynx presence and regular reproduction in a large part of the entire population range, including its outskirts . Thus, fluctuations in the number of floaters may be related to what takes place outside of the protected study area. In particular, increases in the number of floaters reaching the study area from 1 year to the following may be related to a local increase in survival of juveniles and yearlings in part of the Bohemian–Bavarian–Austrian population range. This may be partially related to a local decrease in poaching, possibly as a result of law enforcement and long-term public relations campaigns aimed at improving accept- ance in the wake of high-profile poaching incidents. These findings suggest that the overall conditions of the Bohemian–Bavarian–Austrian population have at least partially improved in recent years and this do not only concern the protected areas lying within this population’s range. This is fundamental for the long term genetic viability of a potential future Central European metapopulation as the conservation of species with large spatial requirements cannot rely on protected areas a lone . Since data on poaching are almost absent, the causes lead- ing to this positive trend of the population remain speculation. However, other potential explanations such as immigration from other populations or variability in reproduction rates within our study area can be excluded since the Bohemian–Bavarian–Austrian is an isolated population . Conclusions High survival probabilities and regular reproduction thanks to local protection indicate that the protected areas act as a stable source of lynx for the wider distribution range of the Bohemian–Bavarian–Austrian population, which is in line with previous studies . Thus our results show the high value of the protected areas for the persistence of the Bohemian–Bavarian–Austrian As proposed for other large carnivores , we conclude that protected areas can be considered as strategic areas for lynx recovery across its range because full protection increases survival probabilities and ensure stable reproduction. However, we must consider that the management of the surrounding landscapes is fundamental for the long term survival of the population, as protected areas in Central Europe are too small to encompass enough animals to maintain healthy populations over the long term. Our study revealed crucial demographic parameters of a lynx population that can be used for improving con- servation and management plans. Finally, our results stress the importance of long-term systematic monitoring as a basis for the understanding of the population dynamics of large carnivore populations and recommend the use of open population SCR methods to achieve such aims. Data availability Data supporting the findings of this study are not included due to conservation concerns. Data may, however, be available from the authors upon reasonable request and with permission of both Administrations of both the BFNP and SNP. Received: 16 February 2021; Accepted: 22 September 2021 References 1. Ripple, W. J. et al. Status and ecological ee ff cts of the world’s largest carnivores. Science 343, 1241484–1241484 (2014). 2. Treves, A. & Karanth, K. U. Human–carnivore conflict and perspectives on carnivore management worldwide. Conserv. Biol. 17, 1491–1499 (2003). 3. Linnell, J. D. C. & Boitani, L. Building biological realism into wolf management policy: The development of the population approach in Europe. Hystrix Ital. J. Mammal. 23, 80–91 (2011). 4. Heurich, M. et al. Illegal hunting as a major driver of the source-sink dynamics of a reintroduced lynx population in Central Europe. Biol. Conserv. 224, 355–365 (2018). 5. Breitenmoser-Würsten, C., Vandel, J.-M., Zimmermann, F. & Breitenmoser, U. Demography of lynx Lynx lynx in the Jura Moun- tains. Wildl. Biol. 13, 381–392 (2007). 6. Clutton-Brock, T. & Sheldon, B. C. Individuals and populations: The role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends Ecol. Evol. 25, 562–573 (2010). 7. O’Connell, A., Nichols, J. D. & Karanth, K. U. Camera Traps in Animal Ecology: Methods and Analyses. (Springer Tokyo, 2011). 8. Noss, A. J. et al. A Camera trapping and radio telemetry study of lowland tapir (Tapirus terrestris) in Bolivian Dry Forests. Tapir Cons. 12, 9 (2003). 9. Karanth, K. U. & Nichols, J. D. Estimation of tiger densities in India using photographic captures and recaptures. Ecology 79, 11 (1998). 10. Satter, C. B., Augustine, B. C., Harmsen, B. J., Foster, R. J. & Kelly, M. J. Sex‐specific population dynamics of ocelots in Belize using open population spatial capture–recapture. Ecosphere 10, e02792 (2019). 11. Silver, S. C. et al. The use of camera traps for estimating jaguar Panthera onca abundance and density using capture/recapture analysis. Oryx 38, 148–154 (2004). 12. Zimmermann, F., Breitenmoser-Würsten, C., Molinari-Jobin, A. & Breitenmoser, U. Optimizing the size of the area surveyed for monitoring a Eurasian lynx (Lynx lynx) population in the Swiss Alps by means of photographic capture–recapture. Integr. Zool. 8, 232–243 (2013). 13. Royle, J. A., Chandler, R. B., Sollmann, R. & Gardner, B. Spatial Capture–Recapture. (Elsevier, 2014). Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 10 Vol:.(1234567890) www.nature.com/scientificreports/ 14. Chandler, R. B. & Clark, J. D. Spatially explicit integrated population models. Methods Ecol. Evol. 5, 1351–1360 (2014). 15. Kaczensky, P . et al. Status, management and distribution of large carnivores—Bear, lynx, wolf and wolverine in Europe (Europe- anCommission, 2013). 16. Magg, N. et al. Habitat availability is not limiting the distribution of the Bohemian–Bavarian lynx Lynx lynx population. Oryx 50, 742–752 (2016). 17. Müller, J. et al. Protected areas shape the spatial distribution of a European lynx population more than 20 years ae ft r reintroduc - tion. Biol. Conserv. 177, 210–217 (2014). 18. Bull, J. K. et al. The effect of reintroductions on the genetic variability in Eurasian lynx populations: The cases of Bohemian–Bavar - ian and Vosges–Palatinian populations. Conserv. Genet. 17, 1229–1234 (2016). 19. Walston, J. et al. Bringing the tiger back from the brink—The six percent solution. PLoS Biol. 8, e1000485 (2010). 20. Schmidt, K., Jędrzejewski, W. & Okarma, H. Spatial organization and social relations in the Eurasian lynx population in Bialowieza Primeval Forest, Poland. Acta Theriol. (Warsz.) 42, 289–312 (1997). 21. Bunnefeld, N., Linnell, J. D. C., Odden, J., van Duijn, M. A. J. & Andersen, R. Risk taking by Eurasian lynx (Lynx lynx) in a human- dominated landscape: Ee ff cts of sex and reproductive status. J. Zool. 270, 31–39 (2006). 22. Gaillard, J.-M., Nilsen, E. B., Odden, J., Andrén, H. & Linnell, J. D. C. One size fits all: Eurasian lynx females share a common optimal litter size. J. Anim. Ecol. 83, 107–115 (2014). 23. Nilsen, E. B., Linnell, J. D. C., Odden, J., Samelius, G. & Andrén, H. Patterns of variation in reproductive parameters in Eurasian lynx (Lynx lynx). Acta Theriol. (Warsz.) 57, 217–223 (2012). 24. O’Brien, T. G., Kinnaird, M. F. & Wibisono, H. T. Crouching tigers, hidden prey: Sumatran tiger and prey populations in a tropical forest landscape. Anim. Conserv. 6, 131–139 (2003). 25. Cailleret, M., Heurich, M. & Bugmann, H. Reduction in browsing intensity may not compensate climate change effects on tree species composition in the Bavarian Forest National Park. For. Ecol. Manag. 328, 179–192 (2014). 26. Heurich, M. et al. Country, cover or protection: What shapes the distribution of red deer and roe deer in the Bohemian Forest Ecosystem?. PLoS ONE 10, e0120960 (2015). 27. van Beeck Calkoen, S. T. S. et al. The blame game: Using eDNA to identify species-specific tree browsing by red deer (Cervus elaphus) and roe deer (Capreolus capreolus) in a temperate forest. For. Ecol. Manag. 451, 117483 (2019). 28. Wöl, fl M. et al. Distribution and status of lynx in the border region between Czech Republic, Germany and Austria. Acta e Th riol. 46, 181–194 (2001). 29. Mináriková, T. et al. Lynx monitoring report for Bohemian–Bavarian–Austrian lynx population for lynx year 2017 (INTERREG Central Europe, 2019). 30. Weingarth, K. et al. First estimation of Eurasian lynx (Lynx lynx) abundance and density using digital cameras and capture–recap- ture techniques in a German national park. Anim. Biodivers. Conserv. 35, 197–207 (2012). 31. Belotti, E. et al. Patterns of lynx predation at the interface between protected areas and multi-use landscapes in Central Europe. PLoS ONE 10, e0138139 (2015). 32. Tobler, M. W. & Powell, G. V. N. Estimating jaguar densities with camera traps: Problems with current designs and recommenda- tions for future studies. Biol. Conserv. 159, 109–118 (2013). 33. Zimmermann, F., Breitenmoser-Würsten, C. & Breitenmoser, U. Natal dispersal of Eurasian lynx ( Lynx lynx ) in Switzerland. J. Zool. 267, 381 (2005). 34. Andrén, H. et al. Survival rates and causes of mortality in Eurasian lynx (Lynx lynx) in multi-use landscapes. Biol. Conserv. 131, 23–32 (2006). 35. Gimenez, O. et al. Spatial density estimates of Eurasian lynx (Lynx lynx) in the French Jura and Vosges Mountains. Ecol. Evol. 9, 11707–11715 (2019). 36. Pesenti, E. & Zimmermann, F. Density estimations of the Eurasian lynx (Lynx lynx) in the Swiss Alps. J. Mammal. 94, 73–81 (2013). 37. Weingarth, K. et al. Hide and seek: Extended camera-trap session lengths and autumn provide best parameters for estimating lynx densities in mountainous areas. Biodivers. Conserv. 24, 2935–2952 (2015). 38. Pollock, K. H. A capture–recapture design robust to unequal probability of capture. J. Wildl. Manag. 46, 752 (1982). 39. Augustine, B. benaug/OpenPopSCR. (2019). https:// github. com/ benaug/ OpenP opSCR. 40. Ergon, T. & Gardner, B. Separating mortality and emigration: Modelling space use, dispersal and survival with robust-design spatial capture–recapture data. Methods Ecol. Evol. 5, 1327–1336 (2014). 41. Schaub, M. & Royle, J. A. Estimating true instead of apparent survival using spatial Cormack–Jolly–Seber models. Methods Ecol. Evol. 5, 1316–1326 (2014). 42. Plummer, M., Best, N., Cowles, K. & Vines, K. CODA: Convergence diagnosis and output analysis for MCMC. R News 6, 7–11 (2005). 43. Brooks, S. P. & Gelman, A. General methods for monitoring convergence of iterative simulations. J. Comput. Graph. Stat. 7, 434–455 (1998). 44. Efford, M. secr 4.1—Spatially explicit capture–recapture in R. (2019). https://cra n.micr os oft.co m/sn apsh ot/2019- 12- 24/ w eb/p acka ges/ secr/ vigne ttes/ secr- overv iew. pdf. 45. Burnham, K. P. & Overton, W. S. Robust estimation of population size when capture probabilities vary among animals. Ecology 60, 927–936 (1979). 46. Burnham, K. P., Anderson, D. R. & Burnham, K. P. Model Selection and Multimodel Inference: A Practical Information–e Th oretic Approach (Springer, 2002). 47. O’Brien, T. G. Abundance, density and relative abundance: A conceptual framework. In Camera Traps in Animal Ecology (eds O’Connell, A. F. et al.) 71–96 (Springer Japan, 2011). https:// doi. org/ 10. 1007/ 978-4- 431- 99495-4_6. 48. Rovero, F. & Zimmermann, F. Camera Trapping for Wildlife Research (Pelagic Publishing Ltd, 2016). 49. Augustine, B. C. et al. Sex-specific population dynamics and demography of capercaillie (Tetrao urogallus L.) in a patchy environ- ment. Popul. Ecol. 62, 80–90 (2020). 50. Duľa, M. et al. Multi-seasonal systematic camera-trapping reveals fluctuating densities and high turnover rates of Carpathian lynx on the western edge of its native range. Sci. Rep. 11, 9236 (2021). 51. Avgan, B., Zimmermann, F., Güntert, M., Arıkan, F. & Breitenmoser, U. The first density estimation of an isolated Eurasian lynx population in southwest Asia. Wildl. Biol. 20, 217–221 (2014). 52. Mengüllüoğlu, D., Ambarlı, H., Berger, A. & Hofer, H. Foraging ecology of Eurasian lynx populations in southwest Asia: Conserva- tion implications for a diet specialist. Ecol. Evol. 8, 9451–9463 (2018). 53. Heurich, M. et al. Activity patterns of Eurasian lynx are modulated by light regime and individual traits over a wide latitudinal range. PLoS ONE 9, e114143 (2014). 54. Jedrzejewski, W . et al. Population dynamics (1869–1994), demography, and home ranges of the lynx in Bialowieza Primeval Forest (Poland and Belarus). Ecography 19, 122–138 (1996). 55. Gardner, B., Sollmann, R., Kumar, N. S., Jathanna, D. & Karanth, K. U. State space and movement specification in open population spatial capture–recapture models. Ecol. Evol. 8, 10336–10344 (2018). 56. López-Bao, J. V. et al. Eurasian lynx t fi ness shows little variation across Scandinavian human-dominated landscapes. Sci. Rep. 9, 8903 (2019). 57. Engleder, T. et al. First breeding record of a 1-year-old female Eurasian lynx. Eur. J. Wildl. Res. 65, 17 (2019). Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 11 Vol.:(0123456789) www.nature.com/scientificreports/ 58. Heurich, M. et al. Selective predation of a stalking predator on ungulate prey. PLoS ONE 11, e0158449 (2016). 59. Andrén, H. & Liberg, O. Large impact of Eurasian lynx predation on roe deer population dynamics. PLoS ONE 10, e0120570 (2015). 60. Elmhagen, B. & Rushton, S. P. Trophic control of mesopredators in terrestrial ecosystems: Top-down or bottom-up?. Ecol. Lett. 10, 197–206 (2007). 61. Wikenros, C. et al. Fear or food—Abundance of red fox in relation to occurrence of lynx and wolf. Sci. Rep. 7, 9059 (2017). 62. Helldin, J. O., Liberg, O. & Glöersen, G. Lynx (Lynx lynx) killing red foxes (Vulpes vulpes) in boreal Sweden? Frequency and population ee ff cts. J. Zool. 270, 657–663 (2006). 63. Sollmann, R., Mohamed, A., Samejima, H. & Wilting, A. Risky business or simple solution—Relative abundance indices from camera-trapping. Biol. Conserv. 159, 405–412 (2013). 64. Linnell, J. D. C., Kaczensky, P., Wotschikowsky, U., Lescureux, N. & Boitani, L. Framing the relationship between people and nature in the context of European conservation: Relationship between people and nature. Conserv. Biol. 29, 978–985 (2015). Acknowledgements The study was funded by the German Academic Exchange Service (DAAD) and allowed by the Administra- tions of both the BFNP and SNP. Financial support was also provided by EU Program Interreg IV (Objective 3 Czech Republic—the Independent State of Bavaria, Grant no. 18) and by the World Wide Fund for Nature Section Germany (Grant no. 533500/53000/10120). We greatly thank Ben Augustine who kindly shared open population SCR models in the R package OpenPopSCR and supported the following analysis as well as helped interpret results. We sincerely thank all the colleagues of both the BFNP and the SNP who facilitated the study. In particular, we thank Nicolò Borgianni who helped with part of the analyses. Author contributions This submission represents original work carried out by the authors. All authors have contributed substantially to the manuscript and have approved the manuscript and agree to its submission. Conceptualisation: S.P., J.P., M.H., methodology: S.P., C.H., J.P., M.H., data collection: M.G., K.W., E.B., L.B., formal analysis and investigation: S.P., C.H., writing—original draft preparation: S.P., J.P., writing—review and editing: J.P., E.B., M.H., funding acquisition: S.P, E.B., L.B., M.H. Funding Open Access funding enabled and organized by Projekt DEAL. Competing interests The authors declare no competing interests. Additional information Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598- 021- 99337-2. Correspondence and requests for materials should be addressed to S.P. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. © The Author(s) 2021 Scientific Reports | (2021) 11:19868 | https://doi.org/10.1038/s41598-021-99337-2 12 Vol:.(1234567890)

Journal

Scientific ReportsSpringer Journals

Published: Oct 6, 2021

There are no references for this article.