Description: The Beaufort Gyre (BG) is the largest liquid freshwater reservoir of the Arctic Ocean. The liquid freshwater content (FWC) significantly increased in the BG in the 2000s during an anticyclonic wind regime and remained at a high level despite a transition to a more cyclonic state in the early 2010s. It is not well understood to what extent the rapid sea ice decline during this period has modified the trend and variability of the BG liquid FWC in the past decade. Our numerical simulations show that about 50% of the liquid freshwater accumulated in the BG in the 2000s can be explained by the sea ice decline caused by the Arctic atmospheric warming. Among this part of the FWC increase, 60% can be attributed to surface freshening associated with the reduction of the net sea ice thermodynamic growth rate, and 40% to changes in ocean circulation, which makes freshwater more accessible to the BG for storage. Thus, the rapid increase of the BG FWC in the 2000s was due to the concurrence of the anticyclonic wind regime and the high freshwater availability. We also find that if the Arctic sea ice had not declined, the liquid FWC in the BG would have shown a stronger decreasing tendency at the beginning of the 2010s owing to the cyclonic wind regime. From our results we argue that changes in sea ice conditions should be adequately taken into account when it comes to understanding and predicting variations of BG liquid FWC in a changing climate.

Description: Sea ice formation is accompanied by the rejection of salt which in nature tends to be mixed vertically by the formation of convective plumes. Here we analyze the influence of a salt plume parameterization (SPP) in an atmosphere-sea ice-ocean model. Two 330 years long simulations have been conducted with the AWI Climate Model. In the reference simulation, the rejected salt in the Arctic Ocean is added to the upper-most ocean layer. This approach is commonly used in climate modelling. In another experiment, employing SPP, the rejected salt is vertically redistributed within the mixed layer based on a power law profile that mimics the penetration of salt plumes. We discuss the effects of this redistribution on the simulated mean state and on atmosphere-ocean linkages associated with the intensity of deep water formation. We find that the salt plume parametrization leads to simultaneous increase of sea ice (volume and concentration) and decrease of sea surface salinity in the Arctic. The SPP considerably alters the interplay between the atmosphere and the ocean in the Nordic Seas. The parameterization modifies the ocean ventilation; however, resulting changes in temperature and salinity largely compensate each other in terms of density so that the overturning circulation is not significantly affected.

Description: Sea ice dynamics determine the drift and deformation of sea ice. Nonlinear physics, usually expressed in a viscous‐plastic rheology, makes the sea ice momentum equations notoriously difficult to solve. At increasing sea ice model resolution the nonlinearities become stronger as linear kinematic features (leads) appear in the solutions. Even the standard elastic‐viscous‐plastic (EVP) solver for sea ice dynamics, which was introduced for computational efficiency, becomes computationally very expensive, when accurate solutions are required, because the numerical stability requires very short, and hence more, subcycling time steps at high resolution. Simple modifications to the EVP solver have been shown to remove the influence of the number of subcycles on the numerical stability. At low resolution appropriate solutions can be obtained with only partial convergence based on a significantly reduced number of subcycles as long as the numerical procedure is kept stable. This previous result is extended to high resolution where linear kinematic features start to appear. The computational cost can be strongly reduced in Arctic Ocean simulations with a grid spacing of 4.5 km by using modified and adaptive EVP versions because fewer subcycles are required to simulate sea ice fields with the same characteristics as with the standard EVP.

Description: Models from phase 5 of the Coupled Model Inter-comparison Project (CMIP5) show substantial biases in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. Here, we analyze the influence of horizontal resolution in a hierarchy of five multi-resolution simulations with the AWI Climate Model (AWI-CM), the climate model used at the Al-fred Wegener Institute, Helmholtz Centre for Polar and Ma-rine Research, which employs a sea ice–ocean model com-ponent formulated on unstructured meshes. The ocean grid sizes considered range from a nominal resolution of ∼ 1◦ (CMIP5 type) up to locally eddy resolving. We show that increasing ocean resolution locally to resolve ocean eddies leads to reductions in deep ocean biases, although these im-provements are not strictly monotonic for the five different ocean grids. A detailed diagnosis of the simulations allows to identify the origins of the biases. We find that two key re-gions at the surface are responsible for the development of the deep bias in the Atlantic Ocean: the northeastern North Atlantic and the region adjacent to the Strait of Gibraltar. Furthermore, the Southern Ocean density structure is equally improved with locally explicitly resolved eddies compared to parameterized eddies. Part of the bias reduction can be traced back towards improved surface biases over outcrop-ping regions, which are in contact with deeper ocean layers along isopycnal surfaces. Our prototype simulations provide guidance for the optimal choice of ocean grids for AWI-CM to be used in the final runs for phase 6 of CMIP (CMIP6) and for the related flagship simulations in the High Resolution Model Intercomparison Project (HighResMIP). Quite remarkably, retaining resolution only in areas of high eddy activity along with excellent scalability characteristics of the unstructured-mesh sea ice–ocean model enables us to per-form the multi-centennial climate simulations needed in a CMIP context at (locally) eddy-resolving resolution with a throughput of 5–6 simulated years per day.

