multimark: an R package for analysis of capture–recapture data consisting of multiple “noninvasive” marks

I describe an open‐source R package, multimark , for estimation of survival and abundance from capture–mark–recapture data consisting of multiple “noninvasive” marks. Noninvasive marks include natural pelt or skin patterns, scars, and genetic markers that enable individual identification in lieu of physical capture. multimark provides a means for combining and jointly analyzing encounter histories from multiple noninvasive sources that otherwise cannot be reliably matched (e.g., left‐ and right‐sided photographs of bilaterally asymmetrical individuals). The package is currently capable of fitting open population Cormack–Jolly–Seber (CJS) and closed population abundance models with up to two mark types using Bayesian Markov chain Monte Carlo (MCMC) methods. multimark can also be used for Bayesian analyses of conventional capture–recapture data consisting of a single‐mark type. Some package features include (1) general model specification using formulas already familiar to most R users, (2) ability to include temporal, behavioral, age, cohort, and individual heterogeneity effects in detection and survival probabilities, (3) improved MCMC algorithm that is computationally faster and more efficient than previously proposed methods, (4) Bayesian multimodel inference using reversible jump MCMC, and (5) data simulation capabilities for power analyses and assessing model performance. I demonstrate use of multimark using left‐ and right‐sided encounter histories for bobcats (Lynx rufus) collected from remote single‐camera stations in southern California. In this example, there is evidence of a behavioral effect (i.e., trap “happy” response) that is otherwise indiscernible using conventional single‐sided analyses. The package will be most useful to ecologists seeking stronger inferences by combining different sources of mark–recapture data that are difficult (or impossible) to reliably reconcile, particularly with the sparse datasets typical of rare or elusive species for which noninvasive sampling techniques are most commonly employed. Addressing deficiencies in currently available software, multimark also provides a user‐friendly interface for performing Bayesian multimodel inference using capture–recapture data consisting of a single conventional mark or multiple noninvasive marks.     

Comparison of photo‐matching algorithms commonly used for photographic capture–recapture studies

Photographic capture–recapture is a valuable tool for obtaining demographic information on wildlife populations due to its noninvasive nature and cost‐effectiveness. Recently, several computer‐aided photo‐matching algorithms have been developed to more efficiently match images of unique individuals in databases with thousands of images. However, the identification accuracy of these algorithms can severely bias estimates of vital rates and population size. Therefore, it is important to understand the performance and limitations of state‐of‐the‐art photo‐matching algorithms prior to implementation in capture–recapture studies involving possibly thousands of images. Here, we compared the performance of four photo‐matching algorithms; Wild‐ID, I3S Pattern+, APHIS, and AmphIdent using multiple amphibian databases of varying image quality. We measured the performance of each algorithm and evaluated the performance in relation to database size and the number of matching images in the database. We found that algorithm performance differed greatly by algorithm and image database, with recognition rates ranging from 100% to 22.6% when limiting the review to the 10 highest ranking images. We found that recognition rate degraded marginally with increased database size and could be improved considerably with a higher number of matching images in the database. In our study, the pixel‐based algorithm of AmphIdent exhibited superior recognition rates compared to the other approaches. We recommend carefully evaluating algorithm performance prior to using it to match a complete database. By choosing a suitable matching algorithm, databases of sizes that are unfeasible to match “by eye” can be easily translated to accurate individual capture histories necessary for robust demographic estimates.     

Validation of photo‐identification as a mark–recapture method in the spotted eagle ray Aetobatus narinari

The spotted eagle ray Aetobatus narinari is characterized by pigmentation patterns that are retained for up to 3·5 years. These pigmentations can be used to identify individuals through photo‐identification. Only one study has validated this technique, but no study has estimated the percentage of correct identification of the rays using this technique. In order to carry out demographic research, a reliable photographic identification technique is needed. To achieve this validation for A. narinari, a double‐mark system was established over 11 months and photographs of the dorsal surface of 191 rays were taken. Three body parts with distinctive natural patterns were analysed (dorsal surface of the cephalic region, dorsal surface of the pectoral fins and dorsal surface of the pelvic fins) in order to determine the body part that could be used to give the highest percentage of correct identification. The dorsal surface of the pectoral fins of A. narinari provides the most accurate photo‐identification to distinguish individuals (88·2%).     

