Use made in marine ecology of methods for estimating demographic parameters from size/frequency data

Macdonald and Pitcher's method of decomposing a sizefrequency histogram into cohorts (mathematical optimization of the fit of the distribution function to the histogram) has been used to estimate the composition of random samples drawn from populations with known cohort structure. The large-sample behaviour of the method is in accordance with the results of asymptotic theory. With sample sizes typical of those used in many ecological studies, good estimates often could not be obtained without imposing constraints upon the estimation procedure, even when the number of age classes in the population was known. If the number of age classes was not known, it was frequently difficult to determine from small samples. When unconstrained solutions were obtainable, confidence limits about estimates were often very wide. Our results and information in the theoretical literature indicate that if the Petersen method (whereby several modes on a size-frequency histogram are taken to represent single age classes and all age classes to be present) does not work, accurate estimates of demographic parameters are unlikely to be obtainable using more rigorous methods. In view of these difficulties, we recommend that an iptimization method, such as that described by Macdonald and Pitcher, be used to estimate demographic parameters. Standard errors of estimates should be reported. Optimization methods give an indication when the data is inadequate to obtain accurate parameter estimates, either by failing to converge or by placing large standard errors about the estimates. Graphical methods do not give a clear warning of this, and should be avoided except where the modes on the size-frequency histogram are very well separated and sample sizes are large. Often, assumptions must be made about population parameters to enable their estimation. This may involve constraining some parameters to particular values, assuming a fixed relationship between cohort mean sizes and their standard deviations, or by assuming that individuals grow according to a von Bertalanffy curve. Any such assumptions need detailed justification in each case.     

Estimation of rare and clustered population mean using stratified adaptive cluster sampling

For many clustered populations, the prior information on an initial stratification exists but the exact pattern of the population concentration may not be predicted. Under this situation, the stratified adaptive cluster sampling (SACS) may provide more efficient estimates than the other conventional sampling designs for the estimation of rare and clustered population parameters. For practical interest, we propose a generalized ratio estimator with the single auxiliary variable under the SACS design. The expressions of approximate bias and mean squared error (MSE) for the proposed estimator are derived. Numerical studies are carried out to compare the performances of the proposed generalized estimator over the usual mean and combined ratio estimators under the conventional stratified random sampling (StRS) using a real population of redwood trees in California and generating an artificial population by the Poisson cluster process. Simulation results show that the proposed class of estimators may provide more efficient results than the other estimators considered in this article for the estimation of highly clumped population.

     

Estimating Population Size of Elusive Animals with DNA from Hunter-Collected Feces: Four Methods for Brown Bears

Abstract:  Noninvasive genetic methods can be used to estimate animal abundances and offer several advantages over conventional methods. Few attempts have been made, however, to evaluate the accuracy and precision of the estimates. We compared four methods of estimating population size based on fecal sampling. Two methods used rarefaction indices and two were based on capture-mark-recapture (CMR) estimators, one combining genetic and field data. Volunteer hunters and others collected 1904 fecal samples over 2 consecutive years in a large area containing a well-studied population of brown bears ( Ursus arctos ). On our 49,000-km² study area in south-central Sweden, population size estimates ranged from 378 to 572 bears in 2001 and 273 to 433 bears in 2002, depending on the method of estimation used. The estimates from the best model in the program MARK appeared to be the most accurate, based on the minimum population size estimate from radio-marked bears in a subsection of our sampling area. In addition, MARK models included heterogeneity and temporal variation in detection probabilities, which appeared to be present in our samples. All methods, though, incorrectly suggested a biased sex ratio, probably because of sex differences in detection probabilities and low overall detection probabilities. The population size of elusive animals can be estimated reliably over large areas with noninvasive genetic methods, but we stress the importance of an adequate and well-distributed sampling effort. In cases of biased sampling, calibration with independent estimates may be necessary. We recommend that this noninvasive genetic approach, using the MARK models, be used in the future in areas where sufficient numbers of volunteers can be mobilized.     

Estimating Population Size with Noninvasive Capture-Mark-Recapture Data

