ALBERT

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 a 10∘ latitude in both hemispheres and also relevant in the very east of the tropical domain. Buoyancy-flux-forced anomalies are correlated with variations driven by wind stress changes, but their effect on the prevailing anomalies and the importance of heat and freshwater fluxes vary locally. In the eastern tropical basin, interannual sea level variability is amplified by anomalous heat fluxes, while the importance of freshwater 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 El Niño–Southern Oscillation (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.

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 process-level 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 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: 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 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.

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 onto 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 sub-mesoscale 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.

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 has been 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 % in 50 years due to a slow-down of the large scale circulation related with the decrease of the intensity of the trade winds. The oxygen levels are mostly constrained by advective processes between 10° N and 10° S while the diffusive processes are dominant poleward of 10°: in a typical PDO positive phase the sluggish equatorial current system provides less oxygen into the eastern equatorial part of the basin while the oxygen transport by diffusive processes significantly decreases south of 10° S. The 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: 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: A hierarchy of global 1/4∘ (ORCA025) and Atlantic Ocean 1/20∘ nested (VIKING20X) ocean–sea-ice models is described. It is shown that the eddy-rich configurations performed in hindcasts of the past 50–60 years under CORE and JRA55-do atmospheric forcings realistically simulate the large-scale horizontal circulation, the distribution of the mesoscale, overflow and convective processes, and the representation of regional current systems in the North and South Atlantic. The representation of the Atlantic Meridional Overturning Circulation (AMOC), and in particular the long-term temporal evolution, strongly depends on numerical choices for the application of freshwater fluxes. The interannual variability of the AMOC instead is highly correlated among the model experiments and also with observations, including the 2010 minimum observed by RAPID at 26.5∘ N. This points to a dominant role of the wind forcing. The ability of the model to represent regional observations in western boundary current (WBC) systems at 53∘ N, 26.5∘ N and 11∘ S is explored. The question is investigated of whether WBC systems are able to represent the AMOC, and in particular whether these WBC systems exhibit similar temporal evolution to that of the zonally integrated AMOC. Apart from the basin-scale measurements at 26.5∘ N, it is shown that in particular the outflow of North Atlantic Deepwater at 53∘ N is a good indicator of the subpolar AMOC trend during the recent decades, once provided in density coordinates. The good reproduction of observed AMOC and WBC trends in the most reasonable simulations indicate that the eddy-rich VIKING20X is capable of representing realistic forcing-related and ocean-intrinsic trends.

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.