Multiscale mushy layer model for Arctic marginal ice zone dynamics

Perhaps the most dynamic component of the Arctic sea ice cover is the marginal ice zone (MIZ), the transitional region between dense pack ice to the north and open ocean to the south. It widens by a factor of four while seasonally migrating more than 1600 km poleward in the Bering–Chukchi Sea sector, impacting climate dynamics, ecological processes, and human accessibility to the Arctic. Here we showcase a transformative mathematical modeling approach to understanding changes in MIZ location and width, focusing on their seasonal cycles as observed by satellites. We view the MIZ as a liquid-solid phase transition region, or mushy zone, on the scale of the Arctic Ocean. Invoking the physics of phase changes, the MIZ is modeled as a dynamic, multiscale composite material layer; this model captures 96% of the annual cycle of MIZ location and 78% of the annual cycle of MIZ width. Temperature in the upper ocean is described by a nonlinear heat equation with effective parameters obtained using homogenization theory for a random medium of ice floes in a sea water host. Observations and simulations together indicate that MIZ location closely tracks the below-ice 273 K isotherm while the width of the MIZ follows vertical heat flux convergence, but with a three-week lag.