We present a series of numerical reconstructions of the Laurentide Ice Sheet during the last deglaciation (18–714Cka) that evaluates the sensitivity of ice-sheet geometry to subglacial sediment deformation. These reconstructions assume that the Laurentide Ice Sheet flowed over extensive areas of water-saturated, deforming sediment (soft beds) corresponding to the St. Lawrence lowland, the Great Lakes region, the western prairies of the U.S. and Canada, and the Hudson Bay and Hudson Strait regions. Sediment rheology is based on a constitutive law that incorporates experimental results from late Wisconsin till deposited by the Laurentide Ice Sheet which suggest only mildly nonlinear viscoplastic behavior. By varying the effective viscosity of till, we produced four reconstructions for the ice sheet during the last glacial maximum 1814Cka, and two reconstructions each of the ice sheet at 14, 13, 12, 11 and 1014Cka. We also produced one reconstruction for 9, 8.4, 8, and 714Cka. Reconstructions that assume a low effective viscosity for all areas of deforming sediment show a multidomed ice sheet with a large bowl-shaped depression over Hudson Bay and thin ice (<1000m above modern sea-level) over the western and southern margins. Those reconstructions that assume a higher effective viscosity of till in the Hudson Bay region than for the western and southern margins also show a multidomed ice sheet but with considerably thicker ice over Hudson Bay and a more northerly position of the central ice divide. These two different geometries may represent ice-sheet orographic changes associated with a Heinrich event. Further increases in effective viscosity of till, approaching the effective viscosity of ice, would result in a high, monolithic ice dome centered over Hudson Bay, reinforcing the notion that a multidomed ice sheet reflects the distribution of substrate geology.Modeled ice-surface geometry at the last glacial maximum shows many of the same general features as previous reconstructions that incorporate deformable beds. Our reconstructions with higher effective till viscosities in Hudson Bay also agree with the ICE-4G reconstructions (Peltier, 1994), which are based on inversion of relative sea-level data, for the early part of the last deglaciation (18–1314Cka), but then depart significantly from ICE-4G beginning at about 1214Cka due to differing assumptions of the history of deglaciation. Modeled ice volume for the last glacial maximum suggests a glacioeustatic change of 50–55m by a soft-bedded Laurentide Ice Sheet; this would increase as the effective viscosity of till increases. Subsequent ice-volume changes through the last deglaciation generally parallel the trend of eustatic rise recorded at Barbados, New Guinea, and Tahiti, but suggest that the Laurentide Ice Sheet was not the source of meltwater pulse 1A.