Roadmap to recovery revealed through the reintroduction of an IUCN Red List species

Reintroductions are powerful tools for tackling biodiversity loss, but the resulting populations can be intrinsically small and vulnerable. It is therefore critical to maximise the number of individuals that are available to contribute to recovery efforts. To address this, we investigated how demographic parameters from a reintroduced population can reveal threats to long-term persistence, inform thresholds for management interventions, and create targets for removing an endangered species from the IUCN Red List. We calculated capture-mark-recapture population estimates for eastern quolls (Dasyurus viverrinus) which had been reintroduced to a fenced reserve in the Australian Capital Territory. We then incorporated the resulting demographic parameters into population viability analyses (PVAs) to estimate probabilities of persistence under several scenarios, including supplementations and harvests (removal of individuals for translocation to other locations). After determining sustainable harvest rates, we then ‘back-cast’ the population size and occupancy area required to remove the species from the IUCN Red List within 10 years. Our demographic results indicated high mean apparent survival (90% ± 5), and PVAs revealed the probability of persistence over a 50-year time horizon was 50.5% with no interventions, 0% when the population was harvested of > 6 individuals, and 100% if harvests ≤ 54 juveniles were combined with an annual supplementation of ten maternal females (with ≤ 6 young each). Based on this model, a total harvest area of 413 km2 and an occupancy area of 437 km2 would be needed to recover the species within 10 years (i.e., 90 similar fenced reserves, not accounting for edge effects). Due to the inherent difficulty in securing large areas for species recovery, we see these ambitious targets as a call to create coordinated and collaborative sanctuary networks where species can be managed as a metapopulation across multiple sites. By taking advantage of a rapid life history and harvesting the ‘doomed surplus’, managers can achieve their stretch goals for species recovery in the long term.


