General Circulation Models (GCMs) are powerful tools for understanding the meteorology and climate of the high latitudes, which are among the most sensitive regions on Earth to global and environmental change. GCMs differ vastly in their numerical treatment of polar regions because of the so-called pole problem (Williamson, 2007). The pole problem refers to numerical instability arising from the convergence of meridian lines into polar singularities on latitude-longitude grids (e.g., Figure 1a, hereafter referred to as lat-lon grids). Depending on the numerics, methods exist to suppress this instability, and lat-lon grids may be advantageous for polar processes by representing structures with finer resolution than elsewhere in the computational domain. With the recent trend toward quasi-uniform unstructured grids, any potential benefits of lat-lon grids in polar regions may be lost (hereafter, quasi-uniform refers to approximately isotropic grids, in contrast to lat-lon grids, which are highly anisotropic due to the polar singularity). In this study, we evaluate a number of grids and dynamical cores (hereafter referred to as dycores) available in the Community Earth System Model, version 2.2 (CESM2.2; Danabasoglu et al., 2020), including new variable-resolution grids (i.e., grids with enhanced resolution over a particular region), to understand their impacts on the simulated Arctic climate. We focus specifically on the climate and surface mass balance of the Greenland Ice Sheet.