Annual populations and spatial density estimates are needed for the threatened Alaska-breeding population of Steller’s eiders (Polysticta stelleri) to assess recovery status and guide recovery actions but these quantities are poorly estimated in this rare species using traditional methods. Therefore, population size and spatial variation in density were estimated across the Arctic Coastal Plain (ACP) and Utqiagvik Triangle (Triangle) survey areas using spatio-temporal generalized additive models and compared to traditional design-based estimates. Compared to design-based estimates, model-based estimates were more precise and less variable between years. Moreover, data sets can be combined to produce common estimates, detection probability can be incorporated to estimate population size, and spatial density maps can be produced at a small scale. Across both data sets and when combined, Steller’s eider populations fluctuated with an approximate 6.52-year period (95% CI: 6.23, 6.76) with populations cycling from modelled posterior lows of 25 to 500 individuals and highs from 230 to >2000 individuals across the ACP. Over the long term, the posterior 25-year geometric mean growth rate was -0.02 (CI: -0.07, 0.02), and shorter term growth rates were less well estimated and fluctuating between positive and near zero. Mean population size over the last 20 years was 405 (CI: 208, 750) across the whole ACP and 214 (CI: 111, 402) in the Triangle area. Density maps showed a high density concentration in the northern part of the Triangle area and lower densities southward and across the ACP. Models fit to just the Triangle area, indicated that eider density has been shifting northward though time, but this pattern was not supported in the much sparser ACP data set, nor was it well estimated when the ACP and Triangle data sets were combined. Given the strong cycle in this population, population size estimates at similar cycle phases should be used for trend estimates.