Introduction
Defaunation in the Anthropocene, driven by human-induced environmental change and destruction, threatens biodiversity, ecosystem function, and human health worldwide (Dirzo et al. 2014). Restoration presents a major challenge for the next century, but to avoid long-term goals being limited by short-term human memory of ecosystems (i.e., 'shifting baseline syndrome ', Pauly 1995;Miller 2005;Manning et al. 2006), robust monitoring and species listings are paramount. The International Union for Conservation of Nature's Red List of Threatened Species (IUCN 2021) provides a powerful tool used by conservationists and researchers to list declining species, and galvanise conservation action and policy change (Rodrigues et al. 2006;Betts et al. 2020). It provides a repository for information related to species range, population size, threats, and conservation actions, and compares these to broadly applicable and standardised criteria to categorise species from 'Critically endangered' to 'Least concern'. These criteria enable managers to quantify the targets required to downlist (reclassify from higher risk to lower risk categories) or delist (remove from the IUCN Red List) a species, and encourages small-scale projects to be unified under long-term visions for species recovery.
Reintroductions are a critical tool used to reverse defaunation and restore ecosystem function (Armstrong and Seddon 2008). While reintroductions aim to establish a viable and self-sustaining population (IUCN 2013), management decisions must always be made in the face of imperfect knowledge about species and ecosystems (Armstrong and Seddon 2008). This uncertainty can be addressed using robust monitoring of demographic parameters including abundance, survival, and reproduction; which can indicate population selfsustainability, density-dependence, carrying capacity, and threats to these processes (e.g., Manning et al. 2019;Parlato et al. 2021). This is particularly crucial when a population transitions through the three phases of reintroduction, from 'establishment' (where post-release effects drive population dynamics) to 'growth' (characterised by high rates of expansion) and, finally, to 'regulation' (where density dependence limits survival and recruitment, Sarrazin 2007). However, many reintroduction studies span short timeframes (i.e., < 5 years, Parlato et al. 2021), limiting their ability to estimate temporal variation in vital rates over the long term (e.g., Leech et al. 2007;Cremona et al. 2017). From the outset, reintroduction programmes should be focused on a long-term vision of species recovery, which necessarily requires robust demographic monitoring to ensure the population reaches the regulation phase (Armstrong and Reynolds 2012;McCarthy et al. 2012;Nichols and Armstrong 2012).
As part of this long-term approach, demographic parameters from a reintroduced population can be built into stochastic population models (i.e., for population viability analyses [PVAs], Lindenmayer et al. 1993). These models can be used to compare alternative management interventions including supplementations (reinforcement of individuals) and harvests (removal of individuals for translocation to other locations), and thus provide crucial insights to guide these management actions. For example, PVAs have revealed that expanding the reintroduction of bearded vultures (Gypaetus barbatus) would dangerously deplete the captive source population (Bustamante 1996), and have highlighted the need to establish a second population of Eurasian beavers (Castor fiber) to ensure their viability in the Netherlands (Nolet and Baveco 1996). Since population growth is often affected by population size, this can result in a trend toward a constant breeding density (i.e., densitydependence, Sibly et al. 2002). Above this density, excess individuals can be considered a 'doomed surplus' (e.g., northern bobwhite Colinus virginianus, Errington 1945, but see Williams et al. 2004). Harvesting these individuals (especially for species with high 1 3 fertility) could offer an ideal opportunity to maximise the number of individuals available for translocation to begin or reinforce other populations (e.g., demonstrated in black-footed ferrets Mustela nigripes, Biggins et al. 2011). Reintroduced populations, however, can be intrinsically small and vulnerable to stochastic effects (Lacy 2000), and therefore the impacts of harvesting from these populations should be simulated before any individuals are removed. PVAs can inform thresholds for harvests of the doomed surplus for translocation by indicating the maximum sustainable number of individuals available to contribute to species recovery.
We used a reintroduced population of eastern quolls (Dasyurus viverrinus) to investigate how demographic monitoring and PVAs can reveal viable interventions (i.e., supplementations and harvests) that can contribute to removing an endangered species from the IUCN Red List. Specifically, we (1) modelled survival and recruitment rates throughout the establishment, growth, and regulation phases of the eastern quoll reintroduction, (2) incorporated these demographic parameters into stochastic population models to reveal long-term viability under different management scenarios (i.e., no intervention, supplementation, harvest, and combinations of these), and (3) determined the contribution of our programme toward species recovery. Since long-term vision is essential for effective restoration, we set an ambitious end point (or 'stretch goal') of eastern quoll species recovery within 10 years, and then retrospectively calculated the area of habitat (henceforth 'area of occupancy') and number of harvests required to achieve this goal (also known as 'back-casting', Manning et al. 2006).
Australia has suffered the highest rate of mammal extinctions of any continent (Woinarski et al. 2015), due in large part to predation by introduced species (e.g., red fox Vulpes vulpes, and feral cat Felis catus, Kinnear et al. 2002;Radford et al. 2018). To circumvent these threatening processes, significant efforts have been made to reintroduce species where introduced predators are absent, such as in conservation-fenced areas (Hayward and Kerley 2009;Moseby et al. 2011;Legge et al. 2018). The benefits of conducting long-term, large-scale experiments under such fenced conditions are increasingly being recognised, and allow researchers to build an understanding of ecological processes which may otherwise be impossible (Hester et al. 2000;Manning et al. 2009).
Our eastern quoll reintroduction took place in this fenced context, and presents a unique model for testing the effect of conservation actions on species recovery because (1) the species is categorised as 'Endangered' under the IUCN Red List (Burbidge and Woinarski 2016), (2) the founding population was small but within the normal range for mammal reintroductions (n = 44), (3) reintroductions to fenced areas provide an ideal opportunity to undertake ecological experiments, (4) the population has been robustly monitored for over five years, and (5) the programme emulates small conservation projects across the globe, creating broad relevance of our outcomes to long-term reintroduction planning.

Ethics statement
Translocations were carried out under licenses from the Tasmanian Department of Primary Industries, Parks, Water and Environment (permits TFA 16025 and 17091, export licences 12818/16 and 13528/17), Victorian Department of Environment, Land, Water and Planning (permit 14505167), and Australian Capital Territory Government (scientific licence 1 3 LT2017959, import licence L120161261). Reintroduction (protocol A2016/02) and monitoring procedures (protocol A2020/40) were approved by The Australian National University Animal Experimentation Ethics Committee.

