@article{AlbrechtWinkelmannLevermann2020,
  author    = {Albrecht, Torsten and Winkelmann, Ricarda and Levermann, Anders},
  title     = {Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)},
  series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union},
  volume    = {14},
  journal   = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union},
  number    = {2},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {1994-0416},
  doi       = {10.5194/tc-14-599-2020},
  pages     = {599 -- 632},
  year      = {2020},
  abstract  = {Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. <br /> We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general. <br /> For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.},
  language  = {en}
}
@article{FeldmannAlbrechtKhroulevetal.2014,
  author    = {Feldmann, J. and Albrecht, Torsten and Khroulev, C. and Pattyn, F. and Levermann, Anders},
  title     = {Resolution-dependent performance of grounding line motion in a shallow model compared with a full-Stokes model according to the MISMIP3d intercomparison},
  series = {Journal of glaciology},
  volume    = {60},
  journal   = {Journal of glaciology},
  number    = {220},
  publisher = {International Glaciological Society},
  address   = {Cambridge},
  issn      = {0022-1430},
  doi       = {10.3189/2014JoG13J093},
  pages     = {353 -- 360},
  year      = {2014},
  abstract  = {Making confident statements about the evolution of an ice-sheet shelf system with a numerical model requires the capability to reproduce the migration of the grounding line. Here we show that the shallow-ice approximation/shallow-shelf approximation hybrid-type Parallel Ice Sheet Model (PISM), with its recent improvements, is capable of modeling the grounding line motion in a perturbed ice-sheet shelf system. The model is set up according to the three-dimensional Marine Ice-Sheet Model Intercomparison Project (MISMIP3d), and simulations are carried out across a broad range of spatial resolutions. Using (1) a linear interpolation of the grounding line with locally interpolated basal friction and (2) an improved driving-stress computation across the grounding line, the reversibility of the grounding line (i.e. its retreat after an advance forced by a local perturbation of basal resistance) is captured by the model even at medium and low resolutions (Delta x > 10 km). The transient model response is qualitatively similar to that of higher-order models but reveals a higher initial sensitivity to perturbations on very short timescales. Our findings support the application of PISM to the Antarctic ice sheet from regional up to continental scales and on relatively low spatial resolutions.},
  language  = {en}
}