CICE (sea ice model)

Output from CICE within a coupled climate model: Averaged 2000-2004 (a) March and (b) September Antarctic sea ice thickness and extent (sea ice with greater than 15% concentration) of five ensemble members from the Community Earth System Model (CESM) large ensemble.[13] The magenta contour is the measured ice edge according to the NOAA Climate Data Record.[14]

Development of CICE began in 1994 by Elizabeth Hunke at Los Alamos National Laboratory (LANL).[12][15] Since its initial release in 1998 following development of the Elastic-Viscous-Plastic (EVP) sea ice rheology within the model,[16] it has been substantially developed by an international community of model users and developers. Enthalpy-conserving thermodynamics and improvements to the sea ice thickness distribution were added to the model between 1998 and 2005.[17][18][19] The first institutional user outside of LANL was Naval Postgraduate School[15] in the late-1990s, where it was subsequently incorporated into the Regional Arctic System Model (RASM) in 2011.[20][21] The National Center for Atmospheric Research (NCAR) was the first to incorporate CICE into a global climate model in 2002,[22] and developers of the NCAR Community Earth System Model (CESM) have continued to contribute to CICE innovations[23][24][25] and have used it to investigate polar variability in Earth's climate system.[13] The United States Navy began using CICE shortly after 2000 for polar research and sea ice forecasting and it continues to do so today.[3][26] Since 2000, CICE development or coupling to oceanic and atmospheric models for weather and climate prediction has occurred at the University of Reading,[27] University College London,[28] the U.K. Met Office Hadley Centre,[29] Environment and Climate Change Canada,[7] the Danish Meteorological Institute,[4] the Commonwealth Science and Industrial Research Organisation,[30] and Beijing Normal University,[8] among other institutions. As a result of model development in the global community of CICE users, the model's computer code now includes a comprehensive saline ice physics and biogeochemistry library that incorporates mushy-layer thermodynamics,[31][32] anisotropic continuum mechanics,[33] Delta-Eddington radiative transfer,[34] melt-pond physics[35][36] and land-fast ice.[37] CICE version 6 is open-source software and was released in 2018 on GitHub.[38]

There are two main physics equations solved using numerical methods in CICE that underpin the model's predictions of sea ice thickness, concentration and velocity, as well as predictions made with many equations not shown here giving, for example, surface albedo, ice salinity, snow cover, divergence, and biogeochemical cycles. The first keystone equation is Newton's second law for sea ice:

${\displaystyle m{\frac {d\mathbf {u} }{dt}}=-mf\mathrm {k} \times \mathbf {u} +\tau _{a}+\tau _{w}-m\mathrm {\hat {g}} \nabla \mu +\nabla \cdot \mathbf {\sigma } }$

where ${\displaystyle m}$ is the mass per unit area of saline ice on the sea surface, ${\displaystyle \mathbf {u} }$ is the drift velocity of the ice, ${\displaystyle f}$ is the Coriolis parameter, ${\displaystyle \mathrm {k} }$ is the upward unit vector normal to the sea surface, ${\displaystyle \tau _{a}}$ and ${\displaystyle \tau _{w}}$ are the wind and water stress on the ice, respectively, ${\displaystyle \mathrm {\hat {g}} }$ is acceleration due to gravity, ${\displaystyle \mu }$ is sea surface height and ${\displaystyle \mathbf {\sigma } }$ is internal ice the two-dimensional stress tensor within the ice.[16] Each of the terms require information about the ice thickness, roughness, and concentration, as well as the state of the atmospheric and oceanic boundary layers. Ice mass per unit area ${\displaystyle m}$ is determined using the second keystone equation in CICE, which describes evolution of the sea ice thickness distribution ${\displaystyle g(h)}$ for different thicknesses ${\displaystyle h}$ spread of the area for which sea ice velocity is calculated above:[18]

${\displaystyle {\frac {dg}{dt}}=\theta +\psi -g(\nabla \cdot \mathbf {u} )}$

where ${\displaystyle \theta }$ is the change in the thickness distribution due to thermodynamic growth and melt, ${\displaystyle \psi }$ is redistribution function due to sea ice mechanics and is associated with internal ice stress ${\displaystyle \mathbf {\sigma } }$, and ${\displaystyle -g(\nabla \cdot \mathbf {u} )}$ describes advection of sea ice in a Lagrangian reference frame.[18][19] From this, ice mass is given by:

${\displaystyle m=\rho \int _{0}^{\infty }h\,g(h)\,dh}$

for density ${\displaystyle \rho }$ of sea ice.[38]

Schematic demonstrating placement of Icepack, in which the thickness distribution ${\displaystyle g(h)}$ is represented (blue), within the MPAS dycore (green) that solves for momentum evolution and horizontal sea ice advection on the E3SM unstructured grid (arrows)

CICE version 6 is coded in FORTRAN90. It is organized into a dynamical core (dycore) and a separate column physics package called Icepack, which is maintained as a CICE submodule on GitHub.[39] The momentum equation and thickness advection described above are time-stepped on a quadrilateral Arakawa B-grid within the dynamical core, while Icepack solves diagnostic and prognostic equations necessary for calculating radiation physics, hydrology, thermodynamics, and vertical biogeochemistry, including terms necessary to calculate ${\displaystyle \tau _{a}}$, ${\displaystyle \tau _{w}}$, ${\displaystyle \mathbf {\sigma } }$, ${\displaystyle \theta }$, and ${\displaystyle \psi }$ defined above. CICE can be run independently, as in the first figure on this page, but is frequently coupled with earth systems models through an external flux coupler, such as the CESM Flux Coupler from NCAR[22] for which results are shown in the second figure for the CESM Large Ensemble. The column physics were separated into Icepack for the version 6 release to permit insertion into earth system models that use their own sea ice dynamical core, including the new DOE Energy Exascale Earth System Model (E3SM),[38][40] which uses an unstructured grid in the sea ice component of the Model for Prediction Across Scales (MPAS),[41][42] as demonstrated in the final figure.

• Sea ice
• Sea ice microbial communities
• Sea ice emissivity modeling
• Sea ice growth processes
• Sea ice concentration
• Sea ice thickness
• Sea ice physics and ecosystem experiment
• Arctic Ocean
• Southern Ocean
• Climate model
• Weather forecasting
• Northern Sea Route
• Northwest Passage
• Antarctica

• CICE Consortium GitHub Information Page
• CICE Consortium Model for Sea-Ice Development
• Icepack: Essential Physics for Sea Ice Models
• Community-Driven Sea Ice Modeling with the CICE Consortium (Witness the Arctic)
• NOAA press release
• Oceans Deeply
• Pacific Standard
• phys.org: Arctic ice model upgrade to benefit polar research, industry and military
• Sea ice:  More than just frozen water (Santa Fe New Mexican)
• Energy Exascale Earth System Model (E3SM)
• Community Earth System Model (CESM)