Sampling variables and their thresholds for the precise estimation of wild felid population density with camera traps and spatial capture–recapture methods

1. Robust monitoring, providing information on population status, is fundamental for successful conservation planning. However, this can be hard to achieve for species that are elusive and occur at low densities, such as felids. These are often keystones of functioning ecosystems and are threatened by habitat loss and human persecution. 2. When elusive species can be individually identified by visible characteristics, for example via camera-trapping, observations of individuals can be used in combination with capture–recapture methods to calculate demographic parameters such as population density. In this context, spatial capture–recapture (SCR) outperforms conventional non-spatial methods, but the precision of results is inherently related to the sampling design, which should therefore be optimised. 3. We focussed on territorial felids in different habitats and investigated how the sampling designs implemented in the field affected the precision of population density estimates. We examined 137 studies that combined camera trapping and SCR methods for density estimation. From these, we collectedspatiotemporal parameters of their sampling designs, monitoring results, such as the number of individuals captured and the number of recaptures, as well as SCR detection parameters. We applied generalised linear mixed-effects models and tree-based regression methods to investigate the influence of variables on the precision of population density estimates and provide numerical thresholds. 4. Our analysis shows that the number of individuals, recapture frequency, and capture probability play the most crucial roles. Surveys yielding over 20 captured individuals that were recaptured on average at least once obtain the most precise population density estimates. 5. Based on our findings, we provide practical guidelines for future SCR studies that apply to all territorial felids. Furthermore, we present a standardised reporting protocol for study transparency and comparability. Our results will improve reporting and reproducibility of SCR studies and aid in setting up optimised sampling designs.publishedVersio

Human disturbance is the most limiting factor driving habitat selection of a large carnivore throughout Continental Europe

Habitat selection is a multi-scale process driven by trade-offs between benefits, such as resource abundance, and disadvantages, such as the avoidance of risk. The latter includes human disturbances, to which large carnivores, with their large spatial requirements, are especially sensitive. We investigated the ecological processes underlying multi-scale habitat selection of a large carnivore, namely Eurasian lynx, across European landscapes characterized by different levels of human modification. Using a unique dataset of 125 lynx from 9 study sites across Europe, we compared used and available locations within landscape and home-range scales using a novel Mixed Effect randomForest approach, while considering environmental predictors as proxies for human disturbances and environmental resources. At the landscape scale, lynx avoided roads and human settlements, while at the home-range scale natural landscape features associated with shelter and prey abundance were more important. The results showed sex was of relatively low variable importance for lynx's general habitat selection behaviour. We found increasingly homogeneous responses across study sites with finer selection scales, suggesting that study site differences determined coarse selection, while utilization of resources at the finer selection scale was broadly universal. Thereby describing lynx's requirement, if not preference, for heterogeneous forests and shelter from human disturbances and implying that regional differences in coarse-scale selection are driven by availability rather than preference. These results provide crucial information for conserving this species in human-dominated landscapes, as well as for the first time, to our knowledge, generalising habitat selection behaviour of a large carnivore species at a continental scale.acceptedVersio

A common statement on anthropogenic hybridization of the European wildcat (Felis silvestris)

Preserving natural genetic diversity and ecological function of wild species is a central goal in conservation biology. As such, anthropogenic hybridization is considered a threat to wild populations, as it can lead to changes in the genetic makeup of wild species and even to the extinction of wild genomes. In European wildcats, the genetic and ecological impacts of gene flow from domestic cats are mostly unknown at the species scale. However, in small and isolated populations, it is known to include genetic swamping of wild genomes. In this context, it is crucial to better understand the dynamics of hybridization across the species range, to inform and implement management measures that maintain the genetic diversity and integrity of the European wildcat. In the present paper, we aim to provide an overview of the current scientific understanding of anthropogenic hybridization in European wildcats, to clarify important aspects regarding the evaluation of hybridization given the available methodologies, and to propose guidelines for management and research priorities

In situ feeding as a new management tool to conserve orphaned Eurasian lynx (lynx lynx)

High human-caused mortality due to wildlife-vehicle-collisions and illegal killing leads to frequent cases of orphaned Eurasian lynx juveniles. Under natural conditions, this would result in starvation of the young. To avoid this, wildlife managers conventionally rear animals in captivity and release them later. However, this measure is an undesirable outcome for species conservation, managers, and animals alike. Increased awareness of Eurasian lynx orphaned by human-caused mortality means managers must often intervene in endangered populations. In this study, we report for the first time a successful case of in situ feeding designed to avoid captivity of two orphaned Eurasian lynx. We exposed 13 roe deer and 7 red deer carcasses in the field to successfully support two orphans to the age of independence and confirm dispersal from the natal range. We present this management approach as a feasible and complimentary tool that can be considered in small or isolated large carnivore populations where every individual counts toward population viability

The boon and bane of boldness: movement syndrome as saviour and sink for population genetic diversity