Abstract:  Estimating population size of elusive and rare species is challenging. The difficulties in catching such species has triggered the use of samples collected noninvasively, such as feces or hair, from which genetic analysis yields data similar to capture-mark-recapture (CMR) data. There are, however, two differences between classical CMR and noninvasive CMR. First, capture and recapture data are gathered over multiple sampling sessions in classical CMR, whereas in noninvasive CMR they can be obtained from a single sampling session. Second, because of genotyping errors and unlike classical CMR, there is no simple relationship between (genetic) marks and individuals in noninvasive CMR. We evaluated, through simulations, the reliability of population size estimates based on noninvasive CMR. For equal sampling efforts, we compared estimates of population size N obtained from accumulation curves, a maximum likelihood, and a Bayesian estimator. For a closed population and without sampling heterogeneity, estimates obtained from noninvasive CMR were as reliable as estimates from classical CMR. The sampling structure (single or multiple session) did not alter the results, the Bayesian estimator in the case of a single sampling session presented the best compromise between low mean squared error and a 95% confidence interval encompassing the parametric value of N in most simulations. Finally, when suitable field and lab protocols were used, genotyping errors did not substantially bias population size estimates (bias < 3.5% in all simulations). The ability to reliably estimate population size from noninvasive samples taken during a single session offers a new and useful technique for the management and conservation of elusive and rare species.     

Effects of sampling error and temporal correlations in population growth on process variance estimators

Estimates of a population’s growth rate and process variance from time-series data are often used to calculate risk metrics such as the probability of quasi-extinction, but temporal correlations in the data from sampling error, intrinsic population factors, or environmental conditions can bias process variance estimators and detrimentally affect risk predictions. It has been claimed (McNamara and Harding, Ecol Lett 7:16–20, 2004) that estimates of the long-term variance that incorporate observed temporal correlations in population growth are unaffected by sampling error; however, no estimation procedures were proposed for time-series data. We develop a suite of such long-term variance estimators, and use simulated data with temporally autocorrelated population growth and sampling error to evaluate their performance. In some cases, we get nearly unbiased long-term variance estimates despite ignoring sampling error, but the utility of these estimators is questionable because of large estimation uncertainty and difficulties in estimating correlation structure in practice. Process variance estimators that ignored temporal correlations generally gave more precise estimates of the variability in population growth and of the probability of quasi-extinction. We also found that the estimation of probability of quasi-extinction was greatly improved when quasi-extinction thresholds were set relatively close to population levels. Because of precision concerns, we recommend using simple models for risk estimates despite potential biases, and limiting inference to quantifying relative risk; e.g., changes in risk over time for a single population or comparative risk among populations.     

Use of the superpopulation approach to estimate breeding population size: an example in asynchronously breeding birds

Many populations of animals are fluid in both space and time, making estimation of numbers difficult. Much attention has been devoted to estimation of bias in detection of animals that are present at the time of survey. However, an equally important problem is estimation of population size when all animals are not present on all survey occasions. Here, we showcase use of the superpopulation approach to capture-recapture modeling for estimating populations where group membership is asynchronous, and where considerable overlap in group membership among sampling occasions may occur. We estimate total population size of long-legged wading bird (Great Egret and White Ibis) breeding colonies from aerial observations of individually identifiable nests at various times in the nesting season. Initiation and termination of nests were analogous to entry and departure from a population. Estimates using the superpopulation approach were 47-382% larger than peak aerial counts of the same colonies. Our results indicate that the use of the superpopulation approach to model nesting asynchrony provides a considerably less biased and more efficient estimate of nesting activity than traditional methods. We suggest that this approach may also be used to derive population estimates in a variety of situations where group membership is fluid.     

Effects of Recent Environmental Change on Accuracy of Inferences of Extinction Status

Correctly classifying a species as extinct or extant is of critical importance if current rates of biodiversity loss are to be accurately quantified. Observing an extinction event is rare, so in many cases extinction status is inferred using methods based on the analysis of records of historic sighting events. The accuracy of such methods is difficult to test. However, results of recent experiments with microcosm communities suggest that the rate at which a population declines to extinction, potentially driven by varying environmental conditions, may alter one's ability accurately to infer extinction status. We tested how the rate of population decline, driven by historic environmental change, alters the accuracy of 6 commonly applied sighting‐based methods used to infer extinction. We used data from small‐scale experimental communities and recorded wild population extirpations. We assessed how accuracy of the different methods was affected by rate of population decline, search effort, and number of sighting events recorded. Rate of population decline and historic population size of the species affected the accuracy of inferred extinction dates; however, faster declines produced more accurate inferred dates of extinction, but only when population sizes were higher. Optimal linear estimation (OLE) offered the most reliable and robust estimates, though no single method performed best in all situations, and it may be appropriate to use a different method if information regarding historic search efforts is available. OLE provided the most accurate estimates of extinction when the number of sighting events used was >10, and future use of this method should take this into account. Data from experimental populations provide added insight into testing techniques to discern wild extirpation events. Care should be taken designing such experiments so that they mirror closely the abundance dynamics of populations affected by real‐world extirpation events. Efectos del Cambio Ambiental Reciente sobre la Precisión de las Inferencias sobre el Estado de Extinción     

