Assessing extinction risk, conservation reliance, and down-listing potential for two endangered Hawaiian waterbirds

. Conservation reliant species complicate the idealized view of endangered species recovery by depending on perpetual management to maintain viable populations, making the goal of de-listing via threat removal potentially unattainable. The Hawaiian Stilt ( Himantopus mexicanus knudseni ) and Hawaiian Coot ( Fulica alai ) are conservation reliant, endangered waterbirds endemic to the Hawaiian Islands. The former is currently under consideration for down-listing to “threatened” under the United States Endangered Species Act. Down-listing criteria for both species include a low probability of extinction according to population viability analysis and a sustained population of over 2000 individuals. We investigated the population viability, equilibrium population size, and degree of conservation reliance of these two species using stochastic individual-based population models and Bayesian state-space analyses on abundance time series. For each species, we evaluated the probability of extinction (one sex remaining) and quasi-extinction (N < 50) by the year 2100, as well as the probability that range-wide equilibrium population size (functional carrying capacity, given current wetland area, configuration, and management) > 2000 birds. We found a low probability of extinction (0% and 11% for stilts and coots, respectively) or quasi-extinction (0% and 16%, respectively). However, sensitivity analysis of simulated populations showed that many vital rates in both species are close to threshold values, across which extinction risk increases dramatically. Cessation of management in major population strongholds would negatively affect several of these parameters, likely resulting in declines. Additionally, the posterior probability estimates of equilibrium population size for the top-performing population models of both species showed that the probability of meeting the population size down-listing criterion was 0.03 for Hawaiian stilts and 0.06 for coots. Our findings suggest that these two species are conservation reliant and not candidates for down-listing


