Winter Arctic Sea Ice Volume Budget Decomposition from Satellite Observations and Model Simulations over the CryoSat-2 period (2010-2019)

In this study, satellite-derived observations of sea ice concentration, drift, and thickness are combined to provide a climatology and inter-annual monthly variability of the sea ice volume budget over the growth season for the CryoSat-2 period Octobre 2010 to April 2019. This allows the first wintertime observational decomposition of the dynamic (advection and divergence) and thermodynamic drivers of ice volume change.

Dynamic and thermodynamic processes will be separated by applying similar methods to Holland and Kwok (2012), which decomposed the governing equation of ice concentration, C:

[1] ∂C/∂t + ∇.(uC) = f_C - r

where u is ice drift motion, f_C is thermodynamic freezing or melting, and r is the concentration change from mass-conserving mechanical ice redistribution processes which convert ice area to ice thickness, such as ridging and rafting. ∂C/∂t represents ice intensification and ∇.(uC) represents ice flux divergence, the dynamic contribution to ice concentration.

Equation 1 can be rearranged and the dynamic contribution, ∇.(uC), expanded to show the contributions from advection, (-u.∇C) and divergence, (-C∇.u), determining four terms:

[2] ∂C/∂t = -u.∇C - C∇.u + f - r

The governing equation of the volume budget is of the same form but combines the thickness and concentration data:

[3] ∂Ch/∂t = -u.∇Ch - Ch∇.u + f_Ch

where h is the thickness of the ice and the resulting product of the two datasets, Ch, is the effective thickness. If Ch were to be multiplied by the grid cell area this would give the volume, V. It is not necessary to take this step because the area remains constant and does not influence the relative values of the terms. However, when deriving an Arctic-wide climatology spatial integration across the grid cell areas is required. Ridging is taken into account by effective thickness change, therefore, it is not included in the calculations.

The method is replicated using model simulations from the Centre for Climate Observation and Modelling (CPOM)-modified Los Alamos sea-ice model (CICE), providing a test of the model’s ability to calculate the volume budgets but also identifying unrealistic growth regimes in the CryoSat-2 observational datasets. Sensitivity to several observational datasets is performed to provide an estimated uncertainty of the budget calculations.

The observational results show ice gain in the central Arctic is dominated by ice freezing with contributions from convergence. Divergence at the coastlines of the Arctic form an ice sink where freezing generates new ice. Advection is shown to drive ice equatorward and induce melting at the ice edge where ice becomes thermodynamically unstable. The dynamic components are found to grow in influence throughout the growth season.

The 2016/17 winter growth season budget shows reduced thermodynamic intensification and stronger dynamic tendencies which may be in response to thin initial ice and an exceptionally warm winter. Compared to the observed volume budget, the CICE model displays similar patterns of thermodynamic freezing, however, dynamic components in the central Arctic are significantly reduced whilst they are over-amplified at the ice edge.