Continuous-time capture-recapture population estimation when capture probabilities vary over time

Closed capture-recapture (CR) estimators have been used extensively to estimate population size. Most closed CR approaches have been developed and evaluated for discrete-time models, but there has been little effort to evaluate their continuous-time counterparts. Continuous-time estimators — developed using maximum likelihood theory by Craig (1953) and Darroch (1958), and martingale theory by Becker (1984) — that allow capture probabilities to vary over time were evaluated using Monte Carlo simulation. Overall, the ML estimators had a smaller MSE. The estimators performed well when model assumptions were upheld, and were somewhat robust to heterogeneity in capture probabilities. However, the estimators were not robust to behavioural effects in the capture probabilities. Time lag effects (periods when animals might be unavailable for immediate recapture) on continuous-time estimates were also investigated and results indicated a positive bias which was greater for smaller populations. There was no gain in performance when using a continuous-time estimator versus a discrete-time estimator on the same simulated data. Usefulness of the continuous-time approach may be limited to study designs where animals are easier to sample using continuous-time methodology.     

Untested assumptions about within-species sample size and missing data in interspecific studies

Phylogenetic comparative studies rely on species-specific data that often contain missing values and/or differ in sample size among species. These phenomena may violate statistical assumptions about the non-random variance component in sampling effort. A major reason why this assumption is often not fulfilled is because the probability of being sampled (i.e., being captured or observed) may depend on species-specific characteristics. Here, we test this assumption by using information on within-species sample sizes and missing data from five independent comparative datasets of European birds. First, we show that the two estimates of data availability (missing values and within-species sample size) are positively correlated and are associated with research effort in general (the number of papers published). Second, we demonstrate biologically meaningful relationships between data availability and phenotypic traits. For example, population size, risk-taking, and habitat specialization independently predicted within-species sample size. The key determinants of missing data were population size and distribution range. However, data availability was not structured by phylogenetic relationships. These results indicate that the accuracy of sampling is repeatable and distributed non-randomly among species, as several species-specific attributes determined the probability of observation. Therefore, data availability seems to be a species-specific trait that can be shaped by ecology, life history, and behavior. Such relationships raise issues about non-random sampling, which requires attention in comparative studies.     

Site-selection bias and apparent population declines in long-term studies

Detecting population declines is a critical task for conservation biology. Logistical difficulties and the spatiotemporal variability of populations make estimation of population declines difficult. For statistical reasons, estimates of population decline may be biased when study sites are chosen based on abundance of the focal species. In this situation, apparent population declines are likely to be detected even if there is no decline. This site-selection bias is mentioned in the literature but is not well known. We used simulations and real population data to examine the effects of site-selection biases on inferences about population trends. We used a left-censoring method to detect population-size patterns consistent with site-selection bias. The site-selection bias is an important consideration for conservation biologists, and we offer suggestions for minimizing or mitigating it in study design and analysis. Article impact statement: Estimates of population declines are biased if studies begin in large populations, and time-series data show a signature of such an effect.     

Estimates of sex- and age-class-specific survival probabilities for a southern Great Barrier Reef green sea turtle population

Sex- and age-class-specific survival probabilities of a southern Great Barrier Reef green sea turtle population were estimated using a capture–mark–recapture (CMR) study and a Cormack–Jolly–Seber (CJS) modelling approach. The CMR history profiles for 954 individual turtles tagged over a 9-year period (1984–1992) were classified into three age classes (adult, subadult, juvenile) based on somatic growth and reproductive traits. Reduced-parameter CJS models, accounting for constant survival and time-specific recapture, fitted best for all age classes. There were no significant sex-specific differences in either survival or recapture probabilities for any age class. Mean annual adult survival was estimated at 0.9482 (95% CI: 0.92–0.98) and was significantly higher than survival for either subadults or juveniles. Mean annual subadult survival was 0.8474 (95% CI: 0.79–0.91), which was not significantly different from mean annual juvenile survival estimated at 0.8804 (95% CI: 0.84–0.93). The time-specific adult recapture probabilities were a function of sampling effort but this was not the case for either juveniles or subadults. The sampling effort effect was accounted for explicitly in the estimation of adult survival and recapture probabilities. These are the first comprehensive sex- and age-class-specific survival and recapture probability estimates for a green sea turtle population derived from a long-term CMR program.Communicated by M.S. Johnson, Crawley     

