CICE (sea ice model)

CICE (/saɪs/) is a computer model that simulates the growth, melt and movement of sea ice. It has been integrated into many coupled climate system models as well as global ocean and weather forecasting models and is often used as a tool in Arctic and Southern Ocean research.^{[1][2][3][4][5][6][7][8][9]} CICE development began in the mid-1990s by the United States Department of Energy (DOE), and it is currently maintained and developed by a group of institutions in North America and Europe known as the CICE Consortium.^[10] Its widespread use in Earth system science in part owes to the importance of sea ice in determining Earth's planetary albedo, the strength of the global thermohaline circulation in the world's oceans, and in providing surface boundary conditions for atmospheric circulation models, since sea ice occupies a significant proportion (4-6%) of Earth's surface.^[11][12] CICE is a type of cryospheric model.

Development

Development of CICE began in 1994 by Elizabeth Hunke at

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]

Code Design

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.

See also

• 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

References