Effective species management and conservation benefit from knowledge of species distribution and status. Surveys to obtain that information often involve replicate sampling, which increases survey effort and costs. We simultaneously modeled species distribution, abundance and spatial correlation, and compared the uncertainty in replicate abundance estimates of the endangered palila (Loxioides bailleui) using hierarchical generalized additive models with a soap film smoother that incorporated random effects for visit. Based on survey coverage and detections, we selected the 2017 point-transect distance sampling survey on Mauna Kea, Hawai'i Island, for our modeling. Our modeling approach allowed us to account for imperfect detections, control the effects of boundary features, and generate visit-specific density surface maps. We found that visit-specific smooths were nearly identical, indicating that little information was gained from a subsequent visit, and that most of the estimator uncertainty was derived from within-visit variability. Scaling back the palila survey to a single visit would halve the survey effort and logistical costs and increase efficiencies in data management and processing. Changing the sampling protocol warrants careful consideration and our findings may help management and regulatory agencies by maximizing efficiency and minimizing costs of surveying protocols, while providing guidelines on how to best collect information critical to species' conservation.

Effective species management and conservation require knowledge of species distribution and status. We used point-transect distance sampling surveys of the endangered palila (Loxioides bailleui), a honeycreeper currently found only on the Island of Hawai'i, USA, to generate robust estimates of total abundance and simultaneously model the distribution, abundance, and spatial correlation of the species as a density surface model (DSM). Point-transect distance sampling is a widely applied method to estimate bird densities accounting for imperfect detection probability. For the DSM we used a generalized additive model framework and soap film smoothers to control the effects of boundary features. This modeling approach allowed us to account for imperfect detection and propagate detection probability uncertainty. We compared the uncertainty in palila abundance estimates using standard point-transect distance sampling to estimates from the DSM. The DSM, accounting for both distance-sampling-derived detection probability variance and the generalized additive model density estimate variance, did not improve population estimator precision; however, it provided insight into the species' distribution, density, and uncertainty. We also applied excursion sets analysis to objectively identify areas where the species occurs in high densities. The 2017 global population of <2,000 individuals was limited to an excursion area of 1,500 ha. Our findings can help management and regulatory agencies by simultaneously mapping a species' distribution and density, improving survey protocols, and providing information important to species conservation.