Description: Previous studies found a negative trend in oxygen concentrations in tropical regions during the last decades. Employing a biogeochemical circulation model, we highlight the importance of wind driven ocean transport associated with the Subtropical-Tropical Cells (STCs) in setting the oxygen levels in the tropical ocean. The observed and simulated slowdown of the STCs by 30 percent from the 1960s to the 1990s caused a decrease in oxygen transport to the tropics. Transport of phosphate was similarly reduced, decreasing export production and respiration. The effects of physical transport and biological consumption partly compensate, damping oxygen interannual and decadal variability. Our results suggest that the observed residual oxygen trend in the tropical Pacific is mainly driven by changes in oxygen transport. Accordingly, the observed recent strengthening of the STCs leads us to expect a pause in the oxygen decrease or even an increase of tropical Pacific oxygen values in the near future.

Description: The response of eddy kinetic energy (EKE) to the strengthening of Southern Hemisphere winds occurring since the 1950s is investigated with a global ocean model having a resolution of 1/12° in the Antarctic Circumpolar Current domain. The simulations expose regional differences in the relative importance of stochastic and wind-related contributions to inter-annual EKE changes. In the Pacific and Indian sectors the model captures the EKE variability observed since 1993 and confirms previous hypotheses of a lagged response to regional wind stress anomalies. Here, the multi-decadal trend in wind stress is reflected in an increase in EKE typically exceeding 5 cm2 sec-2 decade-1. In the western Atlantic EKE variability is mostly stochastic, is weakly correlated with wind fluctuations, and its multi-decadal trends are close to zero. The non-uniform distribution of wind-related changes in the eddy activity could affect the regional patterns of ocean circulation and biogeochemical responses to future climate change.

Description: The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs. OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations

Description: A global ocean model with 1/12∘ horizontal resolution is used to assess the seasonal cycle of surface Eddy Kinetic Energy (EKE). The model reproduces the salient features of the observed mean surface EKE, including amplitude and phase of its seasonal cycle in most parts of the ocean. In all subtropical gyres of the Pacific and Atlantic, EKE peaks in summer down to a depth of ∼350 m, below which the seasonal cycle is weak. Investigation of the possible driving mechanisms reveals the seasonal changes in the thermal interactions with the atmosphere to be the most likely cause of the summer maximum of EKE. The development of the seasonal thermocline in spring and summer is accompanied by stronger mesoscale variations in the horizontal temperature gradients near the surface which corresponds, by thermal wind balance, to an intensification of mesoscale velocity anomalies towards the surface.

Description: The repercussions of surface ocean currents for the near-surface wind and the air-sea momentum flux are investigated in two versions of a global climate model with eddying ocean. The focus is on the effect of mesoscale ocean current features at scales of less than 150 km, by considering high-pass filtered, monthly-mean model output fields. We find a clear signature of a mesoscale oceanic imprint in the wind fields over the energetic areas of the oceans, particularly along the extensions of the western boundary currents and the Antarctic Circumpolar Current. These areas are characterized by a positive correlation between mesoscale perturbations in the curl of the surface currents and the wind curl. The coupling coefficients are spatially non-uniform and show a pronounced seasonal cycle. The positive feedback of mesoscale current features on the near-surface wind acts in opposition to their damping effect on the wind stress. A tentative incorporation of this feedback in the surface stress formulation of an eddy-permitting global ocean-only model leads to a gain in the kinetic energy of up to 10 %, suggesting a fundamental shortcoming of present ocean model configurations.

