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: We present a new framework for global ocean- sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean-sea-ice models (JRA55-do).We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean-ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean-sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80% of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP- 2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP- 2. For example, the sea surface temperatures of the OMIP- 2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating processlevel responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean-sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

Description: The Atlantic meridional overturning circulation (AMOC) represents the zonally integrated stream function of meridional volume transport in the Atlantic Basin. The AMOC plays an important role in transporting heat meridionally in the climate system. Observations suggest a heat transport by the AMOC of 1.3 PW at 26°N—a latitude which is close to where the Atlantic northward heat transport is thought to reach its maximum. This shapes the climate of the North Atlantic region as we know it today. In recent years there has been significant progress both in our ability to observe the AMOC in nature and to simulate it in numerical models. Most previous modeling investigations of the AMOC and its impact on climate have relied on models with horizontal resolution that does not resolve ocean mesoscale eddies and the dynamics of the Gulf Stream/North Atlantic Current system. As a result of recent increases in computing power, models are now being run that are able to represent mesoscale ocean dynamics and the circulation features that rely on them. The aim of this review is to describe new insights into the AMOC provided by high-resolution models. Furthermore, we will describe how high-resolution model simulations can help resolve outstanding challenges in our understanding of the AMOC.

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: Changes in the Atlantic Meridional Overturning Circulation (AMOC) represent a crucial component of Northern Hemisphere climate variability. In modelling studies decadal overturning variability has been attributed to the intensity of deep winter convection in the Labrador Sea. This linkage is challenged by transport observations at sections across the subpolar gyre. Here we report simulations with an eddy-rich ocean model which captures the observed concentration of downwelling in the northeastern Atlantic and the negligible impact of interannual variations in Labrador Sea convection during the last decade. However, the exceptionally cold winters in the Labrador Sea during the first half of the 1990s induced a positive AMOC anomaly of more than 20%, mainly by augmenting the downwelling in the northeastern North Atlantic. The remote effect of excessive Labrador Sea buoyancy forcing is related to rapid spreading of mid-depth density anomalies into the Irminger Sea and their entrainment into the deep boundary current off Greenland.

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: Marine heatwaves along the coast ofWestern Australia, referred to as Ningaloo Niño, have had dramatic impacts on the ecosystem in the recent decade. A number of local and remote forcing mechanisms have been put forward, however little is known about the depth structure of such temperature extremes. Utilizing an eddy-active global Ocean General Circulation Model, Ningaloo Niño and the corresponding cold Ningaloo Niña events are investigated between 1958-2016, with focus on their depth structure. The relative roles of buoyancy and wind forcing are inferred from sensitivity experiments. Composites reveal a strong symmetry between cold and warm events in their vertical structure and associated large-scale spatial patterns. Temperature anomalies are largest at the surface, where buoyancy forcing is dominant and extend down to 300m depth (or deeper), with wind forcing being the main driver. Large-scale subsurface anomalies arise from a vertical modulation of the thermocline, extending from the western Pacific into the tropical eastern Indian Ocean. The strongest Ningaloo Niños in 2000 and 2011 are unprecedented compound events, where long-lasting high temperatures are accompanied by extreme freshening, which emerges in association with La Niñas, more common and persistent during the negative phase of the Interdecadal Pacific Oscillation. It is shown that Ningaloo Niños during La Nina phases have a distinctively deeper reach and are associated with a strengthening of the Leeuwin Current, while events during El Niño are limited to the surface layer temperatures, likely driven by local atmosphere-ocean feedbacks, without a clear imprint on salinity and velocity.

