Description: Decadal variabilities in Indian Ocean subsurface ocean heat content (OHC; 50–300 m) since the 1950s are examined using ocean reanalyses. This study elaborates on how Pacific variability modulates the Indian Ocean on decadal time scales through both oceanic and atmospheric pathways. High correlations between OHC and thermocline depth variations across the entire Indian Ocean Basin suggest that OHC variability is primarily driven by thermocline fluctuations. The spatial pattern of the leading mode of decadal Indian Ocean OHC variability closely matches the regression pattern of OHC on the interdecadal Pacific oscillation (IPO), emphasizing the role of the Pacific Ocean in determining Indian Ocean OHC decadal variability. Further analyses identify different mechanisms by which the Pacific influences the eastern and western Indian Ocean. IPO-related anomalies from the Pacific propagate mainly through oceanic pathways in the Maritime Continent to impact the eastern Indian Ocean. By contrast, in the western Indian Ocean, the IPO induces wind-driven Ekman pumping in the central Indian Ocean via the atmospheric bridge, which in turn modifies conditions in the southwestern Indian Ocean via westward-propagating Rossby waves. To confirm this, a linear Rossby wave model is forced with wind stresses and eastern boundary conditions based on reanalyses. This linear model skillfully reproduces observed sea surface height anomalies and highlights both the oceanic connection in the eastern Indian Ocean and the role of wind-driven Ekman pumping in the west. These findings are also reproduced by OGCM hindcast experiments forced by interannual atmospheric boundary conditions applied only over the Pacific and Indian Oceans, respectively.

Description: Benthic storms are important for both the energy budget of the ocean and for sediment resuspension and transport. Using 30 years of output from a high-resolution model of the North Atlantic, it is found that most of the benthic storms in the model occur near the western boundary in association with the Gulf Stream and the North Atlantic Current, in regions that are generally co-located with the peak near-bottom eddy kinetic energy. A common feature are meander troughs in the near-surface jets that are accompanied by deep low pressure anomalies spinning up deep cyclones with near-bottom velocities of up to more than 0.5 m/s. A case study of one of these events shows the importance of both baroclinic and barotropic instability of the jet, with energy being extracted from the jet in the upstream part of the meander trough and partly returned to the jet in the downstream part of the meander trough. This motivates examining the 30-year time mean of the energy transfer from the (annual mean) background flow into the eddy kinetic energy. This quantity is shown to be co-located well with the region in which benthic storms and large increases in deep cyclonic relative vorticity occur most frequently, suggesting an important role for mixed barotropic-baroclinic instability driven cyclogenesis in generating benthic storms throughout the model simulation. Regions of largest energy transfer and most frequent benthic storms are found to be the Gulf Stream west of the New England Seamounts and the North Atlantic Current near Flemish Cap.

Description: We provide a time series of Agulhas leakage anomalies over the last 20-years from satellite altimetry. Until now, measuring the interannual variability of Indo-Atlantic exchange has been the major barrier in the investigation of the dynamics and large scale impact of Agulhas leakage. We compute the difference of transport between the Agulhas Current and Agulhas Return Current, which allows us to deduce Agulhas leakage. The main difficulty is to separate the Agulhas Return Current from the southern limb of the subtropical "supergyre" south of Africa. For this purpose, an algorithm that uses absolute dynamic topography data is developed. The algorithm is applied to a state-of-the-art ocean model. The comparison with a Lagrangian method to measure the leakage allows us to validate the new method. An important result is that it is possible to measure Agulhas leakage in this model using the velocity field along a section that crosses both the Agulhas Current and the Agulhas Return Current. In the model a good correlation is found between measuring leakage using the full depth velocities and using only the surface geostrophic velocities. This allows us to extend the method to along-track absolute dynamic topography from satellites. It is shown that the accuracy of the mean dynamic topography does not allow us to determine the mean leakage but that leakage anomalies can be accurately computed.

