Recently, multidecadal variability in the Southern Ocean has been found in a strongly eddying global ocean circulation model. In this paper, we study the Lorenz energy cycle of this so‐called Southern Ocean Mode (SOM). The Lorenz energy cycle analysis provides details on the energy pathways associated with the SOM. It shows that ocean eddies and the baroclinic energy pathway together with variations in the kinetic energy input by the wind are crucial aspects of the variability. It is also shown how convective mixing, which is induced by the SOM in particular in the Weddell Gyre, is responsible for the large‐scale multidecadal variability in Antarctic Bottom Water and Atlantic Meridional Overturning Circulation.

Abstract. We investigate the freshwater budget of the Atlantic and Arctic oceans in coupled climate change simulations with the Community Earth System Model and compare a strongly eddying setup with 0.1∘ ocean grid spacing to a non-eddying 1∘ configuration typical of Coupled Model Intercomparison Project phase 6 (CMIP6) models. Details of this budget are important to understand the evolution of the Atlantic Meridional Overturning Circulation (AMOC) under climate change. We find that the slowdown of the AMOC in the year 2100 under the increasing CO2 concentrations of the Representative Concentration Pathway 8.5 (RCP8.5) scenario is almost identical between both simulations. Also, the surface freshwater fluxes are similar in their mean and trend under climate change in both simulations. While the basin-scale total freshwater transport is similar between the simulations, significant local differences exist. The high-ocean-resolution simulation exhibits significantly reduced ocean state biases, notably in the salt distribution, due to an improved circulation. Mesoscale eddies contribute considerably to the freshwater and salt transport, in particular at the boundaries of the subtropical and subpolar gyres. Both simulations start in the single equilibrium AMOC regime according to a commonly used AMOC stability indicator and evolve towards the multiple equilibrium regime under climate change, but only the high-resolution simulation enters it due to the reduced biases in the freshwater budget.

Abstract. Climate variability on multidecadal timescales appears to be organized in pronounced patterns with clear expressions in sea surface temperature, such as the Atlantic Multidecadal Variability and the Pacific Decadal Oscillation. These patterns are now well studied both in observations and global climate models and are important in the attribution of climate change. Results from CMIP5 models have indicated large biases in these patterns with consequences for ocean heat storage variability and the global mean surface temperature. In this paper, we use two multi-century Community Earth System Model simulations at coarse (1∘) and fine (0.1∘) ocean model horizontal grid spacing to study the effects of the representation of mesoscale ocean flows on major patterns of multidecadal variability. We find that resolving mesoscale ocean flows both improves the characteristics of the modes of variability with respect to observations and increases the amplitude of the heat content variability in the individual ocean basins. In the strongly eddying model, multidecadal variability increases compared to sub-decadal variability. This shift of spectral power is seen in sea surface temperature indices, basin-scale surface heat fluxes, and the global mean surface temperature. This implies that the current CMIP6 model generation, which predominantly does not resolve the ocean mesoscale, may systematically underestimate multidecadal variability.

Abstract. A major source of uncertainty in future sea-level projections is the ocean-driven basal melt of Antarctic ice shelves. Whereas ice sheet models require a kilometer-scale resolution to realistically resolve ice shelf stability and grounding line migration, global or regional 3D ocean models are computationally too expensive to produce basal melt forcing fields at this resolution. To bridge this resolution gap, we introduce the 2D numerical model LADDIE (one-Layer Antarctic model for Dynamical Downscaling of Ice–ocean Exchanges) which allows for the computationally efficient modeling of basal melt rates. The model is flexible, and can be forced with output from coarse 3D ocean models or with vertical profiles of offshore temperature and salinity. In this study, we describe the model equations and numerics. To illustrate and validate the model performance, we apply the model to two test cases: the small Crosson-Dotson Ice Shelf in the warm Amundsen Sea region, and the large Filchner-Ronne Ice Shelf in the cold Weddell Sea. At ice-shelf wide scales, LADDIE reproduces observed patterns of basal melt and freezing that are also well reproduced by 3D ocean models. At scales of 0.5–5 km, which are unresolved by 3D ocean models and poorly constrained by observations, LADDIE produces plausible basal melt patterns. Most significantly, the simulated basal melt patterns are physically consistent with the applied ice shelf topography. These patterns are governed by the topographic steering and Coriolis deflection of meltwater flows, two processes that are poorly represented in basal melt parameterisations. The kilometer-scale melt patterns simulated by LADDIE include enhanced melt rates in basal channels, in some shear margins, and nearby grounding lines. As these regions are critical for ice shelf stability, we conclude that LADDIE can provide detailed basal melt patterns at the essential resolution that ice sheet models require. The physical consistency between the applied geometry and the simulated basal melt fields indicates that LADDIE can play a valuable role in the development of coupled ice–ocean modeling.

Abstract. While a global acceleration of sea-level rise (SLR) during the 20th century is now established, locally acceleration is more difficult to detect because additional processes play a role which sometimes mask the acceleration. Here we study the rate of SLR along the coast of the Netherlands from six tide gauge records, covering the period 1890–2000. We focus on the influence of the wind field and the nodal tide variations on the local sea-level trend. We use four generalised additive models, including different predictive variables, and a parametric bootstrap method to compute the sea-level trend. From the sea-level trend, we obtain the continuous evolution of the rate of SLR and its uncertainty over the observational period through differentiation. Accounting for the nodal cycle only or both the nodal cycle and the wind influence on sea level reduces the standard error in the estimation of the rate of SLR. Moreover, accounting for both the nodal and wind influence changes the estimated rate of SLR, unmasking an acceleration of SLR that started in the 1960s. Our best-fitting statistical model yields a rate of SLR of about 1.8 [1.4–2.3] mm/yr in 1900–1919 and 1.5 [1.1–1.8] mm/yr in 1940–1959 compared to 3.0 [2.4–3.5] mm/yr over 2000–2019. If, apart from tidal, wind effects and fluctuations, sea level would have increased at a constant rate, then the probability (the p-value) of finding a rate difference between 1940–1959 and 2000–2019 of at least our estimate is smaller than 1 %. Our findings can be interpreted as an unequivocal sign of the acceleration of current SLR along the Dutch coast since the 1960s. This aligns with global SLR observations and expectations based on a physical understanding of SLR related to global warming. A small but significant part of the long-term sea-level trend is due to wind forcing related to a strengthening and northward shift of the jet stream. Additionally, we detect a multidecadal mode of sea-level variability forced by the wind with an amplitude of around 1 cm. We argue that it is related to multi-decadal sea surface temperature variations in the North Atlantic, similar to the Atlantic Multidecadal Variability.