Articles | Volume 8, issue 6
Model description paper 02 Jun 2015
Model description paper | 02 Jun 2015
Mass-conserving subglacial hydrology in the Parallel Ice Sheet Model version 0.6
E. Bueler and W. van Pelt
Related subject area
Florian Herla, Simon Horton, Patrick Mair, and Pascal Haegeli
Geosci. Model Dev., 14, 239–258,Short summary
The adoption of snowpack models in support of avalanche forecasting has been limited. To promote their operational application, we present a numerical method for processing multivariate snow stratigraphy profiles of mixed data types. Our algorithm enables applications like dynamical grouping and summarizing of model simulations, model evaluation, and data assimilation. By emulating the human analysis process, our approach will allow forecasters to familiarly interact with snowpack simulations.
David Pollard and Robert M. DeConto
Geosci. Model Dev., 13, 6481–6500,Short summary
Buttressing by floating ice shelves at ice-sheet grounding lines is an important process that affects ice retreat and whether structural failure occurs in deep bathymetry. Here, we use a simple algorithm to better represent 2-D grounding-line curvature in an ice-sheet model. Along with other enhancements, this improves the performance in idealized-fjord intercomparisons and enables better diagnosis of potential structural failure at future retreating Antarctic grounding lines.
Fuyuki Saito, Takashi Obase, and Ayako Abe-Ouchi
Geosci. Model Dev., 13, 5875–5896,Short summary
The present study introduces the rational function-based constrained interpolation profile (RCIP) method for use in 1 d dating computations in ice sheets and demonstrates the performance of the scheme. Comparisons are examined among the RCIP schemes and the first- and second-order upwind schemes. The results show that, in particular, the RCIP scheme preserves the pattern of input histories, in terms of the profile of internal annual layer thickness, better than the other schemes.
Tobias Sauter, Anselm Arndt, and Christoph Schneider
Geosci. Model Dev., 13, 5645–5662,Short summary
Glacial changes play a key role from a socioeconomic, political, and scientific point of view. Here, we present the open-source coupled snowpack and ice surface energy and mass balance model, which provides a lean, flexible, and user-friendly framework for modeling distributed snow and glacier mass changes. The model provides a suitable platform for sensitivity, detection, and attribution analyses for glacier changes and a tool for quantifying inherent uncertainties.
Alex West, Mat Collins, and Ed Blockley
Geosci. Model Dev., 13, 4845–4868,Short summary
This study calculates sea ice energy fluxes from data produced by ice mass balance buoys (devices measuring ice elevation and temperature). It is shown how the resulting dataset can be used to evaluate a coupled climate model (HadGEM2-ES), with biases in the energy fluxes seen to be consistent with biases in the sea ice state and surface radiation. This method has potential to improve sea ice model evaluation, so as to better understand spread in model simulations of sea ice state.
Eduardo Moreno-Chamarro, Pablo Ortega, and François Massonnet
Geosci. Model Dev., 13, 4773–4787,Short summary
Climate models need to capture sea ice complexity to represent it realistically. Here we assess how distributing sea ice in discrete thickness categories impacts how sea ice variability is simulated in the NEMO3.6–LIM3 model. Simulations and satellite observations are compared by using k-means clustering of sea ice concentration in winter and summer between 1979 and 2014 at both poles. Little improvements in the modeled sea ice lead us to recommend using the standard number of five categories.
Ilkka S. O. Matero, Lauren J. Gregoire, and Ruza F. Ivanovic
Geosci. Model Dev., 13, 4555–4577,Short summary
The Northern Hemisphere cooled by several degrees for a century 8000 years ago due to the collapse of an ice sheet in North America that released large amounts of meltwater into the North Atlantic and slowed down its circulation. We numerically model the ice sheet to understand its evolution during this event. Our results match data thanks to good ice dynamics but depend mostly on surface melt and snowfall. Further work will help us understand how past and future ice melt affects climate.
Martin Rückamp, Angelika Humbert, Thomas Kleiner, Mathieu Morlighem, and Helene Seroussi
Geosci. Model Dev., 13, 4491–4501,Short summary
We present enthalpy formulations within the Ice-Sheet and Sea-Level System model that show better performance than earlier implementations. A first experiment indicates that the treatment of discontinuous conductivities of the solid–fluid system with a geometric mean produce accurate results when applied to coarse vertical resolutions. In a second experiment, we propose a novel stabilization formulation that avoids the problem of thin elements. This method provides accurate and stable results.
C. Max Stevens, Vincent Verjans, Jessica M. D. Lundin, Emma C. Kahle, Annika N. Horlings, Brita I. Horlings, and Edwin D. Waddington
Geosci. Model Dev., 13, 4355–4377,Short summary
Understanding processes in snow (firn), including compaction and airflow, is important for calculating how much mass the ice sheets are losing and for interpreting climate records from ice cores. We have developed the open-source Community Firn Model to simulate these processes. We used it to compare 13 different firn compaction equations and found that they do not agree within 10 %. We also show that including firn compaction in a firn-air model improves the match with data from ice cores.
Rupert Gladstone, Benjamin Galton-Fenzi, David Gwyther, Qin Zhou, Tore Hattermann, Chen Zhao, Lenneke Jong, Yuwei Xia, Xiaoran Guo, Konstantinos Petrakopoulos, Thomas Zwinger, Daniel Shapero, and John Moore
Geosci. Model Dev. Discuss.,
Revised manuscript accepted for GMDShort summary
Retreat of the Antarctic Ice Sheet, and hence its contribution to sea level rise, is highly sensitive to melting of its floating ice shelves. This melt is caused by warm ocean currents coming into contact with the ice. Computer models used for future ice sheet projections are not able to realistically evolve these melt rates. We describe a new coupling framework to enable ice sheet and ocean computer models to interact, allowing projection of the evolution of melt and its impact on sea level.
Bertrand Cluzet, Matthieu Lafaysse, Emmanuel Cosme, Clément Albergel, Louis-François Meunier, and Marie Dumont
Geosci. Model Dev. Discuss.,
Revised manuscript accepted for GMDShort summary
In mountainous areas, snowpack models suffer from large errors and observations are scarce, a challenge for data assimilation. Here, we develop two variants of the Particle Filter in order to propagate the information content of observations over a complex topography. By adjusting observation errors and exploiting background correlation patterns, these variants demonstrate the potential for partial observations of snow depth and surface reflectance to improve model accuracy in a whole domain.
Alexander Robinson, Jorge Alvarez-Solas, Marisa Montoya, Heiko Goelzer, Ralf Greve, and Catherine Ritz
Geosci. Model Dev., 13, 2805–2823,Short summary
Here we describe Yelmo v1.0, an intuitive and state-of-the-art hybrid ice sheet model. The model design and physics are described, and benchmark simulations are provided to validate its performance. Yelmo is a versatile ice sheet model that can be applied to a wide variety of problems.
Shihe Ren, Xi Liang, Qizhen Sun, Hao Yu, L. Bruno Tremblay, Xiaoping Mai, Fu Zhao, Ming Li, Na Liu, Zhikun Chen, and Yunfei Zhang
Geosci. Model Dev. Discuss.,
Preprint under review for GMDShort summary
Sea ice plays a crucial role in global energy and water budge. To get a better simulation of seaice, we coupled seaice model with atmospheric and ocean model to form a fully coupled system. The seaice simulation results of coupled system demonstrated two-way coupled model has better performance in terms of sea ice, especially in summer. This indicates that sea ice-ocean-atmosphere interaction takes a crucial role in controlling Arctic summertime sea ice distribution.
Shiming Xu, Jialiang Ma, Lu Zhou, Yan Zhang, Jiping Liu, and Bin Wang
Geosci. Model Dev. Discuss.,
Revised manuscript accepted for GMDShort summary
A multi-resolution tripolar grid hierarchy is constructed and integrated in CESM (version 1.2.1). The resolution range includes 0.45-deg, 0.15-deg, and 0.05-deg. Based on atmospherically forced sea ice experiments, we show that the model simulates reasonable sea ice kinematics and scaling properties. Besides, landfast ice thickness can be systematically shifted due to non-convergent solutions to EVP rheology. This work serve as a framework of multi-scale modeling of ocean and sea ice with CESM.
Ahmad Jan, Ethan T. Coon, and Scott L. Painter
Geosci. Model Dev., 13, 2259–2276,Short summary
Computer simulations are important tools for understanding the response of Arctic permafrost to a warming climate. To build confidence in an emerging class of permafrost simulators, we evaluated the Advanced Terrestrial Simulator against field observations from a frozen tundra site near Utqiaġvik (formerly Barrow), Alaska. The 3-year simulations agree well with observations of snow depth, summer water table, soil temperature at multiple locations, and spatially averaged evaporation.
Liz C. Logan, Sri Hari Krishna Narayanan, Ralf Greve, and Patrick Heimbach
Geosci. Model Dev., 13, 1845–1864,Short summary
A new capability has been developed for the ice sheet model SICOPOLIS (SImulation COde for POLythermal Ice Sheets) that enables the generation of derivative code, such as tangent linear or adjoint models, by means of algorithmic differentiation. It relies on the source transformation algorithmic (AD) differentiation tool OpenAD. The reverse mode of AD provides the adjoint model, SICOPOLIS-AD, which may be applied for comprehensive sensitivity analyses as well as gradient-based optimization.
Jean-François Lemieux and Frédéric Dupont
Geosci. Model Dev., 13, 1763–1769,Short summary
Sea ice dynamics plays an important role in shaping the sea cover in polar regions. Winds and ocean currents exert large stresses on the sea ice cover. This can lead to the formation of long cracks and ridges, which strongly impact the exchange of heat, momentum and moisture between the atmosphere and the ocean. It is therefore crucial for a sea ice model to be able to represent these features. This article describes how internal sea ice stresses should be diagnosed from model simulations.
Ludovic Räss, Aleksandar Licul, Frédéric Herman, Yury Y. Podladchikov, and Jenny Suckale
Geosci. Model Dev., 13, 955–976,Short summary
Accurate predictions of future sea level rise require numerical models that predict rapidly deforming ice. Localised ice deformation can be captured numerically only with high temporal and spatial resolution. This paper’s goal is to propose a parallel FastICE solver for modelling ice deformation. Our model is particularly useful for improving our process-based understanding of localised ice deformation. Our solver reaches a parallel efficiency of 99 % on GPU-based supercomputers.
Nander Wever, Leonard Rossmann, Nina Maaß, Katherine C. Leonard, Lars Kaleschke, Marcel Nicolaus, and Michael Lehning
Geosci. Model Dev., 13, 99–119,Short summary
Sea ice is an important component of the global climate system. The presence of a snow layer covering sea ice can impact ice mass and energy budgets. The detailed, physics-based, multi-layer snow model SNOWPACK was modified to simulate the snow–sea-ice system, providing simulations of the snow microstructure, water percolation and flooding, and superimposed ice formation. The model is applied to in situ measurements from snow and ice mass-balance buoys installed in the Antarctic Weddell Sea.
Christiaan T. van Dalum, Willem Jan van de Berg, Quentin Libois, Ghislain Picard, and Michiel R. van den Broeke
Geosci. Model Dev., 12, 5157–5175,Short summary
Climate models are often limited to relatively simple snow albedo schemes. Therefore, we have developed the SNOWBAL module to couple a climate model with a physically based wavelength dependent snow albedo model. Using SNOWBAL v1.2 to couple the snow albedo model TARTES with the regional climate model RACMO2 indicates a potential performance gain for the Greenland ice sheet.
François Massonnet, Antoine Barthélemy, Koffi Worou, Thierry Fichefet, Martin Vancoppenolle, Clément Rousset, and Eduardo Moreno-Chamarro
Geosci. Model Dev., 12, 3745–3758,Short summary
Sea ice thickness varies considerably on spatial scales of several meters. However, contemporary climate models cannot resolve such scales yet. This is why sea ice models used in climate models include an ice thickness distribution (ITD) to account for this unresolved variability. Here, we explore with the ocean–sea ice model NEMO3.6-LIM3 the sensitivity of simulated mean Arctic and Antarctic sea ice states to the way the ITD is discretized.
Shahbaz Memon, Dorothée Vallot, Thomas Zwinger, Jan Åström, Helmut Neukirchen, Morris Riedel, and Matthias Book
Geosci. Model Dev., 12, 3001–3015,Short summary
Scientific workflows enable complex scientific computational scenarios, which include data intensive scenarios, parametric executions, and interactive simulations. In this article, we applied the UNICORE workflow management system to automate a formerly hard-coded coupling of a glacier flow model and a calving model, which contain many tasks and dependencies, ranging from pre-processing and data management to repetitive executions on heterogeneous high-performance computing (HPC) resources.
Sébastien Le clec'h, Aurélien Quiquet, Sylvie Charbit, Christophe Dumas, Masa Kageyama, and Catherine Ritz
Geosci. Model Dev., 12, 2481–2499,Short summary
To provide reliable projections of the ice-sheet contribution to future sea-level rise, ice sheet models must be able to simulate the observed ice sheet present-day state. Using a low computational iterative minimisation procedure, based on the adjustment of the basal drag coefficient, we rapidly minimise the errors between the simulated and the observed Greenland ice thickness and ice velocity, and we succeed in stabilising the simulated Greenland ice sheet state under present-day conditions.
Lionel Favier, Nicolas C. Jourdain, Adrian Jenkins, Nacho Merino, Gaël Durand, Olivier Gagliardini, Fabien Gillet-Chaulet, and Pierre Mathiot
Geosci. Model Dev., 12, 2255–2283,Short summary
The melting at the base of floating ice shelves is the main driver of the Antarctic ice sheet current retreat. Here, we use an ideal set-up to assess a wide range of melting parameterisations depending on oceanic properties with regard to a new ocean–ice-sheet coupled model, published here for the first time. A parameterisation that depends quadratically on thermal forcing in both a local and a non-local way yields the best results and needs to be further assessed with more realistic set-ups.
Hakase Hayashida, James R. Christian, Amber M. Holdsworth, Xianmin Hu, Adam H. Monahan, Eric Mortenson, Paul G. Myers, Olivier G. J. Riche, Tessa Sou, and Nadja S. Steiner
Geosci. Model Dev., 12, 1965–1990,Short summary
Ice algae, the primary producer in sea ice, play a fundamental role in shaping marine ecosystems and biogeochemical cycling of key elements in polar regions. In this study, we developed a process-based numerical model component representing sea-ice biogeochemistry for a sea ice–ocean coupled general circulation model. The model developed can be used to simulate the projected changes in sea-ice ecosystems and biogeochemistry in response to on-going rapid decline of the Arctic.
Katherine J. Evans, Joseph H. Kennedy, Dan Lu, Mary M. Forrester, Stephen Price, Jeremy Fyke, Andrew R. Bennett, Matthew J. Hoffman, Irina Tezaur, Charles S. Zender, and Miren Vizcaíno
Geosci. Model Dev., 12, 1067–1086,Short summary
A robust validation of ice sheet models is presented using LIVVkit, version 2.1. It targets ice sheet and coupled Earth system models, and handles datasets and operations that require high-performance computing and storage. We apply LIVVkit to a Greenland ice sheet simulation to show the degree to which it captures the surface mass balance. LIVVkit identifies a positive bias due to insufficient melting compared to observations that is focused largely around Greenland's southwest region.
Jeremy C. Ely, Chris D. Clark, David Small, and Richard C. A. Hindmarsh
Geosci. Model Dev., 12, 933–953,Short summary
During the last 2.6 million years, the Earth's climate has cycled between cold glacials and warm interglacials, causing the growth and retreat of ice sheets. These ice sheets can be independently reconstructed using numerical models or from dated evidence that they leave behind (e.g. sediments, boulders). Here, we present a tool for comparing numerical model simulations with dated ice-sheet material. We demonstrate the utility of this tool by applying it to the last British–Irish ice sheet.
Fabien Maussion, Anton Butenko, Nicolas Champollion, Matthias Dusch, Julia Eis, Kévin Fourteau, Philipp Gregor, Alexander H. Jarosch, Johannes Landmann, Felix Oesterle, Beatriz Recinos, Timo Rothenpieler, Anouk Vlug, Christian T. Wild, and Ben Marzeion
Geosci. Model Dev., 12, 909–931,Short summary
Mountain glaciers are one of the few remaining subsystems of the global climate system for which no globally applicable community-driven model exists. Here we present the Open Global Glacier Model (OGGM; www.oggm.org), developed to provide a modular and open-source numerical model framework for simulating past and future change of any glacier in the world.
William H. Lipscomb, Stephen F. Price, Matthew J. Hoffman, Gunter R. Leguy, Andrew R. Bennett, Sarah L. Bradley, Katherine J. Evans, Jeremy G. Fyke, Joseph H. Kennedy, Mauro Perego, Douglas M. Ranken, William J. Sacks, Andrew G. Salinger, Lauren J. Vargo, and Patrick H. Worley
Geosci. Model Dev., 12, 387–424,Short summary
This paper describes the Community Ice Sheet Model (CISM) version 2.1. CISM solves equations for ice flow, heat conduction, surface melting, and other processes such as basal sliding and iceberg calving. It can be used for ice-sheet-only simulations or as the ice sheet component of the Community Earth System Model. Model solutions have been verified for standard test problems. CISM can efficiently simulate the whole Greenland ice sheet, with results that are broadly consistent with observations.
Thiago Dias dos Santos, Mathieu Morlighem, Hélène Seroussi, Philippe Remy Bernard Devloo, and Jefferson Cardia Simões
Geosci. Model Dev., 12, 215–232,Short summary
The reduction of numerical errors in ice sheet modeling increases the results' accuracy reliability. We improve numerical accuracy by better capturing grounding line dynamics, while maintaining a low computational cost. We implement an adaptive mesh refinement (AMR) technique in the Ice Sheet System Model and compare AMR simulations with uniformly refined meshes. Our results show that the computational time with AMR is significantly shorter than for uniformly refined meshes for a given accuracy.
David Pollard, Robert M. DeConto, and Richard B. Alley
Geosci. Model Dev., 11, 5149–5172,Short summary
Around the margins of ice sheets in contact with the ocean, calving of icebergs can generate large amounts of floating ice debris called "mélange". In major Greenland fjords, mélange significantly slows down ice flow from upstream. Our study applies numerical models to past and possible future episodes of rapid Antarctic Ice Sheet retreat. We find that, due to larger spatial scales, Antarctic mélange does not significantly impede flow or slow ice retreat and associated sea level rise.
Gerhard Krinner, Chris Derksen, Richard Essery, Mark Flanner, Stefan Hagemann, Martyn Clark, Alex Hall, Helmut Rott, Claire Brutel-Vuilmet, Hyungjun Kim, Cécile B. Ménard, Lawrence Mudryk, Chad Thackeray, Libo Wang, Gabriele Arduini, Gianpaolo Balsamo, Paul Bartlett, Julia Boike, Aaron Boone, Frédérique Chéruy, Jeanne Colin, Matthias Cuntz, Yongjiu Dai, Bertrand Decharme, Jeff Derry, Agnès Ducharne, Emanuel Dutra, Xing Fang, Charles Fierz, Josephine Ghattas, Yeugeniy Gusev, Vanessa Haverd, Anna Kontu, Matthieu Lafaysse, Rachel Law, Dave Lawrence, Weiping Li, Thomas Marke, Danny Marks, Martin Ménégoz, Olga Nasonova, Tomoko Nitta, Masashi Niwano, John Pomeroy, Mark S. Raleigh, Gerd Schaedler, Vladimir Semenov, Tanya G. Smirnova, Tobias Stacke, Ulrich Strasser, Sean Svenson, Dmitry Turkov, Tao Wang, Nander Wever, Hua Yuan, Wenyan Zhou, and Dan Zhu
Geosci. Model Dev., 11, 5027–5049,Short summary
This paper provides an overview of a coordinated international experiment to determine the strengths and weaknesses in how climate models treat snow. The models will be assessed at point locations using high-quality reference measurements and globally using satellite-derived datasets. How well climate models simulate snow-related processes is important because changing snow cover is an important part of the global climate system and provides an important freshwater resource for human use.
Aurélien Quiquet, Christophe Dumas, Catherine Ritz, Vincent Peyaud, and Didier M. Roche
Geosci. Model Dev., 11, 5003–5025,Short summary
This paper presents the GRISLI (Grenoble ice sheet and land ice) model in its newest revision. We present the recent model improvements from its original version (Ritz et al., 2001), together with a discussion of the model performance in reproducing the present-day Antarctic ice sheet geometry and the grounding line advances and retreats during the last 400 000 years. We show that GRISLI is a computationally cheap model, able to reproduce the large-scale behaviour of ice sheets.
Gary D. Clow
Geosci. Model Dev., 11, 4889–4908,Short summary
CVPM is a modular heat-transfer modeling system designed for scientific and engineering studies in permafrost terrain, and as an educational tool. CVPM implements the heat-transfer equations in both Cartesian and cylindrical coordinates. To accommodate a diversity of geologic settings, a variety of materials can be specified within the model domain. CVPM can be used over a broad range of depth, temperature, porosity, water saturation, and solute conditions on either Earth or Mars.
Eef C. H. van Dongen, Nina Kirchner, Martin B. van Gijzen, Roderik S. W. van de Wal, Thomas Zwinger, Gong Cheng, Per Lötstedt, and Lina von Sydow
Geosci. Model Dev., 11, 4563–4576,Short summary
Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve. Therefore, approximations to the FS equations are used, especially when modeling an ice sheet on long time spans. Here, we report a combination of an approximation with the FS equations that allows simulating the dynamics of ice sheets over long time spans without introducing artifacts caused by application of approximations in parts of the domain where they are not valid.
Alek A. Petty, Melinda Webster, Linette Boisvert, and Thorsten Markus
Geosci. Model Dev., 11, 4577–4602,
Matthew J. Hoffman, Mauro Perego, Stephen F. Price, William H. Lipscomb, Tong Zhang, Douglas Jacobsen, Irina Tezaur, Andrew G. Salinger, Raymond Tuminaro, and Luca Bertagna
Geosci. Model Dev., 11, 3747–3780,Short summary
MPAS-Albany Land Ice (MALI) is a new variable-resolution land ice model that uses unstructured grids on a plane or sphere. MALI is built for Earth system modeling on high-performance computing platforms using existing software libraries. MALI simulates the evolution of ice thickness, velocity, and temperature, and it includes schemes for simulating iceberg calving and the flow of water beneath ice sheets and its effect on ice sliding. The model is demonstrated for the Antarctic ice sheet.
Mark Kavanagh and Lev Tarasov
Geosci. Model Dev., 11, 3497–3513,Short summary
We present and validate BrAHMs (BAsal Hydrology Model): a new physically based basal hydrology model, which captures the two main types of subglacial drainage systems (high-pressure distributed systems and low-pressure channelized systems). BrAHMs is designed for continental glacial cycle scale contexts, for which computational speed is essential. This speed is accomplished, in part, by numerical methods novel to basal hydrology contexts.
Aleah Sommers, Harihar Rajaram, and Mathieu Morlighem
Geosci. Model Dev., 11, 2955–2974,Short summary
Meltwater drainage beneath glaciers and ice sheets influences how fast they move and is complicated and constantly changing. Most models distinguish between
slowdrainage with different equations for each system. The SHAKTI model allows for the ice–water drainage arrangement to transition naturally between different types of flow. This model can be used to understand how drainage affects glacier speeds and the associated ice loss to further inform predictions of sea level rise.
Ghislain Picard, Melody Sandells, and Henning Löwe
Geosci. Model Dev., 11, 2763–2788,Short summary
The Snow Microwave Radiative Transfer (SMRT) is a novel model developed to calculate how microwaves are scattered and emitted by snow. The model is built from separate, interconnecting modules to make it easy to compare different aspects of the theory. SMRT is the first model to allow a choice of how to represent the microstructure of the snow, which is extremely important, and has been used to unite multiple previous studies. This model will ultimately be used to observe snow from space.
Lihui Luo, Zhongqiong Zhang, Wei Ma, Shuhua Yi, and Yanli Zhuang
Geosci. Model Dev., 11, 2475–2491,Short summary
Based on the current situation of permafrost modeling in the Qinghai–Tibet Plateau (QTP), a software PIC was developed to evaluate the temporal–spatial change trends of permafrost, which allows us to automatically compute permafrost indices with daily weather and atmospheric forcing datasets. The main features include computing, visualization, and statistics. The software will serve engineering applications and can be used to assess the impact of climate change on permafrost over the QTP.
Alexandra Touzeau, Amaëlle Landais, Samuel Morin, Laurent Arnaud, and Ghislain Picard
Geosci. Model Dev., 11, 2393–2418,Short summary
We introduced a new module of water vapor diffusion into the snowpack model Crocus. Vapor transport locally modifies the density of snow layers, possibly influencing compaction. It also affects the original isotopic signature of snow layers. We also introduced water isotopes (𝛿18O) in the model. Over 10 years, the modeled attenuation of isotopic variations due to vapor diffusion is 7–18 % lower than the observations. Thus, other processes are required to explain the total attenuation.
Joshua K. Cuzzone, Mathieu Morlighem, Eric Larour, Nicole Schlegel, and Helene Seroussi
Geosci. Model Dev., 11, 1683–1694,Short summary
This paper details the implementation of higher-order vertical finite elements in the Ice Sheet System Model (ISSM). When using higher-order vertical finite elements, fewer vertical layers are needed to accurately capture the thermal structure in an ice sheet versus a conventional linear vertical interpolation, therefore greatly improving model runtime speeds, particularly in higher-order stress balance ice sheet models. The implications for paleoclimate ice sheet simulations are discussed.
Kaitlin A. Naughten, Katrin J. Meissner, Benjamin K. Galton-Fenzi, Matthew H. England, Ralph Timmermann, Hartmut H. Hellmer, Tore Hattermann, and Jens B. Debernard
Geosci. Model Dev., 11, 1257–1292,Short summary
MetROMS and FESOM are two ocean/sea-ice models which resolve Antarctic ice-shelf cavities and consider thermodynamics at the ice-shelf base. We simulate the period 1992–2016 with both models, and with two options for resolution in FESOM, and compare output from the three simulations. Ice-shelf melt rates, sub-ice-shelf circulation, continental shelf water masses, and sea-ice processes are compared and evaluated against available observations.
Jeff K. Ridley, Edward W. Blockley, Ann B. Keen, Jamie G. L. Rae, Alex E. West, and David Schroeder
Geosci. Model Dev., 11, 713–723,Short summary
The sea ice component of the Met Office coupled climate model, HadGEM3-GC3.1, is presented and evaluated. We determine that the mean state of the sea ice is well reproduced for the Arctic; however, a warm sea surface temperature bias over the Southern Ocean results in a low Antarctic sea ice cover.
Eric Larour, Daniel Cheng, Gilberto Perez, Justin Quinn, Mathieu Morlighem, Bao Duong, Lan Nguyen, Kit Petrie, Silva Harounian, Daria Halkides, and Wayne Hayes
Geosci. Model Dev., 10, 4393–4403,Short summary
This work presents a new way of carrying out simulations using the C++ based Ice Sheet System Model (ISSM) within a web page. This allows for a new generation of websites that can rely on the entire code of a climate model, without compromising or simplifying the physics implemented in such a model. We believe this approach will enable better education/outreach websites as well as improve access to complex climate models without compromising their integrity.
Christopher J. L. D'Amboise, Karsten Müller, Laurent Oxarango, Samuel Morin, and Thomas V. Schuler
Geosci. Model Dev., 10, 3547–3566,Short summary
We present a new water percolation routine added to the Crocus model. The new routine is physically based, describing motion of water through a layered snowpack considering capillary-driven and gravity flow. We tested the routine on two data sets. Wet-snow layers were able to reach higher saturations than the empirical routine. Meaningful applicability is limited until new and better parameterizations of water retention are developed, and feedbacks are adjusted to handle higher saturations.
Mathieu Barrere, Florent Domine, Bertrand Decharme, Samuel Morin, Vincent Vionnet, and Matthieu Lafaysse
Geosci. Model Dev., 10, 3461–3479,Short summary
Global warming projections still suffer from a limited representation of the permafrost–carbon feedback. This study assesses the capacity of snow-soil coupled models to simulate the permafrost thermal regime at Bylot Island, a high Arctic site. Significant flaws are found in the description of Arctic snow properties, resulting in erroneous heat transfers between the soil and the snow in simulations. Improved snow schemes are needed to accurately predict the future of permafrost.
Per Pemberton, Ulrike Löptien, Robinson Hordoir, Anders Höglund, Semjon Schimanke, Lars Axell, and Jari Haapala
Geosci. Model Dev., 10, 3105–3123,Short summary
The Baltic Sea is seasonally ice covered with intense wintertime ship traffic and a sensitive ecosystem. Understanding the sea-ice pack is important for climate effect studies and forecasting. A NEMO-LIM3.6-based model setup for the North Sea/Baltic Sea is introduced, including a method for ice in the coastal zone. We evaluate different sea-ice parameters and overall find that the model agrees well with the observation though deformed ice is more challenging to capture.
Eleanor J. Burke, Sarah E. Chadburn, and Altug Ekici
Geosci. Model Dev., 10, 959–975,Short summary
There is a large amount of relatively inert organic carbon locked into permafrost soils. In a warming climate the permafrost will thaw and this organic carbon will become vulnerable to decomposition. This process is not typically included within Earth system models (ESMs). This paper describes the development of a vertically resolved soil organic carbon decomposition model which, in the future, can be included within the UKESM to quantify the response of the climate to permafrost carbon loss.
Andrews, L. C., Catania, G. A., Hoffman, M. J., Gulley, J. D., Lüthi, M. P., Ryser, C. , Hawley, R. L., and Neumann, T. A.: Direct observations of evolving subglacial drainage beneath the Greenland Ice Sheet, Nature, 514, 80–83, https://doi.org/10.1038/nature13796, 2014.
Ascher, U. and Petzold, L.: Computer Methods for Ordinary Differential Equations and Differential-algebraic Equations, SIAM Press, Philadelphia, PA, 1998.
Aschwanden, A., Bueler, E., Khroulev, C., and Blatter, H.: An enthalpy formulation for glaciers and ice sheets, J. Glaciol., 58, 441–457, https://doi.org/10.3189/2012JoG11J088, 2012.
Aschwanden, A., Adhalgeirsdóttir, G., and Khroulev, C.: Hindcasting to measure ice sheet model sensitivity to initial states, The Cryosphere, 7, 1083–1093, https://doi.org/10.5194/tc-7-1083-2013, 2013.
Balay, S., Abhyankar, S., Adams, M., Brown, J., Brune, P., Buschelman, K., Eijkhout, V., Gropp, W., Kaushik, D., Knepley, M., Curfman McInnes, L., Rupp, K., Smith, B., and Zhang, H.: PETSc Users Manual, Tech. Rep. ANL-95/11 – Revision 3.2, Argonne National Laboratory, 2011.
Bartholomaus, T. C., Anderson, R. S., and Anderson, S. P.: Response of glacier basal motion to transient water storage, Nature Geosci., 1, 33–37, https://doi.org/10.1038/ngeo.2007.52, 2008.
Bartholomaus, T. C., Anderson, R. S., and Anderson, S. P.: Growth and collapse of the distributed subglacial hydrologic system of Kennicott Glacier, Alaska, USA, and its effects on basal motion, J. Glaciol., 57, 985–1002, 2011.
Bindschadler, R., Nowicki, S., Abe-Ouchi, A., Aschwanden, A., Choi, H., Fastook, J., Granzow, G., Greve, R., Gutowski, G., Herzfeld, U., Jackson, C., Johnson, J., Khroulev, C., Levermann, A., Lipscomb, W. H., Martin, M. A., Morlighem, M., Parizek, B. R., Pollard, D., Price, S. F., Ren, D., Saito, F., Sato, T., Seddik, H., Seroussi, H., Takahashi, K., Walker, R., and Wang, W. L.: Ice-sheet model sensitivities to environmental forcing and their use in projecting future sea-level (The SeaRISE Project), J. Glaciol, 59, 195–224, 2013.
Bueler, E.: Correspondence: Extensions of the lumped subglacial-englacial hydrology model of Bartholomaus, et al. (2011), submitted to J. Glaciol., 2014.
Bueler, E. and Brown, J.: Shallow shelf approximation as a "sliding law" in a thermodynamically coupled ice sheet model, J. Geophys. Res., 114, f03008, https://doi.org/10.1029/2008JF001179, 2009.
Bueler, E., Lingle, C. S., Kallen-Brown, J. A., Covey, D. N., and Bowman, L. N.: Exact solutions and numerical verification for isothermal ice sheets, J. Glaciol., 51, 291–306, 2005.
Clarke, G. K.: Hydraulics of subglacial outburst floods: new insights from the Spring-Hutter formulation, J. Glaciol., 49, 299–313, https://doi.org/10.3189/172756503781830728, 2003.
Clarke, G. K. C.: Subglacial processes, Annu. Rev. Earth Planet. Sci., 33, 247–276, https://doi.org/10.1146/annurev.earth.33.092203.122621, 2005.
Creyts, T. and Schoof, C.: Drainage through subglacial water sheets, J. Geophys. Res., 114, https://doi.org/10.1029/2008JF001215, 2009.
Cuffey, K. M. and Paterson, W. S. B.: The Physics of Glaciers, Elsevier, 4th edn., Burlington, MA, USA, 2010.
Das, S. B., Joughin, I., Behn, M. D., Howat, I. M., King, M. A., Lizarralde, D., and Bhatia, M. P.: Fracture Propagation to the Base of the Greenland Ice Sheet During Supraglacial Lake Drainage, Science, 320, 778–781, https://doi.org/10.1126/science.1153360, 2008.
de Fleurian, B., Gagliardini, O., Zwinger, T., Durand, G., Le Meur, E., Mair, D., and Råback, P.: A double continuum hydrological model for glacier applications, The Cryosphere, 8, 137–153, https://doi.org/10.5194/tc-8-137-2014, 2014.
Ettema J., van den Broeke, M. R., van Meijgaard, E., van de Berg, W. J., Bamber, J. L., Box, J. E., and Bales, R. C.: Higher surface mass balance of the Greenland ice sheet revealed by high-resolution climate modeling, Geophys. Res. Lett., 36, L12501, https://doi.org/10.1029/2009GL038110, 2009.
Flowers, G. E. and Clarke, G. K. C.: A multicomponent coupled model of glacier hydrology 1. Theory and synthetic examples, J. Geophys. Res., 107, 2287, https://doi.org/10.1029/2001JB001122, 2002a.
Flowers, G. E. and Clarke, G. K. C.: A multicomponent coupled model of glacier hydrology 2. Application to Trapridge Glacier, Yukon, Canada, J. Geophys. Res., 107, 2288, https://doi.org/10.1029/2001JB001124, 2002b.
Fountain, A., Jacobel, R. W., Schlichting, R., and Jansson, P.: Fractures as the main pathways of water flow in temperate glaciers, Nature 433, 618–621, https://doi.org/https://doi.org/10.1038/nature03296, 2005.
Goeller, S.: Antarctic Subglacial Hydrology: Interactions of subglacial lakes, basal water flow, and ice dynamics, PhD Dissertation, Universität Bremen, 2014.
Goeller, S., Thoma, M., Grosfeld, K., and Miller, H.: A balanced water layer concept for subglacial hydrology in large-scale ice sheet models, The Cryosphere, 7, 1095–1106, https://doi.org/10.5194/tc-7-1095-2013, 2013.
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Advances in Geophysical and Environmental Mechanics and Mathematics, Springer, Berlin, Germany, 2009.
Harper, J., Bradford, J., Humphrey, N., and Meierbachtol, T.: Vertical extension of the subglacial drainage system into basal crevasses, Nature, 467, 579–582, https://doi.org/10.1038/nature09398, 2010.
Hewitt, I. J.: Modelling distributed and channelized subglacial drainage: the spacing of channels, J. Glaciol., 57, 302–314, 2011.
Hewitt, I. J.: Seasonal changes in ice sheet motion due to melt water lubrication, Earth Planet. Sci. Lett., 371–372, 16–25, https://doi.org/10.1016/j.epsl.2013.04.022, 2013.
Hewitt, I. J., Schoof, C., and Werder, M. A.: Flotation and free surface flow in a model for subglacial drainage. Part II: Channel flow, J. Fluid Mech., 702, 157–188, 2012.
Hoffman, M. J. and Price, S.: Feedbacks between coupled subglacial hydrology and glacier dynamics, J. Geophys. Res. Earth Surf., 119, 414–436, https://doi.org/10.1002/2013JF002943, 2014.
Hooke, R., Hanson, B., Iverson, N. R., Jansson, P., and Fischer, U. H.: Rheology of till beneath Störglaciaren, Sweden, J. Glaciol., 43, 172–179, 1997.
Hundsdorfer, W. and Verwer, J. G.: Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations, Springer Series in Computational Mathematics, Springer, Berlin, Germany, 2003.
Huybrechts, P., Payne, T., and The EISMINT Intercomparison Group: The EISMINT benchmarks for testing ice-sheet models, Ann. Glaciol., 23, 1–12, 1996.
Johnson, J. and Fastook, J. L.: Northern Hemisphere glaciation and its sensitivity to basal melt water, Quat. Int., 95, 65–74, 2002.
Joughin, I., Smith, B. E., Howat, I. M., Scambos, T., and Moon, T.: Greenland flow variability from ice-sheet-wide velocity mapping, J. Glaciol., 56, 415–430, 2010.
Kamb, B.: Glacier surge mechanism based on linked cavity configuration of the basal water conduit system, J. Geophys. Res., 92, 9083–9100, 1987.
Kamb, B.: Rheological nonlinearity and flow instability in the deforming bed mechanism of ice stream motion, J. Geophys. Res.-Sol. Ea., 96, 16585–16595, 1991.
Kinderlehrer, D. and Stampacchia, G.: An Introduction to Variational Inequalities and their Applications, Pure and Applied Mathematics, Academic Press, New York, USA, 1980.
Le Brocq, A., Payne, A., Siegert, M., and Alley, R.: A subglacial water-flow model for West Antarctica, J. Glaciol., 55, 879–888, https://doi.org/10.3189/002214309790152564, 2009.
LeVeque, R. J.: Finite Volume Methods for Hyperbolic Problems, Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, United Kingdom, 2002.
Lingle, C. S. and Brown, T. J.: A subglacial aquifer bed model and water pressure-dependent basal sliding relationship for a West Antarctic ice stream, in: Dynamics of the West Antarctic Ice Sheet, edited by: der Veen, C. J. V. and Oerlemans, J., D. Reidel, Dodrecht, the Netherlands, 1987.
Livingstone, S. J., Clark, C. D., Woodward, J., and Kingslake, J.: Potential subglacial lake locations and meltwater drainage pathways beneath the Antarctic and Greenland ice sheets, The Cryosphere, 7, 1721–1740, https://doi.org/10.5194/tc-7-1721-2013, 2013.
Mahaffy, M. W.: A three-dimensional numerical model of ice sheets: tests on the Barnes Ice Cap, Northwest Territories, J. Geophys. Res., 81, 1059–1066, 1976.
Martin, M. A., Winkelmann, R., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., and Levermann, A.: The Potsdam Parallel Ice Sheet Model (PISM-PIK) – Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet, The Cryosphere, 5, 727–740, https://doi.org/10.5194/tc-5-727-2011, 2011.
Morton, K. W. and Mayers, D. F.: Numerical Solutions of Partial Differential Equations: An Introduction, Cambridge University Press, 2nd edn., Cambridge, United Kingdom, 2005.
Nye, J. F.: Water flow in glaciers: Jökulhlaups, tunnels and veins, J. Glaciol., 17, 181–207, 1976.
PISM authors: PISM, a Parallel Ice Sheet Model: User's Manual, http://www.pism-docs.org, last access: 12 May 2015.
Rooney, S. T., Blankenship, D. D., Alley, R. B., and Bentley, C. R.: Till beneath ice stream B: 2. structure and continuity, J. Geophys. Res.-Sol. Ea., 92, 8913–8920, 1987.
Schoof, C.: The effect of cavitation on glacier sliding, P. R. Soc. A, 461, 609–627, https://doi.org/10.1098/rspa.2004.1350, 2005.
Schoof, C.: A variational approach to ice stream flow, J. Fluid Mech., 556, 227–251, 2006a.
Schoof, C.: Variational methods for glacier flow over plastic till, J. Fluid Mech., 555, 299–320, 2006b.
Schoof, C.: Cavitation on deformable glacier beds, SIAM J. Appl. Math., 67, 1633–1653, 2007.
Schoof, C.: Coulomb friction and other sliding laws in a higher order glacier flow model, Math. Models Methods Appl. Sci. (M3AS), 20, 157–189, https://doi.org/10.1142/S0218202510004180, 2010a.
Schoof, C.: Ice sheet acceleration driven by melt supply variability, Nature, 468, 803–806, 2010b.
Schoof, C., Hewitt, I. J., and Werder, M. A.: Flotation and free surface flow in a model for subglacial drainage. Part I: Distributed drainage, J. Fluid Mech., 702, 126–156, 2012.
Shreve, R.: Movement of water in glaciers, J. Glaciol, 11, 205–214, 1972.
Siegert, M., Le Brocq, A., and Payne, A.: Hydrological connections between Antarctic subglacial lakes, the flow of water beneath the East Antarctic Ice Sheet and implications for sedimentary processes, in: Glacial Sedimentary Processes and Products, edited by: Hambrey, M. J., Christoffersen, P., Glasser, N. F., and Hubbard, B., 3–10, Wiley-Blackwell, Malden, MA, USA, 2007.
Truffer, M. and Harrison, W.: In situ measurements of till deformation and water pressure, J. Glaciol., 52, 175–182, 2006.
Truffer, M., Echelmeyer, K., and Harrison, W.: Glacier motion dominated by processes deep in underlying till, J. Glaciol., 46, 213–221, 2000.
Truffer, M., Echelmeyer, K., and Harrison, W.: Implications of till deformation on glacier dynamics, J. Glaciol., 47, 123–134, https://doi.org/10.3189/172756501781832449, 2001.
Tulaczyk, S., Kamb, W. B., and Engelhardt, H. F.: Basal mechanics of Ice Stream B, West Antarctica 1. Till mechanics, J. Geophys. Res., 105, 463–481, 2000a.
Tulaczyk, S., Kamb, W. B., and Engelhardt, H. F.: Basal mechanics of Ice Stream B, West Antarctica 2. Undrained plastic bed model, J. Geophys. Res., 105, 483–494, 2000b.
Van der Wel, N., Christoffersen, P., and Bougamont, M.: The influence of subglacial hydrology on the flow of Kamb Ice Stream, West Antarctica, J. Geophys. Res.-Earth, 118, 1–14, https://doi.org/10.1029/2012JF002570, 2013.
Van Pelt, W.: Modelling the dynamics and boundary processes of Svalbard glaciers, PhD thesis, Institute for Marine and Atmospheric Research Utrecht (IMAU), The Netherlands, 2013.
Vázquez, J. L.: The Porous Medium Equation, Oxford Mathematical Monographs, The Clarendon Press Oxford University Press, Oxford, 2007.
Walder, J. S.: Stability of sheet flow of water beneath temperate glaciers and implications for glacier surging, J. Glaciol., 28, 273–293, 1982.
Werder, M., Hewitt, I., Schoof, C., and Flowers, G.: Modeling channelized and distributed subglacial drainage in two dimensions, J Geophys. Res.-Earth, 118, 2140–2158, 2013.
Wesseling, P.: Principles of Computational Fluid Dynamics, Springer-Verlag, Berlin, Germany, 2001.
Winkelmann, R., Martin, M. A., Haseloff, M., Albrecht, T., Bueler, E., Khroulev, C., and Levermann, A.: The Potsdam Parallel Ice Sheet Model (PISM-PIK) – Part 1: Model description, The Cryosphere, 5, 715–726, https://doi.org/10.5194/tc-5-715-2011, 2011.
We describe and test a parallel 2-D subglacial hydrology model. Such a hydrology model describes the motion of liquid water, and its pressure, underneath ice sheets and glaciers. Model goals include conserving the mass of water and generating an estimate of the ice basal shear stress. This hydrology model describes the subglacier as a layer of till and linked cavities, but it does not include conduits.
We describe and test a parallel 2-D subglacial hydrology model. Such a hydrology model...