Description: Historical hydrographic data (1940s–2010) show a distinct cross-slope difference of the lower halocline water (LHW) over the Laptev Sea continental margins. Over the slope, the LHW is on average warmer and saltier by 0.2°C and 0.5 psu, respectively, relative to the off-slope LHW. The LHW temperature time series constructed from the on-slope historical records are related to the temperature of the Atlantic Water (AW) boundary current transporting warm water from the North Atlantic Ocean. In contrast, the on-slope LHW salinity is linked to the sea ice and wind forcing over the potential upstream source region in the Barents and northern Kara Seas, as also indicated by hydrodynamic model results. Over the Laptev Sea continental margin, saltier LHW favors weaker salinity stratification that, in turn, contributes to enhanced vertical mixing with underlying AW.

Description: Despite the efforts of the modelling community to improve the representation of the sea surface temperature (SST) over the South Eastern Tropical Atlantic, warm biases still persist. In this work we use four different configurations of the fully-coupled AWI Climate Model (AWI-CM) which allow us to gain physics-based insight into the role of the oceanic and atmospheric resolutions of the model in the regional distribution of the SST. Our results show that a sole refinement of the oceanic resolution reduces warm biases further than a single increase of the atmospheric component. An increased oceanic resolution is required (i) to simulate properly the Agulhas Current and its associated rings; (ii) to reinforce the northward-flowing Benguela Current and (iii) to intensify coastal upwelling. The best results are obtained when both resolutions are refined. However, even in that case, warm biases persist, reflecting that some processes and feedbacks are still not optimally resolved. Our results indicate that overheating is not due to insufficient upwelling, but rather due to upwelling of waters which are warmer than observations as a result of an erroneous representation of the vertical distribution of temperature. Errors in the representation of the vertical temperature profile are the consequence of a warm bias in the simulated climate state.

Description: Climate simulations for the North Atlantic and Europe for recent and future conditions simulated with the regionally coupled ROM model are analyzed and compared to the results from the MPI‐ESM. The ROM simulations also include a biogeochemistry and ocean tides. For recent climate conditions, ROM generally improves the simulations compared to the driving model MPI‐ESM. Reduced oceanic biases in the Northern Atlantic are found, as well as a better simulation of the atmospheric circulation, notably storm tracks and blocking. Regarding future climate projections for the 21st century following the RCP 4.5 and 8.5 scenarios, MPI‐ESM and ROM largely agree qualitatively on the climate change signal over Europe. However, many important differences are identified. For example, ROM shows an SST cooling in the Subpolar Gyre which is not present in MPI‐ESM. Under the RCP8.5 scenario, ROM Arctic sea ice cover is thinner and reaches the seasonally ice‐free state by 2055, well before MPI‐ESM. This shows the decisive importance of higher ocean resolution and regional coupling for determining the regional responses to global warming trends. Regarding biogeochemistry, both ROM and MPI‐ESM simulate a widespread decline in winter nutrient concentration in the North Atlantic of up to ~35%. On the other hand, the phytoplankton spring bloom in the Arctic and in the North‐Western Atlantic starts earlier and the yearly primary production is enhanced in the Arctic in the late 21st century. These results clearly demonstrate the added value of ROM to determine more detailed and more reliable climate projections at the regional scale.

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 climate over the Arctic has undergone changes in recent decades. In order to evaluate the coupled response of the Arctic system to external and internal forcing, our study focuses on the estimation of regional climate variability and its dependence on large-scale atmospheric and regional ocean circulations. A global ocean–sea ice model with regionally high horizontal resolution is coupled to an atmospheric regional model and global terrestrial hydrology model. This way of coupling divides the global ocean model setup into two different domains: one coupled, where the ocean and the atmosphere are interacting, and one uncoupled, where the ocean model is driven by prescribed atmospheric forcing and runs in a so-called stand-alone mode. Therefore, selecting a specific area for the regional atmosphere implies that the ocean–atmosphere system can develop ‘freely’ in that area, whereas for the rest of the global ocean, the circulation is driven by prescribed atmospheric forcing without any feedbacks. Five different coupled setups are chosen for ensemble simulations. The choice of the coupled domains was done to estimate the influences of the Subtropical Atlantic, Eurasian and North Pacific regions on northern North Atlantic and Arctic climate. Our simulations show that the regional coupled ocean–atmosphere model is sensitive to the choice of the modelled area. The different model configurations reproduce differently both the mean climate and its variability. Only two out of five model setups were able to reproduce the Arctic climate as observed under recent climate conditions (ERA-40 Reanalysis). Evidence is found that the main source of uncertainty for Arctic climate variability and its predictability is the North Pacific. The prescription of North Pacific conditions in the regional model leads to significant correlation with observations, even if the whole North Atlantic is within the coupled model domain. However, the inclusion of the North Pacific area into the coupled system drastically changes the Arctic climate variability to a point where the Arctic Oscillation becomes an ‘internal mode’ of variability and correlations of year-to-year variability with observational data vanish. In line with previous studies, our simulations provide evidence that Arctic sea ice export is mainly due to ‘internal variability’ within the Arctic region. We conclude that the choice of model domains should be based on physical knowledge of the atmospheric and oceanic processes and not on ‘geographic’ reasons. This is particularly the case for areas like the Arctic, which has very complex feedbacks between components of the regional climate system.

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: We propose to make the damping time scale, which governs the decay of pseudo-elastic waves in the Elastic Viscous Plastic (EVP) sea-ice solvers, independent of the external time step and large enough to warrant numerical stability for a moderate number of internal time steps. A necessary condition is that the forcing on sea ice varies slowly on the damping time scale, in which case an EVP solution may still approach a Viscous Plastic one, but on a time scale longer than a single external time step. In this case, the EVP method becomes very close to the recently proposed modified EVP (mEVP) method in terms of stability and simulated behavior. In a simple test case dealing with sea ice breaking under the forcing of a moving cyclone, the EVP method with an enlarged damping time scale can simulate linear kinematic features which are very similar to those from the traditional EVP implementation, although a much smaller number of internal time steps is used. There is more difference in sea-ice thickness and linear kinematic features simulated in a realistic Arctic configuration between using the traditional and our suggested choices of EVP damping time scales, but it is minor considering model uncertainties associated with choices of many other parameters in sea-ice models.