Genotyping validates photo‐identification by the head scale pattern in a large population of the European adder (Vipera berus)

Capture‐mark‐recapture procedures are a basic tool in population studies and require that individual animals are correctly identified throughout their lifetime. A method that has become more and more popular uses photographic records of natural markings, such as pigmentation pattern and scalation configuration. As with any other marking tool, the validity of the photographic identification technique should be evaluated thoroughly. Here, we report on a large‐scale double‐marking study in which European adders (Vipera berus) were identified by both microsatellite genetic markers and by the pattern of head scalation. Samples that were successfully genotyped for all nine loci yielded 624 unique genotypes, which matched on a one‐to‐one basis with the individual assignments based on the head scalation pattern. Thus, adders considered as different individuals by their genotypes were also identified as different individuals by their head scalation pattern, and vice versa. Overall, variation in the numbers, shape, and arrangement of the head scales enabled us to distinguish among 3200+ photographed individual snakes. Adders that were repeatedly sequenced genetically over intervals of 2–3 years showed no indication whatsoever for a change in the head scale pattern. Photographic records of 900+ adders that were recaptured over periods of up to 12 years showed a very detailed and precise match of the head scale characteristics. These natural marks are thus robust over time and do not change during an individual's lifetime. With very low frequency (0.3%), we detected small changes in scalation that were readily discernible and could be attributed to physical injury or infection. Our study provides a conclusive validation for the use of photo‐identification by head scale patterns in the European adder.     

Applying computer‐aided photo‐identification to messy datasets: a case study of Thornicroft's giraffe (Giraffa camelopardalis thornicrofti)

Digital photography enables researchers to rapidly compile large quantities of data from individually identifiable animals, and computer software improves the management of such large datasets while aiding the identification process. Wild‐ID software has performed well with uniform datasets controlling for angle and portion of the animal photographed; however, few datasets are collected under such controlled conditions. We examined the effectiveness of Wild‐ID in identifying individual Thornicroft's giraffe from a dataset of photographs (n = 552) collected opportunistically in the Luangwa Valley, Zambia from March to October 2009. We assessed the programme's accuracy in correctly identifying individuals and the effect of five image quality factors on identification success: blurriness, background type and complexity, amount of sky and the presence of other giraffe. The programme correctly identified individuals in 71.6% of photographs. Background complexity was the only significant variable affecting identification success and removing background imagery reduced identification error by 52.8% (from 28.4 to 13.4%). Our results indicate higher levels of error than previously reported for Wild‐ID. However, they also suggest the programme is an effective tool for quickly identifying individuals in large field datasets, especially if photograph backgrounds are removed beforehand and postanalysis visual verification is performed.     

A newt does not change its spots: using pattern mapping for the identification of individuals in large populations of newt species

The correct identification of individuals is a requirement of capture-mark-recapture (CMR) methods, and it is commonly achieved by applying artificial marks or by mutilation of study-animals. An alternative, non-invasive method to identify individuals is to utilize the patterns of their natural body markings. However, the use of pattern mapping is not yet widespread, mainly because it is considered time consuming, particularly in large populations and/or long-term CMR studies. Here we explore the use of pattern mapping for the identification of adult individuals in the alpine (Ichthyosaura alpestris) and smooth (Lissotriton vulgaris) newts (Amphibia, Salamandridae), using the freely available, open-source software Wild-ID. Our photographic datasets comprised nearly 4000 captured animals’ images, taken during a 3-year period. The spot patterns of individual newts of both species did not change through time, and were sufficiently varied to allow their individual identification, even in the larger datasets. The pattern-recognition algorithm of Wild-ID was highly successful in identifying individual newts in both species. Our findings indicate that pattern mapping can be successfully employed for the identification of individuals in large populations of a broad range of animals that exhibit natural markings. The significance of pattern-mapping is accentuated in CMR studies that aim in obtaining long-term information on the demography and population dynamics of species of conservation interest, such as many amphibians facing population declines.     

