Preprints
https://doi.org/10.5194/egusphere-2024-426
https://doi.org/10.5194/egusphere-2024-426
13 Mar 2024
 | 13 Mar 2024
Status: this preprint has been withdrawn by the authors.

A Novel Transformation of the Ice Sheet Stokes Equations and Some of its Properties and Applications

John K. Dukowicz

Abstract. A full-Stokes model provides the most accurate but also the most expensive representation of ice sheet dynamics. The Blatter-Pattyn model is a widely used less expensive approximation that is valid for ice sheets characterized by a small aspect ratio. Here we introduce a novel transformation of the Stokes equations into a form that closely resembles the Blatter-Pattyn equations. The transformed exact Stokes equations only differ from the approximate Blatter-Pattyn equations by a few additional terms, while their variational formulations differ only by the presence of a single term in each horizontal direction (one term in 2D and two terms in 3D). Specifically, the variational formulations differ only by the absence (or the neglect) of the vertical velocity in the second invariant of the strain rate tensor in the Blatter-Pattyn model when compared to the Stokes case. Here we make use of the new transformation in two different ways. First, we consider incorporating the transformed equations into a code that can be very easily converted from a Stokes to a Blatter-Pattyn model, and vice-versa, simply by switching these terms on or off. This may be generalized so that the Stokes model is switched on adaptively only where the Blatter-Pattyn model loses accuracy, hopefully retaining most of the accuracy of the Stokes model but at a lower cost. Second, the key role played by the vertical velocity in converting the transformed Stokes model into the Blatter-Pattyn model motivates new approximations that improve on the Blatter-Pattyn model, heretofore the best approximate ice sheet model. These applications require the use of a grid that enables the discrete continuity equation to be invertible for the vertical velocity in terms of the horizontal velocity components. Examples of such grids, such as the first-order P1-E0 grid and the second-order P2-E1 grid are given in both 2D and 3D. It should be noted, however, that the transformed Stokes model has the same type of gravity forcing as the Blatter-Pattyn model, i.e., determined by the slope of the ice sheet's upper surface, thereby forgoing some of the grid-generality of the traditional formulation of the Stokes model. This is not a serious disadvantage, however, since in practice it has not impaired the widespread use of the Blatter-Pattyn model.

This preprint has been withdrawn.

John K. Dukowicz

Interactive discussion

Status: closed

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • CEC1: 'Comment on egusphere-2024-426', Juan Antonio Añel, 28 Mar 2024
    • AC1: 'Reply on CEC1', John Dukowicz, 28 Mar 2024
      • CEC2: 'Reply on AC1', Juan Antonio Añel, 28 Mar 2024
        • AC2: 'Reply on CEC2', John Dukowicz, 29 Mar 2024
          • CEC3: 'Reply on AC2', Juan Antonio Añel, 29 Mar 2024
            • AC3: 'Reply on CEC3', John Dukowicz, 29 Mar 2024
              • CEC4: 'Reply on AC3', Juan Antonio Añel, 29 Mar 2024
                • AC4: 'Reply on CEC4', John Dukowicz, 01 Apr 2024
  • AC5: 'Comment on egusphere-2024-426', John Dukowicz, 03 Apr 2024

Interactive discussion

Status: closed

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • CEC1: 'Comment on egusphere-2024-426', Juan Antonio Añel, 28 Mar 2024
    • AC1: 'Reply on CEC1', John Dukowicz, 28 Mar 2024
      • CEC2: 'Reply on AC1', Juan Antonio Añel, 28 Mar 2024
        • AC2: 'Reply on CEC2', John Dukowicz, 29 Mar 2024
          • CEC3: 'Reply on AC2', Juan Antonio Añel, 29 Mar 2024
            • AC3: 'Reply on CEC3', John Dukowicz, 29 Mar 2024
              • CEC4: 'Reply on AC3', Juan Antonio Añel, 29 Mar 2024
                • AC4: 'Reply on CEC4', John Dukowicz, 01 Apr 2024
  • AC5: 'Comment on egusphere-2024-426', John Dukowicz, 03 Apr 2024
John K. Dukowicz
John K. Dukowicz

Viewed

Total article views: 218 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
152 42 24 218 6 3
  • HTML: 152
  • PDF: 42
  • XML: 24
  • Total: 218
  • BibTeX: 6
  • EndNote: 3
Views and downloads (calculated since 13 Mar 2024)
Cumulative views and downloads (calculated since 13 Mar 2024)

Viewed (geographical distribution)

Total article views: 243 (including HTML, PDF, and XML) Thereof 243 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 13 Apr 2024
Download

This preprint has been withdrawn.

Short summary
A novel transformation of the rather difficult Stokes ice sheet equations is presented that expands the scope of traditional methods. The new formulation is closely related to a widely used Stokes approximation, the Blatter-Pattyn model, so that an ice sheet model may be easily switched between the two, allowing for adaptive applications. The new formulation also facilitates new approximations that improve on the Blatter-Pattyn model, heretofore the best approximate ice sheet model.