Estimation of population parameters from aerial counts of North American mallards: a cautionary tale.

The accuracy of population estimates strongly interferes with our ability to obtain unbiased estimates of population parameters based on analyses of time series of population fluctuations. Here we use long-term data on fluctuations in the size of Mallard populations collected as part of the May Breeding Waterfowl Survey covering a large section of North America. We assume a log-linear model of density dependence and use a hierarchical Bayesian state-space approach in which all parameters are assumed to be realizations from a common underlying distribution. Thus, parameters for different populations are not allowed to vary independently of each other. We then simulated independent time series of aerial counts, using the estimated parameters and adding various levels of observation error. These simulations showed that the estimates of stochastic population growth rate and strength of density dependence were biased even when moderate sampling errors were present. In contrast, the estimates of the environmental stochasticity and the carrying capacity were unbiased even for short time series and large observation error. Our results underline the importance of reducing the magnitude of sampling error in the design of large-scale monitoring programs of population fluctuations.     

Using citizen science butterfly counts to predict species population trends

Citizen scientists are increasingly engaged in gathering biodiversity information, but trade‐offs are often required between public engagement goals and reliable data collection. We compared population estimates for 18 widespread butterfly species derived from the first 4 years (2011–2014) of a short‐duration citizen science project (Big Butterfly Count [BBC]) with those from long‐running, standardized monitoring data collected by experienced observers (U.K. Butterfly Monitoring Scheme [UKBMS]). BBC data are gathered during an annual 3‐week period, whereas UKBMS sampling takes place over 6 months each year. An initial comparison with UKBMS data restricted to the 3‐week BBC period revealed that species population changes were significantly correlated between the 2 sources. The short‐duration sampling season rendered BBC counts susceptible to bias caused by interannual phenological variation in the timing of species’ flight periods. The BBC counts were positively related to butterfly phenology and sampling effort. Annual estimates of species abundance and population trends predicted from models including BBC data and weather covariates as a proxy for phenology correlated significantly with those derived from UKBMS data. Overall, citizen science data obtained using a simple sampling protocol produced comparable estimates of butterfly species abundance to data collected through standardized monitoring methods. Although caution is urged in extrapolating from this U.K. study of a small number of common, conspicuous insects, we found that mass‐participation citizen science can simultaneously contribute to public engagement and biodiversity monitoring. Mass‐participation citizen science is not an adequate replacement for standardized biodiversity monitoring but may extend and complement it (e.g., through sampling different land‐use types), as well as serving to reconnect an increasingly urban human population with nature.     

Pheromone Monitoring of Rare and Threatened Insects: Exploiting a Pheromone–Kairomone System to Estimate Prey and Predator Abundance

Abstract: Pheromone‐based monitoring is a promising new method for assessing the conservation status of many threatened insect species. We examined the versatility and usefulness of pheromone‐based monitoring by integrating a pheromone–kairomone trapping system and pitfall trapping system in the monitoring of two saproxylic beetles, the hermit beetle Osmoderma eremita (Coleoptera: Scarabaeidae) and its predator Elater ferrugineus (Coleoptera: Elateridae), which live inside hollow trees. We performed mark–recapture studies of both species with unbaited pitfall traps in oak hollows combined with pheromone‐baited funnel traps suspended from oak branches to intercept dispersing individuals. For O. eremita, the integrated trapping system showed that the population in the study sites may be considerably higher than estimates based on extrapolation from pitfall trapping alone (approximately 3400 vs. 1100 or 1800 individuals, respectively). Recaptures between odor‐baited funnel traps showed that males and females had similar dispersal rates, but estimating the number of dispersing individuals was problematic due to declining recapture probability between subsequent capture events. Our conservative estimate, assuming a linear decrease in capture probability, suggested that around 1900 individuals, or at least half of the O. eremita population, may perform flights from their natal host trees, representing higher dispersal rates than previous estimates. E. ferrugineus was rarely caught in pitfall traps. One hundred thirty‐nine individuals, likely almost exclusively females, were caught in odor‐baited funnel traps with approximately 4% recapture probability. If recapture probability over consecutive capture events follows that of O. eremita, this would correspond to a total population size of 2500–3000 individuals of the predator; similar to its supposed prey O. eremita. Our results demonstrate that pheromone‐based monitoring is a valuable tool in the study of species or life‐history stages that would otherwise be inaccessible.     

