Influence of microhabitat on Honduran Emerald ( Amazilia luciae ) abundance in tropical dry forest remnants

Understanding the ecology of at-risk species is the foundation for developing strategies to counteract continued population declines. The Honduran Emerald (Amazilia luciae) is an endemic hummingbird that inhabits tropical dry forest in Honduras. Remaining populations of this endangered species are restricted to habitat fragments located in landscapes dominated by agricultural activities. The conservation of this species is dependent on efforts to protect existing dry forest remnants and to restore additional areas. However, limited information exists regarding Honduran Emerald habitat use on which to base management decisions. We conducted a study to identify point-level habitat features that are important to Honduran Emerald abundance. In 2014–2015, we conducted avian and vegetation surveys at 174 points located within dry forest remnants. We constructed hierarchical multinomial mixture models using time-removal based point count data collected in three sampling periods. Honduran Emerald detection was influenced negatively by Julian date. Local abundance was positively correlated with shrub-sapling density and cacti structural diversity. Mean local abundance per point was (0.39 ± 0.2 SD) and the number of individuals estimated to have used the combined area of the survey points (48.7 ha) was 68 (CI 95% 45-110). The shrub-sapling density and cacti structural diversity levels identified from our abundance models are characteristic to dry forest communities in this region. As such, the need for efforts to protect remaining dry forest remnants cannot be overstated. Future research should evaluate how the microhabitat features we found to be important to Honduran Emerald abundance influence fitness parameters such as survival and reproduction. Studies that examine how reciprocal effects of pollinator species declines and plant species diversity influences the long-term stability of tropical dry forest ecosystems are also warranted. Influence du micro-habitat sur l'abondance de l'Ariane de Lucy (Amazilia luciae) dans les vestiges de forêts tropicales sèches RÉSUMÉ. Comprendre l'écologie des espèces en péril est la base de l'élaboration de stratégies visant à contrer le déclin continu des populations. L'Ariane de Lucy (Amazilia luciae) est un colibri endémique qui vit dans la forêt tropicale sèche du Honduras. Les populations restantes de cette espèce en danger sont limitées à des fragments d'habitat situés dans des paysages dominés par des activités agricoles. La conservation de cette espèce dépend des efforts déployés pour protéger les vestiges de forêt sèche existants et pour restaurer d'autres zones. Cependant, peu d'informations sont disponibles quant à l'utilisation de l'habitat de l'Ariane de Lucy sur lesquelles fonder les décisions de gestion. Nous avons mené une étude pour identifier les caractéristiques de l'habitat au niveau ponctuel qui sont importantes pour l'abondance de l'Ariane de Lucy. En 2014-2015, nous avons effectué des relevés d'oiseaux et de végétation sur 174 points situés dans des vestiges de forêts sèches. Nous avons construit des modèles de mélanges multinomiaux hiérarchiques à l'aide des données des dénombrements ponctuels, basées sur l'élimination du temps, collectées au cours de trois périodes d'échantillonnage. La détection de l'Ariane de Lucy était négativement corrélée au jour Julien. L'abondance locale était positivement corrélée à la densité d'arbustes et à la diversité structurelle des cactus. L'abondance locale moyenne par point était de 0,39 ± 0,2 DS et le nombre d'individus ayant utilisé la superficie totale des points d'étude (48,7 ha) était de 68 (IC 95% 45-110). La densité d'arbustes et les niveaux de diversité structurelle des cactus identifiés à partir de nos modèles d'abondance sont caractéristiques des communautés de forêts sèches de cette région. En conséquent, on ne saurait trop insister sur la nécessité de déployer des efforts pour protéger les vestiges restants des forêts sèches. Des recherches futures devront évaluer dans quelle mesure les caractéristiques du micro-habitat que nous avons trouvées importantes pour l'abondance de l'Ariane de Lucy ont une incidence sur les paramètres de condition physique tels que la survie et la reproduction. Des études qui examineraient comment les effets réciproques du déclin des espèces pollinisatrices et de la diversité des espèces végétales influencent la stabilité à long terme des écosystèmes des forêts tropicales sèches sont également justifiée


