Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.

Description: The skill of numerical Lagrangian drifter trajectories in three numerical models is assessed by comparing these numerically obtained paths to the trajectories of drifting buoys in the real ocean. The skill assessment is performed using the two-sample Kolmogorov-Smirnov statistical test. To demonstrate the assessment procedure, it is applied to three different models of the Agulhas region. The test can either be performed using crossing positions of one-dimensional sections in order to test model performance in specific locations, or using the total two-dimensional data set of trajectories. The test yields four quantities: a binary decision of model skill, a confidence level which can be used as a measure of goodness-of-fit of the model, a test statistic which can be used to determine the sensitivity of the confidence level, and cumulative distribution functions that aid in the qualitative analysis. The ordering of models by their confidence levels is the same as the ordering based on the qualitative analysis, which suggests that the method is suited for model validation. Only one of the three models, a 1/10 degree two-way nested regional ocean model, might have skill in the Agulhas region. The other two models, a 1/2 degree global model and a 1/8 degree assimilative model, might have skill only on some sections in the region.

Description: Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations. The atmospheric forcing from [Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research] was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.

Description: The greater Agulhas Current system has several components with high mesoscale turbulence. The phytoplankton distribution in the southwest Indian Ocean reflects this activity. We have used a regional eddy-permitting, coupled physical–biological model to study the physical–biological interactions and to address the main processes responsible for phytoplankton distribution in three different biogeochemical provinces: the southwest Subtropical Indian Gyre (SWSIG), the subtropical convergence zone (SCZ) and the subantarctic waters (SAW) south of South Africa. The biological model with four compartments (Nitrate–Phytoplankton–Zooplankton–Detritus) adequately reproduces the observed field of chlorophyll a. The phase of the strong modelled seasonality in the SWSIG is opposite to that of the SCZ that forms the southern boundary of the subtropical gyre. Phytoplankton concentrations are governed by the source-minus-sink terms, which are one order of magnitude greater than the dynamical diffusion and advection terms. North of 35°S, in the SWSIG, phytoplankton growth is limited by nutrients supply throughout the year. However, deeper stratification, enhanced cross-frontal transport and higher detritus remineralization explain the simulated higher concentrations of phytoplankton found in winter in the SWSIG. The region between 35° and 40°S constitutes a transition zone between the SCZ and the oligotrophic subtropical province. Horizontal advection is the main process bringing nutrients for phytoplankton growth. The front at 34°S represents a dynamical barrier to an extension further to the north of this advection of nutrients. Within the SCZ, primary production is high during spring and summer. This high productivity depletes the nutrient standing stock built up during winter time. In winter, nutrients supply in the convergence zone is indeed large, but the deep mixing removes phytoplankton from the euphotic zone and inhibits photosynthesis, yielding lower surface chlorophyll a concentrations. Waters south of the Subantarctic Front have a summer biomass close to that of frontal waters and higher than for subtropical waters. However, these simulated concentrations are slightly higher than the observed ones suggesting that limitation by iron and/or silica may play a role

Description: Warming of the North Atlantic region in climate history often was associated with massive melting of the Greenland Ice Sheet. To identify the meltwater’s impacts and isolate these from internal variability and other global warming factors, we run single-forcing simulations including small ensembles using three complex climate models differing only in their ocean components. In 200-year long pre-industrial climate simulations, we identify robust consequences of abruptly increasing Greenland runoff by 0.05 Sv: sea-level rise of 44±10 cm, subpolar North Atlantic surface cooling of 0.7˚C and a moderate AMOC decline of 1.1–2.0 Sv. The latter two emerge in under three decades—and reverse on the same timescale after the perturbation ends in year 100. The ocean translates the step-change perturbation into a multi-decadal to centennial signature in the deep overturning circulation. In all simulations, internal variability creates notable uncertainty in estimating trends, time of emergence and duration of the response.

Description: The Makassar Strait, the main passageway of the Indonesian Throughflow (ITF), is an important component of Indo-Pacific climate through its inter-basin redistribution of heat and freshwater. Observational studies suggest that wind-driven freshwater advection from the marginal seas into the Makassar Strait modulates the strait's surface transport. However, direct observations are too short (〈15 years) to resolve variability on decadal timescales. Here we use a series of global ocean simulations to assess the advected freshwater contributions to ITF transport across a range of timescales. The simulated seasonal and interannual freshwater dynamics are consistent with previous studies. On decadal timescales, we find that wind-driven advection of South China Sea (SCS) waters into the Makassar Strait modulates upper-ocean ITF transport. Atmospheric circulation changes associated with Pacific decadal variability appear to drive this mechanism via Pacific lower-latitude western boundary current interactions that affect the SCS circulation. Key Points: - A global ocean model is used to show how freshwater impacts the decadal variability of transport through the main Indonesian Throughflow pathway - Wind-driven advection of South China Sea freshwater induces an upstream pressure gradient that reduces transport - Freshwater input is modulated by atmospheric circulation changes associated with Pacific decadal variability