Description: Freshwater in the Arctic Ocean is one of the key climate components. It is not well understood how the capability of the Arctic Ocean to store freshwater will develop when freshwater supplies increase in a warming climate. By using numerical experiments, we find that this capability varies with the Arctic sea ice decline nonmonotonically, with the largest capability at intermediate strength of sea ice decline. Through enhancing the ocean surface stress, sea ice decline not only accumulates freshwater toward the Amerasian Basin but also tends to reduce the amount of freshwater in both the Eurasian and Amerasian basins by increasing the occupation of Atlantic-origin water in the upper ocean. An increase in river runoff modulates the counterbalance of the two competing effects, leading to the nonmonotonic changes of the Arctic freshwater storage capability in a warming climate.

Description: Ocean heat transport through the Barents Sea Opening (BSO) has strong impacts on the Barents Sea ice extent and the climate. In this paper we quantified the contributions from different atmospheric forcing components to the trend and interannual variability of the BSO heat transport. Ocean‐ice model simulations were conducted in which the interannual variation of atmospheric forcing was maintained only in or outside the Arctic in two different simulations. The sum of their BSO heat transport anomalies reasonably replicated the trend and variability from a hindcast simulation. The upward trend of the BSO heat transport mainly stems from the increasing ocean temperature in the subpolar North Atlantic. For the interannual variability, the local wind and upstream forcing are similarly important. The location of the Atlantic Water boundary current in the Nordic Seas, influenced by the cyclonic atmospheric circulation, is crucial in determining part of the BSO inflow variability.

Description: A new global climate model setup using FESOM2.0 for the sea ice‐ocean component and ECHAM6.3 for the atmosphere and land surface has been developed. Replacing FESOM1.4 by FESOM2.0 promises a higher efficiency of the new climate setup compared to its predecessor. The new setup allows for long‐term climate integrations using a locally eddy‐resolving ocean. Here it is evaluated in terms of (1) the mean state and long‐term drift under preindustrial climate conditions, (2) the fidelity in simulating the historical warming, and (3) differences between coarse and eddy‐resolving ocean configurations. The results show that the realism of the new climate setup is overall within the range of existing models. In terms of oceanic temperatures, the historical warming signal is of smaller amplitude than the model drift in case of a relatively short spin‐up. However, it is argued that the strategy of “de‐drifting” climate runs after the short spin‐up, proposed by the HighResMIP protocol, allows one to isolate the warming signal. Moreover, the eddy‐permitting/resolving ocean setup shows notable improvements regarding the simulation of oceanic surface temperatures, in particular in the Southern Ocean.

Description: Simulating Arctic Ocean mesoscale eddies in ocean circulation models presents a great challenge because of their small size. This study employs an unstructured‐mesh ocean‐sea ice model to conduct a decadal‐scale global simulation with a 1‐km Arctic. It provides a basinwide overview of Arctic eddy energetics. Increasing model resolution from 4 to 1 km increases Arctic eddy kinetic energy (EKE) and total kinetic energy (TKE) by about 40% and 15%, respectively. EKE is the highest along main currents over topography slopes, where strong conversion from available potential energy to EKE takes place. It is high in halocline with a maximum typically centered in the depth range of 70–110 m, and in the Atlantic Water layer of the Eurasian Basin as well. The seasonal variability of EKE along the continental slopes of southern Canada and eastern Eurasian basins is similar, stronger in fall and weaker in spring.

Description: The response of the oceanic circulation to mesoscale atmospheric forcing is studied by comparing integrations of a global sea ice-ocean model with high-resolution European Centre for Medium-Range Weather Forecasts analysis data (0.4◦) to those with the same forcing coarse grained to a resolution typically employed in climate models and atmospheric reanalyses (1.8◦). It is shown that the representation of mesoscale features in atmospheric forcing fields leads to an increase in the strength of the wind-driven gyres in the North Atlantic and North Pacific regions of about 5–10% of its mean value. An increase of similar magnitude is found for the Atlantic meridional overturning circulation. From the results of this study it is argued that small-scale atmospheric phenomena such as fronts, mesoscale cyclones, and topographic jets play an important role in driving the mean oceanic circulation.