Description: Strong regional sea-level trends, mainly related to basin-wide wind stress anomalies, have been observed in the western tropical Pacific over the last 3 decades. Analyses of regional sea level in the densely populated regions of the neighbouring Australasian Mediterranean Sea (AMS; also called tropical Asian seas) are hindered by its complex topography and respective studies are sparse. We used a series of global eddy-permitting ocean models, including a high-resolution configuration that resolves the AMS with 120∘ horizontal resolution, forced by a comprehensive atmospheric forcing product over 1958–2016 to characterize the patterns and magnitude of decadal sea-level variability in the AMS. The nature of this variability is elucidated further by sensitivity experiments with interannual variability restricted to either the momentum or buoyancy fluxes, building on an experiment employing a repeated-year forcing without interannual variability in all forcing components. Our results suggest that decadal fluctuations of the El Niño–Southern Oscillation (ENSO) account for over 80 % of the variability in all deep basins of the region, except for the central South China Sea (SCS). Changes related to the Pacific Decadal Oscillation (PDO) are most pronounced in the shallow Arafura and Timor seas and in the central SCS. On average, buoyancy fluxes account for less than 10 % of decadal SSH variability, but this ratio is highly variable over time and can reach values of up to 50 %. In particular, our results suggest that buoyancy flux forcing amplifies the dominant wind-stress-driven anomalies related to ENSO cycles. Intrinsic variability is mostly negligible except in the SCS, where it accounts for 25 % of the total decadal SSH variability.

Description: Enhanced Southern Ocean ventilation in recent decades has been suggested to be a relevant modulator of the observed changes in ocean heat and carbon uptake. This study focuses on the Southern Ocean midlatitude ventilation changes from the 1960s to the 2010s. A global 1/4° configuration of the NEMO–Louvain-la-Neuve sea ice model, version 2 (LIM2), including the inert tracer CFC-12 (a proxy of ocean ventilation) is forced with the CORE, phase II (CORE-II), and JRA-55 driving ocean (JRA55-do) atmospheric reanalyses. Sensitivity experiments, where the variability of wind stress and/or the buoyancy forcing is suppressed on interannual time scales, are used to unravel the mechanisms driving ventilation changes. Ventilation changes are estimated by comparing CFC-12 interior inventories among the different experiments. All simulations suggest a multidecadal fluctuation of Southern Ocean ventilation, with a decrease until the 1980s–90s and a subsequent increase. This evolution is related to the atmospheric forcing and is caused by the (often counteracting) effects of wind stress and buoyancy forcing. Until the 1980s, increased buoyancy gains caused the ventilation decrease, whereas the subsequent ventilation increase was driven by strengthened wind stress causing deeper mixed layers and a stronger meridional overturning circulation. Wind stress emerges as the main driver of ventilation changes, even though buoyancy forcing modulates its trend and decadal variability. The three Southern Ocean basins take up CFC-12 in distinct density intervals but overall respond similarly to the atmospheric forcing. This study suggests that Southern Ocean ventilation is expected to increase as long as the effect of increasing Southern Hemisphere wind stress overwhelms that of increased stratification.

Description: Regional anomalies of steric sea level are either due to redistribution of heat and freshwater anomalies or due to ocean-atmosphere buoyancy fluxes. Interannual to decadal variability in sea level across the tropical Pacific is mainly due to steric variations driven by wind stress anomalies. The importance of air--sea buoyancy fluxes is less clear. We use a global, eddy permitting ocean model and a series of sensitivity experiments with quasi-climatological momentum and buoyancy fluxes to identify the contribution of buoyancy fluxes for interannual to decadal sea level variability in the tropical Pacific. We find their contribution on interannual timescales to be strongest in the central tropical Pacific at around 10° latitude in both hemispheres and also relevant in the very east of the tropical domain. Buoyancy flux forced anomalies are in phase with variations driven by wind stress changes but their effect on the prevailing anomalies and the importance of heat and fresh water fluxes vary locally. In the eastern tropical basin interannual sea level variability is amplified by anomalous heat fluxes, while the importance of fresh water fluxes is small and neither has any impact on decadal timescales. In the western tropical Pacific the variability on interannual and decadal timescales is dampened by both, heat and freshwater fluxes. The mechanism involves westward propagating Rossby waves that are triggered during ENSO events by anomalous buoyancy fluxes in the central tropical Pacific and counteract the prevailing sea level anomalies once they reach the western part of the basin.