Recent studies, using data from the OSNAP observational campaign and from numerical ocean models, suggest that surface buoyancy losses over the Iceland Basin and the Irminger Sea may, in contradiction to the established consensus, be more significant than those over the Labrador Sea, and that these former regions are in fact the dominant sites for formation of upper North Atlantic Deep Water), with the Labrador Sea acting mainly as a region of further densification as the dense waters flow around the gyre. Here we present a set of hindcast integrations of a global 1/4° NEMO ocean configuration from 1958 until nearly the present day, forced with three standard surface forcing datasets. We use the surface-forced streamfunction, estimated from surface buoyancy fluxes, along with the overturning streamfunction, similarly defined in potential density space, to investigate the causal link between surface forcing and decadal variability in the strength of the Atlantic meridional overturning circulation (AMOC). A scalar metric based on the surface forced streamfunction, evaluated in critical density and latitude classes, and accumulated in time, is found to be a good predictor of changes in the overturning strength, and the surface heat loss from the Irminger Sea is confirmed to be the dominant mechanism for decadal AMOC variability. We use the streamfunctions to demonstrate that the watermasses in the simulations are transformed to higher densities as they propagate around the subpolar gyre from their formation locations in the north-east Atlantic and the Irminger Sea, consistent with the picture emerging from observations.

Large amplitude oscillations in the meridional overturning circulation (MOC) have been found near the equator in all major ocean basins in the NEMO ocean general circulation model. With periods of 3-15 days and amplitudes of ~±100 Sv in the Pacific, these oscillations have been shown to correspond to zonally integrated equatorially trapped waves forced by winds within 10° N/S of the equator. Observations of dynamic height from the Tropical Atmosphere Ocean (TAO) mooring array in the equatorial Pacific also exhibit spectral peaks consistent with the dispersion relation for equatorially trapped waves. Here, we revisit the TAO observations to confirm that the amplitude of the oscillations is consistent with the simulations, supporting the modelled large amplitude MOC oscillations. We also show that the zonal structure of the frequency spectrum in both observations and simulations is predicted by changes in the baroclinic wave speed with variation in stratification across the ocean basin.