Articles | Volume 17, issue 17
https://doi.org/10.5194/gmd-17-6529-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-17-6529-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance
Danish Meteorological Institute, Sankt Kjelds Plads 11, 2100 Copenhagen, Denmark
Jacob Poulsen
Intel Corporation, 2200 Mission College Blvd., Santa Clara, CA, 95054-1549, USA
Mads Hvid Ribergaard
Danish Meteorological Institute, Sankt Kjelds Plads 11, 2100 Copenhagen, Denmark
Ruchira Sasanka
Intel Corporation, 2200 Mission College Blvd., Santa Clara, CA, 95054-1549, USA
Anthony P. Craig
Science and Technology Corporation, Seattle, WA, 98105, USA
Elizabeth C. Hunke
MS-B216, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
Stefan Rethmeier
Danish Meteorological Institute, Sankt Kjelds Plads 11, 2100 Copenhagen, Denmark
Related authors
Imke Sievers, Henriette Skourup, and Till A. S. Rasmussen
The Cryosphere, 18, 5985–6004, https://doi.org/10.5194/tc-18-5985-2024, https://doi.org/10.5194/tc-18-5985-2024, 2024
Short summary
Short summary
To derive sea ice thickness (SIT) from satellite freeboard (FB) observations, assumptions about snow thickness, snow density, sea ice density and water density are needed. These parameters are impossible to observe alongside FB, so many existing products use empirical values. In this study, modeled values are used instead. The modeled values and otherwise commonly used empirical values are evaluated against in situ observations. In a further analysis, the influence on SIT is quantified.
Jean-François Lemieux, William H. Lipscomb, Anthony Craig, David A. Bailey, Elizabeth C. Hunke, Philippe Blain, Till A. S. Rasmussen, Mats Bentsen, Frédéric Dupont, David Hebert, and Richard Allard
Geosci. Model Dev., 17, 6703–6724, https://doi.org/10.5194/gmd-17-6703-2024, https://doi.org/10.5194/gmd-17-6703-2024, 2024
Short summary
Short summary
We present the latest version of the CICE model. It solves equations that describe the dynamics and the growth and melt of sea ice. To do so, the domain is divided into grid cells and variables are positioned at specific locations in the cells. A new implementation (C-grid) is presented, with the velocity located on cell edges. Compared to the previous B-grid, the C-grid allows for a natural coupling with some oceanic and atmospheric models. It also allows for ice transport in narrow channels.
Imke Sievers, Till A. S. Rasmussen, and Lars Stenseng
The Cryosphere, 17, 3721–3738, https://doi.org/10.5194/tc-17-3721-2023, https://doi.org/10.5194/tc-17-3721-2023, 2023
Short summary
Short summary
The satellite CryoSat-2 measures freeboard (FB), which is used to derive sea ice thickness (SIT) under the assumption of hydrostatic balance. This SIT comes with large uncertainties due to errors in the observed FB, sea ice density, snow density and snow thickness. This study presents a new method to derive SIT by assimilating the FB into the sea ice model, evaluates the resulting SIT against in situ observations and compares the results to the CryoSat-2-derived SIT without FB assimilation.
Imke Sievers, Andrea M. U. Gierisch, Till A. S. Rasmussen, Robinson Hordoir, and Lars Stenseng
The Cryosphere Discuss., https://doi.org/10.5194/tc-2022-84, https://doi.org/10.5194/tc-2022-84, 2022
Preprint withdrawn
Short summary
Short summary
To predict Arctic sea ice models are used. Many ice models exists. They all are skill full, but give different results. Often this differences result from forcing as for example air temperature. Other differences result from the way the physical equations are solved in the model. In this study two commonly used models are compared under equal forcing, to find out how much the models differ under similar external forcing. The results are compared to observations and to eachother.
Imke Sievers, Henriette Skourup, and Till A. S. Rasmussen
The Cryosphere, 18, 5985–6004, https://doi.org/10.5194/tc-18-5985-2024, https://doi.org/10.5194/tc-18-5985-2024, 2024
Short summary
Short summary
To derive sea ice thickness (SIT) from satellite freeboard (FB) observations, assumptions about snow thickness, snow density, sea ice density and water density are needed. These parameters are impossible to observe alongside FB, so many existing products use empirical values. In this study, modeled values are used instead. The modeled values and otherwise commonly used empirical values are evaluated against in situ observations. In a further analysis, the influence on SIT is quantified.
Kavya Sivaraj, Kurt Solander, Charles Abolt, and Elizabeth Hunke
EGUsphere, https://doi.org/10.5194/egusphere-2024-3315, https://doi.org/10.5194/egusphere-2024-3315, 2024
Short summary
Short summary
Melt ponds are seasonal water bodies whose presence affect the rate of Arctic sea ice loss by increasing the absorption of solar radiation. Despite their importance, large-scale observational datasets of Melt Pond Fraction (MPF) are inadequate due to low-resolution sensors and spectral misclassifications caused by different ice types. Our novel ML-based workflow overcomes these limitations by leveraging morphological operators, resulting in an improved Sentinel-2-based mean MPF of 11% from 20%.
Jean-François Lemieux, William H. Lipscomb, Anthony Craig, David A. Bailey, Elizabeth C. Hunke, Philippe Blain, Till A. S. Rasmussen, Mats Bentsen, Frédéric Dupont, David Hebert, and Richard Allard
Geosci. Model Dev., 17, 6703–6724, https://doi.org/10.5194/gmd-17-6703-2024, https://doi.org/10.5194/gmd-17-6703-2024, 2024
Short summary
Short summary
We present the latest version of the CICE model. It solves equations that describe the dynamics and the growth and melt of sea ice. To do so, the domain is divided into grid cells and variables are positioned at specific locations in the cells. A new implementation (C-grid) is presented, with the velocity located on cell edges. Compared to the previous B-grid, the C-grid allows for a natural coupling with some oceanic and atmospheric models. It also allows for ice transport in narrow channels.
Imke Sievers, Till A. S. Rasmussen, and Lars Stenseng
The Cryosphere, 17, 3721–3738, https://doi.org/10.5194/tc-17-3721-2023, https://doi.org/10.5194/tc-17-3721-2023, 2023
Short summary
Short summary
The satellite CryoSat-2 measures freeboard (FB), which is used to derive sea ice thickness (SIT) under the assumption of hydrostatic balance. This SIT comes with large uncertainties due to errors in the observed FB, sea ice density, snow density and snow thickness. This study presents a new method to derive SIT by assimilating the FB into the sea ice model, evaluates the resulting SIT against in situ observations and compares the results to the CryoSat-2-derived SIT without FB assimilation.
Eva Friis Møller, Asbjørn Christensen, Janus Larsen, Kenneth D. Mankoff, Mads Hvid Ribergaard, Mikael Sejr, Philip Wallhead, and Marie Maar
Ocean Sci., 19, 403–420, https://doi.org/10.5194/os-19-403-2023, https://doi.org/10.5194/os-19-403-2023, 2023
Short summary
Short summary
Melt from the Greenland ice sheet and sea ice both influence light and nutrient availability in the Arctic coastal ocean. We use a 3D coupled hydrodynamic–biogeochemical model to evaluate the relative importance of these processes for timing, distribution, and magnitude of phytoplankton production in Disko Bay, west Greenland. Our study indicates that decreasing sea ice and more freshwater discharge can work synergistically and increase primary productivity of the coastal ocean around Greenland.
Younjoo J. Lee, Wieslaw Maslowski, John J. Cassano, Jaclyn Clement Kinney, Anthony P. Craig, Samy Kamal, Robert Osinski, Mark W. Seefeldt, Julienne Stroeve, and Hailong Wang
The Cryosphere, 17, 233–253, https://doi.org/10.5194/tc-17-233-2023, https://doi.org/10.5194/tc-17-233-2023, 2023
Short summary
Short summary
During 1979–2020, four winter polynyas occurred in December 1986 and February 2011, 2017, and 2018 north of Greenland. Instead of ice melting due to the anomalous warm air intrusion, the extreme wind forcing resulted in greater ice transport offshore. Based on the two ensemble runs, representing a 1980s thicker ice vs. a 2010s thinner ice, a dominant cause of these winter polynyas stems from internal variability of atmospheric forcing rather than from the forced response to a warming climate.
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, https://doi.org/10.5194/gmd-15-3721-2022, https://doi.org/10.5194/gmd-15-3721-2022, 2022
Short summary
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.
Imke Sievers, Andrea M. U. Gierisch, Till A. S. Rasmussen, Robinson Hordoir, and Lars Stenseng
The Cryosphere Discuss., https://doi.org/10.5194/tc-2022-84, https://doi.org/10.5194/tc-2022-84, 2022
Preprint withdrawn
Short summary
Short summary
To predict Arctic sea ice models are used. Many ice models exists. They all are skill full, but give different results. Often this differences result from forcing as for example air temperature. Other differences result from the way the physical equations are solved in the model. In this study two commonly used models are compared under equal forcing, to find out how much the models differ under similar external forcing. The results are compared to observations and to eachother.
Milena Veneziani, Wieslaw Maslowski, Younjoo J. Lee, Gennaro D'Angelo, Robert Osinski, Mark R. Petersen, Wilbert Weijer, Anthony P. Craig, John D. Wolfe, Darin Comeau, and Adrian K. Turner
Geosci. Model Dev., 15, 3133–3160, https://doi.org/10.5194/gmd-15-3133-2022, https://doi.org/10.5194/gmd-15-3133-2022, 2022
Short summary
Short summary
We present an Earth system model (ESM) simulation, E3SM-Arctic-OSI, with a refined grid to better resolve the Arctic ocean and sea-ice system and low spatial resolution elsewhere. The configuration satisfactorily represents many aspects of the Arctic system and its interactions with the sub-Arctic, while keeping computational costs at a fraction of those necessary for global high-resolution ESMs. E3SM-Arctic can thus be an efficient tool to study Arctic processes on climate-relevant timescales.
Cited articles
Amdahl, G. M.: Validity of the Single Processor Approach to Achieving Large Scale Computing Capabilities, in: Proceedings of the April 18–20, 1967, Spring Joint Computer Conference, AFIPS '67 (Spring), p. 483–485, Association for Computing Machinery, New York, NY, USA, ISBN 9781450378956, https://doi.org/10.1145/1465482.1465560, 1967. a
Arakawa, A. and Lamb, V. R.: Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model, in: General Circulation Models of the Atmosphere, edited by: Chang, J., vol. 17 of Methods in Computational Physics: Advances in Research and Applications, pp. 173–265, Elsevier, https://doi.org/10.1016/B978-0-12-460817-7.50009-4, 1977. a
Bouchat, A., Hutter, N., Chanut, J., Dupont, F., Dukhovskoy, D., Garric, G., Lee, Y. J., Lemieux, J.-F., Lique, C., Losch, M., Maslowski, W., Myers, P. G., Ólason, E., Rampal, P., Rasmussen, T., Talandier, C., Tremblay, B., and Wang, Q.: Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields, J. Geophys. Res.-Oceans, 127, e2021JC017667, https://doi.org/10.1029/2021JC017667, 2022. a
Bouillon, S., Fichefet, T., Legat, V., and Madec, G.: The elastic–viscous–plastic method revisited, Ocean Model., 71, 2–12, https://doi.org/10.1016/j.ocemod.2013.05.013, 2013. a
Brunke, M. A., Cassano, J. J., Dawson, N., DuVivier, A. K., Gutowski Jr., W. J., Hamman, J., Maslowski, W., Nijssen, B., Reeves Eyre, J. E. J., Renteria, J. C., Roberts, A., and Zeng, X.: Evaluation of the atmosphere–land–ocean–sea ice interface processes in the Regional Arctic System Model version 1 (RASM1) using local and globally gridded observations, Geosci. Model Dev., 11, 4817–4841, https://doi.org/10.5194/gmd-11-4817-2018, 2018. a
CICE Consortium: CICE gx1 WOA Forcing Data – 2020.03.20 (2020.03.20), Zenodo [data set], https://doi.org/10.5281/zenodo.3728600, 2020. a
CICE Consortium: CICE gx1 Grid and Initial Condition Data – 2021.08.16 (2021.08.16), Zenodo [data set], https://doi.org/10.5281/zenodo.5208241, 2021. a
CICE Consortium: CICE gx1 JRA55do Forcing Data by year – 2023.07.03 (2023.07.03), Zenodo [data set], https://doi.org/10.5281/zenodo.8118062, 2023. a
Craig, A. P., Vertenstein, M., and Jacob, R.: A new flexible coupler for earth system modeling developed for CCSM4 and CESM1, The International Journal of High Performance Computing Applications, 26, 31–42, https://doi.org/10.1177/1094342011428141, 2012. a
Craig, A. P., Mickelson, S. A., Hunke, E. C., and Bailey, D. A.: Improved parallel performance of the CICE model in CESM1, The International Journal of High Performance Computing Applications, 29, 154–165, https://doi.org/10.1177/1094342014548771, 2015. a, b
Hibler, W. D. I.: A Dynamic Thermodynamic Sea Ice Model, J. Phys. Oceanogr., 9, 815–846, https://doi.org/10.1175/1520-0485(1979)009<0815:ADTSIM>2.0.CO;2, 1979. a
Hunke, E., Allard, R., Blain, P., Blockley, E., Feltham, D., Fichefet, T., Garric, G., Grumbine, R., Lemieux, J.-F., Rasmussen, T., Ribergaard, M., Roberts, A., Schweiger, A., Tietsche, S., Tremblay, B., Vancoppenolle, M., and Zhang, J.: Should Sea-Ice Modeling Tools Designed for Climate Research Be Used for Short-Term Forecasting?, Curr. Clim. Change Rep., 6, 121–136, https://doi.org/10.1007/s40641-020-00162-y, 2020. a
Hunke, E., Allard, R., Bailey, D. A., Blain, P., Craig, A., Dupont, F., DuVivier, A., Grumbine, R., Hebert, D., Holland, M., Jeffery, N., Lemieux, J.-F., Osinski, R., Poulsen, J., Steketee, A., Rasmussen, T., Ribergaard, M., Roach, L., Roberts, A., Turner, M., Winton, M., and Worthen, D.: CICE-Consortium/CICE: CICE Version 6.5.1, Zenodo [code], https://doi.org/10.5281/zenodo.11223920, 2024. a, b, c
Hunke, E. C. and Dukowicz, J. K.: An elastic-viscous-plastic model for sea ice dynamics, J. Phys. Oceanogr., 27, 1849–1867, 1997. a
Kimmritz, M., Danilov, S., and Losch, M.: The adaptive EVP method for solving the sea ice momentum equation, Ocean Model., 101, 59–67, https://doi.org/10.1016/j.ocemod.2016.03.004, 2016. a
Koldunov, N. V., Danilov, S., Sidorenko, D., Hutter, N., Losch, M., Goessling, H., Rakowsky, N., Scholz, P., Sein, D., Wang, Q., and Jung, T.: Fast EVP Solutions in a High-Resolution Sea Ice Model, J. Adv. Model. Earth Sy., 11, 1269–1284, https://doi.org/10.1029/2018MS001485, 2019. a
Levesque, J. and Vose, A.: Programming for Hybrid Multi/Manycore MPP systems, CRC Press, Taylor and Francis Inc, ISBN 978-1-4398-7371-7, 2017. a
Lynch, P.: The Emergence of Numerical Weather Prediction, Cambridge University Press, ISBN 0521857295 9780521857291, 2006. a
Madec, G., Bell, M., Blaker, A., Bricaud, C., Bruciaferri, D., Castrillo, M., Calvert, D., Chanut, J., Clementi, E., Coward, A., Epicoco, I., Éthé, C., Ganderton, J., Harle, J., Hutchinson, K., Iovino, D., Lea, D., Lovato, T., Martin, M., Martin, N., Mele, F., Martins, D., Masson, S., Mathiot, P., Mele, F., Mocavero, S., Müller, S., Nurser, A. G., Paronuzzi, S., Peltier, M., Person, R., Rousset, C., Rynders, S., Samson, G., Téchené, S., Vancoppenolle, M., and Wilson, C.: NEMO Ocean Engine Reference Manual, Zenodo [data set], https://doi.org/10.5281/zenodo.8167700, 2023. a
Mattson, T. G., Sanders, B. A., and Massingill, B.: Patterns for parallel programming, Addison-Wesley, Boston, ISBN 0321228111 9780321228116, 2005. a
McCalpin, J. D.: Memory Bandwidth and Machine Balance in Current High Performance Computers, IEEE Computer Society Technical Committee on Computer Architecture (TCCA) Newsletter, pp. 19–25, 1995. a
OpenMP Architecture Review Board: OpenMP Application Program Interface Version 5.2, https://www.openmp.org/wp-content/uploads/OpenMP-API-Specification-5-2.pdf (last access: 27 August 2024), 2021. a
Ponsoni, L., Ribergaard, M. H., Nielsen-Englyst, P., Wulf, T., Buus-Hinkler, J., Kreiner, M. B., and Rasmussen, T. A. S.: Greenlandic sea ice products with a focus on an updated operational forecast system, Front. Marine Sci., 10, 979782, https://doi.org/10.3389/fmars.2023.979782, 2023. a, b
Rasmussen, T. A. S., Poulsen, J., and Ribergaard, M. H.: Input data for 1d EVP model, Zenodo [data set], https://doi.org/10.5281/zenodo.11248366, 2024a. a
Rasmussen, T. A. S., Poulsen, J., Ribergaard, M. H., and Rethmeier, S.: dmidk/cice-evp1d: Unit test refactorization of EVP solver CICE, Zenodo [code], https://doi.org/10.5281/zenodo.10782548, 2024b. a
Raul Torres, R. F. and Teruel, X.: A Novel Set of Directives for Multi-device Programming with OpenMP, IEEE International Parallel and Distributed Processing Symposium Workshops (IPDPSW), Lyon, France, 2022, 401–410, https://doi.org/10.1109/IPDPSW55747.2022.00075, 2022. a
Roberts, A. F., Hunke, E. C., Allard, R., Bailey, D. A., Craig, A. P., Lemieux, J.-F., and Turner, M. D.: Quality control for community-based sea-ice model development, Philos. T. Roy. Soc. A, 376, 20170344, https://doi.org/10.1098/rsta.2017.0344, 2018. a
Wadhams, P.: ON SEA ICE. Weeks. 2010, University of Alaska Press, Fairbanks, 664 pp., ISBN 978-1-60223-079-8, Polar Record, 49, e6, https://doi.org/10.1017/S0032247412000502, 2013. a
Short summary
Earth system models (ESMs) today strive for better quality based on improved resolutions and improved physics. A limiting factor is the supercomputers at hand and how best to utilize them. This study focuses on the refactorization of one part of a sea ice model (CICE), namely the dynamics. It shows that the performance can be significantly improved, which means that one can either run the same simulations much cheaper or advance the system according to what is needed.
Earth system models (ESMs) today strive for better quality based on improved resolutions and...