Studies of glacial isostatic adjustment (GIA) based on spherically symmetric viscoelastic Earth models have argued that the rate of change of the degree 2 zonal harmonic of the Earth’s geopotential, or J˙2, provides an important constraint on mean viscosity in the deep mantle (Mitrovica and Peltier, 1993; Nakada et al., 2015; Lau et al., 2016). To refine this constraint, we compute Fr´echet kernels using an adjoint methodology that reveal the sensitivity of the datum to 3D variations in mantle viscosity. We demonstrate that the mantle sensitivity of the datum is largely limited to the region below the ancient Laurentide ice sheet that covered Canada and significant portions of the northeastern United States at Last Glacial Maximum (LGM). In the bottom half of the lower mantle, this region of maximum sensitivity lies outside the location of Large Low Shear Velocity Provinces (LLSVPs) imaged from seismic tomographic studies. Thus, if the low shear velocity of these provinces originates from thermal effects, previous inferences of viscosity based upon the J˙ 2 datum are likely higher than the actual mean viscosity of the lower mantle.