Description: Deep convection and associated deep water formation are key processes for climate variability, since they impact the oceanic uptake of heat and trace gases and alter the structure and strength of the global overturning circulation. For long, deep convection in the subpolar North Atlantic was thought to be confined to the central Labrador Sea in the western subpolar gyre (SPG). However, there is increasing observational evidence that deep convection also has occurred in the eastern SPG south of Cape Farewell and in the Irminger Sea, in particular, in 2015–2018. Here we assess this recent event in the context of the temporal evolution of spatial deep convection patterns in the SPG since the mid-twentieth century, using realistic eddy-rich ocean model simulations. These reveal a large interannual variability with changing contributions of the eastern SPG to the total deep convection volume. Notably, in the late 1980s to early 1990s, the period with highest deep convection intensity in the Labrador Sea related to a persistent positive phase of the North Atlantic Oscillation, the relative contribution of the eastern SPG was small. In contrast, in 2015–2018, deep convection occurred with an unprecedented large relative contribution of the eastern SPG. This is partly linked to a smaller north-westward extent of deep convection in the Labrador Sea compared to previous periods of intensified deep convection, and may be a first fingerprint of freshening trends in the Labrador Sea potentially associated with enhanced Greenland melting and the oceanic advection of the 2012–2016 eastern North Atlantic fresh anomaly.

Description: While forced ocean hindcast simulations are useful for a wide range of applications, a key limitation is their inability to simulate ocean-atmosphere feedbacks. As a consequence, they need to rely on artificial choices such as sea surface salinity restoring and other corrections affecting the surface freshwater fluxes. Fully coupled models overcome these limitations, but lack the correct timing of variability due to weaker observational constraints. This leads to a mismatch between forced and coupled models on interannual to decadal timescales. A possibility to combine the advantages of both modelling strategies is to apply a partial coupling (PCPL), i.e. replacing the surface winds stress in the ocean component by wind stress derived from reanalysis. To identify the capabilities, limitations and possible use cases of partial coupling, we perform a fully coupled, two partially coupled and an ocean-only experiment using an all-Atlantic nested ocean configuration at eddying resolution in a global climate model. We show that the correct timing of Atlantic Meridional Overturning Circulation (AMOC) variability in PCPL experiments is robust on timescales below 5 years. Mid-latitude wind stress curl changes contribute to decadal AMOC variability, but North Atlantic buoyancy fluxes are not significantly altered by incorporating reanalysed wind stress anomalies, limiting the success of PCPL on this timescale. Long term trends of the AMOC in PCPL mode are consistent with fully coupled model experiments under historic atmospheric boundary conditions, suggesting that a partially coupled model is still able to simulate the important ocean-atmosphere feedbacks necessary to maintain a stable AMOC.

Description: Because new observations have revealed that the Labrador Sea is not the primary source for waters in the lower limb of the Atlantic Meridional Overturning Circulation (AMOC) during the Overturning in the Subpolar North Atlantic Programme (OSNAP) period, it seems timely to re-examine the traditional interpretation of pathways and property variability for the AMOC lower limb from the subpolar gyre to 26.5 degrees N. In order to better understand these connections, Lagrangian experiments were conducted within an eddy-rich ocean model to track upper North Atlantic Deep Water (uNADW), defined by density, between the OSNAP line and 26.5 degrees N as well as within the Labrador Sea. The experiments reveal that 77% of uNADW at 26.5 degrees N is directly advected from the OSNAP West section along the boundary current and interior pathways west of the Mid-Atlantic Ridge. More precisely, the Labrador Sea is a main gateway for uNADW sourced from the Irminger Sea, while particles connecting OSNAP East to 26.5 degrees N are exclusively advected from the Iceland Basin and Rockall Trough along the eastern flank of the Mid-Atlantic Ridge. Although the pathways between OSNAP West and 26.5 degrees N are only associated with a net formation of 1.1 Sv into the uNADW layer, they show large density changes within the layer. Similarly, as the particles transit through the Labrador Sea, they undergo substantial freshening and cooling that contributes to further densification within the uNADW layer. Key Points: - The large majority of upper North Atlantic Deep Water (uNADW) sourced from the Irminger Sea transits through the Labrador Sea before reaching 26.5°N - Interior pathways along the eastern flank of the Mid-Atlantic Ridge connect the Iceland Basin and Rockall Trough to 26.5°N - Though uNADW is mainly sourced in the eastern subpolar gyre, its transit in the Labrador Sea is associated with further property changes