Near-field hydro-isostasy: the implementation of a revised sea-level equation

G. A. Milne,J. L. Davis,J. X. Mitrovica

doi:10.1046/j.1365-246x.1999.00971.x

G. A. Milne, J. L. Davis + Show 1 more

Open Access

https://doi.org/10.1046/j.1365-246x.1999.00971.x

Copy DOI

Abstract

Summary We describe how the existing sea-level equation incorrectly predicts the change in sea level (and thus the ocean load) in ice-covered, subgeoidal geographic regions during periods of deglaciation. We go on to present a new sea-level equation that overcomes this problem and we describe how this equation can be solved in a gravitationally self-consistent manner by employing a well-known spectral technique. Application of the new theory to predict relative sea-level (rsl) histories and present-day, 3-D, solid surface deformation rates in northeastern Canada (based on a single earth model characterized by a lithospheric thickness of 100 km and upper and lower mantle viscosities of 5 × 1020 and 5 × 1021 Pa s, respectively) demonstrates that a significant error is introduced when the original theory is employed to predict the oceanic component of the surface load. Predictions of rsl curves show a discrepancy of ~ 40 per cent at sites where data have been obtained and employed to constrain models of earth viscosity structure and ice-sheet histories. This error will significantly bias estimates of mantle viscosity structure and ice thicknesses that are based on the original theory. In addition, predictions of 3-D deformation rates differ by up to 25 per cent in some regions and so future applications that employ these data to constrain models of the glacial isostatic adjustment process should adopt the improved sea-level theory. In contrast, estimates of inverse decay times from the predicted rsl curves are insensitive (to within the observational error) to the improvement in the surface load introduced by the new theory. Thus, viscosity structure inferences based on this parametrization and the original sea-level equation are unaffected by the error in the ocean load. Finally, the new theory predicts a eustatic (i.e. globally uniform) rise in sea level over the postglacial period that is ~ 11 m lower than that determined via the original theory. Therefore, estimates of the global ice budget at the last glacial maximum based on far-field rsl data and the original sea-level theory will be too small by ~ 10 per cent.

Highlights

The response of the Earth to the surface ice^water mass redistribution induced by the Late Pleistocene glacial cycles is manifest in a number of geophysical observables: variations in sea level, anomalies in the geopotential, secular variation of the Earth's rotation vector and 3-D deformation of the solid surface
We describe how the existing sea-level equation incorrectly predicts the change in sea level in ice-covered, subgeoidal geographic regions during periods of deglaciation
Application of the new theory to predict relative sea-level histories and present-day, 3-D, solid surface deformation rates in northeastern Canada demonstrates that a signi¢cant error is introduced when the original theory is employed to predict the oceanic component of the surface load

Summary

This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Other Posted Material, as set forth at http:// Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-42, Cambridge, MA 02138, USA Department of Physics, University of Toronto, 60 St. George Street, Toronto, Ontario, M5S 1A7, Canada Key words: crustal deformation, glacial rebound, sea level.

IN T ROD UCTION

SUMMA RY

Findings

REFERE NCE S