Description: The Agulhas Current plays a crucial role in the thermohaline circulation through its leakage into the South Atlantic. Under both past and present climates, the trade winds and westerlies could have the ability to modulate the amount of Indian-Atlantic inflow. Compelling arguments have been put forward suggesting that trade winds alone have little impact on the magnitude of Agulhas leakage. Here, employing three ocean models for robust analysis – a global coarse resolution, a regional eddy-permitting and a nested high-resolution eddy-resolving configuration – and systematically altering the position and intensity of the westerly wind belt in a series of sensitivity experiments, it is shown that the westerlies, in particular their intensity, control the leakage. Leakage responds proportionally to the westerlies intensity up to a certain point. Beyond this, through the adjustment of the large-scale circulation, energetic interactions occur between the Agulhas Return Current and the Antarctic Circumpolar Current that result in a state where leakage no longer increases. This adjustment takes place within 1 to 2 decades. Contrary to previous assertions, our results further show that an equatorward (poleward) shift in westerlies increases (decreases) leakage. This occurs due to the redistribution of momentum input by the winds. It is concluded that the reported present-day leakage increase could therefore reflect an unadjusted oceanic response mainly to the strengthening westerlies over the last few decades.

Description: Some studies of ocean climate model experiments suggest that regional changes in dynamic sea level could provide a valuable indicator of trends in the strength of the Atlantic meridional overturning circulation (MOC). This paper describes the use of a sequence of global ocean–ice model experiments to show that the diagnosed patterns of sea surface height (SSH) anomalies associated with changes in the MOC in the North Atlantic (NA) depend critically on the time scales of interest. Model hindcast simulations for 1958–2004 reproduce the observed pattern of SSH variability with extrema occurring along the Gulf Stream (GS) and in the subpolar gyre (SPG), but they also show that the pattern is primarily related to the wind-driven variability of MOC and gyre circulation on interannual time scales; it is reflected also in the leading EOF of SSH variability over the NA Ocean, as described in previous studies. The pattern, however, is not useful as a “fingerprint” of longer-term changes in the MOC: as shown with a companion experiment, a multidecadal, gradual decline in the MOC [of 5 Sv (1 Sv ≡ 106 m3 s−1) over 5 decades] induces a much broader, basin-scale SSH rise over the mid-to-high-latitude NA, with amplitudes of 20 cm. The detectability of such a trend is low along the GS since low-frequency SSH changes are effectively masked here by strong variability on shorter time scales. More favorable signal-to-noise ratios are found in the SPG and the eastern NA, where a MOC trend of 0.1 Sv yr−1 would leave a significant imprint in SSH already after about 20 years.

Description: The relationship between the Agulhas Current and the Agulhas leakage is not well understood. Here, this is investigated using two basin-scale and two global ocean models, of incrementally increasing resolution. The response of the Agulhas Current is evaluated under a series of sensitivity experiments, in which idealised anomalies, designed to geometrically modulate zonal trade wind stress, are applied across the Indian Ocean basin. The imposed wind stress changes exceed ±2 standard deviations from the annual mean trade winds and, in the case of intensification, are partially representative of recently observed trends. The Agulhas leakage is quantified using complimentary techniques based on Lagrangian virtual floats and Eulerian passive tracer flux. As resolution increases, model behavior converges and the sensitivity of the leakage to Agulhas Current transport anomalies is reduced. In the two eddy-resolving configurations tested, the leakage is insensitive to changes in Agulhas Current transport at 32°S, though substantial eddy kinetic energy anomalies are evident. Consistent with observations, the position of the retroflection remains stable. The decoupling of Agulhas Current variability from the Agulhas leakage suggests that, while correlations between the two may exist, they may not have a clear dynamical basis. It is suggested that present and future Agulhas leakage proxies be considered in the context of potentially transient forcing regimes.