INTRODUCTION
Effective conservation of biodiversity in landscapes impacted by human activities is aided by the degree to which biodiversity is monitored and studied in the face of habitat degradation and loss (Collen et al. 2013).Ecologists have often considered species abundance as one of the fundamental population metrics to evaluate anthropogenic influences on wildlife (Apps et al. 2004, Heikkinen et al. 2004, Marsh and Trenham 2008).For example, abundance has been used in population-specific studies to inform management and conservation issues of rare or endangered species (Nichols andMackenzie 2004, Kéry andRoyle 2016).When species abundance can be evaluated in the context of habitat structure and composition, studies can offer insight to those habitat conditions that are important for population recovery.Moreover, results from such studies have the potential to inform habitat management and other conservation efforts (Iwata et al. 2003, Heikkinen et al. 2004, Segura and Arturi 2012).
Tropical dry forests are considered a globally threatened ecosystem, with an estimated 60% of their original extent lost and the rest highly threatened by fragmentation (Sánchez-Azofeifa et al. 2013, Blackie et al. 2014).These forests may not match tropical rainforests in their high level of plant diversity, but they hold a significant amount of plant endemism and are characterized by their high floristic turnover across their distribution (Gentry 1995, Banda et al. 2016).Areas that support tropical dry forests have been historically selected by humans for settlement (Gómez-Pompa et al. 1987, Bush et al. 1992), and today, many landscapes that support tropical dry forests remain important for human communities (Blackie et al. 2014).The area of tropical dry forest in the Agalta Valley, Honduras, has been drastically reduced over the past several decades as the valley continues to be transformed into a working landscape driven by dairy production (Banco Interamericano de Desarrollo 2009).Continued encroachment of human land uses on dry forest remnants may jeopardize the persistence of many species that specialize on this plant community.
Among the many species that have been affected by tropical dry forest loss is the Honduran Emerald, Amazilia luciae (Fig. 1), an endangered hummingbird that is endemic to tropical dry forests of Honduras (Anderson et al. 2013, BirdLife International 2016).The species is now restricted to only five locations, i.e., departments, throughout Honduras that support varying amounts of dry forest habitat (Fig. 2A; eBird 2017).The status of the Honduran Emerald in the Agalta Valley is especially tenuous because the species is predicted to become locally extinct within the next decade unless considerable conservation measures are employed (Anderson et al. 2013).
To effectively address conservation needs of the Honduran Emerald it is important to quantify population parameters, i.e., population size (Anderson et al. 2010), and to link such parameters with empirically accessed habitat conditions.Despite first being described by Lawrence (1867), little is known about Honduran Emerald ecology (Anderson et al. 2010).Most of what is known about the species is the result of anecdotal observations, descriptive accounts, contemporary expeditions, and environmental impact assessment reports (Howell and Webb 1989, Thorn et al. 2000, Anderson and Hyman 2007, Anderson et al. 2010, Espinal and Mora 2012, INGTELSIG 2013).To date, no studies have quantified the influence of point-level vegetation on Honduran Emerald abundance.Several studies in avian ecology have previously demonstrated the utility of habitat-species abundance relationships to highlight the role of vegetation structure and composition for species conservation (Pardini et al. 2005, Gómez-Montes andBayly 2010) and to assess the population status of species in a given locality (Sillett et al. 2012, Rivera-Milán et al. 2015).Therefore, blueprints exist to study how Honduran Emerald local abundance is driven by vegetation characteristics in tropical dry forest remnants.Furthermore, floristic communities that shape tropical dry forests across the Americas are considered irreplaceable with areas in need of conservation throughout the continent (Banda et al. 2016).As such, understanding the interrelationships of the species within this ecosystem will increase our knowledge of Honduran Emerald ecology and of the impact that anthropogenic threats will have on these species (Sánchez-Azofeifa et al. 2013).Our aim was to evaluate the relationship between Honduran Emerald abundance and microhabitat features of tropical dry forest remnants in the Agalta Valley.Our use of the term "microhabitat" is similar in context to Cahill and Matthysen (2007) or Renfrew and Ribic (2008) and refers to vegetation structure and species composition features that are quantified at the survey point level.Study objectives were to (i) determine mean detection probability and local abundance of Honduran Emerald, (ii) identify how microhabitat features influence local abundance of Honduran Emerald, and (iii) summarize the microhabitat features of dry forest remnants in which we detected Honduran Emeralds.

