Large-eddy simulations were employed over eight different satellite-sensed sea ice patterns to analyze how the overlying atmospheric boundary layer (ABL) dynamics and thermodynamics differ based on the heterogeneity of the surface. In Suite 1, four simulations were conducted for three different sea ice patterns, resulting in a total of 12 simulations. In Suite 2, two simulations were conducted for eight different sea ice patterns, resulting in a total of 16 simulations. Mean time-averaged variables presented include the geostrophic wind directions u, v, and w, potential temperature theta, and heat and momentum fluxes.