Estimating the complexity of bird song by using capture-recapture approaches from community ecology

Repertoire size, the number of unique song or syllable types in the repertoire, is a widely used measure of song complexity in birds, but it is difficult to calculate this exactly in species with large repertoires. A new method of repertoire size estimation applies species richness estimation procedures from community ecology, but such capture-recapture approaches have not been much tested. Here, we establish standardized sampling schemes and estimation procedures using capture-recapture models for syllable repertoires from 18 bird species, and suggest how these may be used to tackle problems of repertoire estimation. Different models, with different assumptions regarding the heterogeneity of the use of syllable types, performed best for different species with different song organizations. For most species, models assuming heterogeneous probability of occurrence of syllables (so-called detection probability) were selected due to the presence of both rare and frequent syllables. Capture-recapture estimates of syllable repertoire size from our small sample did not differ significantly from previous estimates using larger samples of count data. However, the enumeration of syllables in 15 songs yielded significantly lower estimates than previous reports. Hence, heterogeneity in detection probability of syllables should be addressed when estimating repertoire size. This is neglected using simple enumeration procedures, but is taken into account when repertoire size is estimated by appropriate capture-recapture models adjusted for species-specific song organization characteristics. We suggest that such approaches, in combination with standardized sampling, should be applied in species with potentially large repertoire size. On the other hand, in species with small repertoire size and homogenous syllable usage, enumerations may be satisfactory. Although researchers often use repertoire size as a measure of song complexity, listeners to songs are unlikely to count entire repertoires and they may rely on other cues, such as syllable detection probability.Communicated by A. Cockburn     

Optimization of capture–recapture monitoring of elusive species illustrated with a threatened grasshopper

Information on population sizes and trends of threatened species is essential for their conservation, but obtaining reliable estimates can be challenging. We devised a method to improve the precision of estimates of population size obtained from capture–recapture studies for species with low capture and recapture probabilities and short seasonal activity, illustrated with population data of an elusive grasshopper (Prionotropis rhodanica). We used data from 5 capture–recapture studies to identify methodological and environmental factors affecting capture and recapture probabilities and estimates of population size. In a simulation, we used the population size and capture and recapture probability estimates obtained from the field studies to identify the minimum number of sampling occasions needed to obtain unbiased and robust estimates of population size. Based on these results we optimized the capture–recapture design, implemented it in 2 additional studies, and compared their precision with those of the nonoptimized studies. Additionally, we simulated scenarios based on thresholds of population size in criteria C and D of the International Union for Conservation of Nature (IUCN) Red List to investigate whether estimates of population size for elusive species can reliably inform red-list assessments. Identifying parameters that affect capture and recapture probabilities (for the grasshopper time since emergence of first adults) and optimizing field protocols based on this information reduced study effort (−6% to −27% sampling occasions) and provided more precise estimates of population size (reduced coefficient of variation) compared with nonoptimized studies. Estimates of population size from the scenarios based on the IUCN thresholds were mostly unbiased and robust (only the combination of very small populations and little study effort produced unreliable estimates), suggesting capture–recapture can be considered reliable for informing red-list assessments. Although capture–recapture remains difficult and costly for elusive species, our optimization procedure can help determine efficient protocols to increase data quality and minimize monitoring effort.     

Will Observation Error and Biases Ruin the Use of Simple Extinction Models?

