Calculation of the sensible and latent heat (turbulent) fluxes is required in order to close the surface energy budget of glaciers and model glacial melt. The aerodynamic roughness length, z0, is a key parameter in the bulk approach to calculating sensible heat flux; yet, z0 is commonly considered simply as a tuning parameter or generalized between surfaces and over time. Spatially and temporally distributed observations of z0 over ice are rare. Both direct (from wind towers and sonic anemometers) and indirect (from microtopographic surveys) measurements of z0 are subject to sensitivities and uncertainties that are often unstated or overlooked. In this study, we present a quantitative evaluation of aerodynamic profile-based and microtopographic methods and their effect on z0 using data collected from Storglaciären and Sydöstra Kaskasatjäkkaglaciären, Tarfala Valley, Arctic Sweden. Aggressive data filters discard most of the wind tower data but still produce realistic z0 values of 1.9 mm and 2 mm. Despite uncertainty introduced by scale and resolution dependence, microtopographic methods produced estimates of z0 comparable to wind tower values and those found on similar surfaces. We conclude that (1) in the absence of direct turbulent flux measurements from sonic anemometers, the profile and microtopographic methods provide realistic z0 values, (2) both 2D and 3D microtopographic methods are dependent on scale, resolution, and the chosen detrending method, and (3) careful calibration of these parameters could enable glacier-wide investigations of z0 from remotely sensed data, including those increasingly available from satellite platforms.