Lagranzhev metod rascheta vozrasta i izotopnogo sostava lda v trekhmernoy modeli Antarkticheskogo ledovogo shita (in Russian), A Lagrangian method of calculation of age and isotopic composition of ice in a 3-D model of the Antarctic ice sheet

Two basic approaches exist for describing advection in the ice sheet model. We proceed with.comparing traditional Eulerian approach with a rather new in glaciological applications Lagrangian one. Particle tracing in the ice sheet is required to deal with such problems as detailed ice dating,establis...

Full description

Bibliographic Details
Main Authors: Rybak, Oleg, Huybrechts, Philippe
Format: Article in Journal/Newspaper
Language:unknown
Published: 2007
Subjects:
Online Access:https://epic.awi.de/id/eprint/19698/
https://epic.awi.de/id/eprint/19698/1/Huy2007c.pdf
https://hdl.handle.net/10013/epic.31597
https://hdl.handle.net/10013/epic.31597.d001
Description
Summary:Two basic approaches exist for describing advection in the ice sheet model. We proceed with.comparing traditional Eulerian approach with a rather new in glaciological applications Lagrangian one. Particle tracing in the ice sheet is required to deal with such problems as detailed ice dating,establishing of isotopical composition, or distribution of any other conservative characteristic, which is advected with the ice flow. Within the frame of the Particle-in-Cell-(PlC) algorithm, a particle's trajectory is constructed by numerical interpolation as it moves through an evolving ice sheet within a regular 3D grid. The algorithm is conceptually straightforward, but demanding in terms of its practical implementation. The main advantage of the algorithm as compared to a pure Eulerian approach is that it is diffusion free, which makes it much more accurate in the lower part of the ice sheet. Application of PIC immediately yields the trajectories of the-particles, which is important for reconstruction of the flow peculiarities and the distribution of transported properties on isochronous surfaces. To implement it in the most effective way it is necessary to follow accurate balance between computational overhead and designed accuracy.