BACKGROUND: Many felid species are of high conservation concern, and with increasing human disturbance the situation is worsening. Small isolated populations are at risk of genetic impoverishment decreasing within-species biodiversity. Movement is known to be a key behavioural trait that shapes both demographic and genetic dynamics and affects population survival. However, we have limited knowledge on how different manifestations of movement behaviour translate to population processes. In this study, we aimed to 1) understand the potential effects of movement behaviour on the genetic diversity of small felid populations in heterogeneous landscapes, while 2) presenting a simulation tool that can help inform conservation practitioners following, or considering, population management actions targeting the risk of genetic impoverishment. METHODS: We developed a spatially explicit individual-based population model including neutral genetic markers for felids and applied this to the example of Eurasian lynx. Using a neutral landscape approach, we simulated reintroductions into a three-patch system, comprising two breeding patches separated by a larger patch of differing landscape heterogeneity, and tested for the effects of various behavioural movement syndromes and founder population sizes. We explored a range of movement syndromes by simulating populations with various movement model parametrisations that range from ‘shy’ to ‘bold’ movement behaviour. RESULTS: We find that movement syndromes can lead to a higher loss of genetic diversity and an increase in between population genetic structure for both “bold” and “shy” movement behaviours, depending on landscape conditions, with larger decreases in genetic diversity and larger increases in genetic differentiation associated with bold movement syndromes, where the first colonisers quickly reproduce and subsequently dominate the gene pool. In addition, we underline the fact that a larger founder population can offset the genetic losses associated with subpopulation isolation and gene pool dominance. CONCLUSIONS: We identified a movement syndrome trade-off for population genetic variation, whereby bold-explorers could be saviours - by connecting populations and promoting panmixia, or sinks - by increasing genetic losses via a ‘founder takes all’ effect, whereas shy-stayers maintain a more gradual genetic drift due to their more cautious behaviour. Simulations should incorporate movement behaviour to provide better projections of long-term population viability and within-species biodiversity, which includes genetic diversity. Simulations incorporating demographics and genetics have great potential for informing conservation management actions, such as population reintroductions or reinforcements. Here, we present such a simulation tool for solitary felids

A voxel matching method for effective leaf area index estimation in temperate deciduous forests from leaf-on and leaf-off airborne LiDAR data

The quantification of leaf area index (LAI) is essential for modeling the interaction between atmosphere and biosphere. The airborne LiDAR has emerged as an effective tool for mapping plant area index (PAI) in a landscape consisting of both woody and leaf materials. However, the discrimination between woody and leaf materials and the estimation of effective LAI (eLAI) have, to date, rarely been studied at landscape scale. We applied a voxel matching algorithm to estimate eLAI of deciduous forests using simulated and field LiDAR data under leaf-on and leaf-off conditions. We classified LiDAR points as either a leaf or a woody hit on leaf-on LiDAR data by matching the point with leaf-off data. We compared the eLAI result of our voxel matching algorithm against the subtraction method, where the leaf-off effective woody area index (eWAI) is subtracted from the effective leaf-on PAI (ePAI). Our results, which were validated against terrestrial LiDAR derived eLAI, showed that the voxel matching method, with an optimal voxel size of 0.1 m, produced an unbiased estimation of terrestrial LiDAR derived eLAI with an R2 of 0.70 and an RMSE of 0.41 (RRMSE: 20.1%). The subtraction method, however, yielded an R2 of 0.62 and an RMSE of 1.02 (RRMSE: 50.1%) with a significant underestimation of 0.94. Reassuringly, the same outcome was observed using a simulated dataset. In addition, we evaluated the performance of 96 LiDAR metrics under leaf-on conditions for eLAI prediction using a statistical model. Based on the importance scores derived from the random forest regression, nine of the 96 leaf-on LiDAR metrics were selected. Cross-validation showed that eLAI could be predicted using these metrics under leaf-on conditions with an R2 of 0.73 and an RMSE of 0.27 (RRMSE: 17.4%). The voxel matching method yielded a slightly lower accuracy (R2: 0.70, RMSE:0.41, RRMSE: 20.1%) than the statistical model. We, therefore, suggest that the voxel matching method offers a new opportunity for the estimating eLAI and other ecological applications that require the classification between leaf and woody materials using airborne LiDAR data. It potentially allows transferability to different sites and flight campaigns

Population and landscape genetic analysis of the Malayan sun bear Helarctos malayanus

Conservation genetics can provide data needed by conservation practitioners for their decisions regarding the management of vulnerable or endangered species, such as the sun bear Helarctos malayanus. Throughout its range, the sun bear is threatened by loss and fragmentation of its habitat and the illegal trade of both live bears and bear parts. Sharply declining population numbers and population sizes, and a lack of natural dispersal between populations all threaten the genetic diversity of the remaining populations of this species. In this first population genetics study of sun bears using microsatellite markers, we analyzed 68 sun bear samples from Cambodia to investigate population structure and genetic diversity. We found evidence for two genetically distinct populations in the West and East of Cambodia. Ongoing or recent gene flow between these populations does not appear sufficient to alleviate loss of diversity in these populations, one of which (West Cambodia) is characterized by significant inbreeding. We were able to assign 85% of sun bears of unknown origin to one of the two populations with high confidence (assignment probability ≥ 85%), providing valuable information for future bear reintroduction programs. Further, our results suggest that developed land (mostly agricultural mosaics) acts as a barrier to gene flow for sun bears in Cambodia. We highlight that regional sun bear conservation action plans should consider promoting population connectivity and enforcing wildlife protection of this threatened species

