Articles | Volume 11, issue 8
Development and technical paper
30 Aug 2018
Development and technical paper | 30 Aug 2018
BrAHMs V1.0: a fast, physically based subglacial hydrology model for continental-scale application
Mark Kavanagh and Lev Tarasov
No articles found.
Matthew Drew and Lev Tarasov
The Cryosphere Discuss.,
Preprint under review for TCShort summary
Interaction of fast flowing regions of continental ice sheets with their beds governs how quickly they slide and therefore flow. The coupling of fast ice to its bed is controlled by the pressure of melt water at its base. It is currently poorly understood how the physical details of these hydrologic systems affect ice speed up. Using numerical models we find, surprizingly, that they largely do not – except for the duration of the surge – suggesting that cheap models are sufficient.
Lev Tarasov and Michael Goldstein
Clim. Past Discuss.,
Revised manuscript not acceptedShort summary
This review: 1. Illustrates how current climate and/or ice sheet model-based inferences about the past tend to have limited interpretable value about the real world given inadequate accounting of uncertainties. 2. Explains Bayesian inference to a non-statistical community. 3. Sketches out tractable Bayesian inference for computationally expensive models in a way that meaningfully accounts for uncertainties. 4. Lays out some steps for the community to move forward.
Ryan Love, Heather J. Andres, Alan Condron, and Lev Tarasov
Clim. Past, 17, 2327–2341,Short summary
Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability and climate transitions. We track the routing of glaciologically constrained freshwater volumes in glacial ocean simulations. Our simulations capture important generally not well-represented small-scale features (boundary currents, eddies). We show that the dilution of freshwater as it is transported to key climate regions reduces the freshening to 20 %–60 %.
Taimaz Bahadory, Lev Tarasov, and Heather Andres
Clim. Past, 17, 397–418,Short summary
We present an ensemble of last glacial inception simulations using a fully coupled ice–climate model for the Northern Hemisphere. The ensemble largely captures inferred ice volume changes within proxy uncertainties. Notable features include an ice bridge across Davis Strait and between Greenland and Iceland. Via an equilibrium climate response experiment, we also demonstrate the potential value of fully coupled ice–climate modelling of last glacial inception to constrain future climate change.
Heiko Goelzer, Sophie Nowicki, Anthony Payne, Eric Larour, Helene Seroussi, William H. Lipscomb, Jonathan Gregory, Ayako Abe-Ouchi, Andrew Shepherd, Erika Simon, Cécile Agosta, Patrick Alexander, Andy Aschwanden, Alice Barthel, Reinhard Calov, Christopher Chambers, Youngmin Choi, Joshua Cuzzone, Christophe Dumas, Tamsin Edwards, Denis Felikson, Xavier Fettweis, Nicholas R. Golledge, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Sebastien Le clec'h, Victoria Lee, Gunter Leguy, Chris Little, Daniel P. Lowry, Mathieu Morlighem, Isabel Nias, Aurelien Quiquet, Martin Rückamp, Nicole-Jeanne Schlegel, Donald A. Slater, Robin S. Smith, Fiamma Straneo, Lev Tarasov, Roderik van de Wal, and Michiel van den Broeke
The Cryosphere, 14, 3071–3096,Short summary
In this paper we use a large ensemble of Greenland ice sheet models forced by six different global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The results for two different greenhouse gas concentration scenarios indicate that the Greenland ice sheet will continue to lose mass until 2100, with contributions to sea-level rise of 90 ± 50 mm and 32 ± 17 mm for the high (RCP8.5) and low (RCP2.6) scenario, respectively.
Heather J. Andres and Lev Tarasov
Clim. Past, 15, 1621–1646,Short summary
Abrupt climate shifts of large magnitudes were common during glacial states, with explanations centred on the oceans. However, winds drive ocean surface currents so shifts in mean wind conditions could also have played a critical role. In a small ensemble of transient deglacial simulations, we find abrupt shifts in both jet stream location and variability over the North Atlantic. We show that the eastern North American ice sheet margin strongly constrains regional jet characteristics.
Laurie Menviel, Emilie Capron, Aline Govin, Andrea Dutton, Lev Tarasov, Ayako Abe-Ouchi, Russell N. Drysdale, Philip L. Gibbard, Lauren Gregoire, Feng He, Ruza F. Ivanovic, Masa Kageyama, Kenji Kawamura, Amaelle Landais, Bette L. Otto-Bliesner, Ikumi Oyabu, Polychronis C. Tzedakis, Eric Wolff, and Xu Zhang
Geosci. Model Dev., 12, 3649–3685,Short summary
As part of the Past Global Changes (PAGES) working group on Quaternary Interglacials, we propose a protocol to perform transient simulations of the penultimate deglaciation for the Paleoclimate Modelling Intercomparison Project (PMIP4). This design includes time-varying changes in orbital forcing, greenhouse gas concentrations, continental ice sheets as well as freshwater input from the disintegration of continental ice sheets. Key paleo-records for model-data comparison are also included.
Taimaz Bahadory and Lev Tarasov
Geosci. Model Dev., 11, 3883–3902,Short summary
We describe a new coupling between the Glacial Systems Model and the LOVECLIM intermediate complexity climate model. The coupling is distinguished from that of previous studies by greater completeness and accuracy, with the intent of capturing the major feedbacks between ice sheets and climate on glacial cycle timescales. The fully coupled model will be used to examine the ice/climate phase space of past glacial cycles.
Laurie Menviel, Emilie Capron, Aline Govin, Andrea Dutton, Lev Tarasov, Ayako Abe-Ouchi, Russell Drysdale, Philip Gibbard, Lauren Gregoire, Feng He, Ruza Ivanovic, Masa Kageyama, Kenji Kawamura, Amaelle Landais, Bette L. Otto-Bliesner, Ikumi Oyabu, Polychronis Tzedakis, Eric Wolff, and Xu Zhang
Clim. Past Discuss.,
Preprint withdrawnShort summary
The penultimate deglaciation (~ 138–128 ka), which represents the transition into the Last Interglacial period, provides a framework to investigate the climate and environmental response to large changes in boundary conditions. Here, as part of the PAGES-PMIP working group on Quaternary Interglacials, we propose a protocol to perform transient simulations of the penultimate deglaciation as well as a selection of paleo records for upcoming model-data comparisons.
Masa Kageyama, Pascale Braconnot, Sandy P. Harrison, Alan M. Haywood, Johann H. Jungclaus, Bette L. Otto-Bliesner, Jean-Yves Peterschmitt, Ayako Abe-Ouchi, Samuel Albani, Patrick J. Bartlein, Chris Brierley, Michel Crucifix, Aisling Dolan, Laura Fernandez-Donado, Hubertus Fischer, Peter O. Hopcroft, Ruza F. Ivanovic, Fabrice Lambert, Daniel J. Lunt, Natalie M. Mahowald, W. Richard Peltier, Steven J. Phipps, Didier M. Roche, Gavin A. Schmidt, Lev Tarasov, Paul J. Valdes, Qiong Zhang, and Tianjun Zhou
Geosci. Model Dev., 11, 1033–1057,Short summary
The Paleoclimate Modelling Intercomparison Project (PMIP) takes advantage of the existence of past climate states radically different from the recent past to test climate models used for climate projections and to better understand these climates. This paper describes the PMIP contribution to CMIP6 (Coupled Model Intercomparison Project, 6th phase) and possible analyses based on PMIP results, as well as on other CMIP6 projects.
Mohammad Hizbul Bahar Arif, Lev Tarasov, and Tristan Hauser
Geosci. Model Dev. Discuss.,
Revised manuscript has not been submittedShort summary
This study is a first step answer to the following question: Can you use emulators (machine learning techniques) to make the output of fast simple climate models (a 2-D energy balance model in this test case) indistinguishable from that of a much more computationally expensive General Circulation climate model (GCM) within the uncertainties of GCMs? Our preliminary test of this concept for large spatio-temporal contexts gives a positive answer.
Masa Kageyama, Samuel Albani, Pascale Braconnot, Sandy P. Harrison, Peter O. Hopcroft, Ruza F. Ivanovic, Fabrice Lambert, Olivier Marti, W. Richard Peltier, Jean-Yves Peterschmitt, Didier M. Roche, Lev Tarasov, Xu Zhang, Esther C. Brady, Alan M. Haywood, Allegra N. LeGrande, Daniel J. Lunt, Natalie M. Mahowald, Uwe Mikolajewicz, Kerim H. Nisancioglu, Bette L. Otto-Bliesner, Hans Renssen, Robert A. Tomas, Qiong Zhang, Ayako Abe-Ouchi, Patrick J. Bartlein, Jian Cao, Qiang Li, Gerrit Lohmann, Rumi Ohgaito, Xiaoxu Shi, Evgeny Volodin, Kohei Yoshida, Xiao Zhang, and Weipeng Zheng
Geosci. Model Dev., 10, 4035–4055,Short summary
The Last Glacial Maximum (LGM, 21000 years ago) is an interval when global ice volume was at a maximum, eustatic sea level close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. This paper describes the implementation of the LGM numerical experiment for the PMIP4-CMIP6 modelling intercomparison projects and the associated sensitivity experiments.
Ruza F. Ivanovic, Lauren J. Gregoire, Masa Kageyama, Didier M. Roche, Paul J. Valdes, Andrea Burke, Rosemarie Drummond, W. Richard Peltier, and Lev Tarasov
Geosci. Model Dev., 9, 2563–2587,Short summary
This manuscript presents the experiment design for the PMIP4 Last Deglaciation Core experiment: a transient simulation of the last deglaciation, 21–9 ka. Specified model boundary conditions include time-varying orbital parameters, greenhouse gases, ice sheets, ice meltwater fluxes and other geographical changes (provided for 26–0 ka). The context of the experiment and the choices for the boundary conditions are explained, along with the future direction of the working group.
André Düsterhus, Alessio Rovere, Anders E. Carlson, Benjamin P. Horton, Volker Klemann, Lev Tarasov, Natasha L. M. Barlow, Tom Bradwell, Jorie Clark, Andrea Dutton, W. Roland Gehrels, Fiona D. Hibbert, Marc P. Hijma, Nicole Khan, Robert E. Kopp, Dorit Sivan, and Torbjörn E. Törnqvist
Clim. Past, 12, 911–921,Short summary
This review/position paper addresses problems in creating new interdisciplinary databases for palaeo-climatological sea-level and ice-sheet data and gives an overview on new advances to tackle them. The focus therein is to define and explain strategies and highlight their importance to allow further progress in these fields. It also offers important insights into the general problem of designing competitive databases which are also applicable to other communities within the palaeo-environment.
A. Abe-Ouchi, F. Saito, M. Kageyama, P. Braconnot, S. P. Harrison, K. Lambeck, B. L. Otto-Bliesner, W. R. Peltier, L. Tarasov, J.-Y. Peterschmitt, and K. Takahashi
Geosci. Model Dev., 8, 3621–3637,Short summary
We describe the creation of boundary conditions related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and altitude of ice-free land, at the Last Glacial Maximum (LGM), for use in LGM experiments conducted as part of the Coupled Modelling Intercomparison Project (CMIP5) and Palaeoclimate Modelling Intercomparison Project (PMIP3). The difference in the ice sheet boundary conditions as well as the climate response to them are discussed.
K. Le Morzadec, L. Tarasov, M. Morlighem, and H. Seroussi
Geosci. Model Dev., 8, 3199–3213,Short summary
A long-term challenge for any model of complex large-scale processes is accounting for the impact of unresolved sub-grid (SG) processes. We quantify the impact of SG mass-balance and ice fluxes on glacial cycle ensemble results for North America. We find no easy solutions to accurately capture these impacts. We show that SG process representation and associated parametric uncertainties can have significant impact on coarse resolution model results for glacial cycle ice sheet evolution.
H. Beltrami, G. S. Matharoo, L. Tarasov, V. Rath, and J. E. Smerdon
Clim. Past, 10, 1693–1706,
R. Briggs, D. Pollard, and L. Tarasov
The Cryosphere, 7, 1949–1970,
Related subject area
CryosphereThe Stochastic Ice-Sheet and Sea-Level System Model v1.0 (StISSM v1.0)Improved representation of the contemporary Greenland ice sheet firn layer by IMAU-FDM v1.2GModeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.Benchmarking the vertically integrated ice-sheet model IMAU-ICE (version 2.0)SnowClim v1.0: high-resolution snow model and data for the western United StatesSnow Multidata Mapping and Modeling (S3M) 5.1: a distributed cryospheric model with dry and wet snow, data assimilation, glacier mass balance, and debris-driven meltMPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshesExplicitly modelling microtopography in permafrost landscapes in a land surface model (JULES vn5.4_microtopography)Improving snow albedo modeling in E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan PlateauGeometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)Mapping high-resolution basal topography of West Antarctica from radar data using non-stationary multiple-point geostatistics (MPS-BedMappingV1)NEMO-Bohai 1.0: a high-resolution ocean and sea ice modelling system for the Bohai Sea, ChinaAn improved regional coupled modeling system for Arctic sea ice simulation and prediction: a case study for 2018WIFF1.0: a hybrid machine-learning-based parameterization of wave-induced sea ice floe fractureThe Whole Antarctic Ocean Model (WAOM v1.0): development and evaluationSNICAR-ADv3: a community tool for modeling spectral snow albedoSTEMMUS-UEB v1.0.0: integrated modeling of snowpack and soil water and energy transfer with three complexity levels of soil physical processesA versatile method for computing optimized snow albedo from spectrally fixed radiative variables: VALHALLA v1.0Ice Algae Model Intercomparison Project phase 2 (IAMIP2)A Gaussian process emulator for simulating ice sheet–climate interactions on a multi-million-year timescale: CLISEMv1.0SITool (v1.0) – a new evaluation tool for large-scale sea ice simulations: application to CMIP6 OMIPfenics_ice 1.0: a framework for quantifying initialization uncertainty for time-dependent ice sheet modelsIntroducing CRYOWRF v1.0: Multiscale atmospheric flow simulations with advanced snow cover modellingDevelopment of adjoint-based ocean state estimation for the Amundsen and Bellingshausen seas and ice shelf cavities using MITgcm–ECCO (66j)Sensitivity of Northern Hemisphere climate to ice–ocean interface heat flux parameterizationsicepack: a new glacier flow modeling package in Python, version 1.0Benefits of sea ice initialization for the interannual-to-decadal climate prediction skill in the Arctic in EC-Earth3Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domainPerformance of MAR (v3.11) in simulating the drifting-snow climate and surface mass balance of Adélie Land, East AntarcticaAssessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18The Utrecht Finite Volume Ice-Sheet Model: UFEMISM (version 1.0)PERICLIMv1.0: a model deriving palaeo-air temperatures from thaw depth in past permafrost regionsAssessing the simulated soil hydrothermal regime of the active layer from the Noah-MP land surface model (v1.1) in the permafrost regions of the Qinghai–Tibet PlateauCrocO_v1.0: a particle filter to assimilate snowpack observations in a spatialised frameworkA fully coupled Arctic sea-ice–ocean–atmosphere model (ArcIOAM v1.0) based on C-Coupler2: model description and preliminary resultsThe Framework For Ice Sheet–Ocean Coupling (FISOC) V1.1Comparison of sea ice kinematics at different resolutions modeled with a grid hierarchy in the Community Earth System Model (version 1.2.1)Snow profile alignment and similarity assessment for aggregating, clustering, and evaluating snowpack model output for avalanche forecastingImprovements in one-dimensional grounding-line parameterizations in an ice-sheet model with lateral variations (PSUICE3D v2.1)Implementation of the RCIP scheme and its performance for 1-D age computations in ice-sheet modelsCOSIPY v1.3 – an open-source coupled snowpack and ice surface energy and mass balance modelUsing Arctic ice mass balance buoys for evaluation of modelled ice energy fluxesImpact of the ice thickness distribution discretization on the sea ice concentration variability in the NEMO3.6–LIM3 global ocean–sea ice modelSimulating the Early Holocene demise of the Laurentide Ice Sheet with BISICLES (public trunk revision 3298)Extended enthalpy formulations in the Ice-sheet and Sea-level System Model (ISSM) version 4.17: discontinuous conductivity and anisotropic streamline upwind Petrov–Galerkin (SUPG) methodThe Community Firn Model (CFM) v1.0Description and validation of the ice-sheet model Yelmo (version 1.0)Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra siteSICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenADOn the calculation of normalized viscous–plastic sea ice stresses
Vincent Verjans, Alexander A. Robel, Helene Seroussi, Lizz Ultee, and Andrew F. Thompson
Geosci. Model Dev., 15, 8269–8293,Short summary
We describe the development of the first large-scale ice sheet model that accounts for stochasticity in a range of processes. Stochasticity allows the impacts of inherently uncertain processes on ice sheets to be represented. This includes climatic uncertainty, as the climate is inherently chaotic. Furthermore, stochastic capabilities also encompass poorly constrained glaciological processes that display strong variability at fine spatiotemporal scales. We present the model and test experiments.
Max Brils, Peter Kuipers Munneke, Willem Jan van de Berg, and Michiel van den Broeke
Geosci. Model Dev., 15, 7121–7138,Short summary
Firn covers the Greenland ice sheet (GrIS) and can temporarily prevent mass loss. Here, we present the latest version of our firn model, IMAU-FDM, with an application to the GrIS. We improved the density of fallen snow, the firn densification rate and the firn's thermal conductivity. This leads to a higher air content and 10 m temperatures. Furthermore we investigate three case studies and find that the updated model shows greater variability and an increased sensitivity in surface elevation.
Océane Hames, Mahdi Jafari, David Nicholas Wagner, Ian Raphael, David Clemens-Sewall, Chris Polashenski, Matthew D. Shupe, Martin Schneebeli, and Michael Lehning
Geosci. Model Dev., 15, 6429–6449,Short summary
This paper presents an Eulerian–Lagrangian snow transport model implemented in the fluid dynamics software OpenFOAM, which we call snowBedFoam 1.0. We apply this model to reproduce snow deposition on a piece of ridged Arctic sea ice, which was produced during the MOSAiC expedition through scan measurements. The model appears to successfully reproduce the enhanced snow accumulation and deposition patterns, although some quantitative uncertainties were shown.
Constantijn J. Berends, Heiko Goelzer, Thomas J. Reerink, Lennert B. Stap, and Roderik S. W. van de Wal
Geosci. Model Dev., 15, 5667–5688,Short summary
The rate at which marine ice sheets such as the West Antarctic ice sheet will retreat in a warming climate and ocean is still uncertain. Numerical ice-sheet models, which solve the physical equations that describe the way glaciers and ice sheets deform and flow, have been substantially improved in recent years. Here we present the results of several years of work on IMAU-ICE, an ice-sheet model of intermediate complexity, which can be used to study ice sheets of both the past and the future.
Abby C. Lute, John Abatzoglou, and Timothy Link
Geosci. Model Dev., 15, 5045–5071,Short summary
We developed a snow model that can be used to quantify snowpack over large areas with a high degree of spatial detail. We ran the model over the western United States, creating a snow and climate dataset for three time periods. Compared to observations of snowpack, the model captured the key aspects of snow across time and space. The model and dataset will be useful in understanding historical and future changes in snowpack, with relevance to water resources, agriculture, and ecosystems.
Francesco Avanzi, Simone Gabellani, Fabio Delogu, Francesco Silvestro, Edoardo Cremonese, Umberto Morra di Cella, Sara Ratto, and Hervé Stevenin
Geosci. Model Dev., 15, 4853–4879,Short summary
Knowing in real time how much snow and glacier ice has accumulated across the landscape has significant implications for water-resource management and flood control. This paper presents a computer model – S3M – allowing scientists and decision makers to predict snow and ice accumulation during winter and the subsequent melt during spring and summer. S3M has been employed for real-world flood forecasting since the early 2000s but is here being made open source for the first time.
Adrian K. Turner, William H. Lipscomb, Elizabeth C. Hunke, Douglas W. Jacobsen, Nicole Jeffery, Darren Engwirda, Todd D. Ringler, and Jonathan D. Wolfe
Geosci. Model Dev., 15, 3721–3751,Short summary
We present the dynamical core of the MPAS-Seaice model, which uses a mesh consisting of a Voronoi tessellation with polygonal cells. Such a mesh allows variable mesh resolution in different parts of the domain and the focusing of computational resources in regions of interest. We describe the velocity solver and tracer transport schemes used and examine errors generated by the model in both idealized and realistic test cases and examine the computational efficiency of the model.
Noah D. Smith, Eleanor J. Burke, Kjetil Schanke Aas, Inge H. J. Althuizen, Julia Boike, Casper Tai Christiansen, Bernd Etzelmüller, Thomas Friborg, Hanna Lee, Heather Rumbold, Rachael H. Turton, Sebastian Westermann, and Sarah E. Chadburn
Geosci. Model Dev., 15, 3603–3639,Short summary
The Arctic has large areas of small mounds that are caused by ice lifting up the soil. Snow blown by wind gathers in hollows next to these mounds, insulating them in winter. The hollows tend to be wetter, and thus the soil absorbs more heat in summer. The warm wet soil in the hollows decomposes, releasing methane. We have made a model of this, and we have tested how it behaves and whether it looks like sites in Scandinavia and Siberia. Sometimes we get more methane than a model without mounds.
Dalei Hao, Gautam Bisht, Cenlin He, Edward Bair, Huilin Huang, Cheng Dang, Karl Rittger, Yu Gu, Hailong Wang, Yun Qian, and L. Ruby Leung
Geosci. Model Dev. Discuss.,Short summary
Snow with the highest albedo of land surface plays a vital role in Earth’s surface energy and water cycles. Our study demonstrates that snow grain shape and mixing state of light absorbing particle-snow play an important role in snow processes, surface energy balance and water cycle. This study offers guidance for improving snow simulations and surface radiative forcing estimates in Earth system models under future climate.
Adrian K. Turner, Kara J. Peterson, and Dan Bolintineanu
Geosci. Model Dev., 15, 1953–1970,Short summary
We developed a technique to remap sea ice tracer quantities between circular discrete element distributions. This is needed for a global discrete element method sea ice model being developed jointly by Los Alamos National Laboratory and Sandia National Laboratories that has the potential to better utilize newer supercomputers with graphics processing units and better represent sea ice dynamics. This new remapping technique ameliorates the effect of element distortion created by sea ice ridging.
Zhen Yin, Chen Zuo, Emma J. MacKie, and Jef Caers
Geosci. Model Dev., 15, 1477–1497,Short summary
We provide a multiple-point geostatistics approach to probabilistically learn from training images to fill large-scale irregular geophysical data gaps. With a repository of global topographic training images, our approach models high-resolution basal topography and quantifies the geospatial uncertainty. It generated high-resolution topographic realizations to investigate the impact of basal topographic uncertainty on critical subglacial hydrological flow patterns associated with ice velocity.
Yu Yan, Wei Gu, Andrea M. U. Gierisch, Yingjun Xu, and Petteri Uotila
Geosci. Model Dev., 15, 1269–1288,Short summary
In this study, we developed NEMO-Bohai, an ocean–ice model for the Bohai Sea, China. This study presented the scientific design and technical choices of the parameterizations for the NEMO-Bohai model. The model was calibrated and evaluated with in situ and satellite observations of ocean and sea ice. NEMO-Bohai is intended to be a valuable tool for long-term ocean and ice simulations and climate change studies.
Chao-Yuan Yang, Jiping Liu, and Dake Chen
Geosci. Model Dev., 15, 1155–1176,Short summary
We present an improved coupled modeling system for Arctic sea ice prediction. We perform Arctic sea ice prediction experiments with improved/updated physical parameterizations, which show better skill in predicting sea ice state as well as atmospheric and oceanic state in the Arctic compared with its predecessor. The improved model also shows extended predictive skill of Arctic sea ice after the summer season. This provides an added value of this prediction system for decision-making.
Christopher Horvat and Lettie A. Roach
Geosci. Model Dev., 15, 803–814,Short summary
Sea ice is a composite of individual pieces, called floes, ranging in horizontal size from meters to kilometers. Variations in sea ice geometry are often forced by ocean waves, a process that is an important target of global climate models as it affects the rate of sea ice melting. Yet directly simulating these interactions is computationally expensive. We present a neural-network-based model of wave–ice fracture that allows models to incorporate their effect without added computational cost.
Ole Richter, David E. Gwyther, Benjamin K. Galton-Fenzi, and Kaitlin A. Naughten
Geosci. Model Dev., 15, 617–647,Short summary
Here we present an improved model of the Antarctic continental shelf ocean and demonstrate that it is capable of reproducing present-day conditions. The improvements are fundamental and regard the inclusion of tides and ocean eddies. We conclude that the model is well suited to gain new insights into processes that are important for Antarctic ice sheet retreat and global ocean changes. Hence, the model will ultimately help to improve projections of sea level rise and climate change.
Mark G. Flanner, Julian B. Arnheim, Joseph M. Cook, Cheng Dang, Cenlin He, Xianglei Huang, Deepak Singh, S. McKenzie Skiles, Chloe A. Whicker, and Charles S. Zender
Geosci. Model Dev., 14, 7673–7704,Short summary
We present the technical formulation and evaluation of a publicly available code and web-based model to simulate the spectral albedo of snow. Our model accounts for numerous features of the snow state and ambient conditions, including the the presence of light-absorbing matter like black and brown carbon, mineral dust, volcanic ash, and snow algae. Carbon dioxide snow, found on Mars, is also represented. The model accurately reproduces spectral measurements of clean and contaminated snow.
Lianyu Yu, Yijian Zeng, and Zhongbo Su
Geosci. Model Dev., 14, 7345–7376,Short summary
We developed an integrated soil–snow–atmosphere model (STEMMUS-UEB) dedicated to the physical description of snow and soil processes with various complexities. With STEMMUS-UEB, we demonstrated that the snowpack affects not only the soil surface moisture conditions (in the liquid and ice phase) and energy-related states (albedo, LE) but also the subsurface soil water and vapor transfer, which contributes to a better understanding of the hydrothermal implications of the snowpack in cold regions.
Florent Veillon, Marie Dumont, Charles Amory, and Mathieu Fructus
Geosci. Model Dev., 14, 7329–7343,Short summary
In climate models, the snow albedo scheme generally calculates only a narrowband or broadband albedo. Therefore, we have developed the VALHALLA method to optimize snow spectral albedo calculations through the determination of spectrally fixed radiative variables. The development of VALHALLA v1.0 with the use of the snow albedo model TARTES and the spectral irradiance model SBDART indicates a considerable reduction in calculation time while maintaining an adequate accuracy of albedo values.
Hakase Hayashida, Meibing Jin, Nadja S. Steiner, Neil C. Swart, Eiji Watanabe, Russell Fiedler, Andrew McC. Hogg, Andrew E. Kiss, Richard J. Matear, and Peter G. Strutton
Geosci. Model Dev., 14, 6847–6861,Short summary
Ice algae are tiny plants like phytoplankton but they grow within sea ice. In polar regions, both phytoplankton and ice algae are the foundation of marine ecosystems and play an important role in taking up carbon dioxide in the atmosphere. However, state-of-the-art climate models typically do not include ice algae, and therefore their role in the climate system remains unclear. This project aims to address this knowledge gap by coordinating a set of experiments using sea-ice–ocean models.
Jonas Van Breedam, Philippe Huybrechts, and Michel Crucifix
Geosci. Model Dev., 14, 6373–6401,Short summary
Ice sheets are an important component of the climate system and interact with the atmosphere through albedo variations and changes in the surface height. On very long timescales, it is impossible to directly couple ice sheet models with climate models and other techniques have to be used. Here we present a novel coupling method between ice sheets and the atmosphere by making use of an emulator to simulate ice sheet–climate interactions for several million years.
Xia Lin, François Massonnet, Thierry Fichefet, and Martin Vancoppenolle
Geosci. Model Dev., 14, 6331–6354,Short summary
This study introduces a new Sea Ice Evaluation Tool (SITool) to evaluate the model skills on the bipolar sea ice simulations by providing performance metrics and diagnostics. SITool is applied to evaluate the CMIP6 OMIP simulations. By changing the atmospheric forcing from CORE-II to JRA55-do data, many aspects of sea ice simulations are improved. SITool will be useful for helping teams managing various versions of a sea ice model or tracking the time evolution of model performance.
Conrad P. Koziol, Joe A. Todd, Daniel N. Goldberg, and James R. Maddison
Geosci. Model Dev., 14, 5843–5861,Short summary
Sea level change due to the loss of ice sheets presents great risk for coastal communities. Models are used to forecast ice loss, but their evolution depends strongly on properties which are hidden from observation and must be inferred from satellite observations. Common methods for doing so do not allow for quantification of the uncertainty inherent or how it will affect forecasts. We provide a framework for quantifying how this
initialization uncertaintyaffects ice loss forecasts.
Varun Sharma, Franziska Gerber, and Michael Lehning
Geosci. Model Dev. Discuss.,
Revised manuscript accepted for GMDShort summary
Most current generation climate and weather models have a relatively simplistic description of snow and snow-atmosphere interaction. One reason for this is the belief that including an advanced snow model would make the simulations too computationally demanding. In this study, we bring together two state-of-the-art models for atmosphere (WRF) and snow-cover (SNOWPACK) and highlight both the feasibility and necessity of such coupled models to explore under-explored phenomena in the cryosphere.
Yoshihiro Nakayama, Dimitris Menemenlis, Ou Wang, Hong Zhang, Ian Fenty, and An T. Nguyen
Geosci. Model Dev., 14, 4909–4924,Short summary
High ice shelf melting in the Amundsen Sea has attracted many observational campaigns in the past decade. One method to combine observations with numerical models is the adjoint method. After 20 iterations, the cost function, defined as a sum of the weighted model–data difference, is reduced by 65 % by adjusting initial conditions, atmospheric forcing, and vertical diffusivity. This study demonstrates adjoint-method optimization with explicit representation of ice shelf cavity circulation.
Xiaoxu Shi, Dirk Notz, Jiping Liu, Hu Yang, and Gerrit Lohmann
Geosci. Model Dev., 14, 4891–4908,Short summary
The ice–ocean heat flux is one of the key elements controlling sea ice changes. It motivates our study, which aims to examine the responses of modeled climate to three ice–ocean heat flux parameterizations, including two old approaches that assume one-way heat transport and a new one describing a double-diffusive ice–ocean heat exchange. The results show pronounced differences in the modeled sea ice, ocean, and atmosphere states for the latter as compared to the former two parameterizations.
Daniel R. Shapero, Jessica A. Badgeley, Andrew O. Hoffman, and Ian R. Joughin
Geosci. Model Dev., 14, 4593–4616,Short summary
This paper describes a new software package called "icepack" for modeling the flow of ice sheets and glaciers. Glaciologists use tools like icepack to better understand how ice sheets flow, what role they have played in shaping Earth's climate, and how much sea level rise we can expect in the coming decades to centuries. The icepack package includes several innovations to help researchers describe and solve interesting glaciological problems and to experiment with the underlying model physics.
Tian Tian, Shuting Yang, Mehdi Pasha Karami, François Massonnet, Tim Kruschke, and Torben Koenigk
Geosci. Model Dev., 14, 4283–4305,Short summary
Three decadal prediction experiments with EC-Earth3 are performed to investigate the impact of ocean, sea ice concentration and thickness initialization, respectively. We find that the persistence of perennial thick ice in the central Arctic can affect the sea ice predictability in its adjacent waters via advection process or wind, despite those regions being seasonally ice free during two recent decades. This has implications for the coming decades as the thinning of Arctic sea ice continues.
Moritz Kreuzer, Ronja Reese, Willem Nicholas Huiskamp, Stefan Petri, Torsten Albrecht, Georg Feulner, and Ricarda Winkelmann
Geosci. Model Dev., 14, 3697–3714,Short summary
We present the technical implementation of a coarse-resolution coupling between an ice sheet model and an ocean model that allows one to simulate ice–ocean interactions at timescales from centuries to millennia. As ice shelf cavities cannot be resolved in the ocean model at coarse resolution, we bridge the gap using an sub-shelf cavity module. It is shown that the framework is computationally efficient, conserves mass and energy, and can produce a stable coupled state under present-day forcing.
Charles Amory, Christoph Kittel, Louis Le Toumelin, Cécile Agosta, Alison Delhasse, Vincent Favier, and Xavier Fettweis
Geosci. Model Dev., 14, 3487–3510,Short summary
This paper presents recent developments in the drifting-snow scheme of the regional climate model MAR and its application to simulate drifting snow and the surface mass balance of Adélie Land in East Antarctica. The model is extensively described and evaluated against a multi-year drifting-snow dataset and surface mass balance estimates available in the area. The model sensitivity to input parameters and improvements over a previously published version are also assessed.
Thiago Dias dos Santos, Mathieu Morlighem, and Hélène Seroussi
Geosci. Model Dev., 14, 2545–2573,Short summary
Numerical models are routinely used to understand the past and future behavior of ice sheets in response to climate evolution. As is always the case with numerical modeling, one needs to minimize biases and numerical artifacts due to the choice of numerical scheme employed in such models. Here, we assess different numerical schemes in time-dependent simulations of ice sheets. We also introduce a new parameterization for the driving stress, the force that drives the ice sheet flow.
Constantijn J. Berends, Heiko Goelzer, and Roderik S. W. van de Wal
Geosci. Model Dev., 14, 2443–2470,Short summary
The largest uncertainty in projections of sea-level rise comes from ice-sheet retreat. To better understand how these ice sheets respond to the changing climate, ice-sheet models are used, which must be able to reproduce both their present and past evolution. We have created a model that is fast enough to simulate an ice sheet at a high resolution over the course of an entire 120 000-year glacial cycle. This allows us to study processes that cannot be captured by lower-resolution models.
Tomáš Uxa, Marek Křížek, and Filip Hrbáček
Geosci. Model Dev., 14, 1865–1884,Short summary
We present a simple model that derives palaeo-air temperature characteristics related to the palaeo-active-layer thickness, which can be recognized using many relict periglacial features found in past permafrost regions. Its evaluation against modern temperature records and an experimental palaeo-air temperature reconstruction showed relatively high model accuracy, which suggests that it could become a useful tool for reconstructing Quaternary palaeo-environments.
Xiangfei Li, Tonghua Wu, Xiaodong Wu, Jie Chen, Xiaofan Zhu, Guojie Hu, Ren Li, Yongping Qiao, Cheng Yang, Junming Hao, Jie Ni, and Wensi Ma
Geosci. Model Dev., 14, 1753–1771,Short summary
In this study, an ensemble simulation of 55296 scheme combinations for at a typical permafrost site on the Qinghai–Tibet Plateau (QTP) was conducted. The general performance of the Noah-MP model for snow cover events (SCEs), soil temperature (ST) and soil liquid water content (SLW) was assessed, and the sensitivities of parameterization schemes at different depths were investigated. We show that Noah-MP tends to overestimate SCEs and underestimate ST and topsoil SLW on the QTP.
Bertrand Cluzet, Matthieu Lafaysse, Emmanuel Cosme, Clément Albergel, Louis-François Meunier, and Marie Dumont
Geosci. Model Dev., 14, 1595–1614,Short summary
In the mountains, the combination of large model error and observation sparseness is a challenge for data assimilation. Here, we develop two variants of the particle filter (PF) in order to propagate the information content of observations into unobserved areas. By adjusting observation errors or exploiting background correlation patterns, we demonstrate the potential for partial observations of snow depth and surface reflectance to improve model accuracy with the PF in an idealised setting.
Shihe Ren, Xi Liang, Qizhen Sun, Hao Yu, L. Bruno Tremblay, Bo Lin, Xiaoping Mai, Fu Zhao, Ming Li, Na Liu, Zhikun Chen, and Yunfei Zhang
Geosci. Model Dev., 14, 1101–1124,Short summary
Sea ice plays a crucial role in global energy and water budgets. To get a better simulation of sea ice, we coupled a sea ice model with an atmospheric and ocean model to form a fully coupled system. The sea ice simulation results of this coupled system demonstrated that a two-way coupled model has better performance in terms of sea ice, especially in summer. This indicates that sea-ice–ocean–atmosphere interaction plays a crucial role in controlling Arctic summertime sea ice distribution.
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., 14, 889–905,Short 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.
Shiming Xu, Jialiang Ma, Lu Zhou, Yan Zhang, Jiping Liu, and Bin Wang
Geosci. Model Dev., 14, 603–628,Short summary
A multi-resolution tripolar grid hierarchy is constructed and integrated in CESM (version 1.2.1). The resolution range includes 0.45, 0.15, and 0.05°. Based on atmospherically forced sea ice experiments, the model simulates reasonable sea ice kinematics and scaling properties. Landfast ice thickness can also be systematically shifted due to non-convergent solutions to an elastic–viscous–plastic (EVP) model. This work is a framework for multi-scale modeling of the ocean and sea ice with CESM.
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.
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.
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.
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. a
Bartholomew, I. D., Nienow, P., Sole, A., Mair, D., Cowton, T., King, M. A., and Palmer, S.: Seasonal variations in Greenland Ice Sheet motion: Inland extent and behaviour at higher elevations, Earth Planet. Sc. Lett., 307, 271–278, 2011. a
Clark, P. U. and Walder, J. S.: Subglacial drainage, eskers, and deforming beds beneath the Laurentide and Eurasian ice sheets, Geol. Soc. Am. Bull., 106, 304–314, 1994. a
de Fleurian, B., Morlighem, M., Seroussi, H., Rignot, E., van den Broeke, M. R., Munneke, P. K., Mouginot, J., Smeets, P. C. J. P., and Tedstone, A. J.: A modeling study of the effect of variability on the effective pressure beneath Russell Glacier, West Greenland, J. Geophys. Res.-Earth, 121, 1834–1848, https://doi.org/10.1002/2016JF003842, 2016. a
Dow, C. F., Kulessa, B., Rutt, I. C., Tsai, V. C., Pimentel, S., Doyle, S. H., van As, D., and Lindbäck, K.: Modeling of subglacial hydrological development following rapid supraglacial lake drainage, J. Geophys. Res.-Earth, 120, 1127–1147, 2015. a
Hoffman, M. J. and Price, S.: Feedbacks between coupled subglacial hydrology and glacier dynamics, J. Geophys. Res.-Earth, 119, 414–436, 2014. a
Kamb, B., Raymond, C. F., Harrison, W. D., engelhardt, H., Echelmeyer, K. A., Humphrey, N., Brugman, M. M., and Pfeffer, T.: Glacier Surge Mechanism: 1982–1983 Surge of Variegated Glacier, Alaska, Science, 227, 469–479, 1985. a
Kavanagh, M.: A fast, physically-based subglacial hydrology model applied to the North American ice complex over the last glacial cycle, Master's thesis, Memorial University of Newfoundland, Canada, 2012. a
Mathews, J. H. and Fink, K. K.: Numerical Methods Using Matlab, chap. 9, 4th edn., Prentice-Hall Inc., 2004. a
Person, M., Bense, V., Cohen, D., and Banerjee, A.: Models of ice-sheet hydrogeologic interactions: a review, Geofluids, 12, 58–78, 2012. a
Röthlisberger, H.: Water Pressure in Intra- and Subglacial Channels, J. Glaciology, 11, 177–203, 1972. a
Shchepetkin, A. F. and McWilliams, J. C.: The regional oceanic modeling system (ROMS):a split-explicit, free-surface, topography-following-coordinate oceanic model, Ocean Model., 9, 347–404, 2005. a
Tarasov, L. and Peltier, W. R.: Greenland glacial history and local geodynamic consequences, Geophys. J. Int., 150, 198–229, 2002. a
Tarasov, L. and Peltier, W. R.: A Geophysically constrained large ensemble analysis of the deglacial history of the North Americal ice-sheet complex, Quaternary Sci. Rev., 23, 359–388, 2004. a
Walder, J. S. and Hallet, B.: Geometry of former subglacial water channels and cavities, J. Glaciol., 23, 335–346, 1979. a
Weertman, J.: General theory of water flow at the base of a glacier or ice sheet, Rev. Geophys., 10, 287–333, 1972. a
Werder, M. A., Hewitt, I. J., Schoof, C. G., and Flowers, G. E.: Modeling channelized and distributed subglacial drainage in two dimensions, J. Geophys. Res.-Earth, 118, 2140–2158, 2013. a
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.
We present and validate BrAHMs (BAsal Hydrology Model): a new physically based basal hydrology...