In the polar regions, the atmospheric boundary layer (ABL) characteristics are strongly influenced by convection over leads, which are elongated channels in the sea ice covered ocean. The effects on the ABL depend on meteorological forcing and lead geometry. In non-convection-resolving models, in which several leads of potentially different characteristics might be present in a single grid cell, such surface characteristics and the corresponding ABL patterns are not resolved. Our main goal is to investigate potential implications for such models when these subgrid-scale patterns are not considered appropriately. We performed non-eddy-resolving microscale simulations over five different domains with leads of different widths separated by 100% sea ice. We also performed coarser-resolved simulations over a domain representing a few grid cells of a regional climate model, wherein leads were not resolved but accounted for via a fractional sea ice cover of 91% in each cell. Domain size and mean sea ice concentration were the same in all simulations. Differences in the domain-averaged ABL profiles and patterns of wind, temperature, and turbulent fluxes indicate a strong impact of both the leads and their geometry. Additional evaluations of different turbulence parameterizations show large effects by both gradient-independent heat transport and vertical entrainment.