Most natural landscapes are characterized by multiscale (often multifractal) topography with well-known scale-invariance properties. For example, the spectral density of landscape elevation fields is often found to have a power-law scaling behaviour (with a -2 slope on a log-log scale) over a wide span of spatial scales, typically ranging from tens of kilometres down to a few metres. Even though the effect of topography on the atmospheric boundary layer (ABL) has been the subject of numerous studies, few have focussed on multiscale topography. In this study, large-eddy simulation (LES) is used to investigate boundary-layer flow over multiscale topography, and guide the development of parametrizations needed to represent the effects of subgrid-scale (SGS) topography in numerical models of ABL flow. Particular emphasis is placed on the formulation of an effective roughness used to account for the increased aerodynamic roughness associated with SGS topography. The LES code uses the scale-dependent Lagrangian dynamic SGS model for the turbulent stresses and a terrain-following coordinate transformation to explicitly resolve the effects of the topography at scales larger than the LES resolution. The terrain used in the simulations is generated using a restricted solid-on-solid landscape evolution model, and it is characterized by a -2 slope of the elevation power spectrum. Results from simulations performed using elevation fields band-pass filtered at different spatial resolutions indicate a clear linear relation between the square of the effective roughness and the variance of elevation.