Description: We examine the mean pathways, transit timescales, and transformation of waters flowing from the Pacific and the marginal seas through the Indian Ocean (IO) on their way toward the South Atlantic within a high-resolution ocean/sea-ice model. The model fields are analyzed from a Lagrangian perspective where water volumes are tracked as they enter the IO. The IO contributes 12.6 Sv to Agulhas leakage, which within the model is 14.1 ± 2.2 Sv, the rest originates from the South Atlantic. The Indonesian Through-flow constitutes about half of the IO contribution, is surface bound, cools and salinificates as it leaves the basin within 10–30 years. Waters entering the IO south of Australia are at intermediate depths and maintain their temperature-salinity properties as they exit the basin within 15–35 years. Of these waters, the contribution from Tasman leakage is 1.4 Sv. The rest stem from recirculation from the frontal regions of the Southern Ocean. The marginal seas export 1.0 Sv into the Atlantic within 15–40 years, and the waters cool and freshen on-route. However, the model's simulation of waters from the Gulfs of Aden and Oman are too light and hence overly influenced by upper ocean circulations. In the Cape Basin, Agulhas leakage is well mixed. On-route, temperature-salinity transformations occur predominantly in the Arabian Sea and within the greater Agulhas Current region. Overall, the IO exports at least 7.9 Sv from the Pacific to the Atlantic, thereby quantifying the strength of the upper cell of the global conveyor belt.

Description: The North Atlantic Current (NAC) is subject to variability on multiannual to decadal time scales, influencing the transport of volume, heat, and freshwater from the subtropical to the eastern subpolar North Atlantic (NA). Current observational time series are either too short or too episodic to study the processes involved. Here we compare the observed continuous NAC transport time series at the western flank of the Mid-Atlantic Ridge (MAR) and repeat hydrographic measurements at the OVIDE line in the eastern Atlantic with the NAC transport and circulation in the high-resolution (1/20°) ocean model configuration VIKING20 (1960–2008). The modeled baroclinic NAC transport relative to 3400 m (24.5 ± 7.1 Sv) at the MAR is only slightly lower than the observed baroclinic mean of 27.4 ± 4.7 Sv from 1993 to 2008, and extends further north by about 0.5°. In the eastern Atlantic, the western NAC (WNAC) carries the bulk of the transport in the model, while transport estimates based on hydrographic measurements from five repeated sections point to a preference for the eastern NAC (ENAC). The model is able to simulate the main features of the subpolar NA, providing confidence to use the model output to analyze the influence of the North Atlantic Oscillation (NAO). Model based velocity composites reveal an enhanced NAC transport across the MAR of up to 6.7 Sv during positive NAO phases. Most of that signal (5.4 Sv) is added to the ENAC transport, while the transport of the WNAC was independent of the NAO.

Description: The North Brazil Current (NBC) constitutes a bottleneck for the mean northward return flow of the Atlantic Meridional Overturning Circulation (AMOC) in the tropical South Atlantic. Previous studies suggested a link between interannual to multidecadal NBC and AMOC transport variability and proposed to use NBC observations as an index for the AMOC. Here we use a set of hindcast, sensitivity, and perturbation experiments performed within a hierarchy of ocean general circulation models to show that decadal to multidecadal buoyancy-forced changes in the basin-scale AMOC transport indeed manifest themselves in the NBC. The relation is, however, masked by a strong interannual to decadal wind-driven gyre variability of the NBC. While questioning the NBC transport as a direct index for the AMOC, the results support its potential merit for an AMOC monitoring system, provided that the wind-driven circulation variability is properly accounted for.

Description: The impact of the positive and negative phases of the Pacific Decadal Oscillation (PDO) on the extension of the poorly oxygenated regions of the eastern Pacific Ocean was assessed using a coupled ocean circulation–biogeochemical model. We show that during a “typical” PDO-positive phase the volume of the suboxic regions expands by 7 % over 50 years due to a slowdown of the large-scale circulation related to the decrease in the intensity of the trade winds. Changes in oxygen levels are mostly controlled by advective processes between 10∘ N and 10∘ S, whereas diffusive processes are dominant poleward of 10∘: in a “typical” PDO-positive phase the sluggish equatorial current system provides less oxygen to the eastern equatorial part of the basin while the oxygen transport by diffusive processes significantly decreases south of 10∘ S. The suboxic region located north of 10∘ N displays less sensitivity to the phase of the PDO as the local upwelling-related processes play a dominant role compared to the large-scale circulation in setting the oxygen concentration. Our study suggests that the prevailing PDO-positive conditions since 1975 may explain a significant part of the current deoxygenation occurring in the eastern Pacific Ocean.

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.