Genotyping validates the efficacy of photographic identification in a capture‐mark‐recapture study based on the head scale patterns of the prairie lizard (Sceloporus consobrinus)

Population studies often incorporate capture‐mark‐recapture (CMR) techniques to gather information on long‐term biological and demographic characteristics. A fundamental requirement for CMR studies is that an individual must be uniquely and permanently marked to ensure reliable reidentification throughout its lifespan. Photographic identification involving automated photographic identification software has become a popular and efficient noninvasive method for identifying individuals based on natural markings. However, few studies have (a) robustly assessed the performance of automated programs by using a double‐marking system or (b) determined their efficacy for long‐term studies by incorporating multi‐year data. Here, we evaluated the performance of the program Interactive Individual Identification System (I³S) by cross‐validating photographic identifications based on the head scale pattern of the prairie lizard (Sceloporus consobrinus) with individual microsatellite genotyping (N = 863). Further, we assessed the efficacy of the program to identify individuals over time by comparing error rates between within‐year and between‐year recaptures. Recaptured lizards were correctly identified by I³S in 94.1% of cases. We estimated a false rejection rate (FRR) of 5.9% and a false acceptance rate (FAR) of 0%. By using I³S, we correctly identified 97.8% of within‐year recaptures (FRR = 2.2%; FAR = 0%) and 91.1% of between‐year recaptures (FRR = 8.9%; FAR = 0%). Misidentifications were primarily due to poor photograph quality (N = 4). However, two misidentifications were caused by indistinct scale configuration due to scale damage (N = 1) and ontogenetic changes in head scalation between capture events (N = 1). We conclude that automated photographic identification based on head scale patterns is a reliable and accurate method for identifying individuals over time. Because many lizard or reptilian species possess variable head squamation, this method has potential for successful application in many species.     

The use of a non‐invasive tool for capture–recapture studies on a seahorse Hippocampus guttulatus population

In this study, the spot pattern in Hippocampus guttulatus was analysed using a computer programme algorithm that allowed individual comparison. This methodology was first tested in a controlled environment using 51 adult and 55 juvenile H. guttulatus. Positive matches were obtained in 86·3 and 83·6% of the adults and juveniles, respectively. In a second experiment, monthly surveys were carried out in five selected locations in the Ria Formosa Lagoon, south Portugal, over the course of a year and a total of 980 photographs were analysed. Photographed H. guttulatus were re‐sighted one to nine times during the course of the survey period with an overall re‐sight record of over 30%. Photo‐identification was therefore shown to be a useful tool for non‐invasive mark–recapture studies that can be successfully used to survey the population abundance of H. guttulatus aged 6 months or older in consecutive years. This could be of great value when considering the assessment of H. guttulatus populations and understanding changes over time.     

How common is common? Rapidly assessing population size and structure of the frog Mantidactylus betsileanus at a site in east‐central Madagascar

