Optimising monitoring efforts for secretive snakes: a comparison of occupancy and N-mixture models for assessment of population status

A fifth of reptiles are Data Deficient; many due to unknown population status. Monitoring snake populations can be demanding due to crypsis and low population densities, with insufficient recaptures for abundance estimation via Capture-Mark-Recapture. Alternatively, binomial N-mixture models enable abundance estimation from count data without individual identification, but have rarely been successfully applied to snake populations. We evaluated the suitability of occupancy and N-mixture methods for monitoring an insular population of grass snakes (Natrix helvetica) and considered covariates influencing detection, occupancy and abundance within remaining habitat. Snakes were elusive, with detectability increasing with survey effort (mean: 0.33 ± 0.06 s.e.m.). The probability of a transect being occupied was moderate (mean per kilometre: 0.44 ± 0.19 s.e.m.) and increased with transect length. Abundance estimates indicate a small threatened population associated to our transects (mean: 39, 95% CI: 20–169). Power analysis indicated that the survey effort required to detect occupancy declines would be prohibitive. Occupancy models fitted well, whereas N-mixture models showed poor fit, provided little extra information over occupancy models and were at greater risk of closure violation. Therefore we suggest occupancy models are more appropriate for monitoring snakes and other elusive species, but that population trends may go undetected.