Long‐Wavelength Gravity Field Constraint on the Lower Mantle Viscosity in North America

The long‐wavelength negative gravity anomaly over Hudson Bay coincides with the area depressed by the Laurentide Ice Sheet during the Last Glacial Maximum, suggesting that it is, at least partly, caused by glacial isostatic adjustment (GIA). Additional contributions to the static gravity field stem from surface dynamic topography, core‐mantle boundary (CMB) topography, and density anomalies in the subsurface. Previous estimates of the contribution of GIA to the gravity anomaly range from 25% to more than 80%. However, these estimates did not include uncertainties in all components that contribute to the gravity field. In this study, we develop a forward model for the gravity anomaly based on density models and dynamic models, investigating uncertainty in all components. We derive lithospheric densities from equilibrium constraints but extend the concept of lithospheric isostasy to a force balance that includes the dynamic models. The largest uncertainty in the predicted gravity anomaly is due to the lower mantle viscosity, uncertainties in the ice history, the crustal model, the lithosphere‐asthenosphere boundary, and the conversion from seismic velocities to density have a smaller effect. A preference for lower mantle viscosities >1022 Pa s is found, in which case at least 60% of the observed long‐wavelength gravity anomaly can be attributed to GIA. This lower bound on the lower mantle viscosity has implications for inferences based on models for mantle convection and GIA.