Population‐level data are urgently needed for amphibians in light of the ongoing amphibian extinction crisis. Studies focused on population dynamics are not only important for rare species but also for common species which shape ecosystems to a greater degree than those that are rare. Some of the greatest global amphibian species diversity is found in Madagascar, yet there are few studies on the ecology of frog species on the island. We carried out a mark‐recapture study on the widespread frog Mantidactylus betsileanus (Mantellidae: Mantellinae: Mantidactylus) at two adjacent rainforest sites in east‐central Madagascar to assess its population size and structure. To do so, we validated and implemented an individual identification protocol using photographs of the ventral patterns of frogs and identified individuals with photographic‐matching software. Using this rapid, non‐invasive survey method, we were able to estimate a density of 26 and 28 frogs per 100 m² at each of the two sites sampled. Our results show the rainforests near the village of Andasibe, Madagascar support remarkably high amphibian abundance, helping illustrate the significant ecological role of frogs in this ecosystem. Further, individual frog markings allowed us to develop more precise estimates than traditional survey methods. This study provides a blueprint to augment existing population studies or develop new monitoring programs in Madagascar and beyond.     

Balancing bias and precision in capture-recapture estimates of abundance

Capture‐recapture estimates of abundance using photographic identification data are sensitive to the quality of photographs used and distinctiveness of individuals in the population. Here analyses are presented for examining the effects of photographic quality and individual animal distinctiveness scores and for objectively selecting a subset of data to use for capture‐recapture analyses using humpback whale (Megaptera novaeangliae) data from a 2‐year study in the North Atlantic. Photographs were evaluated for their level of quality and whales for their level of individual distinctiveness. Photographic quality scores had a 0.21 probability of changing by a single‐quality level, and there were no changes by two or more levels. Individual distinctiveness scores were not independent of photographic quality scores. Estimates of abundance decreased as poor‐quality photographs were removed. An appropriate balance between precision and bias in abundance estimates was achieved by removing the lowest‐quality photographs and those of incompletely photographed flukes given our assumptions about the true population abundance. A simulation of the selection process implied that, if the estimates are negatively biased by heterogeneity, the increase in bias produced by decreasing the sample size is not more than 2%. Capture frequencies were independent of individual distinctiveness scores.     

Mark–recapture cloning: a straightforward and cost‐effective cloning method for population genetics of single‐copy nuclear DNA sequences in diploids

We describe a simple protocol to reduce the number of cloning reactions of nuclear DNA sequences in population genetic studies of diploid organisms. Cloning is a necessary step to obtain correct haplotypes in such organisms, and, while traditional methods are efficient at cloning together many genes of a single individual, population geneticists rather need to clone the same locus in many individuals. Our method consists of marking individual sequences during the polymerase chain reaction (PCR) using 5′‐tailed primers with small polynucleotide tags. PCR products are mixed together before the cloning reaction and clones are sequenced with universal plasmid primers. The individual from which a sequence comes from is identified by the tag sequences upstream of each initial primer. We called our protocol mark–recapture (MR) cloning. We present results from 57 experiments of MR cloning conducted in four distinct laboratories using nuclear loci of various lengths in different invertebrate species. Rate of capture (proportion of individuals for which one or more sequences were retrieved) and multiple capture (proportion of individuals for which two or more sequences were retrieved) empirically obtained are described. We estimated that MR cloning allowed reducing costs by up to 70% when compared to conventional individual‐based cloning. However, we recommend to adjust the mark:recapture ratio in order to obtain multiple sequences from the same individual and circumvent inherent technical artefacts of PCR, cloning and sequencing. We argue that MR cloning is a valid and reliable high‐throughput method, providing the number of sequences exceeds the number of individuals initially amplified.     

Population size estimates based on the frequency of genetically assigned parent–offspring pairs within a subsample

Estimating population density as precise as possible is a key premise for managing wild animal species. This can be a challenging task if the species in question is elusive or, due to high quantities, hard to count. We present a new, mathematically derived estimator for population size, where the estimation is based solely on the frequency of genetically assigned parent–offspring pairs within a subsample of an ungulate population. By use of molecular markers like microsatellites, the number of these parent–offspring pairs can be determined. The study's aim was to clarify whether a classical capture–mark–recapture (CMR) method can be adapted or extended by this genetic element to a genetic‐based capture–mark–recapture (g‐CMR). We numerically validate the presented estimator (and corresponding variance estimates) and provide the R‐code for the computation of estimates of population size including confidence intervals. The presented method provides a new framework to precisely estimate population size based on the genetic analysis of a one‐time subsample. This is especially of value where traditional CMR methods or other DNA‐based (fecal or hair) capture–recapture methods fail or are too difficult to apply. The DNA source used is basically irrelevant, but in the present case the sampling of an annual hunting bag is to serve as data basis. In addition to the high quality of muscle tissue samples, hunting bags provide additional and essential information for wildlife management practices, such as age, weight, or sex. In cases where a g‐CMR method is ecologically and hunting‐wise appropriate, it enables a wide applicability, also through its species‐independent use.     