INTRODUCTION
An idealized view of endangered species management is that vulnerable species are investigated, their threats identified, and management is enacted to eliminate or mitigate the threats such that the species is no longer endangered. This scenario is central to the application of many endangered species laws, such as the U.S. Endangered Species Act (Goble 2009, Neel et al. 2012, but it presumes that the threats to a species can be eliminated. For some species, however, threat elimination is not possible and continued management is required to prevent species decline toward extinction (Scott et al. 2005, Scott et al. 2010Goble et al. 2012). Such conservation reliance (also known as conservation dependence) is increasingly widespread (Averill-Murray et al. 2012, Wiens and Gardali 2013, Wilsey et al. 2014, Brown et al. 2017, Connor et al. 2019, Scott et al. 2020, and has been defined for the IUCN Green Status of Species assessments as the shortterm (10 years or 3 generations, whichever is longer) reduction in recovery that would occur if management of a population were to cease (Akçakaya et al. 2018).
Conservation reliance is often described as a binary condition, but more realistically exists on a continuum. Reed et al. (2012) categorized 33 endangered Hawaiian bird species into five groups within this spectrum based on their degree of conservation reliance. These groups ranged from those requiring "modest, often intermittent" or minimal management, to those species that were considered "possibly too late" to recover; only one species was not considered conservation reliant. Developing a system to quantify the degree of conservation reliance would improve prioritization of management actions within and across species following delisting (Bottrill et al. 2008, Brown et al. 2017; this is one function of the IUCN Green Status of Species assessments (Akçakaya et al. 2018). Most assessments of population viability focus on metrics such as the likelihood of persistence for some targeted period of time (Doak et al. 1994, Walters et al. 2002, Traill et al. 2007, with sensitivity analyses identifying combinations of vital rates that would allow the most rapid population growth (Beissinger and Westphal 1998, Mills et al. 2001, Fefferman and Reed 2006, Aiello-Lammens and Akçakaya 2017. The difference between the current value for a vital rate, if the population has a high probability of persistence, and the point at which a population's probability of extinction increases to high (>0.8) levels, making the likelihood of persistence low, may be particularly useful for evaluating conservation reliance (e.g., Mortensen andReed 2016, van Rees and. This is one measure of a population's resistance (Pimm 1984) to extinction, and proximity to this point is relevant to management needs prior to and following delisting. Similarly, Hutchings (2015) identifies population sizes at which species have impaired positive responses to management efforts (an Alleeeffects threshold), and de Silva and Leimgruber (2019) identified vital rate values to look for a tipping point between positive and negative population growth. The effort required to keep vital rates above the threshold between increasing and decreasing could be a useful metric of the degree of conservation reliance.
The catastrophic loss of endemic island species is a global problem largely driven by introduced species, overharvest, and habitat loss (e.g., Milberg andTyrberg 1993, Wood et al. 2017). Although efforts to remove introduced species have grown, and are increasingly successful (e.g., Russell and Holmes 2015), in many cases a return to historical conditions is not possible (Scott et al. 2020). The Hawaiian avifauna is an excellent system to explore species loss and recovery of endangered species in light of these problems (Boyer 2008, Pratt et al. 2009) because it has lost 71 of 113 endemic species, and most that remain are considered globally threatened (Leonard 2008). The extant Hawaiian waterbirds, by contrast, are generally considered to be conservation success stories because management has not only slowed declines but resulted in population increases (Reed et al. 2011, Underwood et al. 2013) and subsequent calls for down-listing some species from endangered to threatened status under the U.S. Endangered Species Act. Nonetheless, concerns remain that any reduced protections could leave these species vulnerable to the same risks that caused the original declines, resulting in the loss of hardearned progress.
Here we use population viability analysis (PVA), perturbation sensitivity analysis, and Bayesian state-space modeling to estimate current extinction risk, the potential to achieve specific downlisting criteria, and degree of conservation reliance for two Hawaiian waterbird taxa with contrasting life histories and levels of scientific understanding (for convenience, we refer to both taxa as "species" following the legal definition embodied in the Endangered Species Act.) The Hawaiian Stilt (Himantopus mexicanus knudseni), which is currently being considered for down-listing to threatened, has a small clutch size, is relatively long-lived, and has a long history of detailed demographic study (Reed et al. 1998a, Pratt and Brisbin 2020. In contrast, the Hawaiian Coot (Fulica alai) has a large clutch size, shorter lifespan, and is poorly studied, with little published information on vital rates (Pratt and Brisbin 2020). Both are considered to be conservation reliant because of our current inability to remove threats caused by exotic invasive predators and plants that make habitat unsuitable (USFWS 2011, Reed et al. 2012, Underwood et al. 2013. For either species to be down-listed from endangered to threatened, three major criteria must be satisfied: (A) a PVA must be conducted to determine the population size necessary for longterm viability, and (B) statewide surveys must show a stable or increasing trend with (C) population sizes above 2000 adult birds for at least 5 consecutive years (USFWS 2011). The last criterion was set in the original recovery plan for the species (USFWS 1978) through expert opinion (Ron Walker, lead author of the original plan, pers. comm. to JMR). Accordingly, the goals of this study were to: 1) quantitatively assess the probability of the statewide populations of either species being sustained at levels at or above the down-listing threshold of 2000 individuals, 2) conduct individual-based PVAs for both species to assess their long-term viability, and 3) use sensitivity analysis of vital rate values within the PVA model to estimate the degree of conservation reliance in these species.

Study system
The Hawaiian Stilt and Hawaiian Coot are year-round residents of coastal wetlands, and they are endangered because of wetland loss, invasion by exotic plants that alter their habitat, and introduced predators and competitors (Coleman 1981, Griffin et al. 1989, USFWS 2011. One result of human activities on the Hawaiian Islands has been a reduction in wetland cover and degradation in habitat quality (Shallenberger 1977, Griffin et al. 1989, van Rees and Reed 2014. Habitat degradation and pressures from anthropogenic impacts like altered hydrology and invasive species require frequent and continuous management action. Management in the form of flooding, mowing, or predator trapping takes place on the scale of days to weeks to adjust to rapidly changing water levels, spreading invasive emergent vegetation, and incursions by invasive predators. Both Hawaiian Stilts and Coots are present on the seven largest Hawaiian Islands (Hawaiʻi, Maui, Molokaʻi, Oʻahu, Kauaiʻi, Lanaʻi, Kahoʻolawe;Pratt 1993, Pratt and, but are absent from smaller islands which have scarce rainfall and little surface water to support wetlands. Higher numbers of both species are reported now than during the 1960s, with annual counts of coots occasionally exceeding 2000 individuals in the last 30 years (Fig. 1;Banko 1987, Engilis and Pratt 1993, Reed et al. 2011 (Reed et al. 2011, USFWS 2011, Camp et al. 2016), so we analyzed long-term time series of statewide abundances using a Bayesian state-space model (Delean et al. 2013, van Rees et al. 2020) to estimate range-wide equilibrium population sizes for the two species and determine if this down-listing criterion is currently achievable. We used the posterior probability distribution of parameter estimates of equilibrium population size to assess whether current habitats within the state of Hawaiʻi are capable of sustaining population levels above the down-listing threshold (goal 1 of this study). This method provides a conservative quantitative metric of the populations' achievement of the down-listing criterion of 2000 individuals by measuring the probability that this abundance could theoretically be achieved. We assessed the long-term viability and conservation reliance of both species' populations using PVA methods outlined in later sections.
We used a Bayesian framework for model-fitting because it has the advantage of producing easily communicable results that can be described using real probabilities derived from the posterior probability distributions. Results presented as actual probabilities are easier to relay to decision-makers for endangered species listings. Bayesian approaches also provide the opportunity to use a priori knowledge about other aspects of a species' biology to guide parameter estimation, such that empirical data on betterstudied parameters can facilitate the estimation of an unknown parameter (Hooten and Hobbs 2015). This approach is broadly applicable to recovery planning where there are sufficient timeseries data, even if vital rate information is missing.
We used Bayesian state-space models to estimate the statewide equilibrium population size for each species with the specific interest of evaluating whether current population dynamics were consistent with an equilibrium population size of at least 2000 individuals. We parameterized these models using time series of abundance data from summer and winter counts of the state of Hawaii's biannual waterbird survey ( Fig. 1), following the approach used by van Rees et al. (2020) for the Hawaiian Stilt. Biannual waterbird survey data consist of both summer (August) and winter (January) counts along consistent routes in all known major habitats for the subspecies (Engilis and Pratt 1993, Camp et al. 2014, Hawaii Division of Forestry and Wildlife 2017. We conducted all analyses separately for winter and summer data. Bayesian state-space models are considered among the best approaches to estimating parameters such as the intrinsic rate of population growth (r max ) and equilibrium population size (sometimes treated as extant carrying capacity, K) using timeseries datasets (Delean et al. 2013). Such analyses explicitly account for detection error which can otherwise bias model outcomes (Dennis et al. 2006, Freckleton et al. 2006, Ahrestani et al. 2013. Van Rees et al. (2020) found that state-space population models did not perform better than non-state-space models for the Hawaiian Stilt, but detection rates may be lower for the Hawaiian Coot, which can occupy more dense habitats than the stilt (Chang 1990). We thus employed both state-space and non-state-space versions of all our nonlinear population models for estimating equilibrium population size to evaluate if fit differed between model types.
For both species, we fit Bayesian state-space and non-state-space versions of five population models (exponential or density independent, Ricker, theta-logistic, Gompertz, and a conditional model; Table A1) as well as one-, two-, and three-year time-lagged versions of the Ricker and conditional models, and compared them using leave-one-out cross-validation (cf., Hooten and Hobbs 2015).
The conditional model was designed to exhibit no densitydependent impacts on population growth (i.e., exponential growth) until a threshold population size was passed, and then Ricker density-dependent dynamics after that (as in van Rees et al. 2020). We selected the level of this threshold by comparing frequentist versions of conditional models for a range of threshold values via AICc and used the threshold value for which AICc was minimized (van Rees et al. 2020). An alternative approach would be to set the population threshold as an additional parameter in the model to be explored via Markovchain Monte Carlo (MCMC) sampling, but models did not converge with this additional parameter. We used the posterior probability distributions of the top (best fit) model for each taxon to estimate the probability that the K parameter (equilibrium population size) was at least 2000. All models were fit using MCMC sampling in JAGS (version 4.3.0; Plummer 2003) using the packages Runjags (Denwood 2016), and jagsUI (Plummer 2018, Kellner et al. 2019) in R (version 3.5.1; R Core Team 2018).
We used vague priors for K (uniform, bounded at 1000-3000 individuals) to achieve an unbiased estimate of carrying capacity under current conditions (equilibrium population size) for both species, and loosely informative priors for observation and process variation (Table A2) to take advantage of prior knowledge of the species' biology from independent data sources. The minimum and maximum values for K were based on observed population fluctuations over the last 30 years. For both taxa, priors for r max were defined by a normal distribution with mean and standard deviation determined by a deterministic demographic model in Vortex (described in Lacy et al. 2017, van Rees et al. 2020, Table  A1). This amounted to running a deterministic version of the model with all parameters set to the maximum observed values for all vital rates. This simulation generated a realistic estimate of a biologically possible maximum growth rate for each species around which to base prior distributions for r max in Bayesian population models. Importantly, this approach served to discourage the search algorithm from using biologically unreasonable values for this parameter.
Priors for process variation in population growth rate and observation error were defined using gamma distributions, and shape parameters were selected to ensure that realistic and observed values for both parameters in both species were within the given distribution, but that larger nonzero values were possible and would be searched by the MCMC algorithm. All priors and initial values used to conduct MCMC sampling are shown in Table A2. We based process variation in growth rate for the stilt on variance in population growth rate from Reed et al. (1998a).
Because we lacked such prior information for coots we used priors that were higher and less informative (more spread) to reflect higher observed apparent variation across time and greater uncertainty for that species. Parameters of prior distributions for observation error (i.e., detection error) for both species were based on plausible detection error from Chang (1990), again with the possibility of much higher errors for coots (Table A2). These were defined using a gamma distribution, which allowed for nonnegative values (i.e., only allowing undercounting), the shape of which was determined using a graphing calculator and values from Chang (1990) for the order of magnitude of likely values. Unlike the Hawaiian Gallinule (Gallinula galeata sandvicensis), detection error in both species is likely to be very low because of their tendency to occupy more open parts of wetland habitats, whereas stilts use open sparsely-vegetated areas and coots are often found in open water.

Goal 2: Baseline population viability analyses
We used an individual-based simulation model conducted in Vortex 10 (Lacy and Pollak 2014) to address goals (2) and (3) of this project: assess the long-term viability of the two focal species and estimate the strength of conservation reliance by testing the effects of changing parameter values on predicted viability. Using Vortex models, we generated population projections for both species until 2100, an approximately 80-year time frame similar to that used in van Rees and  for the Hawaiian Gallinule. Extensive inter-island movements by marked Hawaiian Stilts and Hawaiian Coots (Udvardy 1960, Banko 1988, Reed et al. 1998b, Dibben-Young 2010, Riggs 2016) and high levels of gene flow in the latter (Sonsthagen et al. 2018) suggest little population subdivision, so we modeled both species as single statewide populations.
We evaluated two persistence criteria: probability of extinction, which was defined as when individuals of only one sex remained, and the quasi-extinction probability (e.g., Holmes et al. 2007), estimated as the probability of the population size dropping below 50 adults. We selected the latter criterion because it is generally considered to be small enough that a population is at significant risk of extinction due to demographic stochasticity (Burgman et al. 1993, Lande 1998. Baseline models were run for 1000 replicates; initial population sizes were set from recent biannual waterbird survey data at 1700 for the stilt and 1280 for the coot, starting with an even sex ratio and stable age distribution. We evaluated viability from these models based on (1) probability of extinction by the year 2100, (2) mean population size by the year 2100, and (3) the mean stochastic growth rate estimated in Vortex. We assumed no correlation between reproduction and survival for either species because nest or chick loss rarely involves adult mortality.
The reproductive systems of both taxa were modeled as long-term monogamous based on current knowledge and observations of the species (Byrd et al. 1985, Reed et al. 1998a. Although cooperative breeding has been observed in the stilt (A. Dibben-Young pers. comm.), this was only in a single small group of birds in a highly habitat-limited context, and we have no evidence that it is a sufficiently widespread phenomenon to affect population dynamics. For Hawaiian Stilts, we set a maximum age of 30 based on a large mark-recapture dataset  and the maximum age of reproduction for both sexes at 30 because there is no evidence of reproductive senescence (Table 2). We took all other reproductive parameter estimates for Hawaiian Stilts from Reed et al. (1998a). For Hawaiian Coots, we used an age of first breeding of 1 given anecdotal observations (M. Silbernagle, pers. comm.) and data from the American Coot Fulica americana (Brisbin and Mowbray 2020) and assigned a maximum lifespan of 9 based on a mark-recapture dataset assembled and analyzed in this study. We assumed a 1:1 sex ratio at hatch for both species.
We followed the approach of van Rees and Reed (2018) for implementing density-dependent population dynamics in models for both species, setting the Allee effect parameter to zero in the absence of evidence for such effects, and setting the slope parameter to 20 so that density-dependent impacts were not introduced until population size was above 80% of the model's population ceiling parameter, called carrying capacity K in Vortex, but distinct from the equilibrium population size estimated in our Bayesian time-series models. We assumed density dependence would arise based on territorial dynamics in Hawaiian Coots (Riggs 2016) and recent evidence from Hawaiian Stilts (van Rees et al. 2020). We assumed no mate monopolization and that all males of breeding age were in the breeding pool.
For Hawaiian Stilts, we used the reproduction dataset from Reed et al. (1998a) to define the distribution of brood sizes (number of hatched young) and number of broods per year. Lacking similar information on the Hawaiian Coot we followed the methods of van Rees and  to generate a distribution of brood frequencies using observed nest-failure rates from Chang (1990) and the binomial equation, assuming up to four breeding attempts per breeding season irrespective of fate (M. Silbernagle, pers. comm.). We specified a distribution of brood sizes based on data from Chang (1990).
We estimated mean and standard deviation of juvenile (age 0-1) and first-year (age 1-2) mortality for stilts using draws from the posterior probability distribution produced by Reed et al. (2015), and modeled percent mortality for each subsequent year until age 30 according to the best fit line for annual adult survival from the same study: 22.4 -(0.40*(A-2)), where A is bird age, and the standard deviation of the posterior probability distribution of adult survival was used to account for uncertainty and natural variability in annual survival.
We estimated mean annual apparent survival of adult Hawaiian Coots using standard Cormack-Jolly-Seber (CJS) models in program MARK (version 8.2, White and Burnham 1999) and mark-resight data consisting of >1000 resightings of 596 individuals over 12 years (A. Dibben-Young, unpublished data). Years were defined by calendar date for computational convenience. We compared the likelihood of a suite of potential survival models using AICc: (1) a null model with constant survival (φ) and detection (p); (2) a model with constant survival, but variable detection among years; (3) a model where survival varied by wetland, but detection was constant; (4) a model where both survival and detection varied among years; two models where survival varied by wetland, and detection either (5) varied by year, or (6) varied by wetland; and four models where detection varied by wetland, year, and their interaction and survival was either (7) constant, or varied by (8) year, (9) wetland, or (10) year, wetland, and their interaction). We implemented "goodness-offit" tests using program U-CARE (Choquet et al. 2009).
Given the limited information on within-year variation in adult survival for Hawaiian Coots, we used an estimate for American Coots, which has nearly identical values for apparent adult survival (Arnold et al. 2016). We estimated first-year survival rates from Chang (1990) and Byrd et al. (1985). Because we have little to no information on variance in juvenile survival for Hawaiian Coots we chose a value that resulted in population dynamics qualitatively similar to those observed in recent decades using the trajectories of simulated populations.
We set the population ceiling parameter K for Hawaiian Stilts as the equilibrium population size estimated by the top time-series population model from van Rees et al. (2020). For Hawaiian Coots, however, there is high inter-annual variation in estimated population size (Engilis and Pratt 1993), so we set this parameter at a higher level than the estimated carrying capacity and used the highest observed population size since 1990 (Table 2). We adjusted the shape parameter B in Vortex's density dependence equation such that density-dependent feedback would be exhibited at population sizes approaching actual equilibrium size. This allowed simulated Hawaiian Coot populations to fluctuate similarly to observed statewide counts over the last several decades.

Goal 3: Conservation reliance using PVA Sensitivity analysis
We assessed the effects of changes to vital rate parameters, a proxy for different hypothetical levels of management, on the outcomes of PVAs for both species to provide a measure of conservation reliance (goal 3). While few data exist to compare vital rates between managed and unmanaged populations of either species, a general examination of the vital rate values linked to elevated extinction risk, and the proximity of those values to present estimates for our target species, provides a means of assessing conservation reliance. Examining this relationship in conjunction with the associated uncertainty in parameter values due to measurement error and natural variation provides further guidance on where critical management-relevant information is missing and where research effort should be directed.
We tested the influence of selected vital rate parameter values on model output using two approaches. First, we employed logistic regression to assess the relative contribution of each vital rate parameter to the probability of extinction. We used the Latin hypercube sampling method in Vortex with 10 iterations each for 1000 combinations of parameter values sampled across the range of observed values (Mortensen and Reed, 2016; Table 1). All parameters were varied simultaneously from random distributions based on the ranges provided (Table 2) and a logistic regression was performed relating these values as covariates to the probability of extinction estimated for each simulation run. This method broadly assessed the relative importance of different vital rate parameters to population viability.
Second, we used perturbation analysis to find vital rate parameter values at which large, abrupt increases in the probability of extinction occurred (Reed et al. 1998a, Reed et al. 2011, to highlight vital rates that might be driving conservation reliance. For each perturbation analysis, we held values of all parameters not being assessed at baseline levels and compared the probability of extinction for each simulation. Perturbation analysis is useful in detecting ranges of parameter values where the probability of extinction shifts from being low to high (or vice-versa), providing minimum vital rate values to be maintained by management. We perturbed values of the mean and standard deviation of juvenile survival, adult survival, nestfailure rate, mean brood size, proportion of females breeding at low population densities, and carrying capacity (values and intervals in Table 1). We also perturbed sex ratio at hatch for coots to examine the potential impacts of a population with 69% males, which has been reported on O'ahu (Riggs 2016).
We used the proximity of empirical parameter estimates for vital rates based on field data, and their margin of uncertainty, to observed parameter values with high (≥80%) probability of extinction from our perturbation analysis as a quantitative test for conservation reliance. Because most of the vital rate information collected for the species in this study is derived from field studies on managed populations, we judged this comparison as a reasonable heuristic for degree of conservation reliance. Specifically, where the confidence interval around observed values for a given vital rate encompassed a high extinction risk value, or was 10% of its present value to that point, and if that value was likely to be affected by cessation of management, we judged a given trait as indicative of conservation reliance. If management were to cease, that vital rate would likely be pushed to a value that caused a dramatic shift in population size.

Goal 1: Estimating equilibrium population density using Bayesian time-series models
Among Bayesian population time-series models, the top performing coot model was the non-state-space version of the conditional model using summer data with a one-year lag ( Table  2). The next ranked model was the non-state-space version of the Ricker model with a two-year lag run on winter data, which also had a large difference in cross-validation value from the third best model. Among the top five Bayesian models assessed for Hawaiian Coots, four of which were non-state-space (did not include a detection parameter), four were density-dependent models, and four used summer count data. In both cases, the probability that carrying capacity is at or above the down-listing criterion is ≤6% for the top-performing model † .
Parameter estimates with 95% confidence intervals are shown, where r max , σ proc , K, and σ obs are the annual population growth rate, variation in population growth rate (i.e., process variation), carrying capacity, and detection error, respectively.
Model Name Cross-Validation Statistic

P(K ≥ 2000)
Parameter estimates (95% CI) The best fit model for both species indicated a probability of equilibrium population size (carrying capacity under current conditions) at or above 2000 of 0.06 or lower ( Table 2). The posterior probability distributions of the top five models for the Hawaiian Coot indicated that the probability of the equilibrium population being above the goal of 2000 was between 0.04-0.11, while the probability for Hawaiian Stilts lay between 0.03-0.26 (Table 2).

Goal 2: Baseline population viability analyses
Our baseline population models yielded low probabilities of extinction (0.0 and 0.11) and quasi-extinction (0.0 and 0.16), for Hawaiian Stilt and Hawaiian Coot, respectively. Both species showed similar mean stochastic population growth rates (0.11 and 0.094, respectively). The Hawaiian Stilt population model resulted in no extinctions. For the small number of simulation runs ending in extinction for the Hawaiian Coot, the mean time to extinction was 53.5 years (95% CI: 28.5-78.5); for quasiextinction, it was 40 years (95% CI: 22-73). The mean ending population size at year 80 for the Hawaiian Stilt was 1588, with the entire 95% confidence interval (1321-1855) below the downlisting goal of 2000 individuals. The mean for coots was 1617, also below the down-listing goal, but the confidence interval (0.0-3639) was much broader, encompassing both extinction and the target population size. The mean stochastic growth rate for Hawaiian Stilts was 0.11 (95% CI: -0.27-0.49), and for Hawaiian Coots, it was 0.09 (95% CI: -1.03-1.21). The maximum deterministic growth rate (r max ) based on the highest observed reproduction parameter values was 0.52 for stilts and 1.2 for coots.

Logistic regression
Sensitivity analysis by logistic regression showed mean adult mortality, percentage of adult females in the breeding pool, mean nest-failure rate, and mean juvenile mortality as the most important parameters affecting probability of extinction for Hawaiian Stilts. In contrast, the top four parameters for Hawaiian Coots all related to survival (Table 3).

Perturbation
For both species, sharp increases between low and high probabilities of extinction were identified for most vital rate values (Figs. 1 and 2; Table 4). In Hawaiian Coots, probability of extinction increased abruptly from 0 to near 1 when annual adult mortality exceeded approximately 60%. Uncertainty around this parameter estimate from observations in managed populations encompassed probability of extinction values up to 1. Probability of extinction also transitioned rapidly from 0 to 1 for mean juvenile mortality rates above 74%; no data on the uncertainty of this estimate were available, but the estimate derived from the literature was close to the point of rapid increase in probability of extinction. Nest failure rates above 20% showed a steady increase in probability of extinction, reaching 0.5 at a 34% failure rate and 1 at around 75%. The range of empirical estimates for nest-failure rate produced probabilities of extinction from 0.25-1, with the mean parameter estimate from field observations between 0.9-1. Probability of extinction for Hawaiian Coots was 1 for mean brood sizes below 2.6 chicks and declined to 0 by 4.4 chicks. Observed variation in this parameter value encompassed a full range of probabilities of extinction from 0-1.  (Chang 1990, A) and a national historical park (Morin 1998, B).
For Hawaiian Coots, the probability of extinction increased linearly with variation in annual adult mortality, reaching 1 at a standard deviation of around 25% variation in annual mortality (Fig. A1). Variation in juvenile annual mortality showed a similar trend, with extinction probability increasing with values around 15% and reaching a maximum between 30% and 35% (Fig. A2). Probabilities of extinction decreased rapidly from 1 to near 0 at K above 250 individuals, showing a persistent small, but nonzero, probability of extinction when K exceeded 1000 individuals (Fig.  A3).
Increasing the proportion of male coots at hatch generated a gradual increase in the probability of extinction, approaching 1 when there were >70% of males at hatch (Fig. A4), a value that is similar to the sex bias observed by Riggs (2016) on O'ahu. Probability of extinction was 1 for populations with percentages of females reproductively active below 50%, and declined for higher values, dropping to near 0 at 95% of adult females in the breeding pool (Fig. A5).
For Hawaiian Stilts, mean annual adult mortality rates above 20% led to increased probabilities of extinction, reaching a value of 1 at ~41% mean annual mortality ( Fig. 3; Table 4). Current uncertainty around mean annual adult mortality rates in managed areas corresponds with a probability of extinction below 0.15. Mean annual juvenile mortality showed an increase in probability of extinction for values above 40%, with a rapid rise in probability of extinction from 55-78%, reaching 1 at 79% (Fig. 3). The larger margin of uncertainty around parameter estimates for this vital rate in managed areas showed that the higher end of this estimate corresponded with probabilities of extinction on the order of 0.35-0.5. Simulated Hawaiian Stilt populations showed a rapid decline in probability of extinction for mean brood sizes averaging greater than 0.8 chicks. Probability of extinction dropped to zero above a mean brood size of 1.5. Uncertainty and natural variation in observed brood size from managed populations included probabilities of extinction from 0-1. Nest-failure rates above 50% drove a rapid increase in probability of extinction, reaching 1 at failure rates beyond 75%. The empirical parameter estimates for this value, based on observations in managed habitats, was below this point of increase in probability of extinction, corresponding with a probability of 0.0-0.5.
The variation in mean annual adult mortality for Hawaiian Stilts showed a sharp increase in quasi-extinction risk after 20%, reaching a maximum by 45% ( Fig. A1; Table 4). Variation in juvenile mortality, in contrast, did not affect probability of quasiextinction at values below 80% (Fig. A2). Probability of extinction dropped rapidly for simulations with carrying capacities above 550 individuals, approaching 0 as carrying capacities exceeded 1600-2000 individuals (Fig. A3). For the percentage of adult females in the breeding pool, probability of extinction decreased rapidly above 25%, slowing and approaching zero at around 75% (Fig. A5).  Fig. 3. Results from perturbation sensitivity analysis for the Hawaiian Stilt, displaying the relationship between probability of quasi-extinction (PE) and four vital rate parameters that are especially relevant to management. Dashed lines indicate the present best estimate of these vital rates based on existing field data, and shaded regions indicate this mean ± 1 standard deviation where this information was available.
Importantly, the increases in extinction risk associated with varying vital rate values for mean brood size, mean annual juvenile mortality, and nest-failure rates for Hawaiian Coots were all near to, or within the range of uncertainty of, current estimates for those parameters (Fig. 2; Table 5). This was also true for juvenile mortality and mean brood size in Hawaiian Stilts, but not for other vital rates for that species (Fig. 3; Table 4).

DISCUSSION
Using the best available data on vital rates and individual-based population models we find that both the Hawaiian Stilt and Hawaiian Coot have a high likelihood of persisting to the end of the century under current habitat and management conditions, and assuming no major loss of habitat, as might be expected from development or sea level rise. However, almost all vital rate data used in population viability models came from populations that were managed to control predators and maintain suitable nesting conditions. Moreover, there are few individuals and no reproductive strongholds of these species in unprotected areas. Although few empirical data exist from unmanaged populations, the precipitous declines of these two species prior to management and widespread field surveys, and their coincidence of population recovery with the establishment of protected areas (Reed et al. 2011), strongly suggest that cessation of management would reduce the vital rates identified in this study as most important to population persistence in both species: nest success and both adult and juvenile survival rates. We anticipate that the vital rates used in our models are unlikely to be maintained with reduced management, and our analyses of conservation reliance through the perturbation of individual vital rate values suggest that risk of extinction could increase dramatically if reduced management resulted in less favorable rates.
Although both mortality and reproductive parameters were important for both coots and stilts, their relative importance differed. For Hawaiian Coots, likelihood of extinction was more sensitive to changes in mean and variation in mortality of adults and juveniles than to any reproductive parameters. In contrast, although adult mortality was the most influential variable in extinction risk for Hawaiian Stilts, reproductive parametersespecially the proportion of females nesting and nest-failure rate -ranked highly. This difference between species aligns with life-history differences, with coots having shorter life spans and a much larger reproductive potential. Management activities that focus on both adult and juvenile survival are thus likely to be most important for coots, while activities that focus on maximizing the number of successful nests, both by ensuring that females breed and that their nests survive, will be more important for stilts. Predator management, through removal and exclusion, will likely benefit all of these variables.
The confidence intervals around vital rate parameter estimates for the stilt are considerably smaller than those for the coot due to more abundant data and lower natural variation, but in the case of adult and juvenile mortality as well as brood size, the boundaries of these intervals contain or abut values associated with a high probability of extinction. The proximity of plausible vital rate values in both species to projected high-extinctionprobability states suggests that even small shifts in some vital rates would significantly increase extinction risk. Even at some managed sites, vital rate estimates based on field observations fall within the range of values associated with high probability of extinction (e.g., Morin 1998). Although we do not have such contrasting data for most vital rates, field observations of these species and published findings for an ecologically related species (Hawaiian Gallinule: Nagata 1983 suggest that reduced management for Hawaiian Stilts and Hawaiian Coots would have harmful impacts and easily push some of these parameters toward high-risk values. This concern is supported by a few studies outside protected areas where vital rates are reduced for the Hawaiian Gallinule, implying that analogous risks would be incurred with the cessation of invasive plant management and predator control. The generalizable nature of our perturbation sensitivity analysis of vital rate values provides a quantitative approach to assessing conservation reliance in data-limited species. Particularly where empirical data comparing managed to unmanaged populations are scarce, this approach helps distill relevant information on current estimates of parameter values, uncertainty in those values due to observer variation and measurement error, and the proximity of plausible values to levels at which long-term population viability declines below acceptable levels. Although this is no replacement for quantitative comparisons using baseline data, controls, and replicated observations of wild populations, it provides a formal and organized way to use existing knowledge to inform decision-making with a quantitative illustration of potential risk, especially in situations where key information is impossible to collect without conducting questionable experiments (e.g., by removing protections to test their effects).
Our assessments of extinction risk are likely to be conservative because, in our simulations, only one vital rate needs to change for the worse for a species to decline to extinction, but reduced management would likely impact multiple rates simultaneously (Griffin et al. 1989). Reduced exotic invasive plant control, for example, would rapidly eliminate available habitat (Cox 1999, Bantilan-Smith et al. 2009), reducing number of nests, nest success, and average brood sizes (via competition for space). Similarly, unmanaged mammalian predators would simultaneously impact nest-failure rate, brood size, juvenile survival, and adult survival (Coleman 1981, Ohashi and Oldenburg 1992, Underwood et al. 2014. Thus, although predicted probability of extinction for both species is low, the evidence suggests that this prognosis is precarious and subject to dramatic change in the absence of management. To view this conclusion in a positive light, these two management strategies simultaneously affect a suite of vital rates, and specific management techniques are not needed to address these rates separately. In the absence of empirical data on the impact of reduced or no management on our focal taxa, data on other species may be informative regarding the potential impacts on vital rates. For example, Schüttler et al. (2009) found nest-failure rates of 80-95% due to unmanaged invasive predators in solitary, ground-nesting waterbirds on an island off Chile. Nest failure rates of this magnitude would lead to a high probability of extinction for both of our modeled taxa, even without additional effects on juvenile and adult survival. A recent review of the efficacy of predator management for bird conservation ranked the removal and control of mammalian predators on islands, generally to protect nesting bird communities, as having the strongest empirical evidence for efficacy and conservation benefit (Sutherland et al. 2019), implying strong support for the need to continue management for the focal taxa in our study.
We did not include catastrophic decreases in vital rates in our viability assessment, although it is known that catastrophes can increase extinction risk (Mangel andTier 1994, Dale et al. 1998 (Coleman 1981, Fletcher et al. 2003, Work et al. 2010), such events are not known to have catastrophic effects on our target species.
Our predictions assume maintenance of an ecological status quo for both species, with both continued management and no largescale systematic environmental changes. We know, however, that Hawaiian wetlands are likely to be strongly affected by climate change and sea level rise (Kane et al. 2015, Htun et al. 2016. Like Hawaiian Gallinules, which are typically found only in areas of low salinity, Hawaiian Coots are found primarily in fresh water, but appear to better tolerate brackish water (Byrd et al. 1985, Engilis andPratt 1993). Although it has not been studied in Hawaii, American Coots show reduced breeding density and reproductive success in brackish water (Kantrud in Allen 1985). In contrast, Hawaiian Stilts use wetlands with a wide variety of salinities, from fresh water to hypersaline (Coleman 1981, Engilis andPratt 1993). Their apparent tolerance for increased salinity may buffer the potentially negative impacts of saltwater intrusion on coastal wetland habitats under sea level rise, given that new habitats may be created, and salinized wetlands may not necessarily become uninhabitable to both species. Shorebirds as a group, however, appear to require access to fresh water, particularly for chicks (Rubega and Robinson 1997), so this should be investigated for Hawaiian Stilt management. Accordingly, the type of simple analyses conducted by van Rees and  to estimate loss in available habitats for Hawaiian Gallinules would not be as informative for these species, for which population responses are likely to be more nuanced, requiring hydrological modeling for changes in habitat abundance and quality. Such research should be considered a high priority and perhaps an additional criterion for down-listing of these species, as sea level rise is a potentially critical factor affecting extinction risk (Underwood et al. 2013, Kane et al. 2015. The degree to which individual wetlands of importance for these species have a clear topographic pathway for inland and upland migration under sea level rise will be a major factor in such analyses. Beyond population viability, we also used a Bayesian hierarchical modeling approach to estimate equilibrium population size, which we used as an assessment of whether the down-listing criterion of 2000 individuals is attainable for either species. Given the low probability that the equilibrium population size is above 2000 individuals in either species, our analysis suggests that more habitats will have to be protected and managed to achieve this goal. Although both species use coastal wetlands, they differ significantly in their microhabitat affiliations, so managing for one species would reduce habitat for the other. Many protected wetlands are not currently managed (e.g., large portions of Kawainui marsh on O'ahu) and they could likely support much higher populations with some investment.
One possibility for supplementing protected wetlands for waterbirds is through agricultural wetlands (van Rees 2018, Harmon et al. 2021). Taro lo'i is a traditional agricultural system that is readily used for nesting or foraging by both species (Shallenberger 1977, USFWS 1978, Stone et al. 1989, Rader 2005, Gee 2007). It is possible that restoration of Indigenous land practices such as taro lo'i, which create large tracts of artificial wetland, could support waterbird populations currently reliant upon managed areas (Winter et al. 2020). Little is known, however, about how to manage taro lo'i to meet the needs of any of the Hawaiian waterbirds while also meeting agricultural goals. The potential benefits of jointly managing taro lo'i for waterbirds and taro harvest suggests the need for investigation (Hawaii DLNR 2005). Ideally, such work would involve replicated, controlled experimental assessments (Macnab 1983), similar to studies of rice field management during the non-growing season to support waterbirds in other parts of the world (e.g., Elphick et al. 2010). A key question would be to determine whether there is a risk that taro fields could become ecological traps (Robertson and Hutto 2006).
A complementary alternative method for creating new waterbird habitats and increasing available habitat statewide is via the construction and ecological maintenance of freshwater and coastal natural infrastructure (Da Silva and Wheeler 2017). The restoration and creation of wetland ecosystems is increasingly recognized as an option for improving water security and climate resilience (Thorslund et al. 2017), and strategically locating mitigation or restoration sites could lead to ecological benefits, including increased waterbird habitats, connectivity between habitats, and the protection of offshore marine ecosystems (van Rees and Reed 2015, van Rees 2018). These benefits are especially important in urban communities like the population centers around and including Honolulu, which will be increasingly threatened by sea level rise and intensifying storm activity .
What do these results mean for delisting or down-listing possibilities for the species, particularly with respect to the concern about their conservation reliance? Given the small possibility of removing all of the threats to Hawaiian Stilt and Hawaiian Coot persistence on the islands, intensive management involving predator control and habitat management is critical for the foreseeable future. In the absence of active management, it is clear that these species are likely to decline rapidly. Conservation reliance is not explicitly addressed by the U.S. Endangered Species Act, so the implications of conservation reliance are not addressed in recovery planning (Scott et al. 2020). Therefore, allowing de-/ down-listing in the presence of a mechanism for ongoing management likely would require changing the Act. If either species is down-listed, current legal requirements for management would not change. Delisting would be another matter because the mandate for management would be lifted. The proximity of several vital rate values to thresholds below which population persistence becomes much less likely suggests that easing management could result in either species becoming endangered once again. Since current estimates of all vital rates come from managed populations, even this precarious status quo is the current best-case scenario. Data for the Hawaiian Stilt are more abundant, replicated, and informative than those for the coot, so the viability assessment for that species engenders more confidence. Despite the predicted high chance of persistence, down-listing the species would require additional conservation actions to ensure that the minimum population size criterion of >2000 individuals is met. Such actions would need to be designed in ways that ensure that other endangered waterbirds are not adversely affected (e.g., Griffin et al. 1989).
Responses to this article can be read online at: https://www.ace-eco.org/issues/responses.php/2147 Author Contributions: CBvR conducted analyses, planned, and wrote the manuscript. CSE collected data for analyses, conducted preliminary analyses, and helped with manuscript writing. JMR collected data for analyses, planned, and co-wrote the manuscript.

Acknowledgments:
This work was supported in part by a grant from the USFWS. We thank Carrie Harrington (USFWS) for insightful questions and discussions about the natural history and management prognosis of these two taxa. We also thank Martha Kawasaki and Arleone Dibben-Young for helpful discussions on coot and stilt behavior and access to waterbird banding and resighting data for Hawaii. We thank Bob Lacy for his advice and guidance with respect to the use of VORTEX software. We thank the subject editor from Avian Conservation and Ecology, Alex Bond, and two anonymous reviewers for their insightful comments that significantly improved this manuscript.    Figure A4: Plots of proportion males at birth (i.e., skew in sex ratio) and probability of quasiextinction and stochastic growth rate r for Hawaiian coots.  represents population at time t, r is the intrinsic population growth rate (rmax), K is carrying capacity, and theta (θ) and alpha () are shape parameters for the theta-logistic and Gompertz equations, respectively.

Model Name Equation
Exponential log( ) ~ log( −1 ) + + ∈  Table A2: Prior probability distributions, initial values, and justifications for all parameters used in Bayesian nonlinear population models for Hawaiian Coots and Hawaiian Stilts. rmax, σproc, σobs K, theta, and α are maximum annual population growth rate, variation in population growth rate (i.e., process variation), detection error (i.e., variation due to observation error), carrying capacity, theta parameter for theta-logistic model, and alpha parameter for the Gompertz model, respectively. Table A2: Prior probability distributions, initial values, and justifications for all parameters used in Bayesian nonlinear population models for Hawaiian Coots and Hawaiian Stilts. r max , σ proc , σ obs K, θ, and α are maximum annual population growth rate, variation in population growth rate (i.e., process variation), detection error (i.e., variation due to observation error), carrying capacity, theta parameter for theta-logistic model, and alpha parameter for the Gompertz model, respectively.