Description: The Indian Ocean has sustained robust surface warming in recent decades, but the role of multi-decadal variability remains unclear. Using ocean model hindcasts, characteristics of low-frequency Indian Ocean temperature variations are explored. Simulated upper-ocean temperature changes across the Indian Ocean in the hindcast are consistent with those recorded in observational products and ocean reanalyses. Indian Ocean temperatures exhibit strong warming trends since the 1950s limited to the surface and south of 30°S, while extensive subsurface cooling occurs over much of the tropical Indian Ocean. Previous work focused on diagnosing causes of these long-term trends in the Indian Ocean over the second half of the 20th Century. Instead, the temporal evolution of Indian Ocean subsurface heat content is shown here to reveal distinct multi-decadal variations associated with the Pacific Decadal Oscillation and the long-term trends are thus interpreted to result from aliasing of the low-frequency variability. Transmission of the multi-decadal signal occurs via an oceanic pathway through the Indonesian Throughflow and is manifest across the Indian Ocean centered along 12°S as westward propagating Rossby waves modulating thermocline and subsurface heat content variations. Resulting low-frequency changes in the eastern Indian Ocean thermocline depth are associated with decadal variations in the frequency of Indian Ocean Dipole (IOD) events, with positive IOD events unusually common in the 1960s and 1990s with a relatively shallow thermocline. In contrast, the deeper thermocline depth in the 1970s and 1980s is associated with frequent negative IOD and rare positive IOD events. Changes in Pacific wind forcing in recent decades and associated rapid increases in Indian Ocean subsurface heat content can thus affect the basin’s leading mode of variability, with implications for regional climate and vulnerable societies in surrounding countries.

Description: The upper tropical Atlantic Ocean has markedly warmed since the 1960s. It has been shown that this warming was not due to local heat fluxes, and that the trade winds that drive the coastal and equatorial upwelling have intensified rather than weakened. Remote forcing might thus have played an important role. Here model experiments are used to investigate the contribution from an increased inflow of warm Indian Ocean water through Agulhas leakage. A high-resolution hindcast experiment with interannually varying forcing for the time period 1948 to 2007, in which Agulhas leakage increases by about 45% from the 1960s to the early 2000s, reproduces the observed warming trend. To tease out the role of Agulhas leakage, a sensitivity experiment designed to only increase Agulhas leakage is used. Compared to a control simulation it shows a pronounced warming in the upper tropical Atlantic Ocean. A Lagrangian trajectory analysis confirms that a significant portion of Agulhas leakage water reaches the upper 300m of the tropical Atlantic Ocean within two decades, and that the tropical Atlantic warming in the sensitivity experiment is mainly due to water of Agulhas origin. Therefore, it is suggested that the increased trade winds since the 1960s favor upwelling of warmer subsurface waters, which in parts originate from the Agulhas, leading to higher SSTs in the tropics

Description: A new Earth system model, the Flexible Ocean and Climate Infrastructure (FOCI), is introduced. A first version of FOCI consists of a global high-top atmosphere (ECHAM6.3) and an ocean model (NEMO3.6) as well as sea ice (LIM2) and land surface model components (JSBACH), which are coupled through the OASIS3-MCT software package. FOCI includes a number of optional modules which can be activated depending on the scientific question of interest. In the atmosphere, interactive stratospheric chemistry can be used (ECHAM6-HAMMOZ) to study, for example, the effects of the ozone hole on the climate system. In the ocean, a biogeochemistry model (MOPS) is available to study the global carbon cycle. A unique feature of FOCI is the ability to explicitly resolve mesoscale ocean eddies in specific regions. This is realized in the ocean through nesting; first examples for the Agulhas Current and the Gulf Stream systems are described here. FOCI therefore bridges the gap between coarse-resolution climate models and global high-resolution weather prediction and ocean-only models. It allows to study the evolution of the climate system on regional and seasonal to (multi-) decadal scales. The development of FOCI resulted from a combination of the long-standing expertise in ocean and climate modeling in several research units and divisions at GEOMAR. FOCI will thus be used to complement and interpret long-term observations in the Atlantic, enhance the process understanding of the role of mesoscale oceanic eddies for large-scale oceanic and atmospheric circulation patterns, study feedback mechanisms with stratospheric processes, estimate future ocean acidification, improve the simulation of the Atlantic Meridional Overturning Circulation changes and their influence on climate, ocean chemistry and biology. In this paper we present both the scientific vision for the development of FOCI as well as some technical details. This includes a first validation of the different model components using several configurations of FOCI. Results show that the model in its basic configuration runs stably under pre-industrial control as well as under historical forcing, and produces a mean climate and variability which compares well with observations, reanalysis products and other climate models. The nested configurations reduce some long-standing biases in climate models and are an important step forward to include the atmospheric response in multi-decadal eddy-rich configurations.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.