Study area
This study took place at Mulligans Flat Woodland Sanctuary (henceforth Mulligans Flat), a 485 ha public nature reserve containing critically endangered box-gum grassy woodland (McIntyre et al. 2010) situated on Ngunnawal and Ngambri Country in the Australian Capital Territory (−35.167, 149.158). Mulligans Flat is part of the Mulligans Flat-Goorooyarroo Woodland Experiment (https:// www. mfgow oodla ndexp erime nt. org. au/) and functions as an 'outdoor laboratory' where restoration techniques are trialed (Manning et al. 2011;Shorthouse et al. 2012). This includes the reintroduction of locally extinct species, such as the eastern bettong (Bettongia gaimardi, Manning et al. 2019), bush stone-curlew (Burhinus grallarius, Rapley 2020), and New Holland mouse (Pseudomys novaehollandiae, Abicair et al. 2020, Smith et al. 2022. To enable these reintroductions, the sanctuary is enclosed by a 11.5 km fence which excludes introduced species including the red fox, feral cat, European rabbit (Oryctolagus cuniculus), and European hare (Lepus europaeus, Shorthouse et al. 2012). While these threatening and destructive species have been eradicated from within the sanctuary, conditions are otherwise like other unfenced woodlands in the region. The fence design includes an overhang on the outside which prevents entry by introduced species (Shorthouse et al. 2012), but does not prevent agile species from climbing out of the sanctuary into the surrounding landscape.

Study species
The eastern quoll ('murunguny' in the Indigenous Ngunnawal language) is a solitary, small-to-medium (0.7-1.9 kg, Jones et al. 2001) marsupial carnivore (family Dasyuridae, Stannard and Old 2013). It is an opportunistic hunter with a diet dominated by invertebrates, but will also depredate birds, reptiles, and mammals, and scavenge on carcasses (Blackhall 1980;Godsell 1983).
The species is promiscuous and facultatively monoestrous, breeding synchronously in early Austral winter (Fletcher 1985), giving birth after 20 days gestation, and weaning ≤ 6 young in spring when food availability is high (Godsell 1983). Eastern quolls are sexually dimorphic, with males being larger (0.9-2 kg) and having larger home ranges (mean 44 ha) than females (0.7-1.1 kg, mean 35 ha, Godsell 1983). Populations reach their highest densities in early summer when the annual cohort of juveniles disperse from their natal dens, and lowest densities during winter largely due to juvenile mortality (Godsell 1983). The combination of a short lifespan (3-4 years) and these seasonal population fluctuations result in high population turnover (Jones et al. 2001).
Historically, eastern quolls were distributed throughout south-eastern Australia until the 1960s. The species disappeared from all but the southern island state of Tasmania due to a combination of habitat destruction, disease, human persecution, and predation by introduced species (particularly red foxes, Jones et al. 2001;Peacock and Abbott 2014). While there is no robust assessment of population size, state-wide spotlighting surveys revealed a 52% decline in sightings at 150 sites across Tasmania between 1999 and 2009 (i.e., 10,400 decline in population size from a 20,000 estimate, Fancourt et al. 2013). In addition, trapping surveys revealed a 61-100% decline at historical 'hotspots' (with disproportionately high eastern quoll densities compared with other parts of Tasmania, i.e., Cradoc, Cradle Mountain, and Buckland) compared with trapping conducted 18-31 years earlier. In response, the species was listed as 'Endangered' by the IUCN (Burbidge and Woinarski 2016) and the Australian Threatened Species Scientific Committee (2015).

Reintroduction
To reintroduce eastern quolls to Mulligans Flat, we adopted an adaptive translocation tactics approach (sensu Batson et al. 2015); involving a series of iterative trials where learnings were used to refine tactics for the following trial. When the first trial in 2016 revealed high male mortality associated with increased conspecific aggression and overdispersal (87.5% mortality, Wilson et al. 2020;Wilson et al. 2021), we selected only female founders for subsequent trials to maximise survival (12.5-23.1% mortality). By translocating maternal females (henceforth mothers) in winter, several were either pregnant or carrying pouch young, allowing us to reintroduce juvenile males and females 'via the pouch'. This also potentially increased genetic diversity because multiple sires can be represented in a litter (B. Brockett unpublished data).

Monitoring
We used a Robust Design capture-mark-recapture (CMR) framework to conduct demographic monitoring of the reintroduced eastern quoll population (Pollock 1982;Kendall and Nichols 2002). We conducted eight CMR primary sessions, each consisting of two trap nights 2-3 days apart (secondary sessions), in Austral summer and autumn each year between 2017 and 2022. Trapping during summer was intended to detect the greatest population density following juvenile dispersal in spring, and trapping during autumn was intended to detect the size of the breeding population. We integrated one night of freefeeding before each trap night to encourage the capture of more 'trap-shy' individuals (Biro 2013).
We standardised an array of 92 trap sites across Mulligans Flat, each placed 25 m from vehicle tracks and approximately 200 m apart (Fig. 1). We used wire cage traps (31 cm × 31 cm × 70 cm) baited with sardines, and for each trapped animal we inserted a microchip for identification, sampled fur, scat, and skin (biopsy for genetic material), and recorded sex, weight, and pouch occupancy (as per Portas et al. 2020).

Demography
We fitted Robust Design Pradel Recruitment Closed Population Estimation models (Kendall et al. 1995(Kendall et al. , 1997Pradel 1996) Kendall and Nichols 2002). We assumed emigration and mortality only occurred between primary sessions and that population growth reflected young eastern quolls recruited to adulthood (and not immigration, since Mulligans Flat is fenced). The Pradel Robust Design allowed us to derive estimates of population size (N) at each primary session from initial capture probability (p), recapture probability (c), apparent survival (φ), and recruitment (f) using a logit link function. Since the length of time between primary and secondary sessions varied (either by design or logistical constraints), we ensured the sampling regime was reflected in the models.
We tested 12 candidate models, which included whether p was equal to c or varied from each other (e.g., due to 'trap shyness'), and whether p and c were constant over time (null), or varied by session, season, trap night, minimum temperature (°C), maximum temperature (°C), and rainfall on the day of trapping (Bureau of Meteorology 2021), and whether there was an additive effect of trap availability (fraction of traps available to eastern quolls after removing traps made unavailable by other species and defective traps), across all individuals and between sexes (Table 1). We assessed models based on Akaike's Information Criterion corrected for small sample sizes (AICc) to predict the final parameter estimates (Burnham and Anderson 2002).
We also investigated whether body weight varied by primary session or estimated population size (N) by fitting linear models. Since the eastern quoll is sexually dimorphic, females and males were modelled separately. Finally, we estimated eastern quoll density by dividing the trapped area (485 ha) by the mean female and male estimates (N) across all sessions excluding autumn 2017 (where the population was still establishing). Demographic analyses were conducted within the R environment (version 4.1.2, R Core Team

Viability
We explored the long-term viability of the reintroduced eastern quoll population over a 50-year time horizon using individual-based models simulating a hypothetical, non-spatially explicit population incorporating demographic parameters from our top-ranking model (determined above), and ecological data available for the species (Table 2). Due to their rapid life history (3-4 year lifespan), we considered only two life stages where individuals > 1 year old were classed as adults, otherwise they were classed as juveniles.
Since eastern quolls, and especially males, have naturally high rates of mortality during the juvenile dispersal (summer) and breeding periods (early winter, Godsell 1983;Wilson et al. 2020), we aimed to determine a threshold for sustainably harvesting (i.e., removing individuals for translocation to other locations) this 'doomed surplus' (animals that would never survive the seasonal bottleneck, Errington 1945). In Mulligans Flat, this process likely manifests as an exodus over the conservation fence due to limited territory. In addition, supplementing mothers (carrying ≤ 6 pouch young) allows managers to translocate 'seven for the price of one' (Wilson et al. 2020). As such, we ordered the simulated events in a year as: (1) setting of annual rates (EV), (2) aging, (3) carrying capacity (K) truncation, (4) breeding (early winter) with a census, (5) supplementation (mid-winter) with a census, (6) growth rate (r) calculation, (7) harvest (late spring) with a census, and (8) mortality (summer) with a census (for definitions see Lacy and Pollak 2021).
We modelled the following annual scenarios based on a stable stage distribution: (a) no intervention, (b) harvest (i.e., removed) of six (1:1 sex ratio) juveniles, (c) supplementation (i.e., reinforcement) with one mother (carrying six young, effectively n = 7), (d) supplementation with one mother (effectively n = 7) and harvest of four juveniles, and (e) supplementation with 10 mothers (effectively n = 70) and harvest of 54 juveniles. We simulated models with 1000 iterations to account for stochasticity in parameter estimates and increase model precision, and did not include catastrophes. Inbreeding depression was included for the 'no interventions' and harvest-only scenarios, but was not included for scenarios involving supplementations because inbreeding effects would likely be negated (Mills and Allendorf 1996). Parameters derived from the PVAs included population growth rate (λ) and probability of persistence (percent). PVAs were conducted using Vortex 10.5.5 (Lacy Table 1 Top and second-ranking predictor models for population estimates (by individual and by sex) based on capture-mark-recapture of reintroduced eastern quolls (Dasyurus viverrinus) at Mulligans Flat Woodland Sanctuary, Australian Capital Territory 'p' indicates initial capture probability, 'c' recapture probability, 'K' the number of parameters, 'AICc' Akaike's information criterion corrected for small sample sizes, and '△AICc' the difference between the top-and the second-ranking models

Recovery
To reveal a roadmap towards eastern quoll species recovery, we summarised the IUCN Red List criteria for the status of 'Critically endangered', 'Endangered', and 'Vulnerable' (IUCN 2021), and recommended actions that would result in delisting the species based on the eastern quoll assessment (i.e., Burbidge and Woinarski 2016). We then incorporated these targets with our PVA results (i.e., sustainable harvests) to calculate the annual contribution of Mulligans Flat towards species recovery (i.e., offsetting the population reductions and increasing geographic range that placed the species in the IUCN Red List 'Vulnerable' category). Finally, we back-casted the number of harvests and the area of occupancy that would be required achieve our stretch goal of species recovery within 10 years (Manning et al. 2006).

Demography
During 1,472 trapping nights, we made 421 eastern quoll captures (155 unique individuals, 101 females, 54 males) at Mulligans Flat over eight trapping sessions between 2017 and 2022. 56.44% of females and 70.37% of males were recaptured at least once. The topranking CMR model included capture probability (p) and recapture probability (c) varying from each other, and trap night (Table 1).
In the first monitoring session in autumn 2017, we caught one founder and seven sanctuary-born eastern quolls (n and N = 8). Between 2018 and 2020 population estimates oscillated, with summer estimates being expectedly greater (coinciding with juvenile dispersal, mean = 37 ± 6.24), and autumn estimates being lower and more variable (between the dispersal and breeding periods, mean = 34.56 ± 14.14, Fig. 2a). Sex ratios were relatively balanced until summer 2018 (1.12 females: 1 male), but by the following summer 2019 session they skewed heavily toward females (2.24 females: 1 male). The autumn 2021 population estimate deviated from this oscillating trend and reached a new peak (N = 51), with a sex ratio that was also skewed toward females (2.13 females: 1 male). The final autumn 2022 session achieved a similar population estimate (N = 47) and sex ratio (2.35 females: 1 male).
Mean body weights (females 0.84 kg ± 0.13, males 1.13 kg ± 0.15) oscillated between a maximum in autumn and minimum in summer; the opposite to the trend observed for population estimates (Fig. 2b). However, body weights were unexpectedly lower in autumn 2021 compared to previous autumns (females 0.74 kg ± 0.04, males 1.02 kg ± 0.08), however this was followed by the greatest body weights observed across the study period in autumn 2022 (females 1.02 kg ± 0.04, males 1.35 kg ± 0.08). When fitting a linear model with body weights (kg) against estimated population size (N), we found a significantly 1 3 negative association (p = 0.034), suggesting the population had become density-dependent six years after the first trial reintroduction.
The mean density of eastern quolls across sessions (excluding autumn 2017 where the population was still establishing) was one female per 19.57 ha (± 2.84), one male per 37.51 ha (± 1.06), and for both sexes, one individual per 12.97 ha (± 3.76). The maximum density of eastern quolls in Mulligans Flat was one individual per 9.53 ha in autumn 2021, indicating a maximum carrying capacity of 51 adults (34 females: 17 males).

Viability
PVAs revealed that with no interventions, the Mulligans Flat population would have a 50.5% probability of persistence over the next 50 years (Fig. 3). This probability fell to 0% when > 6 juveniles were harvested from the population annually, but rose to 100% if the population was supplemented with at least one mother annually (carrying ≤ 6 young, effectively n = 7). Further, four juveniles could be sustainably (100% probability of persistence) harvested from the population annually as long as one mother (effectively n = 7) was also supplemented into the population annually. Finally, 54 juveniles could be sustainably harvested if ten mothers (effectively n = 70) were supplemented annually. The deterministic annual population growth rate (r) was 0.5063 across all scenarios, and probability of persistence was sensitive to carrying capacity (K) and the number of individuals supplemented and/or harvested.

Recovery
The estimated eastern quoll population size in Tasmania was ~ 20,000 prior to its 52% decline (Fancourt et al. 2013). To avoid meeting the criteria for IUCN Red List 'Vulnerable' category, the estimated current population of ~ 9600 would need to be increased to   1, IUCN 2001) for the categories of 'Critically endangered', 'Endangered' (orange), 'Vulnerable', and 'Near-threatened' or 'Least concern' (green). Criteria that were achieved by an "any of" clause are indicated in yellow. * was based on Fancourt et al. (2013) 1 3 14,200 (i.e., a < 30% decline from the original population size, Table 3). This provides us with a stretch goal of producing 4600 eastern quolls, equivalent to 460 individuals per year for 10 years, to achieve species recovery. Since we can sustainably harvest the Mulligans Flat population of 54 individuals per year across its 4.85 km 2 area, this implies that to harvest 460 juveniles, an area of 41.31 km 2 would be needed (this area would need to have conditions comparable to Mulligans Flat, i.e., without introduced predators). Finally, a total of 437.45 km 2 (4600 individuals / Mulligans Flat density [10.52 individuals per km 2 ]) would be required to sustain these individuals beyond 10 years (i.e., 90 sanctuaries similar to Mulligans Flat, not accounting for edge effects).
In addition, eastern quolls need to occur in > 10 'locations' ("geographically or ecologically distinct areas in which a single threatening event can rapidly affect all individuals of the taxon present", IUCN 2001) to avoid meeting the criteria for the IUCN Red List 'Vulnerable' category (Table 3). We interpreted a 'single threatening event' as being spatially and temporally explicit (e.g., flood) as opposed to a threatening process which is spatially and temporally dynamic (e.g., disease, climate change). Burbidge and Woinarski (2016) indicated there were two locations: Tasmania and Bruny Island. We suggest that Mulligans Flat and Mt Rothwell should be also considered locations because they support self-sustaining, density-dependent, and geographically distinct populations since their reintroductions in 2016 and 2003, respectively. Thus, an additional seven locations with similar levels of introduced predator mitigation will be required. Finally, while not quantified, Burbidge and Woinarski (2016) reported a decline in the eastern quoll's extent of occurrence (EOO) and area of occupancy (AOO). These criteria could be addressed by increasing the number of eastern quoll locations on mainland Australia.

Discussion
We have demonstrated how incorporating a reintroduced population's demographic parameters into PVAs to determine the viability of harvesting a 'doomed surplus' can contribute measurably to species recovery. We also highlighted how back-casting can reveal the pathway to removing an endangered species from the IUCN Red List within 10 years (i.e., by reversing population declines and increasing the number of locations). While our targets may appear daunting, a stretch goal, by definition, must be ambitious enough to inspire the creativity and innovation to achieve long-term outcomes that currently seem impossible (Manning et al. 2006). Progress toward large-scale conservation and restoration requires such innovation to prevent the 'locking-in' of the current shifting baseline ).

Demography
The reintroduced eastern quoll population at Mulligans Flat grew rapidly in the absence of introduced predators, despite a limited number of founders (n = 44). After autumn 2017, the population oscillated between maximums in summer and minimums in autumn for the next 3 years (Fig. 2a), depicting a transition from the establishment phase (2017) to the growth phase (2018-2019), and finally, to the regulation phase (2020-2022, sensu Sarrazin 2007). The combination of a rapid life history (2-3 year lifespan), high reproductive success (100% female breeding, Godsell 1983), and high fertility (≤ 6 progeny per year) created an ideal scenario for producing an insurance population on mainland Australia.
Mean body weights oscillated between maximums in autumn and minimums in winter (Fig. 2b) and had a significantly negative association with population estimates (p = 0.034). This indicates density-dependence, a fundamental objective of any reintroduction, and mechanistically could have resulted from conspecific competition for prey or territory (e.g., dens). Interestingly, the autumn 2021 and 2022 sessions deviated from the oscillating trends in both population estimates and mean body weights. Australia suffered longterm rainfall deficiencies between 2017 and 2020, with conditions easing in March 2021 (Bureau of Meteorology 2021). The severe climatic conditions likely affected available prey for eastern quolls (similarly observed in Tasmania by Fancourt et al. 2018), thereby limiting the carrying capacity of Mulligans Flat in these years. In the 2021 session the population estimates reached a new maximum (N = 51) but mean body weights were lower (0.82 kg ± 0.04) than in previous autumn sessions (mean 1.033 kg ± 0.11). This suggests that as population size increased, resources (e.g., prey items) may have become limiting, resulting in lower body weights.
Interestingly, mean apparent survival was similar for females (91% ± 6) and males (89% ± 6) in Mulligans Flat, whereas historical survival rates in Tasmania differed between the sexes (females 63%, males 25%, Godsell 1983). Sex ratios were relatively balanced until summer 2018 (1.12 females: 1 male), after which it became heavily skewed toward females (2.24 females: 1 male in summer 2019). Our population estimates indicate a maximum carrying capacity of 34 females, each with longer lifespans (2-3 years) than males (1-2 years). While females have smaller and overlapping home ranges (Wilson et al. 2020), there may not be enough territory to support all the juveniles they produced (34 females produce ≤ 204 young per year, Godsell 1983); explaining high adult survival but low recruitment. There is no evidence to suggest eastern quolls reduce their fecundity with density-dependence (unlike reintroduced eastern bettongs, Manning et al. 2019). Rather, if favourable environmental conditions increased the carrying capacity of Mulligans Flat, territory could become limiting by lowering the fecundity of females that cannot secure natal territory (as observed in bobcats Lynx rufus, Knick 1990).
Estimated recruitment was low but consistent (females 7 ± 5%, males 7 ± 7%), suggesting that outcompeted individuals either died in Mulligans Flat or emigrated over the conservation fence into the surrounding landscape (note that we estimated apparent, rather than true, survival because we could not distinguish mortality from emigration, Williams et al. 2002). Such overdispersal is a problem in reintroductions, where individuals disperse away from the recipient site and do not contribute to population establishment (Richards and Short 2003). Serendipitously, some of these migrants have colonised the adjoining Goorooyarroo Nature Reserve which now also has an conservation fence (S. Stratford pers comms), thereby founding a new population (the Mulligans Flat fence only allows one-way passage to Goorooyarroo). Additionally, this spill-over or halo effect (Tanentzap and Lloyd 2017) could be used to colonise the landscape 'beyond-the-fence' if introduced predators can be maintained below the tolerance levels of eastern quolls (sensu Evans et al. 2021).
It has been suggested that a population must contain at least 1000 individuals to maintain "adequate adaptive potential … in the face of environmental change" (Willi et al. 2006;Weeks et al. 2011). Due to its limited size this target is not feasible for Mulligans Flat, but reinforcement translocations using the one-migrant-per-generation method (Mills and Allendorf 1996) could negate the effects of small population size (i.e., inbreeding depression, Weeks et al. 2011). It is encouraging that 100% probability of persistence was achieved by supplementing the population with one new mother annually, highlighting the need to increase gene flow between isolated populations. As such, we recommend that Mulligans Flat and other reintroduction locations be treated as a metapopulation; 1 3 translocating individuals between them to promote in situ genetic diversity (Weeks et al. 2011;Frankham 2015).

Viability
Despite reaching the regulation phase, the Mulligans Flat eastern quoll population is inherently small and vulnerable to demographic stochasticity (Caughley 1994), and self-sustainability does not necessarily translate to long-term persistence (Seddon 1999). There is a 50.5% likelihood of persistence over the next 50 years (Fig. 3), emphasising the importance of ongoing management interventions. Annual supplementation of one mother could stabilise the population over 50 years, though we note two assumptions: ecological conditions will remain similar to the most recent monitoring sessions (autumn 2021 and 2022), and the population will maintain similar vital rates. Though we assumed high juvenile mortality (females 64.17% ± 19.92, males 64.93% ± 19.87, based on Godsell 1983), similar studies have demonstrated how increases in juvenile mortality can trigger comprehensive mortality and recruitment failure (e.g., in northern quolls Dasyurus hallucatus, Cremona et al. 2017, and African lions Panthera leo, Barthold et al. 2016). For example, after persisting for two decades, reintroduced Arabian oryx (Oryx leucoryx) suffered a poaching epidemic which rendered the population non-viable (Stanley Price 1989). While we present more than the minimum 5 years of vital rates required to identify temporal parameters driving variation (Gelman and Hill 2006;Parlato et al. 2021), each additional year of monitoring will improve our inferences; highlighting the value of long-term datasets that inform longterm goals.

Recovery
We explored how our programme could contribute measurably to eastern quoll species recovery. For the species to be delisted, we need to produce ≥ 4,600 eastern quolls to raise the total number from 9,600 (after 52% decline from 20,000 estimate in Tasmania, Fancourt et al. 2013) to the 14,200 (≤ 71% of 20,000 estimate) required to avoid meeting the IUCN Red List criteria for the 'Vulnerable' category (i.e., ≥ 30% decline). In addition, the eastern quoll must occur in > 10 locations ("geographically or ecologically distinct areas in which a single threatening event can rapidly affect all individuals of the taxon present", IUCN 2001). While there are currently two recognised locations (Tasmania and Bruny Island, Burbidge and Woinarski 2016), based on the definition for 'location' we suggest that Mulligans Flat and Mt Rothwell should also be recognised. Further, we posit that Tasmania (2,320 km 2 occupancy area) may represent more than one location, since significant genetic structure with consistent regional differentiation related to geographic distance has been found between populations (Cardoso et al. 2014). This geographical and/or behavioural separation suggests a single threatening event is unlikely to endanger all eastern quolls in Tasmania.
To achieve a stretch goal of species recovery within 10 years, a 41 km 2 harvest area with conditions comparable to Mulligans Flat (i.e., without introduced predators) will need to be harvested of 460 individuals annually, and a 437 km 2 occupancy area would be required to sustain these individuals beyond 10 years (i.e., 90 sanctuaries similar to Mulligans Flat, not accounting for edge effects). Due to the inherent difficulty in securing large areas for species recovery, we see these ambitious targets as a call to create a coordinated and collaborative sanctuary network where the eastern quoll, and other species, can be managed as a metapopulation across multiple sites (e.g., South Eastern Australia Sanctuary Operations Network or 'SEASON', Sharp 2021). In such a network, decisions regarding management interventions should be based on robust monitoring of the populations' demographics and genetic composition. This network would buffer against the demographic and genetic perils facing isolated populations, and limit edge effects associated with small occupancy areas (McGregor et al. 2020). To identify appropriate areas for future reintroduction sites, we recommend incorporating eastern quoll occurrence data from established mainland populations, such as Mulligans Flat, into broad-scale habitat modelling across the species' former range (e.g., maximum entropy species distribution modelling).
While fenced sanctuaries have produced insurance populations of at least 38 species that are susceptible to introduced predators, they are limited in area and capacity to expand, and maintaining them comes at a cost (Ringma et al. 2017;Legge et al. 2018). Reintroducing species 'beyond-the-fence' where introduced predators are actively managed (and adaptively calibrated) to remain below species' tolerance levels (sensu Evans et al. 2021) is the next frontier to establish viable, self-sustaining populations and return ecological functions to our increasingly defaunated landscape (James and Eldridge 2007). Maintaining such a 'Goldilocks zone' of tolerance (the 'just right' predation level needed to drive selection for predator-resistant traits, Evans et al. 2021) in the area surrounding a conservation fence could deliver a great return on investment by protecting migrants and aiding in their establishment 'beyond-the-fence' (i.e., spill-over or halo effect, Tanentzap and Lloyd 2017). Finally, to prevent the "locking-in" of the degraded shifting baseline (where native species vulnerable to introduced predators are accepted as permanently absent from the wild), we must explore innovative solutions to drive or enable adaptive evolution of threatened species and introduced predators alike (i.e., 'coexistence conservation', Evans et al. 2022).

Conclusion
Here we demonstrated how demographic parameters from a reintroduced population can inform management interventions and create targets for delisting an endangered species. We also highlighted the value of conducting ecological experiments within conservationfenced sanctuaries. If treated as 'outdoor laboratories', these areas provide unique opportunities to measure vital rates in free-ranging, endangered species when it would otherwise be difficult or impossible (Hester et al. 2000;Manning et al. 2009). Ironically, these populations may be better understood than the extant populations (e.g., there are no equivalent robust eastern quoll population estimates for Tasmania, Burbidge and Woinarski 2016) for which knowledge of their population dynamics would greatly assist conservation efforts (Ashbrook et al. 2016). The current extinction crisis highlights the need for managers to use evidence and collaboration to orient decisions and contribute lasting progress towards species recovery.