# Spatial density estimates of Eurasian lynx (Lynx lynx) in the French Jura and Vosges Mountains

Olivier Gimenez, Sylvain Gatti, Christophe Duchamp, Estelle Germain, Alain

Laurent, Fridolin Zimmermann, Eric Marboutin

► To cite this version:

Olivier Gimenez, Sylvain Gatti, Christophe Duchamp, Estelle Germain, Alain Laurent, et al.. Spatial density estimates of Eurasian lynx (Lynx lynx) in the French Jura and Vosges Mountains. Ecology and Evolution, 2019, 10.1002/ece3.5668. hal-02329835

HAL Id: hal-02329835

https://hal.science/hal-02329835

Submitted on 4 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ORIGINAL RESEARCH

Ecology and Evolution
Open Access

WILEY

# Spatial density estimates of Eurasian lynx (Lynx lynx) in the French Jura and Vosges Mountains

Olivier Gimenez<fnref n="1" /> ID | Sylvain Gatti<fnref n="2" /> | Christophe Duchamp<fnref n="3" /> | Estelle Germain<fnref n="4" /> |
Alain Laurent<fnref n="2" /> | Fridolin Zimmermann<fnref n="5" /> ID | Eric Marboutin<fnref n="2" />

## Correspondence

Olivier Gimenez, CEFE, CNRS, EPHE, IRD, Univ Montpellier, Univ Paul Valéry Montpellier 3, Montpellier, France.
Email: olivier.gimenez@cefe.cnrs.fr

## Funding information

Centre National de la Recherche Scientifique; European Union within the framework of the Operational Program FEDER-FSE "Lorraine et Massif des Vosges 2014-2020": Direction Régionale de l'Environnement, l'Aménagement et le Logement; Région Grand Est; Zoo d'Amnéville; Fondation d'entreprise UEM; Fondation Nature & Découvertes; Fondation Le Pal Nature; Chasseur d'images magazine; ONCFS; CIL&B; MTES (ITTECOP); FRB; "Mission pour l'Interdisciplinarité" through the "Osez l'Interdisciplinarité" initiative; "Commissariat à l'Aménagement du Massif des Vosges" for the FNADT ("Fonds National d'Aménagement et de Développement du Territoire")

## Abstract

Obtaining estimates of animal population density is a key step in providing sound conservation and management strategies for wildlife. For many large carnivores however, estimating density is difficult because these species are elusive and wide-ranging. Here, we focus on providing the first density estimates of the Eurasian lynx (Lynx lynx) in the French Jura and Vosges mountains. We sampled a total of 413 camera trapping sites (with two cameras per site) between January 2011 and April 2016 in seven study areas across seven counties of the French Jura and Vosges mountains. We obtained 592 lynx detections over 19,035 trap days in the Jura mountains and 0 detection over 6,804 trap days in the Vosges mountains. Based on coat patterns, we identified a total number of 92 unique individuals from photographs, including 16 females, 13 males, and 63 individuals of unknown sex. Using spatial capture-recapture (SCR) models, we estimated abundance in the study areas between 5 (SE = 0.1) and 29 (0.2) lynx and density between 0.24 (SE = 0.02) and 0.91 (SE = 0.03) lynx per 100 km<fnref n="2" />. We also provide a comparison with nonspatial density estimates and discuss the observed discrepancies. Our study is yet another example of the advantage of combining SCR methods and noninvasive sampling techniques to estimate density for elusive and wide-ranging species, like large carnivores. While the estimated densities in the French Jura mountains are comparable to other lynx populations in Europe, the fact that we detected no lynx in the Vosges mountains is alarming. Connectivity should be encouraged between the French Jura mountains, the Vosges mountains, and the Palatinate Forest in Germany where a reintroduction program is currently ongoing. Our density estimates will help in setting a baseline conservation status for the lynx population in France.

## KEYWORDS

camera trapping, large carnivores, noninvasive sampling, photo identification, spatially explicit capture-recapture models

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2019 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.

<footnote n="1">¹CEFE, CNRS, EPHE, IRD, Univ Montpellier, Univ Paul Valéry Montpellier 3, Montpellier, France</footnote>

<footnote n="2">²Office National de la Chasse et de la Faune Sauvage, Gières, France</footnote>

<footnote n="3">³Office National de la Chasse et de la Faune Sauvage, Gap, France</footnote>

<footnote n="4">⁴Centre de Recherche et d'Observation sur les Carnivores (CROC), Lucy, France</footnote>

<footnote n="5">⁵KORA, Muri, Switzerland</footnote>
# 1 | INTRODUCTION

Obtaining estimates of animal population density is a key step in providing sound conservation and management strategies for wildlife (Williams, Nichols, & Conroy, 2002). For many large carnivores however, estimating density is difficult because these species are elusive and wide-ranging, resulting in low detection rates (Obbard, Howe, & Kyle, 2010). To deal with these issues, noninvasive techniques, such as camera trapping and DNA sampling, are increasingly used (Kelly, Betsch, Wultsch, Mesa, & Mills, 2012). These noninvasive techniques generate data that can be analyzed with capture-recapture methods to estimate densities (Royle, Chandler, Sollmann, & Gardner, 2014).

Standard capture-recapture models for closed populations (Otis, Burnham, White, & Anderson, 1978) have long been used to estimate animal abundance and density, including many large carnivores (Gerber, Ivan, & Burnham, 2014; Mumma, Zieminski, Fuller, Mahoney, & Waits, 2015). However, when converting abundance into density, density estimates are highly sensitive to the size of user-defined area assumed to reflect the effective sampling area (White, Anderson, Burnham, & Otis, 1982). In addition, individual heterogeneity in detection due to spatial variation in the distance of home ranges to the sampling devices may lead to biased density estimates (Otis et al., 1978). Spatial capture-recapture (SCR) models deal with these issues by explicitly incorporating spatial locations of detections (Borchers, 2012; Borchers & Efford, 2008; Efford, 2004; Royle & Young, 2008), and they are increasingly used to estimate densities of large carnivores (Alexander, Gopalaswamy, Shi, & Face, 2015; Broekhuis & Gopalaswamy, 2016; Goldberg et al., 2015; López-Bao et al., 2018; Pesenti & Zimmermann, 2013; Stetz, Mitchell, & Kendall, 2018).

