We present new inferences of the radial profile of mantle viscosity that simultaneously fit long-wavelength free-air gravity harmonics associated with mantle convection and a large set of decay times estimated from the postglacial uplift of sites within previously glaciated regions (Hudson Bay, Arctic Canada, and Fennoscandia). The relative sea level variation at these latter sites is constrained by age-height pairs obtained by geological survey, rather than the subjective trends which are commonly used in glacial isostatic adjustment (GIA) studies. Our viscosity inferences are generated using two approaches. First, we adopt a relative viscosity profile which is known to provide a good fit to the free-air gravity harmonics and determine an absolute scaling which yields a best fit to the GIA decay time constraints. Second, we perform an iterative, nonlinear, joint inversion of the two data sets. In both cases our inferred profiles are characterized by a significant increase of viscosity (~2 orders of magnitude), with depth, to values of ~1022 Pa s in the bottom half of the lower mantle. The new viscosity profiles are shown to satisfy constraints based on the postglacial uplift of both Fennoscandia (the classic Haskell <1935> number) and Hudson Bay which have commonly been invoked to argue for an isoviscous mantle. Furthermore, the models are used to predict a set of long-wavelength signatures of the GIA process. These include predictions of GIA-induced variations in (1) the length-of-day over the late Holocene period; (2) the Earth's precession constant and obliquity over the last 2.6 Myr; and (3) the present-day zonal harmonics of the geopotential, J˙l(l≤7). The predictions (1) and (3) bound the late Holocene (and ongoing) mass flux between the large polar ice sheets (Greenland and Antarctic) and the global oceans to small values (≤0.4 mm/yr equivalent eustatic sea level rise).¿ 1997 American Geophysical Union