Abstract
Abstract. Ice-sheet age computations are formulated using an Eulerian advection equation, and there are many schemes that can be used to solve them numerically. Typically, these differ in numerical characteristics such as stability, accuracy, and diffusivity. Furthermore, although various methods have been presented for ice-sheet age computations, the constrained interpolation profile method and its variants have not been examined in this context. The present study introduces one of its variants, a rational function-based constrained interpolation profile (RCIP) scheme, to one-dimensional ice age computation and demonstrates its performance levels via comparisons with those obtained from first- and second-order upwind schemes. Our results show that the RCIP scheme preserves the pattern of input surface mass balance histories in terms of the vertical profile of internal annual layer thickness better than the other schemes.
Highlights
Core samples extracted from ice sheets can provide an archive of past climate history data, and a major issue for researchers attempting to utilize ice-core properties is defining the age of ice along the depth of the ice sheet
The experiments in the present study were limited to 1-D computations under summits, we believe the characteristics of the RCIP schemes have been presented sufficiently to allow evaluations of their performance levels
The computed vertical profiles of the annual layer thicknesses produced by RCIP schemes follow the expected depth profiles more reasonably than the other methods
Summary
Core samples extracted from ice sheets can provide an archive of past climate history data, and a major issue for researchers attempting to utilize ice-core properties is defining the age of ice along the depth of the ice sheet. In Rybak and Huybrechts (2003), the authors compared the Lagrangian and Eulerian schemes for simulated Antarctic ice sheets under various schematic steadystate conditions and analytical solutions, as well as under different 3-D velocity fields. They concluded that the Lagrangian method produced less error than an Euler approach, the difference was small over a large part of the domain. The present study introduces a CIP method variant named the rational function-based constrained interpolation profile (RCIP) method (Xiao et al, 1996) for use in 1-D ice age computations and demonstrates the performance of the scheme
Sign-in/Register to view highlights & summary Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