Here, we focus on the threatened Eurasian lynx (Lynx lynx) in the French Jura and Vosges mountains (see Chapron et al., 2014 for a map of its distribution in Europe; see also https://www.lcie.org/Large-carnivores/Eurasian-lynx for recent updates). As in many regions of western Europe (Breitenmoser, 1998), lynx was extirpated from France between the 17th and 20th centuries due to habitat degradation, persecution by humans and decrease in prey availability (Vandel & Stahl, 2005). Shortly, after their initial reintroduction in Switzerland in the 1970s (Breitenmoser, Breitenmoser-Würsten, & Capt, 1998), lynx naturally increased their range and started recolonizing France by repopulating forests on the French side of the Jura (Vandel & Stahl, 2005). Reintroductions also occurred in the French Vosges mountains between 1983 and 1993 with the perspective of establishing a population there (Vandel, Stahl, Herrenschmidt, & Marboutin, 2006). The species is listed as endangered in the IUCN Red list and is of conservation concern in France due to habitat fragmentation, poaching, and collisions with cars and trains. Currently, the French population of lynx is restricted to three mountain ranges: the Vosges in northeastern France, the Jura, and the Alps, with little connectivity between them most likely due to human-made linear infrastructures. While the Northern Alps are slowly being recoloned with lynx mostly coming from the Jura (Marboutin et al., 2012), the Jura holds the bulk of the French lynx population. In contrast, the lynx presence in the Vosges mountains remained stable following

the reintroductions and has been continuously decreasing since 2005 (Laurent et al., 2012).

Despite their conservation status, little information on abundance and density of lynx in France exist. In this study, we used SCR and standard capture-recapture models to provide the first estimate of lynx abundance and density using camera trap surveys implemented in the French Jura and Vosges mountains from 2011 to 2016. Based on these results, we discuss research and management priorities for the effective conservation of lynx in France.

# 2 | METHODS

## 2.1 | Ethics statement

We used noninvasive methods for data collection, which did not involve manipulation or handling of any living organism. Therefore, approval from an animal ethics committee was not required. Cameras were set on public or private forests with the permission of local authorities or local owners, respectively. We advertised the study and the presence of camera traps to the local stakeholders and the public visiting the areas. In agreement with French legislation, we deleted photographs permitting the identification of goods or people.

## 2.2 | Study area and sampling design

The study area encompassed three counties of the French Jura mountains, namely Ain, Doubs and Jura, and four counties of the Vosges mountains, namely Vosges, Haut-Rhin, Bas-Rhin, and Moselle (Figure 1). Elevation ranged from 163 to 1,718 m above sea level in the Jura mountains, and from 104 to 1,422 m in the Vosges mountains. The human population density was 88 per km² in the Jura mountains and 170 per km² in the Vosges mountains. The Jura mountains were 50% forest on average (Breitenmoser et al., 2007), and the Vosges mountains were 70% forest on average (DREAL Grand Est, 2018). The rest of the area was permanent pastures, able land, and human settlements. Sampling occurred over 6 years, between January 2011 and April 2016, mostly in winter and spring, with surveys lasting between 2 and 4 months. We considered two study areas in 2011, 2014, and 2015, three study areas in 2013, and one study area in 2012 and 2016 through camera trapping (Figure 1).

We divided each study area into a grid of 2.7 × 2.7 km cells applying a systematic design where one out of two cells was sampled (Zimmermann, Breitenmoser-Würsten, Molinari-Jobin, & Breitenmoser, 2013), hence ensuring that at least one camera trap was set in each potential lynx home range (between 100 and 250 km², see Breitenmoser-Würsten et al., 2007). To maximize detectability, we set (nonbaited) camera traps in forested habitats, based on previous signs of lynx presence and on local knowledge, at optimal locations where landscape and terrain features were likely to channel lynx movements on more predictable paths (on forest roads, hiking trails, and to a lesser extent on game paths; Blanc, Marboutin, Gatti, & Gimenez, 2013). Camera was settled within a single session continuously during 60 days between February and beginning of March with little variation between sites.
**FIGURE 1** Map of the study area in the French Jura and Vosges mountains. The study area encompassed seven counties (Ain, Jura and Doubs in the Jura mountains and Vosges, Haut-Rhin, Bas-Rhin and Moselle in the Vosges mountains) that were monitored through 413 camera trapping sites (298 in the Jura mountains and 115 in the Vosges mountains; two camera traps were set per site), each within a 2.7 × 2.7 km cell. The inset map represents the French counties (gray borders), the counties that were considered in the study (black borders), the Jura mountains (green shaded area) and the Vosges mountains (red shaded area)

<figure><img src="image_3.png" /></figure>

At each trapping location, we set two Xenon white flash camera traps (models: Capture, Ambush and Attack; Cuddeback) with passive infrared trigger mechanisms to photograph both flanks of an animal. We checked camera traps weekly to change memory cards, batteries and to remove fresh snow after heavy snowfall. Based on unique coat patterns, we identified individual lynx on photographs (Zimmermann & Foresti, 2016). The recognition of individual was computer-induced, not fully automated. We used the Extract-compare® software that compares the lynx spot pattern with a library of previously extracted pattern and proposes potential matches according to a score (http://conservationresearch.org.uk/Home/ExtractCompare). The observer can confirm the lynx identification or not and browse through the highest-ranking proposed matches. The final decision is made by the observer based on an additional visual examination of the entire photograph set for this particular lynx. Pictures for which no match was found with the software were visually checked against our entire photograph library. Only when the match was undeniable was the individual recorded as a match, otherwise it was recorded as a new individual. All captures that did not fit automated or associated visual confirmation with no doubt, because of a poor picture quality (e.g., blurry, overexposed), were classified as “unconfirmed” and excluded from the analyses. We recorded the date, time, sex whenever possible, and location of each photographic capture of a lynx. During the time of year our study took place, juvenile lynx (<1 year old) can still be with their mother (Zimmermann, Breitenmoser-Würsten, & Breitenmoser, 2005). In our analysis, we retained only independent lynx, that is, adult lynx

or emancipated individuals based on physical characteristics or previous knowledge of their age or status (from photographic evidence). We defined a capture occasion as 5 successive trap nights (Blanc et al., 2013), dissociating trapping events from individual photograph to avoid pseudo-replications.

## 2.3 | Spatial capture-recapture analyses

We used spatial capture-recapture (SCR) models to estimate lynx densities (Royle et al., 2014). In contrast with standard (nonspatial) capture-recapture models, SCR models use the spatial locations of captures to infer the activity center (or home range) of each individual. We assumed that individual encounters are Bernoulli random variables with individual- and trap-specific detection probabilities. More precisely, the detection probability $p_{ij}$ of an individual $i$ at trap $j$ is assumed to decrease as the distance ($d_{ij}$) from its activity center increases according to a detection function. We used the half-normal detection function, $p_{ij} = p_0 \exp(-d_{ij}^2/(2\sigma^2))$, where $p_0$ is the probability of detecting an individual when the trap is located exactly at its center of activity and $\sigma$ is the spatial scale (or movement) parameter that controls the shape of the detection function. For one of the two study areas in the French Jura mountains in years 2011 and 2013, we detected only a few individuals (see the columns Doubs in Table 1). To increase the effective sample size, we combined the data from these two sampling areas using common detection and spatial parameters for both areas, while estimating density
<table><caption>TABLE 1 Main characteristics and results of the lynx camera trap survey carried out in (a) the French Jura mountains and (b) the French Vosges mountains</caption><thead><tr><th>(a) Year/County</th><th>2011/Doubs</th><th>2011/Jura</th><th>2012/Jura & Doubs</th><th>2013/Doubs</th><th>2013/Ain & Jura</th><th>2014/Ain</th><th>2015/Ain</th></tr></thead><tbody><tr><td>Period of trap activity</td><td>January-April</td><td>February-April</td><td>February-April</td><td>February-April</td><td>February-April</td><td>February-April</td><td>February-May</td></tr><tr><td>Number of active camera traps</td><td>48</td><td>66</td><td>148</td><td>44</td><td>142</td><td>118</td><td>30</td></tr><tr><td>Number of trapping days (average/area)</td><td>63</td><td>59</td><td>69</td><td>63</td><td>58</td><td>59</td><td>99</td></tr><tr><td>Number of capture occasions<sup>i</sup></td><td>15</td><td>15</td><td>17</td><td>14</td><td>13</td><td>13</td><td>21</td></tr><tr><td>Number of detections</td><td>22</td><td>42</td><td>130</td><td>25</td><td>117</td><td>158</td><td>38</td></tr><tr><td>Number of detected individuals</td><td>4</td><td>9</td><td>21</td><td>6</td><td>19</td><td>23</td><td>10</td></tr><tr><td>Number of females, unknown, males</td><td>1, 1, 2</td><td>1, 7, 1</td><td>2, 14, 5</td><td>1, 4, 1</td><td>2, 13, 4</td><td>4, 16, 3</td><td>2, 8, 0</td></tr><tr><td>Number of detections/ind: mean, min, max</td><td>3, 2, 4</td><td>2, 8, 1, 6</td><td>2, 5, 1, 10</td><td>2, 7, 1, 6</td><td>3, 6, 1, 11</td><td>3, 3, 1, 9</td><td>2, 2, 1, 5</td></tr></tbody></table><table><thead><tr><th>(b) Year/County</th><th>2013/Haut-Rhin & Vosges</th><th>2014/Bas-Rhin & Moselle</th><th>2015/Bas-Rhin & Moselle</th><th>2016/Bas & Haut-Rhin</th></tr></thead><tbody><tr><td>Period of trap activity</td><td>December-January</td><td>February-April</td><td>February-April</td><td>February-April</td></tr><tr><td>Number of active traps</td><td>60</td><td>50</td><td>60</td><td>60</td></tr><tr><td>Number of trapping days (average/area)</td><td>52</td><td>59</td><td>57</td><td>59</td></tr><tr><td>Number of detections</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>

separately (e.g., Rocha, Sollmann, Ramalho, Ilha, & Tan, 2016). We defined a state-space, that is, the area encompassing all potential activity centers of the observed individuals, by building a grid that buffered outermost camera trap locations by 15 km (corresponding to at least $2σ$; Royle et al., 2014) with a resolution of 1.5 km (or pixels of area 2.25 km²). We fitted SCR models in the maximum likelihood framework using the R package oSCR (Sutherland, Muñoz, Miller, & Grant, 2016; Sutherland, Royle, & Linden, 2019).

For comparison, we also estimated abundance using standard (nonspatial) capture-recapture models (Otis et al., 1978). We dropped the spatial information and considered only the detections and nondetections for each individual. We considered two models, M0 in which the detection probability is the same for all individuals, and Mh in which the detection probability varies among individuals. We fitted standard models in the maximum likelihood framework using the R package Rcapture (Baillargeon & Rivest, 2007). We estimated density as the ratio of estimated abundance over an effective trapping area (ETA). ETA was estimated by adding a buffer to the trapping area equal to the mean maximum distance moved (MMDM) or half of it (HMMDM). We calculated the MMDM by averaging the maximum distances between capture locations for all individuals detected at more than one site.

## 3 | RESULTS

We collected data from 413 camera trapping sites (two camera traps were set per site) resulting in 25,839 trap days (Table 1). In total, we identified 92 lynx over 532 detection events in the Jura mountains,

including 16 females, 13 males, and 63 individuals of unknown sex. The number of detections per individual was 2.6 on average and varied from 1 to 11. In contrast, we collected no lynx photograph in the Vosges mountains; therefore, we did not proceed with analyses for this area.

For the Jura mountains, abundance estimates were similar whether we used spatial or nonspatial models, although always slightly higher for the former. Estimated abundance among study areas varied between 5 ($SE = 0.1$) and 29 (0.2) lynx in the spatial analyses, between 4 (0.7) and 23 (0.7) with model M0, and between 5 (1.7) and 28 (3.6) with model Mh. Estimated density varied between 0.24 (0.02) and 0.91 (0.03) lynx per 100 km² in the spatial analyses (Table 2). In the nonspatial analyses, the density varied between 0.31 (0.05) and 0.78 (0.02) lynx per 100 km² under model M0 and between 0.34 (0.06) and 0.95 (0.12) under model Mh when the MMDM was used. When we used HMMDM, the density varied between 0.57 (0.10) and 1.46 (0.16) lynx per 100 km² under model M0 and between 0.67 (0.12) and 1.43 (0.16) under model Mh.

From the spatial analyses, we used the model estimates to produce density surfaces within the state-space (Figure 2). The density per pixel of area 2.25 km² ranged from 0 to 0.20 individuals in the Jura mountains.

## 4 | DISCUSSION

By using camera trap sampling and SCR models, we provided the first multi-site density estimates for lynx that will help in setting a baseline conservation status for the French lynx population. The

<footnote>iA capture occasion is defined as 5 successive trap days.</footnote>
<table><caption>**TABLE 2** Lynx abundance and density estimates obtained from spatial and nonspatial capture-recapture analyses of camera trapping data collected in the French Jura mountains</caption><thead><tr><th>Year/County</th><th>2011/Doubs</th><th>2011/Jura</th><th>2012/Jura-Doubs</th><th>2013/Doubs</th><th>2013/Ain-Jura</th><th>2014-Ain</th><th>2015-Ain</th></tr></thead><tbody><tr><td>SCR abundance (SE)</td><td>5 (0.1)</td><td>12 (0.1)</td><td>29 (0.2)</td><td>7 (0.1)</td><td>21 (0.1)</td><td>29 (0.2)</td><td>12 (0.1)</td></tr><tr><td>SCR density (SE)</td><td>0.24 (0.02)</td><td>0.44 (0.02)</td><td>0.67 (0.02)</td><td>0.36 (0.02)</td><td>0.54 (0.02)</td><td>0.91 (0.03)</td><td>0.64 (0.03)</td></tr><tr><td>p<sub>0</sub> logit scale (SE)</td><td>-2.94 (0.24)</td><td></td><td>-2.01 (0.20)</td><td>-2.57 (0.20)</td><td></td><td>-2.34 (0.19)</td><td>-3.01 (0.42)</td></tr><tr><td>σ log scale (SE)</td><td>8.89 (0.14)</td><td></td><td>8.54 (0.08)</td><td>8.95 (0.06)</td><td></td><td>8.80 (0.07)</td><td>8.97 (0.19)</td></tr><tr><td>M0 abundance (SE)</td><td>4 (0.7)</td><td>9 (0.7)</td><td>21 (0.6)</td><td>6 (0.3)</td><td>19 (0.8)</td><td>23 (0.7)</td><td>11 (1.2)</td></tr><tr><td>Mh abundance (SE)</td><td>5 (1.7)</td><td>10 (1.8)</td><td>25 (2.8)</td><td>7 (1.2)</td><td>25 (4.1)</td><td>28 (3.6)</td><td>11 (1.2)</td></tr><tr><td>MMDM (km)</td><td>9.1</td><td>16.2</td><td>8.9</td><td>9.1</td><td>18.2</td><td>13.6</td><td>12.1</td></tr><tr><td>ETA with MMDM (km<sup>2</sup>)</td><td>1,991</td><td>2,930</td><td>3,089</td><td>1,171</td><td>4,954</td><td>2,936</td><td>1,549</td></tr><tr><td>M0 density MMDM (SE)</td><td>0.31 (0.05)</td><td>0.31 (0.02)</td><td>0.68 (0.02)</td><td>0.51 (0.02)</td><td>0.38 (0.02)</td><td>0.78 (0.02)</td><td>0.71 (0.08)</td></tr><tr><td>Mh density MMDM (SE)</td><td>0.39 (0.13)</td><td>0.34 (0.06)</td><td>0.81 (0.09)</td><td>0.60 (0.10)</td><td>0.50 (0.08)</td><td>0.95 (0.12)</td><td>0.70 (0.08)</td></tr><tr><td>ETA with HMMDM (km<sup>2</sup>)</td><td>697</td><td>1,491</td><td>2,111</td><td>659</td><td>2,673</td><td>1,668</td><td>753</td></tr><tr><td>M0 density HMMDM (SE)</td><td>0.57 (0.10)</td><td>0.60 (0.05)</td><td>0.99 (0.03)</td><td>0.91 (0.05)</td><td>0.71 (0.03)</td><td>1.38 (0.04)</td><td>1.46 (0.16)</td></tr><tr><td>Mh density HMMDM (SE)</td><td>0.72 (0.24)</td><td>0.67 (0.12)</td><td>1.18 (0.13)</td><td>1.06 (0.18)</td><td>0.93 (0.15)</td><td>1.68 (0.21)</td><td>1.43 (0.16)</td></tr></tbody></table>

Note: Densities are provided in number of lynx per 100 km². For 2011 and 2013, parameters of the spatial capture-recapture model ($p_0$ and $σ$) are common to both areas in each year.

Abbreviations: ETA, effective trapping area; MMDM, half mean maximum distance moved; M0, the (nonspatial) capture-recapture model with homogeneous detection probability; Mh, the (nonspatial) capture-recapture model with heterogeneous detection probability; MMDM, mean maximum distance moved; SCR, spatial capture-recapture; SE, standard error.

multi-site dimension of our study allows exploring variability in the density estimates across landscapes. Our study is yet another example of the potential of combining SCR methods and noninvasive sampling techniques to estimate abundance and density for elusive and wide-ranging species, like large carnivores (Alexander et al., 2015; Broekhuis & Gopalaswamy, 2016; Goldberg et al., 2015; López-Bao et al., 2018; Pesenti & Zimmermann, 2013; Stetz et al., 2018).

When examining densities across study areas in the French Jura mountains, we found spatial variation between the three counties, with Doubs area having the lowest densities, Ain the highest densities, and Jura intermediate densities. Our density estimates

were of similar magnitude to other lynx populations in Europe: 1.47 and 1.38 lynx/100 km² in the Northwestern Swiss Alps (Pesenti & Zimmermann, 2013), 0.58 (Štiavnica mountains) and 0.81 individuals/100 km² (Velká Fatra National Park) in Slovakia (Kubala et al., 2017) and 0.9 individuals/100 km² in the Bavarian Forest National Park in Germany (Weingarth et al., 2012).

While Kubala et al. (2017) and Pesenti and Zimmermann (2013) used SCR models, (Weingarth et al., 2012) used standard capture-recapture models with MMDM to estimate densities, which makes them difficult to compare (Gerber, Karpanty, & Kelly, 2012). Indeed, in other carnivore studies, the use of MMDM also produced similar

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Lynx (Lynx lynx) density maps in the French Jura mountains. The density scale is in lynx per 2.25 km² (pixel resolution is 1,500 m × 1,500 m). We obtained the estimated abundance in each map by summing up the densities in each pixel altogether. Yellow is for low densities, green for medium densities, and blue for high densities; the density scales are specific to each map. Note that the interpretation of these plots as density maps is subject to caution (see the vignette “secr-densitysurface” of the SECR R package; Efford, 2019)</figcaption></figure>
density estimates to SCR models (Pesenti & Zimmermann, 2013), while in others, including ours, the SCR estimates were closer to the MMDM estimates (Obbard et al., 2010) or intermediate between the MMDM and HMMDM estimates (Reppucci, Gardner, & Lucherini, 2011). When looking at reference values for densities across the distribution range of the species, it may be biologically meaningful to use the MMDM density estimate as a reference as it covers the whole potential of animal movements. On the other hand, because SCR models make space explicit whereas standard model-based densities are sensitive to the definition of the effective sampling area, we recommend the use of SCR models to estimate lynx densities.

Our lynx density estimates might suffer from potential sources of bias that need to be discussed. First, the period of sampling is important to account for when setting up camera trap surveys (Weingarth et al., 2015). We conducted our survey outside the dispersal period, during the lynx mating season (February–March mostly). We did so to avoid capturing transient individuals and to increase detectability because of high lynx activity and relatively reduced human activities (Zimmermann & Foresti, 2016). However, some individuals might have moved in and out of the study areas, especially males who cover greater distances during the mating season. Whereas the presence of nonresident individuals can affect the calculation of (H) MMDM, and in turn density estimated with standard capture-recapture models, SCR density estimates were found to be robust to the presence of transient individuals (Royle, Fuller, & Sutherland, 2016). Second, males have larger home ranges than females (Pesenti & Zimmermann, 2013), which leads to heterogeneity in the SCR model parameter estimates. Because there were too few males and females identified and lots of individuals with unknown sex, sex-specific SCR analyses (Sollmann et al., 2011) produced unreliable abundance and density estimates (results not shown). If detection heterogeneity is ignored in capture-recapture models, abundance is underestimated (Cubaynes et al., 2010), therefore our density estimates are probably biased low and should be considered as a conservative metric. The determination of sex could be improved by (a) combining the photographic surveys with genetic surveys, (b) conducting deterministic surveys over several years (e.g., Pesenti & Zimmermann, 2013), (c) conducting an opportunistic camera trapping survey all over years and setting camera traps at fresh lynx kills, (d) setting infrared flash camera traps capable of taking burst of images in rapid sequence at marking sites regularly used by the lynx (e.g., Vogt, Zimmermann, Kölliker, & Breitenmoser, 2014).

Last, we did not detect any individuals in the Vosges mountains, even though the sampling effort was similar to that implemented in the Jura mountains (Table 1). Despite the release of 21 lynx between 1983 and 1993 (out of which only 10 survived; Vandel et al., 2006), connectivity with the Jura mountains was and still is difficult to emerge because of artificial habitat fragmentation (highways and high-speed train railways). In recent years, very few signs of presence of lynx have been collected in the Vosges mountains through opportunistic monitoring (mostly direct observations, more rarely footprints or hairs). In 2018, the regular presence area of lynx was of 400 km² with signs of presence both in the north, in the center,

and south of the massif. Currently, only two males are identified with photographs. One came from the Palatinate Forest in Germany (from a reintroduction program) and installed his home range in the Hautes-Vosges in 2017. The second one came from the Jura mountains from where he dispersed in 2015. There are also some punctual incursions of lynx from Palatinate forest in the north Vosges. Overall, our findings are likely to be representative of the current critical situation of the lynx in the Vosges mountains.

We envision several perspectives to our work. First, while density estimates are of primary interest for conservation, understanding the mechanisms underlying trends in abundance is required to make sound conservation decisions (Williams et al., 2002). SCR models have been extended to open populations (Gardner, Reppucci, Lucherini, & Royle, 2010) and can be used to estimate demographic parameters (survival, reproduction) of large carnivores (Whittington & Sawaya, 2015). Unfortunately, because of logistic constraints, we could not sample the same areas over several years, which precludes a standard application of these models. A solution may lie in the combination of the data we collected through systematic camera trap surveys with additional data in the SCR framework, such as occupancy data (Blanc, Marboutin, Gatti, Zimmermann, & Gimenez, 2014) or opportunistic camera trap data (Tenan, Pedrini, Bragalanti, Groff, & Sutherland, 2017). Second, in addition to traffic-induced mortality and conflicts with human activities, the expansion of lynx populations is limited by habitat fragmentation (Kramer-Schadt, Revilla, Wiegand, & Breitenmoser, 2004), hence the need to assess connectivity with other populations (Zimmermann & Breitenmoser, 2007). SCR models can be used to quantify landscape connectivity by replacing the Euclidean distance between camera traps and home range centers by the least-cost path (Royle, Chandler, Gazenski, & Graves, 2013; Sutherland, Fuller, & Royle, 2015). For lynx, this will require setting up traps across a gradient of habitat types, not only forested habitats, so that resistance to movement can be estimated.

In conclusion, our lynx density estimates for the French Jura mountains complement nicely the estimates recently provided for the Northwestern Swiss Alps (Pesenti & Zimmermann, 2013). The use of camera trapping coupled with SCR models in both France and Switzerland was the result of a cooperation between the two countries with the perspective of a transboundary monitoring (Gervasi et al., 2016; Vitkalova et al., 2018). This approach would prove useful to accurately estimate densities in other areas where habitats and prey availability might differ, and overall lynx detectability varies. Also, collecting and adding movement data from GPS-collared lynx would be useful (Linden, Sirén, & Pekins, 2018; Tenan et al., 2017) to try and infer the connections between subpopulations.

The case can be made for monitoring the return of the lynx in the French Alps. Indeed, small-scale camera trapping surveys and opportunistic observations are currently active and producing signs of lynx presence. However, the lack of a coordinated and intensive sampling effort prevents us from being able to estimate abundance and density and inferring trends.

In contrast, the situation in the Vosges mountains is alarming with no individuals detected over the study period. Because the
Vosges mountains are located between the French Jura mountains and the Palatinate Forest in Germany where a reintroduction program is ongoing (program LIFE13 NAT/DE/000755), the lynx colonization in the Vosges mountains remains possible both by the north and the south. Incidentally, two cases of lynx dispersal in the Vosges mountains from neighboring mountains have been recently observed (Hurstel & Laurent, 2016; program LIFE13 NAT/DE/000755). To ensure the detection of lynx in the Vosges mountains, we recommend reinforcing collaborative monitoring by involving all field stakeholders and enhancing communication on the species signs of presence.

In this context, obtaining accurate and comparable lynx densities will be crucial to closely monitor population trends at the national scale and inform management policies for the effective conservation of the Eurasian lynx in France.

## ACKNOWLEDGMENTS

We thank the staff from the French National Game and Wildlife Agency (ONCFS), the CROC, the Forest National Agency (ONF), the "Directions Départementales des Territoires," the "Fédérations Départementales des Chasseurs", the Regional Natural Parks, the environmental protection associations and all the volunteers from the "Réseau Loup Lynx" who collected the photographs during the camera trapping session. OG was funded by CNRS and the "Mission pour l'Interdisciplinarité" through the "Osez l'Interdisciplinarité" initiative. CD, EM, and SG were funded by ONCFS. CEFE and CROC were funded by CIL&B, MTES (ITTECOP), and FRB through the research program ERC-Lynx. CROC was funded by the European Union within the framework of the Operational Program FEDER-FSE "Lorraine et Massif des Vosges 2014-2020," the "Commissariat à l'Aménagement du Massif des Vosges" for the FNADT ("Fonds National d'Aménagement et de Développement du Territoire"), the DREAL Grand Est ("Direction Régionale de l'Environnement, l'Aménagement et le Logement"), the "Région Grand Est," the "Zoo d'Amnéville," the "Fondation d'entreprise UEM," the "Fondation Nature & Découvertes," the "Fondation Le Pal Nature," and the *Chasseur d'images* magazine. This study could not have been conducted without authorizations from and agreements with municipalities, environmental managers, and owners.

## CONFLICT OF INTEREST

None declared.

## AUTHOR CONTRIBUTIONS

OG wrote the paper and all co-authors commented on the manuscript. OG and SG analyzed the data. AL, CD, EG, EM, and SG coordinated the study designs, the data collection, and interpretation, with help from FZ for setting the experimental design in the Jura mountains.

## DATA AVAILABILITY STATEMENT

The Eurasian lynx is an endangered species with high conservation stakes. Interactions with human activities are problematic and lead to poaching and anthropogenic pressures. Providing accurate information on lynx locations can be detrimental to the conservation status of the species. As a consequence, the original data could not be shared.

## ORCID

Olivier Gimenez ID https://orcid.org/0000-0001-7001-5142

Fridolin Zimmermann ID https://orcid.org/0000-0003-3861-246X

## REFERENCES

- Alexander, J. S., Gopalaswamy, A. M., Shi, K., & Face, R. P. (2015). Value: Towards robust estimates of snow leopard densities. *PLoS ONE*, *10*, e0134815. https://doi.org/10.1371/journal.pone.0134815

- Baillargeon, S., & Rivest, L.-P. (2007). Rcapture: Loglinear models for capture-recapture in R. *Journal of Statistical Software*, *19*, 1–31.

- Blanc, L., Marboutin, E., Gatti, S., & Gimenez, O. (2013). Abundance of rare and elusive species: Empirical investigation of closed versus spatially explicit capture-recapture models with lynx as a case study. *The Journal of Wildlife Management*, *77*, 372–378. https://doi.org/10.1002/jwmg.453

- Blanc, L., Marboutin, E., Gatti, S., Zimmermann, F., & Gimenez, O. (2014). Improving abundance estimation by combining capture-recapture and occupancy data: Example with a large carnivore. *Journal of Applied Ecology*, *51*, 1733–1739. https://doi.org/10.1111/1365-2664.12319

- Borchers, D. (2012). A non-technical overview of spatially explicit capture-recapture models. *Journal of Ornithology*, *152*, 435–444. https://doi.org/10.1007/s10336-010-0583-z

- Borchers, D. L., & Efford, M. G. (2008). Spatially explicit maximum likelihood methods for capture-recapture studies. *Biometrics*, *64*, 377–385. https://doi.org/10.1111/j.1541-0420.2007.00927.x

- Breitenmoser, U. (1998). Large predators in the Alps: The fall and rise of man's competitors. *Biological Conservation*, *83*, 279–289. https://doi.org/10.1016/S0006-3207(97)00084-0

- Breitenmoser, U., Breitenmoser-Würsten, C., & Capt, S. (1998). Introduction and present status of the lynx (Lynx lynx) in Switzerland. *Hystrix*, the *Italian Journal of Mammalogy*, *10*. https://doi.org/10.4404/hystrix-10.1-4118

- Breitenmoser, U., Breitenmoser-Würsten, C., Capt, S., Molinari-Jobin, A., Molinari, P., & Zimmermann, F. (2007). Conservation of the lynx Lynx lynx in the Swiss Jura Mountains. *Wildlife Biology*, *13*, 340–355. https://doi.org/10.2981/0909-6396(2007)13[340:COTLL]2.0.CO;2

- Breitenmoser-Würsten, C., Zimmermann, F., Stahl, P., Vandel, J.-M., Molinari-Jobin, A., Molinari, P., ... Breitenmoser, U. (2007). Spatial and Social stability of a Eurasian lynx Lynx lynx population: An assessment of 10 years of observation in the Jura Mountains. *Wildlife Biology*, *13*, 365–380. https://doi.org/10.2981/0909-6396(2007)13[365:SASSO]A2.0.CO;2

- Broekhuis, F., & Gopalaswamy, A. M. (2016). Counting cats: Spatially explicit population estimates of cheetah (Acinonyx jubatus) using unstructured sampling data. *PLoS ONE*, *11*, e0153875. https://doi.org/10.1371/journal.pone.0153875

- Chapron, G., Kaczensky, P., Linnell, J. D. C., von Arx, M., Huber, D., Andren, H., ... Boitani, L. (2014). Recovery of large carnivores in
Europe's modern human-dominated landscapes. *Science*, 346, 1517-1519. https://doi.org/10.1126/science.1257553

Cubaynes, S., Pradel, R., Choquet, R., Duchamp, C., Gaillard, J.-M., Lebreton, J.-D., ... Gimenez, O. (2010). Importance of accounting for detection heterogeneity when estimating abundance: The case of French wolves. *Conservation Biology*, 24, 621–626. https://doi.org/10.1111/j.1523-1739.2009.01431.x

da Rocha, D. G., Sollmann, R., Ramalho, E. E., Ilha, R., & Tan, C. K. W. (2016). Ocelot (*Leopardus pardalis*) density in Central Amazonia. *PLoS ONE*, 11, e0154624. https://doi.org/10.1371/journal.pone.0154624

DREAL Grand Est (2018). *Diagnostic partagé du massif des Vosges pour une vision transversale du territoire et de ses enjeux*. Retrieved from https://www.grandest.fr/piece-n07-rapport-du-sraddet-1sur3-diagnostic

Efford, M. (2004). Density estimation in live-trapping studies. *Oikos*, 106, 598–610. https://doi.org/10.1111/j.0030-1299.2004.13043.x

Efford, M. (2019). *secr: Spatially explicit capture-recapture models*. R package version 3.2.1. [Internet]. Retrieved from https://CRAN.R-proje ct.org/package=secr

Gardner, B., Reppucci, J., Lucherini, M., & Royle, J. A. (2010). Spatially explicit inference for open populations: Estimating demographic pa rameters from camera-trap studies. *Ecology*, 91, 3376–3383. https:// doi.org/10.1890/09-0804.1

Gerber, B. D., Ivan, J. S., & Burnham, K. P. (2014). Estimating the abun dance of rare and elusive carnivores from photographic-sampling data when the population size is very small. *Population Ecology*, 56, 463–470. https://doi.org/10.1007/s10144-014-0431-8

Gerber, B. D., Karpanty, S. M., & Kelly, M. J. (2012). Evaluating the po tential biases in carnivore capture-recapture studies associated with the use of lure and varying density estimation techniques using pho tographic-sampling data of the Malagasy civet. *Population Ecology*, 54, 43–54. https://doi.org/10.1007/s10144-011-0276-3

Gervasi, V., Brøseth, H., Gimenez, O., Nilsen, E. B., Odden, J., Flagstad, Ø., & Linnell, J. D. C. (2016). Sharing data improves monitoring of trans-boundary populations: The case of wolverines in central Scandinavia. *Wildlife Biology*, 22, 95–106. https://doi.org/10.2981/wlb.00142

Goldberg, J. F., Tempa, T., Norbu, N., Hebblewhite, M., Mills, L. S., Wangchuk, T. R., & Lukacs, P. (2015). Examining temporal sample scale and model choice with spatial capture-recapture models in the common leopard *Panthera pardus*. *PLoS ONE*, 10, e0140757. https:// doi.org/10.1371/journal.pone.0140757

Hurstel, A., & Laurent, A. (2016). Première preuve de dispersion du Lynx d'Eurasie (*Lynx lynx*) du Jura vers les Vosges. *Ciconia*, 40, 1-6.

Kelly, M., Betsch, J., Wultsch, C., Mesa, B., & Mills, L. (2012). Noninvasive sampling for carnivores. In L. Boitani & R. A. Powell (Eds.), *Carnivore ecology and conservation: A handbook of techniques* (pp. 47-69). Oxford, UK: Oxford University Press.

Kramer-Schadt, S., Revilla, E., Wiegand, T., & Breitenmoser, U. (2004). Fragmented landscapes, road mortality and patch con nectivity: Modelling influences on the dispersal of Eurasian lynx. *Journal of Applied Ecology*, 41, 711–723. https://doi.org/10.1111/j.0021-8901.2004.00933.x

Kubala, J., Smolko, P., Zimmermann, F., Rigg, R., Tám, B., Il'ko, T., ... Breitenmoser, U. (2017). Robust monitoring of the Eurasian lynx Lynx lynx in the Slovak Carpathians reveals lower numbers than officially reported. *Oryx*, 53, 548–556. https://doi.org/10.1017/S003060531 700076X

Laurent, A., Léger, F., Briaudet, P. E., Léonard, Y., Bataille, A., & Goujon, G. (2012). Evolution récente (2008-2010) de la population de Lynx en France. *Faune Sauvage*, 294, 38-39.

Linden, D. W., Sirén, A. P. K., & Pekins, P. J. (2018). Integrating telem etry data into spatial capture-recapture modifies inferences on multi-scale resource selection. *Ecosphere*, 9, e02203. https://doi. org/10.1002/ecs2.2203

López-Bao, J. V., Godinho, R., Pacheco, C., Lema, F. J., García, E., Llaneza, L., ... Jiménez, J. (2018). Toward reliable population estimates of

wolves by combining spatial capture-recapture models and non- invasive DNA monitoring. *Scientific Reports*, 8, 2177. https://doi. org/10.1038/s41598-018-20675-9

Marboutin, E., Duchamp, C., Moris, P., Briaudet, P.-E., Léonard, Y., & Catusse, M. (2012). Survey of the Lynx distribution in the French Alps: 2005-2009 update. *Acta Biologica Slovenica*, 55, 29-34.

Mumma, M. A., Zieminski, C., Fuller, T. K., Mahoney, S. P., & Waits, L. P. (2015). Evaluating noninvasive genetic sampling techniques to es timate large carnivore abundance. *Molecular Ecology Resources*, 15, 1133–1144. https://doi.org/10.1111/1755-0998.12390

Obbard, M. E., Howe, E. J., & Kyle, C. J. (2010). Empirical comparison of density estimators for large carnivores. *Journal of Applied Ecology*, 47, 76–84. https://doi.org/10.1111/j.1365-2664.2009.01758.x

Otis, D. L., Burnham, K. P., White, G. C., & Anderson, D. R. (1978). Statistical inference from capture data on closed animal populations. *Wildlife Monographs*, 76, 3–135. https://doi.org/10.2307/2287873

Pesenti, E., & Zimmermann, F. (2013). Density estimations of the Eurasian lynx (*Lynx lynx*) in the Swiss Alps. *Journal of Mammalogy*, 94, 73-81. https://doi.org/10.1644/11-MAMM-A-322.1

Reppucci, J., Gardner, B., & Lucherini, M. (2011). Estimating detection and density of the Andean cat in the high Andes. *Journal of Mammalogy*, 92, 140–147. https://doi.org/10.1644/10-MAMM-A-053.1

Royle, J. A., Chandler, R. B., Gazenski, K. D., & Graves, T. A. (2013). Spatial capture-recapture models for jointly estimating population den sity and landscape connectivity. *Ecology*, 94, 287-294. https://doi. org/10.1890/12-0413.1

Royle, J. A., Chandler, R. B., Sollmann, R., & Gardner, B. (2014). *Spatial capture-recapture* [Internet]. New York, NY: Academic Press. https:// doi.org/10.1016/C2012-0-0122-7

Royle, J. A., Fuller, A. K., & Sutherland, C. (2016). Spatial capture-recapture models allowing Markovian transience or dispersal. *Population Ecology*, 58, 53–62. https://doi.org/10.1007/s10144-015-0524-z

Royle, J. A., & Young, K. V. (2008). A hierarchical model for spa tial capture-recapture data. *Ecology*, 89, 2281–2289. https://doi. org/10.1890/07-0601.1

Sollmann, R., Furtado, M. M., Gardner, B., Hofer, H., Jácome, A. T. A., Tôrres, N. M., & Silveira, L. (2011). Improving density estimates for elusive carnivores: Accounting for sex-specific detection and movements using spatial capture-recapture models for jaguars in central Brazil. *Biological Conservation*, 144, 1017-1024. https://doi.org/10.1016/j.biocen.2010.12.011

Stetz, J. B., Mitchell, M. S., & Kendall, K. C. (2018). Using spatially-explicit capture-recapture models to explain variation in seasonal density patterns of sympatric ursids. *Ecography*, 42, 237-248. https://doi.org/10.1111/ecog.03556

Sutherland, C., Fuller, A. K., & Royle, J. A. (2015). Modelling non-Euclidean movement and landscape connectivity in highly structured ecologi cal networks. *Methods in Ecology and Evolution*, 6, 169-177. https:// doi.org/10.1111/2041-210X.12316

Sutherland, C., Muñoz, D. J., Miller, D. A. W., & Grant, E. H. C. (2016). Spatial capture-recapture: A promising method for analyzing data collected using artificial cover objects. *Herpetologica*, 72, 6-12. https://doi.org/10.1655/HERPETOLOGICA-D-15-00027

Sutherland, C., Royle, J. A., & Linden, D. W. (2019) oSCR: A spatial capture-recapture R package for inference about spatial ecologi cal processes. *Ecography*, 42, 1459-1469. https://doi.org/10.1111/ ecog.04551

Tenan, S., Pedrini, P., Bragalanti, N., Groff, C., & Sutherland, C. (2017). Data integration for inference about spatial processes: A model-based approach to test and account for data inconsis tency. *PLoS ONE*, 12, e0185588. https://doi.org/10.1371/journ al.pone.0185588

Vandel, J.-M., & Stahl, P. (2005). Distribution trend of the Eurasian lynx Lynx lynx populations in France. *Mammalia*, 69, 145-158. https://doi.org/10.1515/mamm.2005.013
Vandel, J.-M., Stahl, P., Herrenschmidt, V., & Marboutin, E. (2006). Reintroduction of the lynx into the Vosges mountain massif: From animal survival and movements to population development. *Biological Conservation*, 131, 370-385. https://doi.org/10.1016/j.biocon.2006.02.012

Vitkalova, A. V., Feng, L., Rybin, A. N., Gerber, B. D., Miquelle, D. G., Wang, T., ... Ge, J. (2018). Transboundary cooperation improves endangered species monitoring and conservation actions: A case study of the global population of Amur leopards. *Conservation Letters*, 11, e12574. https://doi.org/10.1111/conl.12574

Vogt, K., Zimmermann, F., Kölliker, M., & Breitenmoser, U. (2014). Scent-marking behaviour and social dynamics in a wild population of Eurasian lynx Lynx lynx. *Behavioural Processes*, 106, 98-106. https://doi.org/10.1016/j.beproc.2014.04.017

Weingarth, K., Heibl, C., Knauer, F., Zimmermann, F., Bufka, L., & Heurich, M. (2012). First estimation of Eurasian lynx (Lynx lynx) abundance and density using digital cameras and capture-recapture techniques in a German national park. *Animal Biodiversity and Conservation*, 35, 197-207.

Weingarth, K., Zeppenfeld, T., Heibl, C., Heurich, M., Bufka, L., Daniszová, K., & Müller, J. (2015). Hide and seek: Extended camera-trap session lengths and autumn provide best parameters for estimating lynx densities in mountainous areas. *Biodiversity and Conservation*, 24, 2935-2952. https://doi.org/10.1007/s10531-015-0986-5

White, G. C., Anderson, D. R., Burnham, K. P., & Otis, D. L. (1982). *Capture-recapture and removal methods for sampling closed populations* (p. 14). Los Alamos, NM: Los Alamos National Laboratory. Retrieved from http://pubs.er.usgs.gov/publication/70120390

Whittington, J., & Sawaya, M. A. (2015). A comparison of grizzly bear demographic parameters estimated from non-spatial and spatial open

population capture-recapture models. *PLoS ONE*, 10, e0134446. https://doi.org/10.1371/journal.pone.0134446

Williams, B. K., Nichols, J. D., & Conroy, M. J. (2002). *Analysis and management of animal populations*. San Diego, CA: Academic Press.

Zimmermann, F., & Breitenmoser, U. (2007). Potential distribution and population size of the Eurasian lynx Lynx lynx in the jura Mountains and possible corridors to adjacent ranges. *Wildlife Biology*, 13, 406-416. https://doi.org/10.2981/0909-6396(2007)13[406:PDAPS0]2.0.CO;2

Zimmermann, F., Breitenmoser-Würsten, C., & Breitenmoser, U. (2005). Natal dispersal of Eurasian lynx (Lynx lynx) in Switzerland. *Journal of Zoology*, 267, 381-395. https://doi.org/10.1017/S0952836905007545

Zimmermann, F., Breitenmoser-Würsten, C., Molinari-Jobin, A., & Breitenmoser, U. (2013). 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. *Integrative Zoology*, 8, 232-243. https://doi.org/10.1111/1749-4877.12017

Zimmermann, F., & Foresti, D. (2016). Capture-recapture methods for density estimation. In F. Rovero & F. Zimmermann (Eds.), *Camera trapping for wildlife research* (pp. 95-141). Exeter, UK: Pelagic Publishing.

**How to cite this article:** Gimenez O, Gatti S, Duchamp C, et al. Spatial density estimates of Eurasian lynx (Lynx lynx) in the French Jura and Vosges Mountains. *Ecol Evol*. 2019;9:11707-11715. https://doi.org/10.1002/ece3.5668