Measuring leaf angle distribution using terrestrial laser scanning in a European Beech forest

Leaf angle distribution (LAD) is an important canopy structure metric. It controls the flux of radiation, carbon and water, and has therefore been used in many radiative transfer, meteorological and hydrological models. However, LAD is too tedious to measure using conventional manual methods. Terrestrial laser scanning (TLS) has recently been proposed to estimate LAD due to its ability to record unprecedented detailed plant 3D structure. However, previous research was restricted to a controlled environment with simple canopy structure. In this research, TLS was used in a natural deciduous European beech forest to estimate LAD. Digital hemispherical photograph (DHP) was also used as a reference. The results demonstrated that both TLS and DHP could capture a variation of LAD in beech plots at different succession stages. Compared to DHP, TLS has the advantage of resolving foliar and woody materials, as well as deriving the 3D distribution of leaf angles

Comparison of terrestrial LiDAR and digital hemispherical photography for estimating leaf angle distribution in European broadleaf beech forests

Leaf inclination plays a crucial role in regulating the radiation, carbon and water fluxes in plant canopies. Accurate measurement of the probability density function of leaf inclination (i.e., the leaf angle distribution or leaf inclination distribution function (LIDF)), is very important for modelling photosynthesis as well as measuring leaf area index. In spite of its importance, in situ measurement of LIDF is very challenging. Both digital hemispherical photography (DHP) and terrestrial LiDAR (TLS) have been used to measure LIDF. However, the consistency and relative accuracy of these two techniques has never been evaluated. In this research, we aimed to evaluate which in situ technique, either DHP or TLS, could measure LIDF more accurately, with respect to both field-based and synthetic datasets. The field-based datasets were collected from 36 natural European beech stands covering a range of forest structures. The synthetic datasets were generated from 44 virtual forest scenes using TLS and DHP simulators. Due to the inability of differing leaf and woody materials in DHP, the average plant inclination angle ( ) from DHP and TLS was selected for LIDF comparison. The average inclination angle ( ) was retrieved from TLS point clouds using a geometric method, and from DHP using three gap fraction inversion methods including the NC method (Norman and Campbell, 1989), as well as the CAN-EYE and Hemisfer software. Results from the field-based datasets showed a significant difference and inconsistency between the average inclination angle ( ) retrieved from TLS and DHP (respectively (40°, 58°) from TLS, (15°, 86°) from DHP NC, (36°, 78°) from DHP CAN-EYE, (0°, 67°) from DHP Hemisfer). Results from the synthetic datasets demonstrated that the accuracy of from TLS was considerably higher than that obtained from DHP (R2: 0.90 > 0.74; RMSE: 5.38° < 13.30°). This study demonstrated that the LIDF estimated from TLS and DHP were not coherent. Based on the experimental results as well as deduction from theoretical arguments, we recommended using TLS when measuring leaf inclination in broadleaf forests, especially for stands with a heterogeneous structure

Improving leaf area index (LAI) estimation by correcting for clumping and woody effects using terrestrial laser scanning

Leaf area index (LAI) has frequently been measured in the field using traditional optical methods such as digital hemispherical photography (DHP). However, in the DHP retrieved LAI, there is always contribution of woody components due to the difficulty in distinguishing woody and foliar materials. In addition, the leaf angle distribution which strongly affects the estimation of LAI is either ignored while using the convergent angle 57.5°, or inversed simultaneously with LAI using multiple directions. Terrestrial laser scanning (TLS) provides a 3-dimensional view of the forest canopy, which we used in this study to improve LAI estimation by directly retrieving leaf angle distribution, and subsequently correcting foliage clumping and woody effects. The leaf angle distribution was retrieved by estimating the angle between the leaf normal vectors and the zenith vectors. The clumping index was obtained by using the gap size distribution method, while the woody contribution was evaluated based on an improved point classification between woody and foliar materials. Finally, the gap fraction derived from TLS was converted to effective LAI, and thence to LAI. The study was conducted for 31 forest plots including deciduous, coniferous and mixed plots in Bavarian Forest National Park. The classification accuracy was improved by approximately 10% using our method. Results showed that the clumping caused an underestimation of LAI ranging from 1.2% to 48.0%, while woody contribution led to an overestimation from 3.0% to 31.9% compared to the improved LAI. The combined error ranged from −46.2% to 32.6% of the leaf area index (LAI) measurements. The error was largely dependent on forest types. The clumping index of coniferous plots on average was lower than that of deciduous plots, whereas deciduous plots had a higher woody-to-total area ratio. The proposed method provides a more accurate estimate of LAI by eliminating clumping and woody effects, as well as the effect of leaf angle distribution