AbstractObservations in polar regions show that sea ice deformations are often narrow linear features. These long bands of deformations are referred to as Linear Kinematic Features (LKFs). Viscous‐plastic sea ice models have the capability to simulate LKFs and more generally sea ice deformations. Moreover, viscous‐plastic models simulate a larger number and more refined LKFs as the spatial resolution is increased. Besides grid spacing, other aspects of a numerical implementation, such as the placement of velocities and the associated degrees of freedom, may impact the formation of simulated LKFs. To explore these effects this study compares numerical solutions of sea ice models with different velocity staggering in a benchmark problem. Discretizations based on A‐,B‐, and C‐grid systems on quadrilateral meshes have similar resolution properties as an approximation with an A‐grid staggering on triangular grids (with the same total number of vertices). CD‐grid approximations with a given grid spacing have properties, specifically the number and length of simulated LKFs, that are qualitatively similar to approximations on conventional Arakawa A‐grid, B‐grid, and C‐grid approaches with half the grid spacing or less, making the CD‐discretization more efficient with respect to grid resolution. One reason for this behavior is the fact that the CD‐grid approach has a higher number of degrees of freedom to discretize the velocity field. The higher effective resolution of the CD‐discretization makes it an attractive alternative to conventional discretizations.