Introduction
One of the major consequences of the current biodiversity crisis is defaunation, the local loss or decline of medium- and large-sized vertebrates (Dirzo et al., Reference Dirzo, Young, Galetti, Ceballos, Isaac and Collen2014). Local extinctions of large frugivores are common (e.g. Canale et al., Reference Canale, Peres, Guidorizzi, Gatto and Kierulff2012), leading to changes in seed dispersal and predation (Asquith et al., Reference Asquith, Terborgh, Arnold and Riveros1999; Galetti et al., Reference Galetti, Bovendorp and Guevara2015), which compromise the structure and diversity of plant communities in tropical forests (Janzen, Reference Janzen1970; Connell, Reference Connell, den Boer and Gradwell1971). Such changes in community structure may ultimately affect global ecosystem services such as carbon storage through the replacement of late-successional, dense wooded species by early-successional, rarer ones (Bello et al., Reference Bello, Galetti, Pizo, Magnago, Rocha and Lima2015).
Reintroductions aim to restore populations in places where they were historically extinct (IUCN/SSC, 2013). By restoring populations of large vertebrates the process of defaunation can be effectively reversed (Oliveira-Santos & Fernandez, Reference Oliveira-Santos and Fernandez2010; Seddon et al., Reference Seddon, Griffiths, Soorae and Armstrong2014). However, reintroduction efforts have a low success rate, and unmonitored releases do little to shed light on the causes of either failure or success (Fischer & Lindenmayer, Reference Fischer and Lindenmayer2000). Because knowledge about reintroduced populations is imperfect, management decisions must be made in the face of uncertainty (Armstrong & Seddon, Reference Armstrong and Seddon2008). Consequently, monitoring should be designed to fill gaps in ecological knowledge to assess reintroduction success and compare the value of alternative management options (Nichols & Armstrong, Reference Nichols, Armstrong, Ewen, Armstrong, Parker and Seddon2012). Population monitoring and the estimation of parameters such as abundance, survival and fecundity are paramount to ascertain the fulfilment of demographic criteria of success (Sarrazin, Reference Sarrazin2007; Converse et al., Reference Converse, Moore and Armstrong2013).
Historically, the main focus of reintroductions has been the recovery of threatened species. Efforts directed at restoring ecological processes are still scarce (Polak & Saltz, Reference Polak and Saltz2011). Nevertheless, it has been argued that the scope of reintroduction biology should be broadened to encompass questions about effects on ecosystem processes (e.g. Lipsey & Child, Reference Lipsey and Child2007; Armstrong & Seddon, Reference Armstrong and Seddon2008). Defaunated areas benefit from initiatives that aim not only to recover extirpated populations but also to bring back interactions and processes lost with them (Oliveira-Santos & Fernandez, Reference Oliveira-Santos and Fernandez2010). Agoutis Dasyprocta spp. are arguably the most efficient dispersers of large seeds in Neotropical forests. Their scatter-hoarding behaviour enhances seed dispersal and recruitment of large-seeded plants (Asquith et al., Reference Asquith, Terborgh, Arnold and Riveros1999; Pires & Galetti, Reference Pires and Galetti2012), including threatened taxa such as the Critically Endangered Vouacapoua americana (Forget, Reference Forget1990; Varty & Guadagnin, Reference Varty and Guadagnin1998). The disperser role Dasyprocta plays goes even further; Dasyprocta punctata individuals are able to locate, steal and relocate seed caches of other individuals, increasing the dispersal distance (Jansen et al., Reference Jansen, Hirsch, Emsens, Zamora-Gutierrez, Wikelski and Kays2012). Furthermore, this competition drives D. punctata to use caches far from the parental trees (Hirsch et al., Reference Hirsch, Kays, Pereira, Jansen and Rejmanek2012), which increases the seed germination rate (Carson et al., Reference Carson, Anderson, Leigh, Schnitzer, Carson and Schnitzer2008). Moreover, Dasyprocta spp. tend to cache large seeds more frequently and over longer distances (Galetti et al., Reference Galetti, Donatti, Steffler, Genini, Bovendorp and Fleury2010) than smaller ones, potentially increasing carbon storage in defaunated areas.
A programme was initiated in 2010 to reintroduce the red-rumped agouti Dasyprocta leporina in Tijuca National Park, an impoverished Atlantic Forest reserve in the city of Rio de Janeiro, Brazil. The goal was to restore the recruitment of large-seeded plants by re-establishing a self-sustaining population of D. leporina (Cid et al., Reference Cid, Figueira, Mello, Pires and Fernandez2014), a large rodent, weighing 3–6 kg, whose diet consists mainly of seeds and fruits (Reis et al., Reference Reis, Peracchi, Pedro and Lima2006). Since 2013 we have monitored the reintroduced population through live capture and camera-trapping, estimating demographic parameters such as abundance, survival rate and recruitment. The monitoring is ongoing but we report here the results up to March 2015, with the main goal of assessing the success of the reintroduction through monitoring of population growth, thereby evaluating the need for further releases or alternative management actions. We also aimed to identify threats to the long-term persistence of the population, as well as to direct future monitoring to provide useful information for its adaptive management.
Study area
The area now occupied by Tijuca National Park (3,953 ha; Fig. 1) was used mostly for coffee farming up to the 19th century; it was then reforested, with the intent of restoring the degraded water supplies of Rio de Janeiro city (Pádua, Reference Pádua2002). Despite the reforestation many animal species are absent, including important seed dispersers such as the lowland tapir Tapirus terrestris and the woolly spider monkey Brachyteles arachnoides.
Prior to the beginning of our project D. leporina, which is categorized as Least Concern on the IUCN Red List (Emmonds & Reid, Reference Emmonds and Reid2016), had not been recorded in the Park for at least 2 decades, despite release efforts in the 1970s (Coimbra-Filho et al., Reference Coimbra-Filho, Aldrighi and Martins1973). Thirty-one adults were released in the Park during 2010–2014, in a mixed regimen of hard and soft releases (see Supplementary Table S1 for information on all individuals released). Previously the animals had lived in semi-captivity at two urban parks in Rio de Janeiro city. Cid et al. (Reference Cid, Figueira, Mello, Pires and Fernandez2014) reported the short-term success of the release effort. The reintroduction re-established lost ecological interactions to some extent, through dispersal and caching of palm seeds (Zucaratto, Reference Zucaratto2013).
Methods
Capture and marking
We conducted seven capture sessions from July 2013 to March 2015 (Supplementary Table S2). All individuals captured were handled by veterinarians and sedated using a ketamine–midazolam solution. Adults were fitted with coloured plastic collars. We also marked the fur with individual codes of dots and bars, using two techniques: freeze-branding and fur bleaching. Freeze-branding consisted of cooling an aluminium branding iron in dry ice and 90% alcohol solution, and then applying it to the agouti's shaved skin for 70 seconds. This procedure kills the melanocytes on the epidermis, making the fur grow permanently white (Hadow, Reference Hadow1972). We believe the combination of the two techniques is useful for camera-trapping studies of mammals because (a) the codes are readable at a distance, (b) they may be identified in black & white photographs (many camera trap models use infrared light instead of flash in low-light environments) and (c) they are complementary: fur bleaching is temporary but immediate, whereas freeze-branding is permanent but requires months to reach readability (Plate 1).
Resighting surveys
Resighting surveys were carried out using camera trapping after each capture session, from November 2013 onwards. We used a single camera trap (Bushnell Outdoor Products, Overland Park, USA; Tigrinus Equipamentos para Pesquisa, Timbó, Brazil) per station during most surveys, because agoutis were marked on both sides. Traps were checked and baited with sweet potato every 6 days for a total of 30 days per survey. Trap failure events were recorded to track the variation in sampling effort (Foster & Harmsen, Reference Foster and Harmsen2012).
We used 21 stations during the first two surveys, and 33 thereafter. The additional stations were deployed within the original grid, increasing resighting effort without changing the area surveyed (Fig. 1). We accounted for effort disparity by considering absent stations as trap failures. Trap spacing was always smaller than the lower confidence limit of home range diameters in the study area (351 m; Cid et al., Reference Cid, Figueira, Mello, Pires and Fernandez2014); thus all individuals in the grid were subject to detection (Foster & Harmsen, Reference Foster and Harmsen2012).
Data analysis
Photographic records were identified at the individual level, or as unmarked when there were no marks on the animal. When there were visible marks but the individual could not be identified with certainty the record was defined as marked but unidentified. Records for which it was not possible to determine the presence or absence of marks were defined as unidentifiable and discarded from the analysis, providing information only on the rate of record loss. Resightings of the same individual were considered to be independent when they were at least 1 hour apart or from different stations. Using different time thresholds for independence did not affect our results (Supplementary Fig. S1). Independent records were pooled for each 6-day interval. Following Armstrong & Reynolds (Reference Armstrong, Reynolds, Ewen, Armstrong, Parker and Seddon2012), we excluded records from the first 45 days after release, to separate survival in the wild from post-release survival.
We used the robust design Poisson-log normal mixed-effects mark–resight model of McClintock & White (Reference McClintock and White2009) to estimate the number of unmarked individuals U and apparent survival ϕ (the complement of the sum of mortality and permanent emigration). Other parameters estimated by this model are α (intercept for mean resighting rate, on the log scale), σ (individual heterogeneity in detection rate, on the log scale), γ″ (probability of temporarily moving out of the study area) and γ′ (probability of remaining outside the study area). The model assumes (a) population closure within each primary sampling interval, (b) no loss of marks within each primary interval, (c) perfect distinction between marked and unmarked individuals, and (d) independently and identically distributed resighting probabilities for marked and unmarked individuals.
Sampling effort and identifiability were not constant throughout the study. We sought to account for this variability with two covariates, which should influence resighting rates: sampling effort (i.e. the sum of all camera trap days over each 6-day interval) and record loss (i.e. the proportion of unidentifiable records over the same period). We used the spatial heterogeneity (variance between camera trap stations) of these variables as a covariate for heterogeneity in detection.
We considered the population to be biologically closed throughout each resighting survey, based on a priori analyses (Supplementary Table S3). As the model estimates six parameters simultaneously, the potential number of models can be large. We therefore used a stepwise approach to model selection. We began modelling α while keeping all other parameters except U constant, and used only the best models (those with ΔAICc < 2) for estimation of a subsequent parameter, then repeated the same process for σ, ϕ, γ′′ and γ′, in that order. Covariates used in parameter modelling were origin, age, sex, 6-day interval, resighting survey, sampling effort, spatial variance of sampling effort, record loss and spatial variance of record loss (Table 1).
1α, intercept on the log scale of the resighting rate; σ, variation in resighting rate as a result of individual heterogeneity; ϕ, apparent survival; γ″, probability of temporarily emigrating out of the study area; γ′, probability of remaining outside of the study area; U, number of unmarked individuals in the population.
2σ was a function of eff var or rloss var in models that had an effect of eff or rloss on α.
3γ″ and γ′ are presented together as they were always modelled as a function of the same covariates.
We used parameter estimates from the model to derive parameters of interest for population management. We used a conservative estimate of $\hat N$ by adding the minimum number of marked individuals known alive to the estimated number of unmarked individuals Û. Recruitment was then calculated as
where $\hat \phi _t^* $ is uncorrected, or realized, survival (this distinction is important because most estimation packages correct transition parameters for a constant period of time). Population growth was measured by the finite rate of increase, or λ. An indicator of reintroduction success is λ > 1 (positive growth) without the aid of additional releases (Ostermann et al., Reference Ostermann, Deforge and Edge2001). The finite rate of increase over a period t − 1 to t is calculated as
As $\hat N_t $ and $\hat \phi _t^* $ are estimated parameters with associated standard errors, we were able to derive standard errors for $\hat B_t $ and $\hat \lambda _t $ through parametric bootstrapping, assuming a log-normal distribution for $\hat N$ and a normal distribution for $\hat \phi ^* $ on the logit scale (Williams et al., Reference Williams, Nichols and Conroy2002). All analyses were conducted in R v.3.3.1 (R Development Core Team, 2016), using the package RMark to interface with MARK (White & Burnham, Reference White and Burnham1999), and the source code for all analyses is available at https://doi.org/10.5281/zenodo.202290.
Results
During 428 live-trapping nights we obtained 38 captures of 18 individuals. The number of captures was highest in August 2014 and lowest in November 2014 and March 2015. Half of the individuals captured were young, and this proportion was higher in August 2014 (57.1%; Supplementary Table S4).
The total sampling effort of resighting surveys was 4,764 trap-days, with a mean of 159 ± 34 per interval (Supplementary Table S5). The trap failure rate was high (19.8%). However, we explicitly incorporated variation in sampling effort when estimating resighting rates, and therefore this variation should not bias our estimates. We obtained 3,719 independent records of agoutis, of which 17.64% were discarded as unidentifiable (Supplementary Table S6).
After running the stepwise model selection for all parameters, eight models had ΔAICc < 2 (Table 2; see Supplementary Tables S7–S10 for a full account of models run). Model-averaged estimates did not differ from the best model, and thus we report here the best-fitting model estimates rather than model-averaged ones. Resighting rate was influenced by effort ($\hat \beta _{eff}^\alpha = 0.003 \pm 0.002$), record loss ($\hat \beta _{rloss}^\alpha = - 3.10 \pm 0.58$), captivity ($\hat \beta _{captv}^\alpha = 0.26 \pm 0.11$) and sex ($\hat \beta _{sex}^\alpha = 0.46 \pm 0.20$). Individual heterogeneity was affected by variance of record loss ($\hat \beta _{rloss\; var}^\sigma = 6.50 \pm 3.53$), and survival was influenced by age ($\hat \beta _{age}^\phi = 2.34 \pm 0.97$). Mean apparent survival $\hat \phi $ was estimated to be 0.92 ± 0.04 for adults and 0.54 ± 0.20 for juveniles. Estimated monthly recruitment between surveys was mostly low, with a peak of 8.71 (95% CI 6.23–11.59) during May–August 2014. Estimated population size of wild individuals was lowest in November 2013 (15.68; 95% CI 12.7–19.58) and highest in August 2014 (41.6; 95% CI 37.38–46.46), declining to 29.8 (95% CI 25.55–34.94) in the last survey, in March 2015 (Fig. 2a). The overall finite rate of increase $\hat \lambda _{1,t} $ for the whole study period was 1.92 (95% CI 1.44–2.48).
* Parameter U is omitted because it was always modelled as survey-dependent parameter.
Discussion
Age was the main driver of variation in survival rates; adults had higher probabilities of surviving than young individuals. High mortality of juveniles has also been reported for D. punctata on Barro Colorado Island, Panama (Smythe, Reference Smythe1978), where it was mainly attributed to mammalian predators, such as South American coatis Nasua narica and ocelots Leopardus pardalis (Smythe, Reference Smythe1978; Silvius & Fragoso, Reference Silvius and Fragoso2003). There are no large native mammals in Tijuca National Park that could prey upon adult agoutis (ICMBio, 2008). The abundant coati is a potential predator of young agoutis but we have never observed a coati chasing an agouti; rather, we observed a released agouti foraging alongside a group of coatis. However, we believe a major driver of agouti mortality is predation by domestic dogs Canis lupus familiaris. We repeatedly observed dogs chasing and killing agoutis in the Park.
Release into the wild is often associated with reduced survival, especially when released individuals come from captivity (Jule et al., Reference Jule, Leaver and Lea2008; Tavecchia et al., Reference Tavecchia, Viedma, Martínez-Abraín, Bartolomé, Gómez and Oro2009). However, we found the opposite: Cid et al. (Reference Cid, Figueira, Mello, Pires and Fernandez2014) reported higher post-release survival rates during 2010–2011 (monthly $\hat \phi $: 0.98) than those we estimated for adult individuals during 2013–2015. Because population size and density had increased since the time of release, an increase in competition for food could have worsened mean individual body condition and increased intraspecific aggression. Additionally, predation by dogs may have increased since the 2010 releases because individuals residing in the Park could have developed a search image as encounters with D. leporina became more frequent (generally speaking, a search image develops when a predator becomes accustomed to a common prey item, and thus becomes more efficient at preying upon it; Ishii & Shimada, Reference Ishii and Shimada2010). Furthermore, as we could estimate only apparent survival ϕ and not true survival S, we cannot distinguish mortality from permanent emigration (Williams et al., Reference Williams, Nichols and Conroy2002). Cid et al. (Reference Cid, Figueira, Mello, Pires and Fernandez2014), however, used a known fate model, which could distinguish these confounding rates. Our sampling area is not geographically closed, and therefore low survival may be partially explained by individuals, especially juveniles, dispersing away from the sampling area.
Juveniles were recorded in all surveys. However, the high number of captures of this age group and recruitment in the middle of the year suggest a peak in recruitment. The only recorded decrease in abundance occurred after this peak. Although this decrease is certainly attributable to higher mortality in young individuals, losses must also have been caused by juvenile dispersal. Besides an increase in population size, spatial expansion is also expected in thriving reintroduced populations (Bar-David et al., Reference Bar-David, Saltz and Dayan2005; Molinari-Jobin et al., Reference Molinari-Jobin, Marboutin, Wölfl, Wölfl, Molinari and Fasel2010). Larger-scale monitoring is needed to assess population expansion throughout the Park and to distinguish mortality from dispersal.
In total, 31 individuals were released since 2010; 11 died or were removed within 12 weeks of release, and therefore we consider the release group size to be 20 (Supplementary Table S1). By March 2015 numbers of wild-born individuals fluctuated around 30. The number of wild-born individuals had thus surpassed the release group size. Most of the increase was as a result of the reproductive success of the wild population, as we detected only four reintroduced individuals in our monitoring. Besides, our population estimates are likely to be conservative because our sampling grid was not geographically closed and the population had probably already expanded beyond it. Unassisted population increase is a basic criterion for the success of reintroduced populations (Armstrong & Seddon, Reference Armstrong and Seddon2008) and represents a critical demographic landmark that many reintroductions fail to achieve (Fischer & Lindenmayer, Reference Fischer and Lindenmayer2000). The reintroduction of agoutis had been considered to be successful in the short term (Cid et al., Reference Cid, Figueira, Mello, Pires and Fernandez2014), and we conclude that it has also been successful in the medium term, entering the population growth phase as described by Sarrazin (Reference Sarrazin2007). The long-term success can be ascertained only by the onset of population regulation (Sarrazin, Reference Sarrazin2007), leading to eventual population stability and viability.
We are not aware of what caused the previous reintroduction effort (Coimbra-Filho et al., Reference Coimbra-Filho, Aldrighi and Martins1973) to fail, as no post-release monitoring was undertaken. The number of individuals released was low (25), and thus the release group was vulnerable to demographic stochasticity. Nevertheless, their numbers were not much smaller than ours. The previous releases (Coimbra-Filho et al., Reference Coimbra-Filho, Aldrighi and Martins1973) were sparse both in space and time; they were spread across three release sites c. 1 km apart and individuals were released at a mean rate of 1.5 per month at the site with most releases. In addition, individuals may have been killed by domestic dogs roaming the Park, or by illegal hunters. Even with a positive trend in growth, the same threats may still be present in the Park, and render the agouti population vulnerable to extirpation. Continuous monitoring of the population and further investigation of the drivers of mortality in all age classes are needed to assess the long-term viability of the population.
Despite positive growth overall, the population is still small and thus vulnerable to demographic stochasticity (Caughley, Reference Caughley1994). Additional releases may increase persistence probability and population growth, and decrease the likelihood of bottleneck-induced endogamic depression. On the other hand, during the last 2 years of monitoring the population increased with little aid from additional releases. It is important to assess when releases no longer significantly increase viability (Schaub et al., Reference Schaub, Zink, Beissmann, Sarrazin and Arlettaz2009). At this point financial and logistical resources should be directed to other management actions (such as predator control) or to monitoring programmes that will provide the most useful information for future decision making (Nichols & Armstrong, Reference Nichols, Armstrong, Ewen, Armstrong, Parker and Seddon2012).
The main goal of the reintroduction of D. leporina in Tijuca National Park is to restore recruitment of large-seeded plants (Cid et al., Reference Cid, Figueira, Mello, Pires and Fernandez2014). Seed burial of Astrocaryum aculeatissimum is conditional on the local presence of D. leporina (Zucaratto, Reference Zucaratto2013); however, the agouti's effect on other plant species has not yet been assessed. Moreover, there is no quantitative information available about the influence of local agouti abundance on the rate of seed burial and dispersal distance. On Barro Colorado Island intraspecific competition in D. punctata increases dispersal distances through cache theft (Jansen et al., Reference Jansen, Hirsch, Emsens, Zamora-Gutierrez, Wikelski and Kays2012), and scatter-hoarding is directed to areas with low conspecific tree densities to avoid such thefts (Hirsch et al., Reference Hirsch, Kays, Pereira, Jansen and Rejmanek2012). As the positive effect of these behaviours on seed dispersal should increase with increased density, monitoring D. leporina density and its effect on large-seeded plant demography could establish a minimum density threshold as a criterion of success for the goal of restoring ecological interactions (Polak & Saltz, Reference Polak and Saltz2011).
Despite the high failure rate of reintroduction attempts, some projects succeed even with a small number of founders (Taylor et al., Reference Taylor, Jamieson and Armstrong2005; Reynolds et al., Reference Reynolds, Hatfield, Laniawe, Vekasy, Klavitter and Berkowitz2012). The reintroduction of D. leporina was successful in the medium term after the release of only 31 individuals. Thus, D. leporina has the potential to achieve high numbers from small release groups. The total annual cost of the reintroduction is USD 6,300, including releases in the wild and post-release and population monitoring (but excluding personnel salaries; Table S11), which is significantly smaller than costs (up to USD 1,000,000) reported for other reintroductions (Fischer & Lindenmayer, Reference Fischer and Lindenmayer2000). Thus the reintroduction of D. leporina is a cost-effective tool for habitat restoration, given its key role in the enhancement of large-seeded plant recruitment, low reintroduction costs and high likelihood of success. Moreover, long-term monitoring of both demographic and ecosystem processes resulting from these efforts can identify threats to population persistence and drivers of reintroduction success, as well as the impact of agoutis on the plant community through seed consumption and dispersal.
Acknowledgements
This study received financial support from Conselho Nacional de Desenvolvimento Científico e Tecnológico, Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro, Coordenação de Aperfeiçoamento de Pessoal de Nível Superior and Fundação Grupo Boticário de Proteção à Natureza. The agoutis released were provided by Fundação Parques e Jardins and the administration of Palacete Modesto Leal. Licenses for capture, transport and release were granted by Instituto Chico Mendes de Conservação da Biodiversidade. We thank the staff of Fundação Jardim Zoológico do Rio de Janeiro and our partners at Fundação Oswaldo Cruz for providing veterinary screening and handling support both ex situ and in situ, as well as animal housing pre-release. We are grateful to B. Araujo, A. Arevalo, A. Landim, A.F. Mello, H. Menezes, M. Rheingantz and many others for help in the field, and to two anonymous reviewers for their comments and insights.
Author contributions
The study was conceived and designed by CFK and FASF. Fieldwork, analysis of photographs and curation of data were organized and conducted by CFK, RS and CK. Data analyses were conducted by CFK, with insights from FASF. All authors contributed to writing the article.
Biographical sketches
Caio Fittipaldi Kenup has worked on the reintroduction of the red-rumped agouti in Tijuca National Park since 2011. He is interested in wildlife management, population modelling and reintroduction biology. Raíssa Sepulvida Alves is interested in wildlife management, conservation and restoring ecological interactions. Catharina Kreischer works on the population ecology of the paca and is interested in environmental conservation and reintroduction of fauna. Fernando Fernandez is a population ecologist who has turned his interests to the conservation and management of populations, paleoecology (late Quaternary extinctions) and fostering public understanding of these subjects.