Field‐readable alphanumeric flags are valuable markers for shorebirds: use of double‐marking to identify cases of misidentification

Implicit assumptions for most mark‐recapture studies are that individuals do not lose their markers and all observed markers are correctly recorded. If these assumptions are violated, e.g., due to loss or extreme wear of markers, estimates of population size and vital rates will be biased. Double‐marking experiments have been widely used to estimate rates of marker loss and adjust for associated bias, and we extended this approach to estimate rates of recording errors. We double‐marked 309 Piping Plovers (Charadrius melodus) with unique combinations of color bands and alphanumeric flags and used multi‐state mark recapture models to estimate the frequency with which plovers were misidentified. Observers were twice as likely to read and report an invalid color‐band combination (2.4% of the time) as an invalid alphanumeric code (1.0%). Observers failed to read matching band combinations or alphanumeric flag codes 4.5% of the time. Unlike previous band resighting studies, use of two resightable markers allowed us to identify when resighting errors resulted in reports of combinations or codes that were valid, but still incorrect; our results suggest this may be a largely unappreciated problem in mark‐resight studies. Field‐readable alphanumeric flags offer a promising auxiliary marker for identifying and potentially adjusting for false‐positive resighting errors that may otherwise bias demographic estimates.     

Capture–recapture models with heterogeneity to study survival senescence in the wild

Detecting senescence in wild populations and estimating its strength raise three challenges. First, in the presence of individual heterogeneity in survival probability, the proportion of high‐survival individuals increases with age. This increase can mask a senescence‐related decrease in survival probability when the probability is estimated at the population level. To accommodate individual heterogeneity we use a mixture model structure (discrete classes of individuals). Second, the study individuals can elude the observers in the field, and their detection rate can be heterogeneous. To account for detectability issues we use capture–mark–recapture (CMR) methodology, mixture models and data that provide information on individuals’ detectability. Last, emigration to non‐monitored sites can bias survival estimates, because it can occur at the end of the individuals’ histories and mimic earlier death. To model emigration we use Markovian transitions to and from an unobservable state. These different model structures are merged together using hidden Markov chain CMR models, or multievent models. Simulation studies illustrate that reliable evidence for survival senescence can be obtained using highly heterogeneous data from non site‐faithful individuals. We then design a tailored application for a dataset from a colony of black‐headed gull Chroicocephalus ridibundus. Survival probabilities do not appear individually variable, but evidence for survival senescence becomes significant only when accounting for other sources of heterogeneity. This result suggests that not accounting for heterogeneity leads to flawed inference and/or that emigration heterogeneity mimics survival heterogeneity and biases senescence estimates.     

Evaluating the interaction of faecal pellet deposition rates and DNA degradation rates to optimize sampling design for DNA‐based mark–recapture analysis of Sonoran pronghorn

