Simulations of the Northern Hemisphere ice sheet at the Last Glacial Maximum (LGM; 21 kyr BP) are performed using a high-resolution atmospheric general circulation model (AGCM) in order to re-evaluate the conventional surface temperature- or elevation-based parameterization. The influence of precipitation change on the steady-state topography of the Laurentide ice sheet at the LGM is estimated using an AGCM with a horizontal resolution of ∼1° and a three-dimensional thermomechanically coupled ice-sheet model. The ice volume estimated by the AGCM simulation is much larger than that indicated by the conventional parameterization. Through sensitivity analysis of the AGCM and ice-sheet model, it is found that the rate of precipitation change depends on the location of the ice sheet, and that the rate of precipitation change due to surface elevation change is higher than the rate unrelated to surface elevation change on the Laurentide ice sheet. The rate of precipitation is also shown to exhibit seasonality and regionality due to effects such as interior desertification and the concentration of storm tracks.