On post-glacial sea level – II. Numerical formulation and comparative results on spherically symmetric models

TL;DR: In this article, Mitrovica and Milne extended their previous work in three ways: first, they derived an efficient, iterative numerical algorithm for solving the generalized sea-level equation, and then they considered a special case of the new sea level theory involving spherically symmetric earth models.

Abstract: SUMMARY Theoretical approaches to computing gravitationally self-consistent sea-level changes in consequence of ice growth and ablation are comprised of two parts. The first is a mapping between variations in global sea level and changes in ocean height (required to define the surface load), and the second is a method for computing global sea-level change arising from an arbitrary surface loading. In Mitrovica & Milne (2003) (Paper I) we described a new, generalized mapping between sea-level change and ocean height that takes exact account of the evolution of shorelines associated with both transgression and regression cycles and time-dependent marine-based ice margins. The theory is valid for any earth model. In this paper we extend our previous work in three ways. First, we derive an efficient, iterative numerical algorithm for solving the generalized sea-level equation. Secondly, we consider a special case of the new sea-level theory involving spherically symmetric earth models. Specifically, we combine our iterative numerical formulation with viscoelastic Love number theory to derive an extended pseudo-spectral algorithm for solving the new sea-level equation. This algorithm represents an extension of earlier methods developed for the fixed-shoreline case to precisely incorporate shoreline migration processes. Finally, using this special case, we quantitatively assess errors incurred in previous efforts to extend the traditional (fixed shoreline) sea-level equation of Farrell & Clark (1976) to treat time-dependent shorelines. We find that the approximations adopted by Johnston (1993) and Milne (1998) to treat transgression and regression at shorelines introduce negligible (∼1 per cent) error into predictions of post-glacial relative sea-level histories. In contrast, the errors associated with the Peltier (1994) sea-level equation are an order of magnitude larger, and comparable to the error incurred using the traditional sea-level theory. Furthermore, our numerical tests verify the high accuracy of the Milne (1998) approximation for treating the influence of grounded, marine-based ice.