Knowledge of population demographics is important for species management but can be challenging in low‐density, wide‐ranging species. Population monitoring of the endangered Sonoran pronghorn (Antilocapra americana sonoriensis) is critical for assessing the success of recovery efforts, and noninvasive DNA sampling (NDS) could be more cost‐effective and less intrusive than traditional methods. We evaluated faecal pellet deposition rates and faecal DNA degradation rates to maximize sampling efficiency for DNA‐based mark–recapture analyses. Deposition data were collected at five watering holes using sampling intervals of 1–7 days and averaged one pellet pile per pronghorn per day. To evaluate nuclear DNA (nDNA) degradation, 20 faecal samples were exposed to local environmental conditions and sampled at eight time points from one to 124 days. Average amplification success rates for six nDNA microsatellite loci were 81% for samples on day one, 63% by day seven, 2% by day 14 and 0% by day 60. We evaluated the efficiency of different sampling intervals (1–10 days) by estimating the number of successful samples, success rate of individual identification and laboratory costs per successful sample. Cost per successful sample increased and success and efficiency declined as the sampling interval increased. Results indicate NDS of faecal pellets is a feasible method for individual identification, population estimation and demographic monitoring of Sonoran pronghorn. We recommend collecting samples >7 days old and estimate that a sampling interval of 4–7 days in summer conditions (i.e. extreme heat and exposure to UV light) will achieve desired sample sizes for mark–recapture analysis while also maximizing efficiency.     

Combining noninvasive genetics and a new mammalian sex‐linked marker provides new tools to investigate population size,structure and individual behaviour: An application to bats

Monitoring wild populations is crucial for their effective management. Noninvasive genetic methods provide robust data from individual free‐ranging animals, which can be used in capture–mark–recapture (CMR) models to estimate demographic parameters without capturing or disturbing them. However, sex‐ and status‐specific behaviour, which may lead to differences in detection probabilities, is rarely considered in monitoring. Here, we investigated population size, sex ratio, sex‐ and status‐related behaviour in 19 Rhinolophus hipposideros maternity colonies (Northern France) with a noninvasive genetic CMR approach (using faeces) combined with parentage assignments. The use of the DDX3X/Y‐Mam sexual marker designed in this study, which shows inter‐ and intrachromosomal length polymorphism across placental mammals, together with eight polymorphic microsatellite markers, produced high‐quality genetic data with limited genotyping errors and allowed us to reliably distinguish different categories of individuals (males, reproductive and nonreproductive females) and to estimate population sizes. We showed that visual counts represent well‐adult female numbers and that population composition in maternity colonies changes dynamically during the summer. Before parturition, colonies mainly harbour pregnant and nonpregnant females with a few visiting males, whereas after parturition, colonies are mainly composed of mothers and their offspring with a few visiting nonmothers and males. Our approach gives deeper insight into sex‐ and status‐specific behaviour, a prerequisite for understanding population dynamics and developing effective monitoring and management strategies. Provided sufficient samples can be obtained, this approach can be readily applied to a wide range of species.     

Combining capture–recapture data and known ages allows estimation of age‐dependent survival rates

In many animal populations, demographic parameters such as survival and recruitment vary markedly with age, as do parameters related to sampling, such as capture probability. Failing to account for such variation can result in biased estimates of population‐level rates. However, estimating age‐dependent survival rates can be challenging because ages of individuals are rarely known unless tagging is done at birth. For many species, it is possible to infer age based on size. In capture–recapture studies of such species, it is possible to use a growth model to infer the age at first capture of individuals. We show how to build estimates of age‐dependent survival into a capture–mark–recapture model based on data obtained in a capture–recapture study. We first show how estimates of age based on length increments closely match those based on definitive aging methods. In simulated analyses, we show that both individual ages and age‐dependent survival rates estimated from simulated data closely match true values. With our approach, we are able to estimate the age‐specific apparent survival rates of Murray and trout cod in the Murray River, Australia. Our model structure provides a flexible framework within which to investigate various aspects of how survival varies with age and will have extensions within a wide range of ecological studies of animals where age can be estimated based on size.     

Computer-Aided Pattern Recognition of Large Reptiles as a Noninvasive Application to Identify Individuals