Abstract: Estimating the risk of extinction for populations of endangered species is an important component of conservation biology. These estimates must be made from data that contain both environmental noise in the year-to-year transitions in population size (so-called "process error"), random errors in sampling, and possible biases in sampling ( both forms of observation errors). To determine how much faith to place in estimated extinction rates, it is important to know how sensitive they are to observation error. We used three simple, commonly employed models of population dynamics to generate simulated population time series. We then combined random observation error or systematic biases with those data, fit models to the time series data, and observed how close the extinction dynamics of the fitted models compared with the dynamics of the underlying models. We found that systematic biases in sampling rarely affected estimates of extinction risk. We also found that even moderate levels of random observation error do not significantly affect extinction estimates except over a small range of process errors, corresponding to the region where extinction risk is most uncertain. With more substantial sampling error, estimates of extinction risk degraded rapidly. Field census techniques for a variety of taxa often involve observation errors within ±32% of actual population sizes. For typical time series used in conservation, therefore, we often may not need to be overly concerned about observation errors as an extra source of imperfection in our estimated extinction rates.     

Estimating the Effective Population Size of Conserved Populations

Accurate estimation of effective population size is important in attempts to conserve small populations of animals or plants. We review the genetic and ecological methods that have been used to estimate effective population size in the past and suggest that, while genetic methods may often be appropriate for the estimation of N _e, and its monitoring, ecological methods have the advantage of providing data that can help predict the effect of a changed environment on N _e. Estimation of N _e, is particularly complex in populations with overlapping generations, and we summarize previous empirical estimates of N _e that used ecological methods in such populations. Since it is often difficult to assess what parameters and assumptions have been used in previous calculations, we suggest a method that provides a good estimate of N _e, makes clear what assumptions are involved, and yet requires a minimum of information. The method is used to analyze data from 14 studies. In 36% (5) of these studies, our estimate is in excellent agreement with the original, and yet we use significantly less information, in 21% (3) the original estimate is markedly lower, in 43% (6) it is markedly higher. Reasons for the discrepancies are suggested. Two of the underestimates involve a failure in the original to account for a long maturation time, and four of life overestimates involve problems in the original with the correction for overlapping generations.     

Can modeling improve estimation of desert tortoise population densities?

The federally listed desert tortoise (Gopherus agassizii) is currently monitored using distance sampling to estimate population densities. Distance sampling, as with many other techniques for estimating population density, assumes that it is possible to quantify the proportion of animals available to be counted in any census. Because desert tortoises spend much of their life in burrows, and the proportion of tortoises in burrows at any time can be extremely variable, this assumption is difficult to meet. This proportion of animals available to be counted is used as a correction factor (g0) in distance sampling and has been estimated from daily censuses of small populations of tortoises (6-12 individuals). These censuses are costly and produce imprecise estimates of go due to small sample sizes. We used data on tortoise activity from a large (N = 150) experimental population to model activity as a function of the biophysical attributes of the environment, but these models did not improve the precision of estimates from the focal populations. Thus, to evaluate how much of the variance in tortoise activity is apparently not predictable, we assessed whether activity on any particular day can predict activity on subsequent days with essentially identical environmental conditions. Tortoise activity was only weakly correlated on consecutive days, indicating that behavior was not repeatable or consistent among days with similar physical environments.     

Temporal Analysis of Genetic Structure to Assess Population Dynamics of Reintroduced Swift Foxes

Reintroductions are increasingly used to reestablish species, but a paucity of long‐term postrelease monitoring has limited understanding of whether and when viable populations subsequently persist. We conducted temporal genetic analyses of reintroduced populations of swift foxes (Vulpes velox) in Canada (Alberta and Saskatchewan) and the United States (Montana). We used samples collected 4 years apart, 17 years from the initiation of the reintroduction, and 3 years after the conclusion of releases. To assess program success, we genotyped 304 hair samples, subsampled from the known range in 2000 and 2001, and 2005 and 2006, at 7 microsatellite loci. We compared diversity, effective population size, and genetic connectivity over time in each population. Diversity remained stable over time and there was evidence of increasing effective population size. We determined population structure in both periods after correcting for differences in sample sizes. The geographic distribution of these populations roughly corresponded with the original release locations, which suggests the release sites had residual effects on the population structure. However, given that both reintroduction sites had similar source populations, habitat fragmentation, due to cropland, may be associated with the population structure we found. Although our results indicate growing, stable populations, future connectivity analyses are warranted to ensure both populations are not subject to negative small‐population effects. Our results demonstrate the importance of multiple sampling years to fully capture population dynamics of reintroduced populations. Análisis Temporal de la Estructura Genética para Evaluar la Dinámica Poblacional de Zorros (Vulpes velox) Reintroducidos