Contribution of buoyancy fluxes to tropical Pacific sea level variability

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

Late 20th Century Indian Ocean Heat Content Gain Masked by Wind Forcing

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

Regional Imprints of Changes in the Atlantic Meridional Overturning Circulation in the Eddy-rich Ocean Model VIKING20X

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, and in particular the long-term temporal evolution, of the Atlantic Meridional Overturning Circulation (AMOC) 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 pointing at a dominant role of the forcing. Regional observations in western boundary current systems at 53° N, 26.5° N and 11° S are explored in respect to their ability to represent the AMOC and to monitor the temporal evolution of the 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, if the latter is 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 in representing realistic forcing-related and ocean-intrinsic trends

Mesoscale perturbations control inter-ocean exchange south of Africa

The quantification of inter-ocean leakage from the South Indian to the South Atlantic Ocean is an important measure for the role of the Agulhas system in the global thermohaline circulation. To explore the specific role of mesoscale variability (such as Agulhas rings and Mozambique eddies) in this process a high-resolution model (based on NEMO-ORCA) for the Agulhas region has been set up. It is nested into a global coarse-resolution model. The high-resolution nest captures all salient features of the greater Agulhas region, including the upstream perturbations of the Agulhas Current and Natal Pulses along the African coast. A comparison of the inter-ocean exchange in the high-resolution nest with its coarse resolution counterpart reveals that the latter significantly over-estimates the amount of water flowing into the Atlantic Ocean, demonstrating the need to explicitly simulate the mesoscale features. A sensitivity experiment that excludes the upstream perturbations revealed no difference in the amount of inter-ocean exchange. However, the realistic representation of Agulhas rings and their drift path into the South Atlantic depends on the simulation of those upstream perturbations

New Antarctic gravity anomaly grid for enhanced geodetic and geophysical studies in Antarctica

Gravity surveying is challenging in Antarctica because of its hostile environment and inaccessibility. Nevertheless, many ground-based, airborne, and shipborne gravity campaigns have been completed by the geophysical and geodetic communities since the 1980s. We present the first modern Antarctic-wide gravity data compilation derived from 13 million data points covering an area of 10 million km2, which corresponds to 73% coverage of the continent. The remove-compute-restore technique was applied for gridding, which facilitated leveling of the different gravity data sets with respect to an Earth gravity model derived from satellite data alone. The resulting free-air and Bouguer gravity anomaly grids of 10 km resolution are publicly available. These grids will enable new high-resolution combined Earth gravity models to be derived and represent a major step forward toward solving the geodetic polar data gap problem. They provide a new tool to investigate continental-scale lithospheric structure and geological evolution of Antarctica

Seasonal variability of eddy kinetic energy in a global high-resolution ocean model

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

The INALT family – a set of high-resolution nests for the Agulhas Current system within global NEMO ocean/sea-ice configurations

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

Oceanic hindcast simulations at high resolution suggest that the Atlantic MOC is bistable

All climate models predict a freshening of the North Atlantic at high latitude that may induce an abrupt change of the Atlantic Meridional Overturning Circulation (hereafter AMOC) if it resides in the bistable regime, where both a strong and a weak state coexist. The latter remains uncertain as there is no consensus among observations and ocean reanalyses, where the AMOC is bistable, versus most climate models that reproduce a mono-stable strong AMOC. A series of four hindcast simulations of the global ocean at 1/12° resolution, which is presently unique, are used to diagnose freshwater transport by the AMOC in the South Atlantic, an indicator of AMOC bistability. In all simulations, the AMOC resides in the bistable regime: it exports freshwater southward in the South Atlantic, implying a positive salt advection feedback that would act to amplify a decreasing trend in subarctic deep water formation as projected in climate scenarios

Evaluation of global ocean–sea-ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2)

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