Study area
During August 2014 to July 2015, we conducted surveys in 35 dry forest remnants (0.5-380 ha) within the Agalta Valley (Fig. 2B).The geographic extent of the Agalta Valley is approximately

Avian surveys
We used QGIS to randomly select point count survey locations within each of the 35 tropical dry forest remnants (Quantum GIS Development Team, https://qgis.org/en/site/).This resulted in the selection of 179 points that were spaced at least 300 m apart (Fig. 2B).Of these we only analyzed data from 174 points because of loss of habitat or access to the points within the study period.We conducted point count surveys between 0600 and 1000 hrs on days with appropriate weather conditions, e.g., no rain (Lynch 1995).We surveyed each point during three periods throughout the year-long study (30 September-9 December; 4 March-28 April; 21 June-27 July).At the beginning of each point count survey, we allowed a 1-minute "settle period" to elapse before starting the survey (Bibby et al. 1992).Each point count consisted of a 10-minute passive listening phase in which all Honduran Emeralds seen or heard were recorded.Listening phases were organized into five 2-minute time-removal intervals (Farnsworth et al. 2002).
Surveys throughout the study period were conducted by two observers (FR and DE).These observers trained together during two days in August 2014, in one of the largest dry forest fragments and conducted trial point counts to learn the point count protocols.Subsequently, the observers conducted the first period of sampling (September-December) together to standardize their ability to record data: detections of Honduran Emeralds and recording of Honduran Emeralds in the adequate intervals.In the subsequent sampling periods (March-April and June-July), each observer worked separately to maximize the points visited per survey.For this reason in these two sampling periods, we considered observer as a potential variable in our models to account for possible differences between them in our ability to detect Honduran Emeralds.

Vegetation surveys
We conducted vegetation surveys at each point where we had conducted avian point counts.In each plot we collected vegetation structure and plant species composition data within nested plots based on Nudds (1977), Bibby et al. (1992), Martin et al. (1997), Nassar et al. (2008), and García-Villacorta (2009).Within a 5-m radius plot we recorded (1) average shrub-sapling height (m) estimated by observers, and (2) shrub-sapling stem count including woody species in the shrub strata < 2 m.Within a 12.6 m radius plot we recorded the following: (1) diameter at breast height (DBH) for woody plants ≥ 3 cm DBH; (2) tree height (m) for each woody plant ≥ 3 cm DBH estimated by observers; (3) tree species; (4) count of trees or cacti with presence of epiphytic bromeliads; (5) number of cacti; (6) cacti height (m) for each individual estimated by observer; (7) count of woody stems < 3 cm DBH; and (8) average horizontal vegetation density measured using a 2-m tall density board consisting of a cloth scroll comprising twenty 20X20 cm grid cells.To quantify this measurement an observer stood at plot center while another observer held the density board 10 m away in each of the four cardinal directions.The observer at plot center recorded the number of grid cells covered ≥ 50% by vegetation.To minimize errors in estimation, all vegetation measurements were collected by the same observers during all surveys.Once vegetation measurements were completed, we collected representative plant samples throughout each plot, preferably from individuals with fruit or flowers.Plants were pressed, transported to the Universidad Nacional Autónoma de Honduras (UNAH), to be identified to species and deposited at the Herbarium Cyril Hardy Nelson Sutherland (TEFH) in UNAH.We deposited a duplicate of the samples at the A. G. Shields Herbarium in Indiana University of Pennsylvania.

Statistical analysis
First, we summarized the vegetation metrics to generate 15 variables (Table 1).These variables reflect the structure and taxonomic composition of microhabitat features of each survey location (Table 1).Prior to analysis, we standardized all microhabitat features to have a mean zero and standard deviation of one.Additionally, we evaluated for correlations using the Pearson coefficient (r ≥ 0.75 or ≤ -0.75).Only two covariates, woody species richness and Shannon Diversity, were highly correlated (Pearson coefficient = 0.85), and we chose to exclude the latter from our analyses.Ultimately, we included 14 variables (Table 1) in the hierarchical models to relate Honduran Emerald abundance to microhabitat features.
We analyzed data from our three Honduran Emerald survey periods separately to limit our inferences to time periods when the population was sampled consistently.We assume that within each survey period, populations are likely to be closed.However, because of the mobility of our study species we interpreted our abundance values as the number of individuals that potentially use the sampling points at a given time, as opposed to absolute number of individuals (Chandler et al. 2011).We screened the raw counts of Honduran Emeralds and included detections that were ≤ 30 meters from the observer.Observations beyond this distance were sparse and may not have captured the microhabitat features that we measured.Honduran Emerald counts within each sampling period were tested for spatial autocorrelation using Moran's I test within the "ape" package in R (Paradis and Schliep 2018).
We analyzed the local abundance patterns with the "unmarked" package in program R for model building (Fiske and Chandler 2011, R Core Team 2017).We followed the multinomial Nmixture modeling framework using the Poisson mixture, incorporating a time-removal specification with the information from the count time-intervals (Farnsworth et al. 2002, Royle 2004a).We fit models in a two phase process.First, we fitted models that only considered the detection probability covariates of Julian date, observer, and vegetation density and kept abundance constant.We used the highest-ranking model from this phase for posterior modeling of abundance.Then, we fit models with the top detection covariate of the previous phase and the vegetation variables that may influence Honduran Emerald abundance (Table 1).To avoid overparameterization in the models, because of low number of detections, we fit single variable models and single variable with quadratic terms for each of the 14 variables.We ranked the models based on the Akaike Information Criteria corrected for small sample sizes (AICc) using the package "AICc modavg" (Hurvich andTsai 1989, Mazerolle 2017).The top model was selected for the estimation of mean abundance per point and local abundance.We evaluated the confidence intervals at the 85% level because we consider them appropriate for exploring the informative value of our proposed variables (Arnold 2010).Additionally, we conducted a goodness of fit test with three fit statistics (Chi-square, FreeTukey, and sums of squares) and 10,000 repetitions following Kéry and Royle (2016).Models up to ΔAICc = 7 were evaluated to identify variables that could be considered in future studies, albeit these are less likely to influence abundance based on our data (Burnham and Anderson 2002).We present the results of the models fit in phase two and the null model in a table.
Finally, we summarize the vegetation variables of used and unused point count locations.For all vegetation variables, we conducted a two-sample t-test and applied a Holm's correction factor to account for multiple testing (Holm 1979).We present test results for both corrected and uncorrected comparisons.The corrected comparisons were used to consider the effect of the differences between used and unused points, whereas the uncorrected comparisons were used to identify microhabitat features that warrant inclusion in future Honduran Emerald studies.

RESULTS
We detected Honduran Emeralds at 86 of 174 point count locations in tropical dry forest remnants during all three survey periods combined.In period one (September-December), we recorded 46 Honduran Emeralds at 41 of 174 (24%) survey points.During period two (March-April), we recorded 67 Honduran Emeralds at 57 of 174 (33%) survey points, and during period three (June-July), 30 Honduran Emeralds at 28 of 174 (16%) survey points.A total of five detections > 30 m were not included in the analyses (four and one for the first and second sampling periods, respectively).We did not find evidence for spatial autocorrelation for Honduran Emerald counts for all sampling periods (Appendix 1).This index combines height and presence of species to reflect the contribution of cacti species richness in each of three levels of a vertical distribution (Pretzsch 2009).Three levels were determined based on the raw individual height distribution for cacti: 0-7.5 m (level III or lower stratum), 7.5-12 m (level II or medium stratum), and 12-15 m (level III or high stratum).This standardized species profile index (Arel) was adapted to cacti.Three levels were selected based on applications in other vegetation ecology studies (Jiménez et al. 2001, Mora-Donjuán et al. 2014).Cacti are one of the characteristic families of tropical dry forest floristic composition (Gentry 1995).

Proportion of Opuntia hondurensis cacti
Proportion of individuals of Opuntia hondurensis relative to all cacti individuals present in the plot.
Indicates the availability of one of the Honduran Emerald common food sources in the height stratum where most individuals of this cactus species were concentrated (0-7.5 m).

Proportion of Pilosocereus leucocephalus cacti
Proportion of individuals of Pilosocereus leucocephalus relative to all cacti individuals present in the plot.
Indicates the availability of one of the Honduran Emerald common food sources in the height stratum where most individuals of this cactus species were concentrated (0-7.5 m).Shrub-sapling density (no.stems/ha) Shrub density is represented by all woody plants under two meters that have solely a shrub growth habit or may have both shrub and tree growth habit.
Shrub strata encompasses plant species that are Honduran Emerald feeding source or nesting substrate.
Vegetation density (%) Percentage of horizontal cover present within two m of height.
A complementary value to complexity of a stand that represents the cover wildlife uses (Nudds 1977).
It may indicate the potential cover that may serve the Honduran Emerald for several purposes such as predation protection or nesting attributes.

Importance Value (IVI) Acacia picachensis
Quantified at the plot level with the two available parameters of height and basal area (Mostacedo and Fredericksen 2000) for Acacia picachensis one of the highest IVI values overall for Agalta as composition microhabitat features.
To incorporate the composition of Agalta's plant community.

Importance Value (IVI) Eugenia hondurensis
Quantified at the plot level with the two available parameters of height and basal area (Mostacedo and Fredericksen 2000) for Eugenia hondurensis one of the highest IVI values overall for Agalta as composition microhabitat features.
To incorporate the composition of Agalta's plant community.

Importance Value (IVI) Erythroxylum areolatum
Quantified at the plot level with the two available parameters of height and basal area (Mostacedo and Fredericksen 2000) for Erythroxylum areolatum one of the highest IVI values overall for Agalta as composition microhabitat features.
To incorporate the composition of Agalta's plant community.

Importance Value (IVI) Casearia nitida
Quantified at the plot level with the two available parameters of height and basal area (Mostacedo and Fredericksen 2000) for Casearia nitida one of the highest IVI values overall for Agalta as composition microhabitat features.
To incorporate the composition of Agalta's plant community.

Epiphyte tree count
Count per plot of the individual woody plants greater than three cm DBH with presence of at least one epiphytic bromeliad.
Indicates the availability of plant species that may be used by Honduran Emeralds as a food resource.

Patch size
Area in ha of the polygon representing the study site that contains the point count location.
Indicates the extent of the property that holds the mosaic of tropical dry forest habitat and cattle farm use areas.Shannon Diversity Index Shannon Diversity Index presented in effective species units.Effective species units (D) is calculated by taking the base of the natural logarithm of the Shannon Diversity Index.D = eH'.
Indicates plant species diversity.
Woody plant species richness Species richness of all individuals ≥ three cm DBH.Indicates plant species richness.The probability of detecting a Honduran Emerald declined over time (Julian date) during the first two sampling periods (Fig. 3).During the second sampling period this decline was best characterized by a polynomial model with a quadratic term.Mean detection probability was the highest in the first sampling period (Table 2).Models constructed using data collected during our third sampling period were unable to estimate detection probability and abundance.We expect this is due to the lower number of detections.As such, we did not attempt further modeling using data collected during the third survey period.During our first sampling period, Honduran Emeralds were most abundant at locations with greater shrub-sapling density.This relationship was described best with a polynomial model including a quadratic term.This positive relationship (Table 2; Fig. 4) appeared to peak around 200 stems/ha, beyond which our model was not able to characterize a relationship.During our second sampling period, Honduran Emeralds were most abundant at locations with greater cacti structural diversity (Table 2; Fig. 4).Cacti structural diversity suggested that cacti at the lower strata between zero and seven meters were more prevalent with mean 22.62 ± 15.17% SD for all points used during sampling period two (N = 57 points).These values of cacti structural diversity represent the used points with mean cacti height of 2.47 ± 1.66 m (range = 1-8.2m), total cacti density of all species of 235.44 ± 89.84 individuals/ha (range = 40-400), and cacti species richness of 2.09 ± 1.1 (range = 0-4 species) where Opuntia hondurensis and Pilosocereus leucocephalus were the most frequent cacti species in the points.The highest local abundance was estimated during the second sampling period, albeit this period had the lowest detection probability (Table 2).Models from both sampling periods were not overdispersed (c-hat = 1.01 sampling period one, and 0.99 for sampling period two) and had good fit according to the sums of squares, Chi-square, and Freeman Tukey statistics (p > 0.05).Our t-test results indicated that the proportion of the cactus Pilosocereus leucocephalus, relative to other cacti at a given plot, was higher at used points (Mean = 0.2; SD = 0.27; Table 3) compared with unused points (Mean = 0.08; SD = 0.18; Table 3).We propose the following variables to be included in future studies: cacti structural diversity, dominance (IVI) of Casearia nitida, Eugenia hondurensis, species richness of woody plants, epiphyte tree count, Holdridge Complexity Index, stems under 3 cm DBH, and patch size.These variables approached the threshold of significance in our t-tests that evaluated sites with and without presence of Honduran Emeralds (Table 3) and were within ∆ AICc=4-7 (Table 4).

DISCUSSION
Our study revealed that shrub-sapling density and cacti structural diversity explained Honduran Emerald local abundance in tropical dry forest remnants of the Agalta Valley.In other parts of the species range, shrubs and cacti species have been highlighted as components of the tropical dry forest used by this species (Thorn et al. 2000, House 2004).In fact, it is well documented that shrub and cacti are ecologically important to many species of tropical hummingbirds for feeding, or perching close to feeding areas (Skutch 1958, Wolf 1964, Wolf and Stiles 1970, Feinsinger 1976, Hainsworth 1977, Snow and Snow 1986, Fraga 1989).In addition, small trees or saplings are considered to have suitable floral resources that are suitable for hummingbirds (Snow and Snow 1986).Honduran Emeralds' use of shrubs and cacti as a floral resource has been reported elsewhere (Anderson et al. 2010).The use of the shrub-sapling strata by Honduran Emeralds as a source of floral resources was observed previously in the Agalta Valley (Mora et al. 2016).Moreover, field observations during our study revealed that the species used floral resources of four shrub and two cacti species: Aphelandra scabra (Acanthaceae), Combretum fruticosum (Combretaceae), Pedilanthus tithymaloides (Euphorbiaceae), Cnidoscolus aconitifolius (Euphorbiaceae), and Opuntia hondurensis (Cactaceae) and Pilosocereus leucocephalus (Cactaceae), respectively.Honduran Emeralds also used shrubsapling and cacti strata as nesting substrate (Rodríguez et al. 2016).
In addition to shrub-sapling density, cacti structural diversity most explained variation in Honduran Emerald abundance across our survey locations.This finding reinforces the previous observations and highlights the importance of cacti species as nesting and feeding sites (Anderson et al. 2010, Rodríguez et al. 2016).Furthermore, when all sampling periods were pooled, the percent cover of the cactus, Pilosocereus leucocephalus, was the only vegetation feature we measured that differed between used and unused points whereby used points had higher values.Our results regarding the influence of shrub-sapling density and cacti structural diversity on Honduran Emerald abundance should be considered in future studies that attempt to evaluate the influence of habitat conditions on the species fitness parameters (Johnson 2007).For example, most Honduran Emerald nests found in dry forest remnants in the Agalta Valley were placed in the shrubsapling strata (Rodríguez et al. 2016).But it remains unknown how shrub-sapling density or vegetation composition affects nest success.Our results in combination with those from studies that evaluate reproductive success and other fitness parameters can guide habitat management recommendations for the species.
We also found that the probability of detecting Honduran Emerald varied within sampling periods and that it was best explained by Julian date.Although we are uncertain why Julian date had a negative influence on Honduran Emerald detection, we consider that detectability was influenced negatively by variability in the detection history for any survey location, e.g., one individual detected in one point versus three in a different point (Royle 2004b, Dénes et al. 2015).In turn, we propose that the variability in detection history within each sampling period may reflect the patterns of hummingbird movement within the heavily fragmented study area that are driven by seasonal variability and availability of resources.That is, the need for individuals to periodically shift activity centers throughout the year in response to shifting floral resources, may translate to variability in detection history for any given survey location.Past studies have indicated how hummingbird detections vary per month, most likely because of flower availability driven by vegetation phenology in these environments (Arizmendi and Ornelas 1990, Murphy and Lugo 1995, Rodríguez-Flores and Wethington 2017).Furthermore, if floral resources are patchily distributed, hummingbird movements could be influenced by the landscape configuration making their availability and detection more variable (Hadley andBetts 2009, Volpe et al. 2014).Future studies that examine Honduran Emerald ecology at the territory or home range scales throughout the entire annual cycle will likely provide insight needed to understand mechanisms that explain the influence of Julian date on detection probability.
Our local abundance estimates extrapolated to the extent of dry forest remnants (2840 Honduran Emerald individuals, CI 95% 1880-4720) are comparable to population estimates for the Agalta Valley presented in the Endangered Species Act ruling for the species (United States Fish and Wildlife Service 2015).However, we suggest a caveat to this observation, which is the observed variability of abundance across our survey points that was influenced by microhabitat features described above.This influence reflects the variation of habitat conditions throughout the study area, thus making the assumption of a uniform distribution of Honduran Emeralds in the dry forest remnants difficult to meet (Kéry and Schmidt 2008).Nevertheless, our population estimate is valuable for understanding the current conservation status of this species in the Agalta Valley.The rate at which dry forest is being converted to cattle farms in the Agalta Valley will result in continued habitat loss that will negatively influence Honduran Emerald local abundance and ultimately its population dynamics.Between 2016 and 2017, two years after our study was conducted, at least 30 ha of dry forest were lost within our study sites based on a global forest cover change database (Hansen et al. 2013).As such, our abundance estimate supports the assertion by others that conservation of existing habitat is needed to prevent further population decline (Anderson et al. 2010, Mora et al. 2016).Such conservation actions should begin immediately and be implemented at biologically meaningful scales.
Our study is the first to account for imperfect detection in abundance estimates for the endangered Honduran Emerald.We provided additional insight on the species abundance patterns by estimating detection probability with a hierarchical modeling approach.In doing so, we avoided the assumption that detection probability is always perfect, and we correct by estimating the variability in detection to obtain an abundance estimate (Kéry and Schmidt 2008).We stress that future studies or monitoring aimed at estimating Honduran Emerald local abundance or population size use a framework that accounts for imperfect detection (Kéry and Schmidt 2008).Indeed, future Honduran Emerald surveys could be conducted when detection probability is expected to be high, i.e., September-December p = 0.23 ± 0.12).However, to relate abundance with reproductive success and assess habitat quality, surveys should be conducted within the breeding season (January-April; Rodríguez et al. 2016) and should include quantification of cacti structural diversity and shrub-sapling strata density.Equally important to future monitoring will be to conduct multiple visits to each survey point within these months, e.g., two-three visits, to counter the factors that contribute to imperfect detection (Royle 2004b).Parallel to the considerations in survey design to model Honduran Emerald abundance, our characterization of the vegetation at points used by Honduran Emeralds can be informative for future studies and conservation actions (Table 4).The microhabitat values we quantified may serve as a guideline for habitat restoration activities (Janzen 1988).Subsequently, different habitat managed treatments could be evaluated to assess their effectiveness for restoring and maintaining Honduran Emerald habitat.
Our research revealed how microhabitat features influence the local abundance of the endangered Honduran Emerald and how the use of modeling techniques for avian counts can be fine-tuned to inform the species status in the Agalta Valley.Similar work should be conducted in the other valleys where the Honduran Emerald is known to occur to obtain insight to the species status across its entire distributional range.Our study framework could be useful in other dry forests in the valleys of Aguán and Telica, which share similar composition and structure to the Agalta Valley (Anderson et al. 2010).Furthermore, because of past detections of the species in ecotones between dry forest and pineoak in the Agalta Valley, a study that compares Honduran Emerald abundance between these two cover types should be conducted (Mora et al. 2016).Until more insight on Honduran Emerald ecology is gained from additional research and monitoring, conservation activities for this species should target dry forest protection and restoration, as well as promote dense shrub-sapling strata and structurally diverse cacti components.Periodic monitoring of Honduran Emerald response to conservation activities is warranted, but in a geography with limited conservation resources we urge that the immediate priority should be to focus available funding on working with landowners to protect existing dry forest and to restore additional areas.
Although our study focused on the relationship between microhabitat features and the abundance of a single species, our findings provide a broad reminder of the implications species declines can potentially have on the intricate codependencies that exist between plant communities and their associated pollinators.For instance, Honduran Emerald local abundance was positively influenced by cacti structural diversity that included Opuntia hondurensis, suggesting that this imperiled hummingbird is likely a key pollinator of this dry forest endemic plant, as seen in other Opuntia sp.pollinated by hummingbirds (Díaz and Cocucci 2003).Researchers have previously highlighted the importance of hummingbirds as "generators of biodiversity" due to their pollination services (Arizmendi 2017).We are unaware of the magnitude of the role the Honduran Emerald plays as a pollinator of Opuntia hondurensis and other dry forest species, but it is well known that reductions in pollinator species diversity can have negative impacts on plant communities (Potts et al. 2010).For example, increasing the functional diversity of both pollinators and plants led to the persistence of more diverse plant communities (Fontaine et al. 2006).Moreover, reductions in pollinator functional diversity has the potential to initiate plant population declines or extinctions (Kearns et al. 1998).With tropical dry forests being floristically distinct across the Americas (Banda et al. 2016), understanding how hummingbirds and dry forest plant communities interact will need to go hand-in-hand to guarantee effective conservation of these tropical dry forests and the species they support.Indeed, studies that examine how reciprocal effects of pollinator species declines and plant species diversity influences the long-term stability of tropical dry forest ecosystems are warranted.Conservation efforts such as the one ongoing in the Agalta Valley, that acknowledge the needs of both human and natural communities, should increase awareness for the importance of plant-animal interactions and the need to work at both microhabitat and landscape scales to most effectively conserve those interspecific relationships.

