Skip to main content Accessibility help
×
Home

A smoothed-particle hydrodynamics model for ice-sheet and ice-shelf dynamics

  • Wenxiao Pan (a1), Alexandre M. Tartakovsky (a1) and Joe J. Monaghan (a2)

Abstract

Mathematical modeling of ice sheets is complicated by the nonlinearity of the governing equations and boundary conditions. Standard grid-based methods require complex front-tracking techniques and have a limited capability to handle large material deformations and abrupt changes in bottom topography. Consequently, numerical methods are usually restricted to shallow ice-sheet and ice-shelf approximations. We propose a new smoothed-particle hydrodynamics (SPH) model for coupled ice-sheet and ice-shelf dynamics. SPH, a fully Lagrangian particle method, is highly scalable and its Lagrangian nature and meshless discretization are well suited to the simulation of free surface flows, large material deformation and material fragmentation. In this paper, we use the SPH model to study ice-sheet/ice-shelf behavior, and the dynamics of the grounding line. The steady-state position of the grounding line obtained from SPH simulations is in good agreement with laboratory observations for a wide range of simulated bedrock slopes and density ratios, similar to those of ice and sea water. The numerical accuracy of the SPH algorithm is verified by simulating the plane-shear flow of two immiscible fluids and the propagation of a highly viscous blob of fluid along a horizontal surface. In the experiment, the ice was represented with a viscous Newtonian fluid. For consistency, in the described SPH model the ice is also modeled as a viscous Newtonian fluid. Typically, ice sheets are modeled as a non-Newtonian fluid, accounting for the changes in the mechanical properties of the ice. Implementation of a non-Newtonian rheology in the SPH model is the subject of our ongoing research.

  • View HTML
    • Send article to Kindle

      To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

      Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

      Find out more about the Kindle Personal Document Service.

      A smoothed-particle hydrodynamics model for ice-sheet and ice-shelf dynamics
      Available formats
      ×

      Send article to Dropbox

      To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

      A smoothed-particle hydrodynamics model for ice-sheet and ice-shelf dynamics
      Available formats
      ×

      Send article to Google Drive

      To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

      A smoothed-particle hydrodynamics model for ice-sheet and ice-shelf dynamics
      Available formats
      ×

Copyright

References

Hide All
Bueler, E and Brown, J (2009) Shallow shelf approximation as a ‘sliding law’ in a thermomechanically coupled ice sheet model. J. Geophys. Res, 114(F3), F03008 (doi: 10.1029/2008JF001179)
Bueler, E, Brown, J and Lingle, C (2007) Exact solutions to the thermomechanically coupled shallow-ice approximation: effective tools for verification. J. Glaciol., 53(182), 499-516
Cleary, PW and Monaghan, JJ (1999) Conduction modelling using smoothed particle hydrodynamics. J. Comput. Phys., 148(1), 227-264
Colagrossi, A and Landrini, M (2003) Numerical simulation of interfacial flows by smoothed particle hydrodynamics. J. Comput. Phys., 191(2), 448-475
Durand, C, Cagliardini, O, Zwinger, T, Le Meur, E and Hindmarsh, RCA (2009) Full Stokes modeling of marine ice sheets: influence of the grid size. Ann. Glaciol., 50(52), 109-114
Cagliardini, O, Cohen, D, Raback, P and Zwinger, T (2007) Finite- element modeling of subglacial cavities and related friction law. J. Geophys. Res., 112(F2), F02027 (doi: 10.1029/2006JF000576)
Goldsby, DL and Kohlsted, DL (2001) Superplastic deformation of ice: experimental observations. J. Geophys. Res., 106(B6), 11 017-11 030
Hosseini, SM, Manzari, MT and Hannani, SK (2007) A fully explicit three-step SPH algorithm for simulation of non-Newtonian fluid flow. Int. J. Num. Meth. Heat Fluid Flow, 17(7), 715-735
Huppert, HE (1982) The propagation of two-dimensional and axisymmetric viscous gravity currents over a rigid horizontal surface. J. Fluid Mech., 121, 43-58
Huybrechts, P, Payne, T and the EISMINT Intercomparison Group (1996) The EISMINT benchmarks for testing ice-sheet models. Ann. Glaciol., 23, 1-12
Marshall, SJ (2005) Recent advances in understanding ice sheet dynamics. Earth Planet. Sci. Lett., 240(2), 191-204
Monaghan, JJ (1992) Smoothed particle hydrodynamics. Annu. Rev. Astron. Astrophys., 30, 543-574
Monaghan, JJ (1994) Simulating free surface flows with SPH. J. Comput. Phys., 110(2), 399-406
Monaghan, JJ (2005) Smoothed particle hydrodynamics. Rep. Prog. Phys., 68(8), 1703-1759
Monaghan, JJ and Kajtar, JB (2009) SPH particle boundary forces for arbitrary boundaries. Comput. Phys. Commun., 180(10), 1811-1820
Monaghan, JJ, Huppert, HE and Worster, MG (2005) Solidification using smoothed particle hydrodynamics. J. Comput. Phys., 206(2), 684-705
Morland, LW and Johnson, IR (1980) Steady motion of ice sheets. J. Glaciol., 25(92), 229-246
Pollard, D and DeConto, RM (2009) Modelling West Antarctic ice sheet growth and collapse through the past five million years. Nature, 458(7236), 329-332
Price, SF, Payne, AJ, Howat, IM and Smith, BE (2011) Committed sea-level rise for the next century from Greenland ice sheet dynamics during the past decade. Proc. Natl Acad. Sci. USA (PNAS), 108(22), 8978-8983
Robison, RAV, Huppert, HE and Worster, MG (2010) Dynamics of viscous grounding lines. J. Fluid Mech., 648, 363-380
Ryan, EM, Tartakovsky, AM and Amon, C (2010) A novel method for modeling Neumann and Robin boundary conditions in smoothed particle hydrodynamics. Comput. Phys. Commun., 181(12), 2008-2023
Seddik, H, Greve, R, Zwinger, T and Placidi, L (2011) A full Stokes ice flow model for the vicinity of Dome Fuji, Antarctica, with induced anisotropy and fabric evolution. Cryosphere, 5(2), 495-508
Seroussi, H and 6 others (2011) Ice flux divergence anomalies on 79north Glacier, Greenland. Geophys. Res. Lett., 38(9), L09501 (doi: 10.1029/2011GL047338)
Shao, S and Lo, EYM (2003) Incompressible SPH method for simulating Newtonian and non-Newtonian flows with a free surface. Adv. Water Resour., 26(7), 787-800
Tartakovsky, AM and Meakin, P (2005) Simulation of unsaturated flow in complex fractures using smoothed particle hydrodynamics. Vadose Zone J., 4(3), 848-855
Tartakovsky, AM, Meakin, P, Scheibe, TD and Eichler West, RM (2007) Simulations of reactive transport and precipitation with smoothed particle hydrodynamics. J. Comput. Phys., 222(2), 654-672
Tartakovsky, AM, Tartakovsky, DM, Scheibe, TD and Meakin, P (2008) Hybrid simulations of reaction-diffusion systems in porous media. SIAM J. Sci. Comput., 30(6), 2799-2816
Tartakovsky, AM, Ferris, KM and Meakin, P (2009) Lagrangian particle model for multiphase flows. Comput. Phys. Commun., 180(10), 1874-1881
Weis, M, Greve, R and Hutter, K (1999) Theory of shallow ice shelves. Contin. Mech. Thermodyn., 11(1), 15-50

Metrics

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed