Numerical investigation of the Arctic ice-ocean boundary layer; implications for air-sea gas fluxes
Abstract. In ice-covered regions it can be challenging to determine air-sea exchange – for heat and momentum, but also for gases like carbon dioxide and methane. The harsh environment and relative data scarcity make it difficult to characterize even the physical properties of the ocean surface. Here, we seek a mechanistic interpretation for the rate of air-sea gas exchange (k) derived from radon-deficits. These require an estimate of the water column history extending 30 days prior to sampling. We used coarse resolution (36 km) regional configuration of the MITgcm with fine near surface vertical spacing (2 m) to evaluate the capability of the model to reproduce conditions prior to sampling. The model is used to estimate sea-ice velocity, concentration and mixed-layer depth experienced by the water column. We then compared the model results to existing field data including satellite, moorings and Ice-tethered profilers. We found that model-derived sea-ice coverage is 88 to 98 % accurate averaged over Beaufort Gyre, sea-ice velocities have 78 % correlation which resulted in 2 km/day error in 30 day trajectory of sea-ice. The model demonstrated the capacity to capture the broad trends in the mixed layer although with a bias and model water velocities showed only 29 % correlation with actual data. Overall, we find the course resolution model to be an inadequate surrogate for sparse data, however the simulation results are a slight improvement over several of the simplifying assumptions that are often made when surface ocean geochemistry, including the use of a constant mixed layer depth and a velocity profile that is purely wind-driven.