The neutral form drag coefficient is an important parameter when estimating surface turbulent fluxes over Arctic sea ice. The form drag caused by surface features (𝑪_{𝒅𝒏,𝒇𝒓}) dominates the total drag in the winter, but long-term pan-Arctic records of 𝑪_{𝒅𝒏,𝒇𝒓} are still lacking for Arctic sea ice. In this study, we first developed an improved surface feature detection algorithm and characterized the surface features (including height and spacing) over Arctic sea ice during the late winter of 2009-2019 using the full-scan laser altimeter data obtained in the Operation IceBridge mission. 𝑪_{𝒅𝒏,𝒇𝒓} was then estimated using an existing parameterization scheme. This was followed by applying a satellite-derived backscatter coefficient (𝝈^𝒐_𝒗𝒗 ) to 𝑪_{𝒅𝒏,𝒇𝒓} regression model to extrapolate, for the first time, 𝑪_{𝒅𝒏,𝒇𝒓} to the pan-Arctic scale for the entire winter season over two decades (from 1999 to 2021). We found that the surface features have a larger height and smaller spacing over multi-year ice (1.15 ± 0.21 m and 142 ± 49 m) than over first-year ice (0.90 ± 0.16 m and 241 ± 129 m). The monthly mean 𝑪_{𝒅𝒏,𝒇𝒓} increases through the winter, from 0.2 × 10 −3 in November to 0.4-0.5 × 10 −3 in April. The central Arctic has the largest 𝑪_{𝒅𝒏,𝒇𝒓} (up to 2 × 10 −3), but experienced a drop of ~50% in the period from 2001/2002 to 2008/2009. The interannual fluctuations in 𝑪_{𝒅𝒏,𝒇𝒓} are strongly linked to the variability of sea ice thickness and deformation, and the latter has become increasingly important for 𝑪_{𝒅𝒏,𝒇𝒓} since 2009.