For large species, the capture and handling of individuals in capture–mark–recapture studies introduces nonhuman animal welfare issues associated with handling, physical marking, and possible wounding due to tag loss. The use of photographic identification for these species offers an alternative and less invasive marking technique. This study investigated the opportunity offered by photo identification to individually mark individuals of a large reptile, the perentie (Varanus giganteus), in Australia and therefore avoid the stress of physically capturing and handling. Photographs submitted by a remotely located community were first validated to confirm whether perenties could be individually identified from their spots electronically using a computer program. Computer-aided selection of unique patterns was found to be appropriate for the identification of individuals and confirmed 38 individuals during the sampling period. The value of this approach is 2-fold: There is a benefit to animal welfare in that handling an animal is not required to capture him or her, thus reducing capture-related stress; and confirmation that photo identification of distinctive patterns of the perentie is valid and offers a useful option to identify individuals of this large species.     

Assessing the performance of index calibration survey methods to monitor populations of wide‐ranging low‐density carnivores

Apex carnivores are wide‐ranging, low‐density, hard to detect, and declining throughout most of their range, making population monitoring both critical and challenging. Rapid and inexpensive index calibration survey (ICS) methods have been developed to monitor large African carnivores. ICS methods assume constant detection probability and a predictable relationship between the index and the actual population of interest. The precision and utility of the resulting estimates from ICS methods have been questioned. We assessed the performance of one ICS method for large carnivores—track counts—with data from two long‐term studies of African lion populations. We conducted Monte Carlo simulation of intersections between transects (road segments) and lion movement paths (from GPS collar data) at varying survey intensities. Then, using the track count method we estimated population size and its confidence limits. We found that estimates either overstate precision or are too imprecise to be meaningful. Overstated precision stemmed from discarding the variance from population estimates when developing the method and from treating the conversion from tracks counts to population density as a back‐transformation, rather than applying the equation for the variance of a linear function. To effectively assess the status of species, the IUCN has set guidelines, and these should be integrated in survey designs. We propose reporting the half relative confidence interval width (HRCIW) as an easily calculable and interpretable measure of precision. We show that track counts do not adhere to IUCN criteria, and we argue that ICS methods for wide‐ranging low‐density species are unlikely to meet those criteria. Established, intensive methods lead to precise estimates, but some new approaches, like short, intensive, (spatial) capture–mark–recapture (CMR/SECR) studies, aided by camera trapping and/or genetic identification of individuals, hold promise. A handbook of best practices in monitoring populations of apex carnivores is strongly recommended.     

Tadpole fingerprinting: Using tail venation patterns to photo-identify tadpole individuals of a threatened frog

Traditional methods for identifying individual amphibians in capture–mark–recapture (CMR) studies have been primarily confined to post-metamorphic stages, using artificial markers that come with a variety of limitations. An alternative that may open CMR studies to earlier life stages involves the use of a species' natural external markers in photo-based identification. In this study, we investigated whether it was possible to distinguish tadpoles of the threatened green and golden bell frog (Litoria aurea) at the individual level based on tail venation patterns. We collected photographs of the tails of captive-raised tadpoles using a smartphone over a 4-week period. This photo-library was used to create an electronic survey where participants were asked to detect matches for query tadpoles from small image pools. We found that most participants agreed on a match for each query, with perfect consensus achieved for most queries (83%). We detected a 14% decline in perfect consensus when participants were asked to match images of tadpoles separated by longer time intervals, suggesting that it is more difficult to visually identify recapture events of L. aurea tadpoles over extended periods due to changes to tail appearance. However, consensus was obtained by participants for all queries, with all matches verified as being correct by the primary researcher. The strength of agreement among participants with no prior experience in matching tadpole tails suggests that there is sufficient inter-individual variation in this feature for individuals to be manually identified. We thus propose that photo-identification is likely to be a valid, non-invasive technique that can be used for short-term studies on tadpole populations that display tail venation. This offers an alternative to artificial markers that may not allow for individual identification, while also opening up tadpole monitoring programmes to citizen scientists who can be recruited online to process image data from home.