Description: A set of experiments utilizing different implementations of the global ORCA-LIM model with horizontal resolutions of 2°, 0.5° and 0.25° is used to investigate tropical and extra-tropical influences on equatorial Pacific SST variability at interannual to decadal time scales. The model experiments use a bulk forcing methodology building on the global forcing data set for 1958 to 2000 developed by Large and Yeager (2004) that is based on a blend of atmospheric reanalysis data and satellite products. Whereas representation of the mean structure and transports of the (sub-) tropical Pacific current fields is much improved with the enhanced horizontal resolution, there is only little difference in the simulation of the interannual variability in the equatorial regime between the 0.5° and 0.25° model versions, with both solutions capturing the observed SST variability in the Niño3-region. The question of remotely forced oceanic contributions to the equatorial variability, in particular, the role of low-frequency changes in the transports of the Subtropical Cells (STCs), is addressed by a sequence of perturbation experiments using different combinations of fluxes. The solutions show the near-surface temperature variability to be governed by wind-driven changes in the Equatorial Undercurrent. The relative contributions of equatorial and off-equatorial atmospheric forcing differ between interannual and longer, (multi-) decadal timescales: for the latter there is a significant impact of changes in the equatorward transport of subtropical thermocline water associated with the lower branches of the STCs, related to variations in the off-equatorial trade winds. A conspicuous feature of the STC variability is that the equatorward transports in the interior and along the western boundary partially compensate each other at both decadal and interannual time scales, with the strongest transport extrema occurring during El Niño episodes. The behaviour is rationalized in terms of a wobbling in the poleward extents of the tropical gyres, which is manifested also in a meridional shifting of the bifurcation latitudes of the North and South Equatorial Current systems.

Description: Predicting the evolution of climate over decadal timescales requires a quantitative understanding of the dynamics that govern the meridional overturning circulation (MOC)1. Comprehensive ocean measurement programmes aiming to monitor MOC variations have been established in the subtropical North Atlantic2, 3 (RAPID, at latitude 26.5° N, and MOVE, at latitude 16° N) and show strong variability on intraseasonal to interannual timescales. Observational evidence of longer-term changes in MOC transport remains scarce, owing to infrequent sampling of transoceanic sections over past decades4, 5. Inferences based on long-term sea surface temperature records, however, supported by model simulations, suggest a variability with an amplitude of plusminus1.5–3 Sv (1 Sv = 106 m3 s-1) on decadal timescales in the subtropics6. Such variability has been attributed to variations of deep water formation in the sub-arctic Atlantic, particularly the renewal rate of Labrador Sea Water7. Here we present results from a model simulation that suggest an additional influence on decadal MOC variability having a Southern Hemisphere origin: dynamic signals originating in the Agulhas leakage region at the southern tip of Africa. These contribute a MOC signal in the tropical and subtropical North Atlantic that is of the same order of magnitude as the northern source. A complete rationalization of observed MOC changes therefore also requires consideration of signals arriving from the south.

Description: Large amounts of the greenhouse gas methane are released from the seabed to the water column1, where it may be consumed by aerobic methanotrophic bacteria2. The size and activity of methanotrophic communities, which determine the amount of methane consumed in the water column, are thought to be mainly controlled by nutrient and redox dynamics3–7. Here, we report repeated measurements of methanotrophic activity and community size at methane seeps west of Svalbard, and relate them to physical water mass properties and modelled ocean currents. We show that cold bottom water, which contained a large number of aerobic methanotrophs, was displaced by warmer water with a considerably smaller methanotrophic community within days. Ocean current simulations using a global ocean/sea-ice model suggest that this water mass exchange is consistent with short-term variations in the meandering West Spitsbergen Current. We conclude that the shift from an offshore to a nearshore position of the current can rapidly and severely reduce methanotrophic activity in the water column. Strong fluctuating currents are common at many methane seep systems globally, and we suggest that they affect methane oxidation in the water column at other sites, too.

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: 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: Observations show a significant intensification of the Southern Hemisphere westerlies, the prevailing winds between the latitudes of 30° and 60° S, over the past decades. A continuation of this intensification trend is projected by climate scenarios for the twenty-first century. The response of the Antarctic Circumpolar Current and the carbon sink in the Southern Ocean to changes in wind stress and surface buoyancy fluxes is under debate. Here we analyse the Argo network of profiling floats and historical oceanographic data to detect coherent hemispheric-scale warming and freshening trends that extend to depths of more than 1,000 m. The warming and freshening is partly related to changes in the properties of the water masses that make up the Antarctic Circumpolar Current, which are consistent with the anthropogenic changes in heat and freshwater fluxes suggested by climate models. However, we detect no increase in the tilt of the surfaces of equal density across the Antarctic Circumpolar Current, in contrast to coarse-resolution model studies. Our results imply that the transport in the Antarctic Circumpolar Current and meridional overturning in the Southern Ocean are insensitive to decadal changes in wind stress.

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.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence on the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding “nests” covering the South Atlantic and the western Indian oceans at 1/10∘ (INALT10) and 1/20∘ (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60∘ resolution (INALT60) has been developed to gain insights into submesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations, indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both the southward-flowing Agulhas Current and the northward-flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.