Fig. 1 .
Fig. 1.Honduran Emerald hummingbird (Amazilia luciae) is a species endemic to Honduras.(A) An area of dry forest dominated by shrub.(B) Honduran Emerald feeding from Aechmea bracteata in Santa Bárbara department, the western range of the species (Photo by: Mayron Mejía, 2018).(C) A Honduran Emerald hovering in its habitat in the department of Cortés (Photo by: Alex Martínez Matute, 2018).

Fig. 3 .
Fig. 3. Honduran Emerald (Amazilia luciae) detection probability varied with Julian date, with a higher mean detection probability during September-December 2014 (left) than March-April 2015 (right).Graphs represent the topranked models for each survey period, and gray lines represent 95% confidence intervals.

Fig. 4 .
Fig. 4. Expected abundance of Honduran Emerald (Amazilia luciae) in the Agalta Valley was most influenced by shrubsapling density during September-December 2014 (left) and most influenced by cacti structural diversity during March-April 2015 (right).Gray lines indicate 95% confidence intervals.

Table 1 .
Structure and composition microhabitat features of the Agalta Valley's tropical dry forest remnants used in the analysis of Honduran Emerald (Amazilia luciae) local abundance.

Table 2 .
Honduran Emerald (Amazilia luciae) top detection probability and abundance models, and abundance estimates for two sampling periods in the Agalta Valley, during 2014-2015 at 174 sampled point locations (48.7 ha).Mean values with SD for detection and abundance per point also shown.A complete AICc table of all fitted models of detection and abundance can be found in Appendices 2 and 3.

Table 3 .
Comparison of vegetation variables between point count locations that were used and unused by Honduran Emeralds (Amazilia luciae) in the Agalta Valley during 2014-2015.

Table 4 .
Honduran Emerald (Amazilia luciae) models of abundance in the Agalta Valley during 2014-2015 for the first (September-December) and second (March-April) sampling periods.We present models up to Δ AICc = 7, which illustrate microhabitat features that may influence Honduran Emerald abundance.Coefficients (β) with SE in parenthesis and 85% confidence intervals (CI) shown.
|Variable with quadratic term