Constraints on mantle viscosity and Laurentide ice sheet evolution from pluvial paleolake shorelines in the western United States

Abstract

The deformation pattern of the paleoshorelines of extinct Lake Bonneville were among the first features to indicate that Earth's interior responds viscoelastically to changes in surface loads (Gilbert, 1885). Here we revisit and extend this classic study of isostatic rebound with updated lake chronologies for Lake Bonneville and Lake Lahontan as well as revised elevation datasets of shoreline features. The first order domal pattern in the shoreline elevations can be explained by rebound associated with the removal of the lake load. We employ an iterative scheme to calculate the viscoelastic lake rebound, which accounts for the deformation of the solid Earth and gravity field, to calculate a lake load that is consistent with the load-deformed paleotopography. We find that the domal deformation requires a regional Earth structure that exhibits a thin elastic thickness of the lithosphere (15–25 km) and low sublithospheric Maxwell viscosity (~1019 Pa s). After correcting for rebound due to the lake load, shoreline feature elevations reveal a statistically significant northward dipping trend. We attribute this trend to continent-scale deformation caused by the ice peripheral bulge of the Laurentide ice sheet, and take advantage of the position of these lakes on the distal flank of the peripheral bulge to provide new insights on mantle viscosity and Laurentide ice sheet reconstructions. We perform ice loading calculations to quantify the deformation of the solid Earth, gravity field, and rotation axis that is caused by the growth and demise of the Laurentide ice sheet. We test three different ice reconstructions paired with a suite of viscosity profiles and confirm that the revealed trend can be explained by deformation associated with the Laurentide ice sheet when low viscosities below the asthenosphere are adopted. We obtain best fits to shoreline data using ice models that do not have the majority of ice in the eastern sectors of the Laurentide ice sheet, with the caveat that this result can be affected by lateral variations in viscosity. We show that pluvial lakes in the western United States can place valuable constraints on the Laurentide ice sheet, the shape of its peripheral bulge, and the underlying mantle viscosity.