Articles | Volume 17, issue 5
Research article
 | Highlight paper
28 Sep 2021
Research article | Highlight paper |  | 28 Sep 2021

A dynamically based method for estimating the Atlantic meridional overturning circulation at 26° N from satellite altimetry

Alejandra Sanchez-Franks, Eleanor Frajka-Williams, Ben I. Moat, and David A. Smeed

The large-scale system of ocean currents that transport warm waters in the upper 1000 m northward and return deeper cooler waters southward is known as the Atlantic meridional overturning circulation (AMOC). Variations in the AMOC have significant repercussions for the climate system; hence, there is a need for long-term monitoring of AMOC fluctuations. Currently the longest record of continuous directly measured AMOC changes is from the RAPID-MOCHA-WBTS programme, initiated in 2004. The RAPID programme and other mooring programmes have revolutionised our understanding of large-scale circulation; however, by design they are constrained to measurements at a single latitude and cannot tell us anything pre-2004.

Nearly global coverage of surface ocean data from satellite altimetry has been available since the launch of the TOPEX/Poseidon satellite in 1992 and has been shown to provide reliable estimates of surface ocean transports on interannual timescales including previous studies that have investigated empirical correlations between sea surface height variability and the overturning circulation. Here we show a direct calculation of ocean circulation from satellite altimetry of the upper mid-ocean transport (UMO), the Gulf Stream transport through the Florida Straits (GS), and the AMOC using a dynamically based method that combines geostrophy with a time mean of the vertical structure of the flow from the 26 N RAPID moorings. The satellite-based transport captures 56 %, 49 %, and 69 % of the UMO, GS, and AMOC transport variability, respectively, from the 26 N RAPID array on interannual (18-month) timescales. Further investigation into the vertical structure of the horizontal transport shows that the first baroclinic mode accounts for 83 % of the interior geostrophic variability, and the combined barotropic and first baroclinic mode representation of dynamic height accounts for 98 % of the variability. Finally, the methods developed here are used to reconstruct the UMO and the AMOC for the time period pre-dating RAPID, 1993 to 2003. The effective implementation of satellite-based method for monitoring the AMOC at 26 N lays down the starting point for monitoring large-scale circulation at all latitudes.

1 Introduction

The Atlantic meridional overturning circulation (AMOC) is the large-scale oceanic circulation comprised of currents that carry warm, shallow water northward and return cold deep water southward. Variations in the AMOC’s strength have a significant impact on the Earth’s climate system (Srokosz et al.2012) as the AMOC's northward transport of warm surface waters releases heat on the order of 1 PW over the North Atlantic (Trenberth and Caron2001), which is key in maintaining the relatively mild winter climate of northwestern Europe (Hall and Bryden1982; Pohlmann et al.2006). On a decadal timescale, the AMOC is identified as the underlying driver of latitudinal shifts in the Gulf Stream path (Sanchez-Franks and Zhang2015), and changes in the AMOC are a key driver of Atlantic multidecadal variability and associated climate variability (Zhang et al.2019; Moat et al.2019). On longer timescales, coupled climate models predict significant weakening of the AMOC, which results in a reduced northward heat transport associated with increased greenhouse gas forcing (Caesar et al.2018).

One of the existing programmes monitoring the AMOC is the RAPID-MOCHA-WBTS (RAPID – Meridional Overturning Circulation and Heatflux Array – Western Boundary Time Series; hereafter RAPID) 26 N mooring array, which has been measuring the meridional volume and heat transport at 26 N since 2004 (Rayner et al.2011; Johns et al.2011; Smeed et al.2014; McCarthy et al.2015; Moat et al.2020a, b). The RAPID programme is an international campaign that uses an array of moored instruments at strategic locations in order to continuously estimate AMOC fluctuations (further details in McCarthy et al.2015). Key to the RAPID strategy is the fact that the interior of the ocean is largely in geostrophic balance; placing moorings in the western and eastern end points of the Atlantic basin, as well as both flanks of the mid-Atlantic ridge, allows fluctuations in the vertical density profiles to be simultaneously measured at the mooring locations. By taking the difference between the end point measurements, the basin-wide interior flow can be estimated (e.g. Cunningham et al.2007; Kanzow et al.2010; Rayner et al.2011). The vertical density profiles estimated from the RAPID moorings are affected by long Rossby waves that start off as westward-travelling perturbations (Johnson and Marshall2002), altering the east–west density gradient driving AMOC variability on seasonal and subannual timescales (Hirschi et al.2007).

Other international efforts to measure the AMOC via mooring array programmes include SAMBA (South Atlantic MOC Basin-wide Array) at 34.5 S (Meinen et al.2018; Kersale et al.2020), TRACOS (Tropical Atlantic Circulation and Overturning) at 11 S (Hummels et al.2015; Herrford et al.2021), and OSNAP (Overturning in the Subpolar North Atlantic Program) in the subpolar North Atlantic (Lozier et al.2019). RAPID and these other mooring array programmes have made step-change advancements in our understanding of the AMOC over the last 15 years (Frajka-Williams et al.2019); however, they are limited to a single line of latitude and alone cannot be used to infer upstream–downstream changes in the larger-scale structure of the AMOC. Further, the large meridional distances between the arrays (e.g. 34.5 S, 11 S, 26 N, 58 N), combined with the relatively short length and/or life span of the programmes (e.g. the oldest is RAPID, initiated in 2004), preclude our ability to understand AMOC connectivity, i.e. its meridional coherence, which is still poorly understood (e.g. Bingham et al.2007).

Satellite data provide high spatial coverage of the global oceans. Satellite altimetry measures the surface topography of the ocean (i.e. sea surface height; SSH), which is an effective tool for measuring large-scale surface circulation and associated dynamics. In particular, westward-propagating Rossby waves, or surface-intensified flow, which is largely of baroclinic structure, are captured by satellite observations (Chelton et al.2007). Previous efforts to use Earth-observing methods have shown capacity to estimate the large-scale ocean circulation. For example, empirical orthogonal function (EOF) analysis of satellite sea level anomaly (SLA) has suggested that altimetry could be a useful indicator for AMOC transports along 40–50 N on interannual timescales (Bingham and Hughes2009). Frajka‐Williams (2015) found that the relationship between satellite SLA and RAPID dynamic height was robust enough to be used to create a proxy for the upper mid-ocean transport, which when combined with Ekman transport (derived from reanalysis products) and Gulf Stream transport (from a submarine telephone cable) accounted for more than 90 % of MOC variability on interannual timescales. The method developed in Frajka‐Williams (2015) was purely statistical. The relationship between dynamic height from RAPID moorings and satellite SSH at 26 N was also found to be robust, accounting for up to 72 % of variance at a station 500 km offshore (Kanzow et al.2009); however, this relationship was found to deteriorate in proximity to the western boundary (Bryden et al.2009; Clément et al.2014).

One of the limitations of satellite altimetry is that it only provides information about the sea surface. The information captured by satellite altimetry at the surface is representative of variability associated with the first baroclinic mode (Hirschi et al.2007). Because of the baroclinic nature of the surface-intensified ocean, Rossby wave theory suggests that some MOC variability can be linked to satellite altimetry on subannual (Hirschi et al.2007) and subannual to interannual (Hirschi et al.2009) timescales. However, satellite altimetry cannot tell us anything about the vertical shape of the horizontal velocity, which forms part of the AMOC structure. To get around this, previous studies have combined satellite altimetry data with vertical structure of the flow derived from in situ datasets in the North Atlantic at 26 N (e.g. Hirschi et al.2007, 2009; Kanzow et al.2009; Szuts et al.2012). In particular, Hirschi et al. (2009) found that combining SSH with baroclinic structure using a numerical model yielded a robust estimate of eastern and western branches of the MOC, but results could not be replicated for reconstructing variability of the basin-wide MOC on subannual to interannual timescales (Hirschi et al.2007, 2009). They found that even small errors in the western and eastern components could overwhelm MOC variability, and only highly accurate estimates of the western and eastern branches would yield a reliable reconstruction of the basin-wide MOC. In a study focused on the impact of eddy dynamics on basin-wide transport, Kanzow et al. (2009) tested whether geostrophic transport could be obtained by combining geostrophic surface flow from SSH with the dominant vertical mode of horizontal velocity; they found that the correlation between the difference in SSH at eastern and western end points of the basin had little skill in capturing the upper mid-ocean transport and suggested that the reason for this was not enough satellite data close to the shelf as well as changes in the vertical structure of the flow moving toward the western boundary. Szuts et al. (2012) similarly found that SSH captured fluctuations in the centre of the basin but struggled to capture the variability closer to the boundaries, where the vertical structure of the flow was not as well described by the dominant baroclinic mode.

The length of the RAPID records has more than doubled since the aforementioned studies that combined geostrophy and vertical structure of horizontal velocity (e.g. Hirschi et al.2009; Kanzow et al.2009; Szuts et al.2012), so the time is right to re-examine the issues (e.g. problems at the boundary, including contribution from vertical modes) and build on the methods developed in these previous studies. Thus, here we combine geostrophic principles with the vertical structure of the horizontal flow to develop a satellite-based MOC transport on interannual timescales. In contrast to earlier efforts by Hirschi et al. (2009) and Kanzow et al. (2009), the longer records now available from RAPID enable us to test the methods for longer timescales (interannual rather than subannual), and we find that the skill is generally higher at interannual timescales than at shorter timescales. Unlike Frajka‐Williams (2015), this method now uses geostrophy (east minus west differences in sea level anomaly) rather than finding the point location where the sea level anomaly has the strongest correlation with the transport variability. This change in the method is a prerequisite to developing similar methods at other latitudes as it reduces (though it does not entirely eliminate) the requirement for in situ data to “train” the method. In addition, while Szuts et al. (2012) and Clément et al. (2014) found that multiple modes are required to explain the variance in the western boundary profiles, we show that on interannual timescales, the contributions from higher modes are reduced and the first mode explains a majority of the variance. In the following section (Sect. 2), a brief overview of data and methods used here is presented. Then an evaluation of the satellite data and parameters from the RAPID moorings follows in Sect. 3. Rossby wave theory and an in-depth analysis of horizontal velocity normal modes from RAPID mooring data are shown in Sect. 4. Sections 5 and 6 investigate construction of the upper mid-ocean, the Gulf Stream, and the AMOC transports from satellite altimetry. Finally, a summary and conclusions are given in Sect. 7.

2 Data and methods

2.1 Satellite and mooring data

For this study, SLA from the Copernicus Marine Environment Monitoring Service (CMEMS:, last access: 10 December 2020; product ID: SEALEVEL_GLO_PHY_L4_REP_OBSERVATIONS_008_047) gridded multimission satellite altimetry is used to create a proxy for AMOC transport. The CMEMS gridded altimetry includes data from the following satellite missions: Envisat, Geosat Follow-On, Jason-1, and TOPEX/Poseidon interleaved. The satellite product has a 1/4 resolution in space and monthly in time, though the resolution is limited by the underlying track spacing between altimeter passes.

The SLA is compared to data from the RAPID 26 N mooring array. The RAPID array is designed to measure the AMOC at 26 N by using strategically placed moorings in three key regions (Rayner et al.2011; McCarthy et al.2015). The first is in the western boundary, the second at the mid-Atlantic ridge, and the third along the eastern boundary (see inset in Fig. 1a). The RAPID moorings provide continuous dynamic height profiles from conductivity, temperature, and depth using MicroCAT CTDs. Temperature and salinity from individual MicroCAT records are vertically interpolated using the climatological profiles of dT/dp and dS/dp to produce a regularly gridded profile of temperature and salinity on a 20 dbar grid (Johns et al.2005; McCarthy et al.2015). Dynamic height (ϕ) is then calculated from the gridded density profiles as

(1) ϕ ( p ) = 4820 p δ ( p ) d p ,

where δ is the specific volume anomaly. The western boundary used in the RAPID AMOC calculation includes moorings WB2, WBH2, and WB3 (concatenated into a single profile and hereafter referred to as West). When comparing with the SLA, data from the RAPID moorings WB2 (26.5 N and 76.75 W) and WB3 (26.5 N and 76.5 W) as well as West are used for the western component, and an amalgamation of data from the eastern boundary moorings between 13.75–24.22 W and 23.7–27.9 N (concatenated into a single profile and hereafter referred to as EB; Chidichimo et al.2010; McCarthy et al.2015) is used for the eastern component (Sects. 3 and 4) (Fig. 1a). Further details on the instrumentation and observational strategy of the RAPID mooring array can be found in McCarthy et al. (2015).

Figure 1Absolute dynamic topography (ADT) (a) mean and (b) standard deviation (SD) across the North Atlantic over the 2004 to 2018 time period. Western (purple diamonds) and eastern (orange diamonds) moorings are indicated along the RAPID array (black line). The inset in panel (a) shows a vertical cross section of the basin with the location and label of RAPID moorings. (c) Time series of the western (purple line) and eastern (orange line) ADT near the locations of the western (WB2, WB3) and eastern moorings (EB), respectively. Units are metres (m).

RAPID MicroCATs record on an hourly basis. These data are then post-processed with a 2 d low-pass filter to remove tidal fluctuations. The resulting data are used to compute transport estimates, which are gridded at 12 h resolution. A 10 d low-pass filter is applied to final transport estimates (Kanzow et al.2007). RAPID data used here include temperature, salinity, pressure, dynamic height, and estimates of the upper mid-ocean and AMOC transports from 2004 to 2018 (Moat et al.2020a), which can be obtained here: (last access: 2 September 2020).

In this study, all 12 h and daily data are monthly averaged; the annual seasonal climatology is removed and then smoothed with an 18-month filter over the available time periods, except where indicated otherwise. The filter is a Gaussian-weighted moving average filter; end points are computed such that when the number of points available is fewer than the window size, the window is resized and/or truncated, and the average is taken from the elements available in the truncated window. Clément et al. (2014) found that Rossby waves and eddies are important contributors to geostrophic transport on timescales of ∼3–8 months, and Kanzow et al. (2009) found that eddies do not affect AMOC on low-frequency (interannual to decadal) timescales; hence, an 18-month filter is deemed appropriate for removing the influence of eddies. Finally, the mean over 2004 to 2018, which is the overlapping time period between RAPID and the SLA, is removed from their respective time series. When computing correlations and regressions, the linear trend from 2004–2018 is removed from the respective time series, and the statistical significance is determined using the effective degrees of freedom calculated as the length of each time series over the integral timescale of decorrelation (Emery and Thomson2001).

2.2 RAPID transport estimates and SLA

RAPID defines the AMOC transport at 26 N as the sum of three components: the Gulf Stream transport (TGS), the Ekman transport (TEK), and the upper mid-ocean transport (TUMO) (McCarthy et al.2015; Kanzow et al.2010, 2007; Johns et al.2011; Rayner et al.2011) such that

(2) T MOC ( t ) = T GS ( t ) + T EK ( t ) + T UMO ( t ) .

The first component, TGS, has been measured by submarine telephone cables in the Florida Straits since 1982 (Baringer and Larsen2001; Meinen et al.2010) (, last access: 7 September 2020). The second component, TEK, is derived from zonal wind stress from ERA5 reanalysis products (Hersbach et al.2020). The third component, TUMO, is the sum of the western boundary wedge transport (TWBW), the hypsometric mass compensation (TEXT), and the internal geostrophic transport (TINT) over the top 1100 m. TWBW is the northward transport measured between the Abaco Island continental shelf and the WB2 mooring. TEXT functions to compensate for TUMO so that the net meridional flow is zero, and TINT is the internal geostrophic (southward) transport.

Through geostrophic balance, the meridional mass transport is proportional to the integrated pressure difference between the eastern (PE) and western (PW) basin end points such that

(3) T ( t ) = 1 f ρ 0 z ref 0 [ P E ( z , t ) - P W ( z , t ) ] d z ,

where f is the Coriolis parameter, ρ0 is reference density, and zref is the reference level at 4740 m. Pressure at z=-h is then related to sea level displacement and the vertical profile of density as follows:

(4) P = ρ 0 g η - - h 0 ρ g d z ,

where g is gravitational acceleration and η is satellite SLA.

Satellite altimetry measures the surface geostrophic velocity, but surface geostrophic velocity does not give us information about the vertical structure of the horizontal velocity; thus, time-varying geostrophic flow is combined with the first mode of horizontal velocity to obtain TUMO (Sects. 4 and 5). The TGS is estimated using satellite altimetry and linear regression (Sect. 5). The TMOC is constructed by adding the satellite-derived TGS and TUMO with the TEK, obtained from ERA5 wind stress, as per Eq. (2) (Sects. 5 and 6). In the following sections, the relationship between satellite data and the vertical structure of the flow from the moorings is used to develop a new method for estimating AMOC transport.

3 Evaluation of the satellite and mooring data

The CMEMS absolute dynamic topography (ADT) is computed as the sum of SLA and mean dynamic topography (MDT). The MDT is a mean estimate of SSH above the geoid over the given reference period 1993 to 2013 (Rio et al., 2018; further details on this data product here:, last access: 24 September 2021). The ADT shows characteristics of the mean state North Atlantic, with negative ADT marking the subpolar gyre in the north ( 48 N; Fig. 1a) and positive ADT delimiting the subtropical gyre, largely occupying the region between 15 and 40 N (Fig. 1a). Key features are observed at 36 N, where positive ADT represents the Gulf Stream as it separates from the coast at Cape Hatteras (Fig. 1a) and flows north as a free-wheeling jet, eventually feeding into the North Atlantic Current on its path northward. The negative ADT around the region of the Grand Banks (48 N) is representative of the equatorward-flowing Labrador Current that supplies the Slope Sea just north of the Gulf Stream (e.g. Petrie and Drinkwater1993; Fratantoni and Pickart2007), influencing its position (Peña-Molino and Joyce2008; Sanchez-Franks et al.2016). The Labrador Current is also part of the surface-flowing limb of the AMOC. The apparent bipolar structure over the North Atlantic is similar to the EOF mode 1 of the North Atlantic, characterised by Zhang (2008) as the AMOC fingerprint. The standard deviation of the ADT shows that most of the variability is contained within the region of the Gulf Stream after it separates from Cape Hatteras (Fig. 1b). This variability is due to the large latitudinal shifts in the Gulf Stream position following changes in the North Atlantic oscillation (e.g. Sanchez-Franks et al.2014; Sanchez-Franks and Zhang2015; Pérez-Hernández and Joyce2014; Bisagni et al.2017; Taylor and Stephens1998). Further details of the AMOC are reviewed in Zhang et al. (2019).

In the region of the RAPID 26 N mooring array (Fig. 1a, b), the time-varying ADT at western and eastern points shows substantially higher mean and fluctuations in the western boundary (0.73 m root mean square – rms) compared to the eastern boundary (0.21 m rms) over the 2004 to 2018 period (Fig. 1c). Variability from the cross-basin pressure gradient is largely driven by the variability along the western boundary on lower-frequency (periods longer than a year) timescales (Frajka‐Williams2015).

3.1 Surface variability

The RAPID programme estimates geostrophic transport from the difference between the basin-wide eastern and western end point pressure fluctuations (Eq. 3). Correspondingly, to construct a TUMO estimate from satellite altimetry, the time-varying east–west SLA difference (Δη) is compared at each grid point with the RAPID TUMO.

Figure 2a shows the correlation between RAPID TUMO and Δη at each grid point; the eastern point remains fixed at the first easternmost SLA grid point at each latitude, and the western point shifts longitudinally in the westward direction. RAPID TUMO and η have had the annual seasonal climatology removed and an 18-month Gaussian filter applied. In general correlations are highest in the western part of the basin in the 26–30 N latitudinal band and in the east for roughly the same latitudes between 25 and 35 W (Fig. 2a). It is interesting to note that the correlations along the RAPID line are lower than those found at higher latitudes of 27–30 N. This could be due to an improvement in the signal-to-noise ratio further north.

Figure 2(a) Correlation map between the east–west difference in SLA (Δη) and the RAPID upper mid-ocean transport (TUMO). Magenta squares in the western and eastern part of basin indicate regions of maximum correlation between Δη and TUMO. Western (purple diamonds) and eastern (orange diamonds) moorings are indicated along the RAPID array (black line). (b) Correlation between SLA and dynamic height from RAPID moorings West (purple triangle line), WB3 (purple circle line), and EB (orange solid line) at each pressure. Dynamic height from RAPID moorings is referenced to SLA at the surface (ϕη). The dashed lines indicate the 95 % significance level for SLA and dynamic height r values per mooring using the two-tailed t test.

The correlation map illustrates the optimal choice of longitude and latitude coordinates for maximum correlation of Δη with RAPID TUMO, r=0.74 (statistically significant at the 95 % level); the easternmost SLA grid point is located at 27.875 N and 13.125 W, and the western SLA grid point is located at 27.875 N and 74.375 W (Fig. 2a). The east and west dynamic height measurements from the moorings make up the interior geostrophic component of the RAPID TUMO (Sect. 2.2); however, the TUMO also includes a contribution from the Antilles Current. For the satellite data to account for as much of the upper mid-ocean transport variability as possible, it is advantageous to use satellite end points as shown in Fig. 2a to calculate Δη, which here appears to better reflect the changes in the meridional mass transport. This is consistent with other studies which have similarly found higher agreement between RAPID TUMO transport and satellite variability north of 26 N (e.g. Frajka‐Williams2015). Further, sensitivity of the correlation between RAPID TUMO and Δη to the choice of filter was tested for 6-month intervals at 6, 12, and 24 months, in addition to the 18-month Gaussian filter (Fig. S1 in the Supplement). The pattern of correlations between RAPID TUMO and Δη across the subtropical North Atlantic remained generally consistent between the four different choices of filter (as described above for the 18-month filter); correlations decreased (for the eastern SLA grid point located at 27.875 N and 13.125 W and the western SLA grid point located at 27.875 N and 74.375 W) and were lowest at r=0.48 (statistically significant at the 95 % level) using 6-month filtering, increasing to r=0.67 using 12-month filtering, and r=0.79 when using the 24-month filter, contrasted with r=0.74 when using the 18-month filter.

3.2 Variability in the vertical

In order to determine how well the SLA can estimate the upper 1000 m circulation, it is useful to assess to what depth the variability from SLA captures subsurface variability. Thus, SLA is compared to the dynamic height (Eq. 1) at every depth from the surface to 1100 dbar from the RAPID moorings in the western and eastern boundary. Figure 2b shows the correlation (r values) between SLA west and east grid points (as indicated in previous section) and dynamic height from the RAPID moorings: West, WB3, and EB. In this section only, the respective RAPID mooring dynamic height has been referenced at the surface to SLA (hereafter ϕη) to get a correlation value of 1 at the surface, which decreases with depth to assess to what depth surface variability is coherent. In the western boundary, SLA at 27.875 N and 74.375 W is compared with both WB3 and West, and the correlation in the top 1100 dbar is found to be greater than 0.79 everywhere (statistically significant at the 95 % level). SLA has a higher correlation with ϕη at mooring West compared to WB3, maintaining a statistically significant correlation coefficient above r=0.88 in the top 1100 dbar. Correlation between SLA and WB3 decreases more abruptly at around 300 dbar. In the eastern boundary, correlation between SLA at 27.875 N, 13.125 W and ϕη at EB is similarly high throughout the top 1100 dbar of the ocean, albeit weaker compared with the western boundary mooring, with correlations from r=1 at the surface decreasing to r=0.77 at a depth of 1100 dbar (Fig. 2b). These results suggest that the variability observed at the sea surface is a good measure and coherent with variability to at least a depth of 1100 dbar. This is in agreement with Clément et al. (2014), who showed that isopycnal displacements at the RAPID western mooring locations agree well with satellite data. Specifically, Clément et al. (2014) did a similar analysis using profiles of isopycnal displacements, instead of dynamic height, from the moorings and satellite. They found statistically significant correlation (r=0.5–0.6) between WB3 and WB2 with SSHA, respectively, in the upper 1000 m over 2004 to 2011. They did not analyse moorings in the eastern boundary. Given that we have referenced the dynamic height to SLA at the surface here and used a longer time period and different filtering, the discrepancies between our studies are reasonable.

For completeness, the SLA is also compared to dynamic height anomaly thickness for the 0–1100 dbar layer (Fig. S2). The dynamic thickness is calculated from the difference in dynamic height (referenced to 4820 dbar) at 0 and 1100 dbar and gives a measure of the shear between those levels (independent from choice of reference level). In both the western and eastern basin, the dynamic height thickness anomaly has a standard deviation smaller (0.30 m2 s−2 at WB3 and 0.07 m2 s−2 at EB) than the SLA multiplied by gravitational acceleration (0.70 m2 s−2 in the west and 0.14 m2 s−2 in the east). This result suggests that the baroclinic structure in the upper 1100 dbar is not overwhelming the SLA signal, and the SLA can be used for transport estimates in the upper 1100 dbar.

4 Rossby wave theory and horizontal velocity modes

4.1 Westward propagation

Wind stress and density fluctuations dominate AMOC variability on seasonal and subannual timescales (Hirschi et al.2007). Wind-driven variability may also play a role on interannual timescales; for example, during the winter of 2009/2010, anomalous wind-driven Ekman transport contributed to interannual AMOC fluctuations (Zhao and Johns2014; Evans et al.2017). However, density variability in the upper 1000 m has also been identified as a leading driver of the AMOC on interannual timescales (Hirschi et al.2007). This density variability is associated with isopycnal perturbations that travel westward as long Rossby waves (Hirschi et al.2007, 2009). The Rossby waves impact the upper mid-ocean component of the AMOC transport through their effect on the east–west density structure of the basin (Hirschi et al.2007; Cabanes et al.2008; Hirschi et al.2009).

The changes in the upper 1000 m of the density field from westward propagation are visible as a sea surface signal and have been shown to be captured by satellite altimetry. Altimeters are most likely to reflect the first baroclinic mode, and by association motion of the main thermocline due to the nature of baroclinic modes, which is surface-intensified (Wunsch1997). In the subtropical North Atlantic, studies have used satellite altimetry to track westward-propagating anomalies (e.g. Hirschi et al.2007; Kanzow et al.2009; Clément et al.2014). A Hovmöller diagram of SLA along the 26 N RAPID line shows how perturbations in the eastern boundary propagate westwards up to  79 W before they reach the Bahamas and abruptly diminish (Fig. 3). Westward propagation is visible at monthly resolution and even to some extent with 18-month smoothing applied to the data (Fig. 3b). Amplitudes in the western boundary are larger compared to the eastern boundary, and they are particularly prominent during 1996–1997, 2003–2004, 2006, and post-2015. Interannual variations are observed in the western basin (between 70 and 80 W), and the strong positive anomalies apparent during 2003–2004 and negative anomalies in 2006–2007 are consistent with Kanzow et al. (2009, their Fig. 5). The speed of these westward-propagating anomalies is typically similar to baroclinic Rossby wave phase speed (Gill1982; Killworth and Blundell2003; Hirschi et al.2007).

Figure 3Hovmöller diagram of sea level anomaly (SLA) at 1/4 resolution in space at 26.625 N over 1993 to 2018 for monthly resolution (a) and 18-month smoothed (b). The annual seasonal climatology has been removed from the SLA. Units are metres (m).


4.2 Modal decomposition

Though satellite altimetry can measure surface geostrophic velocities, it cannot be used to infer the vertical structure of the flow. To understand the contribution of the flow's vertical structure to the meridional mass transport and how that impacts satellite-derived transport, pressure modes derived from RAPID mooring data are examined.

The vertical structure of the flow is assessed using normal-mode decomposition, assuming a flat-bottomed and motionless ocean (Gill1982), using data from the RAPID moorings at the western (West, WB3) and eastern (EB) boundary. The normal-mode decomposition for the given mode number n=0,1,2 is determined from the Sturm–Liouville equation:


where N(z) is the Brunt–Vaisala frequency and cn is the phase/modal speed of the waves. The boundary conditions are defined as Gn=0 at z=0,-H.

The buoyancy frequency (N) is estimated from temperature and salinity profiles at moorings West, WB3, and EB, which are converted into density profiles, ρ(z), and averaged over the 2004 to 2018 time period:

(7) N ( z ) = - g ρ 0 ρ ( z ) z .

Normalisation is then computed to satisfy Kronecker delta, δij (rendering the resulting pressure modes dimensionless), such that

(8) - H 0 F i F j d z = H δ i j ,

where H is 4820 dbar, which is the bottom reference depth used in the AMOC calculation.

The buoyancy frequency and the pressure modes, Fn (as per Eqs. 5 and 6), are estimated at West, WB3, and EB (Fig. 4). The stratification profiles are characteristic of buoyancy profiles in the North Atlantic (e.g. Szuts et al.2012; Clément et al.2014): in the west, the buoyancy frequency shows a sharp maximum at ∼80dbar and a second maxima around 600 to 800 dbar before decreasing to background stratification levels around 1100 dbar. These two peaks indicate strong stratification linked to the seasonal and main pycnocline (Siegel et al.1999). Moving eastward across the basin, at WB3, the upper surface stratification decreases and the peak deepens to 100 dbar; the second peak also deepens and has slightly lower stratification than at mooring West. The stratification minima apparent in West and WB3 at ∼300–400 dbar is indicative of 18 C water (i.e. subtropical mode water), characteristic of the western subtropical North Atlantic. On the other side of the basin, in the eastern boundary, the stratification at EB is decreased though the subsurface maximum remains at 100 dbar.

Figure 4The buoyancy frequency (N2; dashed purple line) and the horizontal velocity normal modes: the barotropic mode (i.e. mode 0; grey line) and the baroclinic normal mode 1 (blue line), mode 2 (orange line), and mode 3 (yellow line) at RAPID moorings (a) West, (b) WB3, and (c) EB. The modes have been normalised to satisfy Eq. (8) and are dimensionless.


Figure 4 also shows the first four pressure modes at West, WB3, and EB. Mode 0 is representative of the barotropic mode, while modes 1 to 3 are the first three baroclinic modes. The first mode shows a zero crossing roughly between 1000 and 1100 dbar in the western boundary (West and WB3), which deepens to ∼1500dbar by the time it reaches the eastern boundary (EB). In general, western moorings show shallower zero crossings and more complex structure, reflective of the peaks in stratification in the west compared to the eastern mooring, which has a deeper zero crossing and smoother structure. The vertical structure of the horizontal velocity shown here is consistent with the characteristic shape of known modes in the western North Atlantic (at this latitude) (e.g. Gill1982; Szuts et al.2012, their Figs. 6.14c and 2, respectively). The shape of the first baroclinic mode in the surface to its first zero crossing (1100–1500 dbar) will be key in informing the satellite-derived estimates of transport, as altimetry has been shown to capture the upper 1000 m intensified structure of the flow (away from the boundaries by ∼50km or more), reflective of the first baroclinic normal mode (Wunsch1997; Szuts et al.2012). The zeroth mode, i.e. the barotropic mode, is important in the deeper ocean (>1000 m), which is less stratified (Wunsch1997; Kanzow et al.2008). However, below the upper 1000 m, changes in the deeper transport are not captured by satellite altimetry (Kanzow et al.2008).

In the following sections, in-depth analysis of the normal modes from moorings and their relationship with satellite altimetry is explored.

4.3 Modal amplitude and variability

To understand the importance of each mode to the total variability observed in dynamic height anomaly, ϕ (Eq. 1, referenced to 4820 dbar), measured by the RAPID moorings, the modal amplitude is analysed and used to construct time-varying ϕ from the normal modes. The amplitudes of the modes, an, are estimated by integrating the product of the normal mode, Fn, and ϕ such that

(9) a n ( t ) = 1 H - H 0 F n ( z ) ϕ ( z , t ) d z ,

where n=0,1,2,3 is the mode number and H is the reference pressure 4820 dbar. Because satellite altimetry has been shown to reflect fluctuations associated with the amplitude of the first baroclinic mode (Wunsch and Stammer1997; Hirschi et al.2009), the amplitude for the first mode, a1, can also be estimated from SLA, η, such that

(10) a ^ ( t ) = s g F 1 ( z = 0 ) η ( t ) ,

where F1(z=0) is the surface value of the first mode, and s is an empirically determined scale factor equal to 0.25.

Because of baroclinic shear, SLA alone does not determine transport in the upper 1100 m of the water column: if we integrated the SLA over the upper 1100 m to obtain geostrophic transport, the SLA would overestimate the transport magnitude compared to dynamic height (e.g. Fig. S3). For this reason, it is necessary to combine the SLA with a function that describes the shape of the geostrophic transport profile. We have used the first baroclinic mode for this function; however, a comparison between the dynamic height from the RAPID moorings and the SLA indicates that further correction is needed. Thus, a scale factor is determined empirically by examining the signal from ϕ at the surface (z=0) at moorings West and EB against η.

Figure 5(a) The root mean square (rms) at every longitude point of the SLA (η) at 26.125 N (black dashed line; latitude of RAPID moorings) and at 27.875 N (black solid line; latitude of maximum correlation between Δη and TUMO – see Fig. 2). Symbols indicate the rms of the dynamic height anomaly (ϕ) at the surface (z=0) at the West (purple triangle) and EB (orange square) moorings. (b) Scatterplot of the east–west difference in η, i.e. Δη, vs. east (EB)–west (West) difference in the ϕ, i.e. Δϕ, and best-fit regression line (solid line). In (b) data are detrended and subsampled once every 6 months. In both plots, η is multiplied by gravitational acceleration. Dashed lines indicate the 95 % confidence interval (units: m2 s−2).


Figure 5a shows the rms of dynamic height and sea level anomaly (η multiplied by gravitational acceleration) at two latitudes across the Atlantic. At the eastern boundary, the rms of gη is about twice (0.18 m2 s−2) that of the dynamic height from the moorings (0.08 m2 s−2) (Fig. 5a). Moving west, rms values of gη peak at 74 and 75 W (rms of gη(27.875 N, 74 W) = 0.70 m2 s−2 and rms of gη(26.125 N, 75 W) =0.59 m2 s−2) before decreasing to the western boundary (rms of gη(27.875 N, 77 W) =0.56 m2 s−2 and rms of gη(26.125 N,77 W) =0.39 m2 s−2). The rapid decrease in sea level variability at the western boundary is due to mesoscale suppression associated with the continental slope (Kanzow et al.2009). The moorings, in contrast, show lower variance in dynamic height at both the eastern and western boundaries (rms of ϕ(West) =0.31 m2 s−2 and rms of ϕ(EB) =0.08 m2 s−2). A scatterplot of the east–west difference in surface ϕ (i.e. Δϕ=ϕEB(z=0)-ϕWest(z=0)) and the equivalent from sea level anomaly (gΔη=g(η(27.875 N, 13.125 W)-η(27.875 N, 74.375 W))) shows general agreement between the two values. The slope of the regression line is ∼0.25, and the intercept goes through the origin (Fig. 5b). Figure 5 suggests that the satellite altimetry alone does not capture the same magnitude of signal observed by the moorings. One of the reasons for this discrepancy may be the proximity of the moorings (e.g. WB2) to land, where variability experiences changes due to coastal processes (Kanzow et al.2009). Kanzow et al. (2009) suggest that near (within 100 km) the western boundary, variability becomes influenced by a combination of barotropic–baroclinic flow over the slope, reduced eddy variability, and coastally trapped waves. The value of the slope of the regression line, i.e. 0.25, is thus used as the scale factor for the satellite where indicated (e.g. Eq. 10).

The modal amplitude, a^, computed from η (Eq. 10) is compared with the modal amplitude, an, computed from the first three baroclinic normal modes (Eq. 9) at West and EB (Fig. 6). The amplitudes of the first mode, a1, and a^ have a correlation of r=0.57 (significant at the 90 % level) at West and a correlation of r=0.67 (significant at the 90 % level) at EB. In the east (west), a^ is constructed using η at 27.875 N and 13.125 W (27.875 N and 74.375 W). Comparison of a1 and a^ time series at West shows a deviation in 2005–2006, when a^ has a larger negative anomaly (0.15 m2 s−2) compared to a1 (0.02 m2 s−2); this is likely linked to the collapse of the WB2 mooring from November 2005 to March 2006, when the WB3 mooring was used to fill the gap in the WB2 and West data. In general, all three modes start off with the same sign (negative) and eventually change sign (positive) with the exception of mode a2, which starts positive. The modes have been defined or set to positive at the surface. The higher modes, a2 and a3, have smaller amplitudes compared to a1.

The modal amplitude (Eq. 9) can also be combined with the pressure modes (Eq. 6) to reconstruct ϕ as follows:

(11) ϕ n * ( z , t ) = a n ( t ) F n ( z ) .

Figure 7 shows the reconstructed ϕ, i.e. ϕ*, using only the first baroclinic mode at West and EB as per Eq. (11) from the surface to the reference pressure 4820 dbar over the 2004 to 2018 period. Using only the first baroclinic mode, ϕn* for n=1, at West shows a general undulating pattern of negative and positive anomalies in the upper 1000 dbar: the anomalies are negative for the years 2007–2008, 2011, and 2014 and positive during 2009, 2013, and from 2015 to 2017, similar to the surface anomalies observed in ϕ at West (Fig. 7a, b). The difference between ϕ and ϕ* at West suggests that although the reconstruction generally captures the intensified structure in the upper 1000 dbar, it underestimates the magnitude of ϕ and does not capture any of the signal below 1000 dbar (Fig. 7c). These differences may be due to the fact that the reconstruction uses only the first baroclinic mode, reflective of changes in the main thermocline, and thus has an oversimplified structure in the upper layer of the ocean compared to the ϕ at West. Further, the baroclinic structure, in general, is more likely to reflect changes in the stratified upper ocean to roughly 1000 m, while in the deeper less stratified ocean, barotropic motion is more important (Kanzow et al.2008). In the eastern boundary, ϕ* is reconstructed using the first baroclinic mode at EB (Fig. 7e), similar to pressure perturbation analysis as per Szuts et al. (2012). The reconstruction captures similar negative and positive anomalies in the upper 1000 dbar as those observed in ϕ at EB (Fig. 7d). The magnitudes of the ϕ anomalies are generally substantially smaller at the eastern part of the basin compared to the west.

Figure 6Modal amplitudes, an, fit to the dynamic height anomaly (ϕ) and the baroclinic modes, Fn, as per Eq. (9) at (a) the West and (b) EB moorings. The first three modal amplitudes a1 (blue line), a2 (orange dashed line), and a3 (yellow dashed line) are compared with first mode amplitude estimated from SLA, a^ (black line), as per Eq. (10) (units: m2 s−2).


Figure 7Dynamic height anomaly (ϕ) from the surface to the reference pressure 4820 dbar over the 2004 to 2018 period from RAPID moorings at West (a) and EB (d). Reconstruction of the dynamic height anomaly (ϕ*) at West (b) and EB (e) using the first baroclinic mode (F1) and the corresponding modal amplitude (a1) as per Eq. (11). The difference between the ϕ* for the first baroclinic mode and ϕ at West (c) and EB (f) (units: m2 s−2).


The contribution of the barotropic mode (mode 0) and higher baroclinic modes (modes 1 and 2) to ϕ* is examined in Fig. 8. The inclusion of the barotropic mode and higher baroclinic modes in ϕ* shows substantially improved patterns of the ϕ variability in the upper 1000 dbar, with reduced differences at every pressure level (Fig. 8) compared to ϕ* using only mode 1 (Fig. 7). To assess whether the inclusion of the barotropic and higher baroclinic modes makes a meaningful contribution to the upper ocean transport, transports from the ϕ* reconstruction using the normal modes and the ϕ at West and EB are presented in the following section.

4.4 Transport anomalies from West and EB moorings

Transport anomalies at West and EB are estimated using ϕ, ϕ*, or η, with the aim of understanding the contribution of the normal modes to ϕ variability and implications for satellite-based transport. Kanzow et al. (2010), Chidichimo et al. (2010), and Szuts et al. (2012) set the precedent for estimating transport from a single mooring (as opposed to a horizontal gradient) to separate the contribution of the western and eastern components of the basin-wide geostrophic transport. Therefore, transport anomalies from ϕ at the West or EB mooring (as indicated) is defined as

(12) T ϕ ( t ) = 1 f - 1100 0 ϕ ( z , t ) d z .

The mode-reconstructed ϕ* transport is calculated using the modes (Eq. 6) and modal amplitude (Eq. 9) such that

(13) T F n ( t ) = Σ T F n ( t ) = Σ 1 f a n ( t ) - 1100 0 F n ( z ) d z .

If the amplitude of the first mode is taken to correspond to η (i.e. a^; Eq. 10), then the satellite-based transport can be constructed using the vertical structure from the first baroclinic mode F1 in the following manner.

(14) T η ( t ) = 1 f a ^ ( t ) - 1100 0 F 1 ( z ) d z

A comparison of the Tϕ (for simplicity, the primes are dropped hereafter) and the reconstructed TF transports shows that TF accounts for 83 % of the variability of Tϕ (r=0.91, statistically significant at the 95 % level) at West and 90 % of the variance of Tϕ (r=0.95, statistically significant at the 95 % level) at EB, where TF is reconstructed using only the first mode derived from the dynamic height profiles (Fig. 9a, d). The satellite-based transport, Tη, on the other hand, shows lower correlation with Tϕ at both West (r=0.60, statistically significant at the 95 % level) and EB (r=0.49, statistically significant at the 95 % level). Though Tη captures the general patterns of variability observed in Tϕ, there are several differences apparent between Tη and Tϕ occurring around 2005/2006, which is when WB2 collapsed and data from WB3 were used instead, and from 2004–2006 at EB. Additionally, the Tη time series generally appears to underestimate the magnitude of the Tϕ at West during 2009–2010, 2012, and 2016–2017 and at EB during 2010–2012 and 2016 (Fig. 9a, d). The results shown here are comparable to a satellite and RAPID mooring comparison featured in Szuts et al. (2012, their Fig. 7), wherein correlations between SSH and geopotential anomalies from WB3 and WB2 were approximately r=0.5 and 0.7, respectively; r=0.75 at EB1 and was close to zero at EBH. In Szuts et al. (2012) a 10 d low-pass filter was used prior to correlation. Contrasted with results here, which include a longer time period, 18-month smoothing, and differences in the location of the SSH data, it is not surprising there are differences to results from Szuts et al. (2012). Szuts et al. (2012) further posited that those differences between moorings and satellites are likely due to (a) loss of small spatial and temporal variability that occurs with satellite products and (b) satellite data include surface-layer processes which moorings do not entirely capture.

Figure 8Same as in Fig. 7, except the contribution of the barotropic mode (n=0) and the first two baroclinic modes (n=1,2) are included in the dynamic height anomaly reconstruction (ϕ*) at West (b) and EB (e).


Figure 9Transport anomalies estimated from dynamic height anomalies, Tϕ (orange line), as per Eq. (12) at West (a, b, c), EB (d, e, f), and EB–West (g, h, i). Transport anomalies estimated from the reconstructed dynamic height anomaly, TFn as per Eq. (13), for mode 1 (yellow line; a, d, g), modes 1 and 2 (green line; b, e, h), and modes 0 and 1 (purple line; c, f, i) at West, EB, and EB–West, respectively. Satellite-based transport, Tη (blue line), as per Eq. (14) is estimated in the (a) west (η at 27.875 N and 74.375 W), (d) east (η at 27.875 N and 13.125 W), and (g) east–west. Correlation between Tϕ and TF is indicated in each plot (units: Sv).


Including the second baroclinic mode as well as the first in the reconstructed transport, TF, shows some improvement in the correlation with Tϕ at West (r=0.98, statistically significant at 99 %) and a slight decrease at EB (r=0.94, statistically significant at 99 %), though the amplitude of TF still underestimates Tϕ (Fig. 9b, e). The addition of the barotropic mode to the first baroclinic mode in TF noticeably improves the amplitude of TF and correlation with Tϕ at West (r=0.99) and EB (r=0.999) (Fig. 9c, f).

The basin-wide geostrophic transport can be constructed for ϕ and ϕ*, respectively, by integrating the east–west difference such that Eqs. (12) and (13) become


where the subscripts E and W denote east (EB) and west (West), respectively. TΔF and TΔϕ have a correlation of r=0.91 (significant at the 99 % level) when TΔF is constructed using only mode 1 (Fig. 9g). The correlation between TΔF and TΔϕ increases to r=0.97 when TΔF is constructed using modes 1 and 2 and to r=0.99 when TΔF is constructed using modes 0 and 1 (Fig. 9h, i). These results suggest that the time-averaged first baroclinic mode accounts for most of the interior geostrophic transport variability, and the combined barotropic and first baroclinic mode accounts for 98 %. The barotropic mode is reflective of changes in the deeper less stratified ocean.

5 Construction of the GS, the UMO, and AMOC transports

The Gulf Stream within the Florida Straits has a mean of 31.2 Sv and is balanced by the UMO and Ekman transports, with means of 18 and 3.74 Sv, respectively, to yield the total mean AMOC transport of 17 Sv (McCarthy et al.2015). The Gulf Stream time series is based on a submarine telephone cable that has been recording data between the Bahamas and Florida at 27 N since 1982 (Baringer and Larsen2001). Principles of geostrophy have been previously used to provide alternative mechanisms for estimating the cable-based TGS using satellite altimetry (Volkov et al.2020) and pressure gauges (Meinen et al.2020).

Here, the east–west difference in η (Δη) in the western end of the basin (i.e. west of 77 W) is compared with Gulf Stream transport anomalies from the submarine cable data (Baringer and Larsen2001) over the RAPID time period (i.e. 2004 to 2018). Maximum correlation between TGS from the cable data and the satellite Δη (r=0.70, statistically significant at the 95 % level) is found when using ηE located at 27.625 N and 77.125 W and ηW located at 27.625 N and 80.125 W (Fig. S4). Similarly, Volkov et al. (2020) used along-track satellite altimetry to infer TGS measured at the level of the Florida Straits; they found that satellite altimetry captures 56 % of variability observed from the TGS estimated from submarine cable records for the 2006–2020 time periods, when the cable-based transport estimates were subsampled at 10 d intervals to coincide with the along-track satellite passes. We use gridded altimetry with a similar objective (to determine the Gulf Stream transport via the Florida Straits from satellite altimetry) but with a focus on the lower-frequency variability rather than an instantaneous comparison with the cable measurements.

Linear regression is used against the cable-based TGS to obtain the regression coefficients (a,b) needed to produce a final estimate for satellite-based GS transport:

(17) T GS * ( t ) = a Δ η + b ,

where a=8.0087, b=0.0026, and the eastern and western SLA points that correspond to Δη are located at 27.625 N and 77.125 W for ηE and at 27.625 N and 80.125 W for ηW. The accuracy of TGS* was determined using a Monte Carlo technique, whereby 90 % of the time series was randomly sampled 10 000 times, performing a linear regression to obtain an upper and lower bound (mean value of ±0.0614 Sv) for the 90 % confidence intervals. The satellite-derived TGS* captures 49 % of the variability of the in situ TGS and generally underestimates the amplitude of variability (Fig. 10a, b).

Figure 10(a) Time series of the Gulf Stream transport (TGS) from submarine cables (orange line; (Baringer and Larsen2001)) and the satellite-derived TGS* estimates as per Eq. (17) (blue line). (b) Scatterplot between east (27.625 N and 77.125 W) minus west (27.625 N and 80.125 W) SLA (Δη) and TGS with the line of best-fit linear regression (note y and x axis limits reduced to improved visibility of data points). (c) RAPID upper mid-ocean transport (TUMO; orange line) and the satellite-based TUMO* transport estimates as per Eq. (18) (blue line). (d) Scatterplot between the TUMO* and TUMO with the line of best-fit linear regression. (e) RAPID MOC transport (TMOC; orange line) and the satellite-based TMOC* transport estimates as per Eq. (19) (blue line). (f) Scatterplot between the TMOC* and TMOC with the line of best-fit linear regression. Blue shading indicates uncertainty estimates using a Monte Carlo technique (a, c) as described in text, and blue shading in (e) is the sum of the confidence intervals from (a), (c). In (b), (d), (f) data are subsampled once every 6 months. Units in Sv.


The satellite-based TUMO* is constructed by combining the satellite-based mode amplitude, a^ (Eq. 10), with the first baroclinic mode, F1(z) (Eq. 6), as follows:

(18) T UMO * ( t ) = 1 f a E ^ ( t ) - 1100 0 F E ( z ) d z - a W ^ ( t ) - 1100 0 F W ( z ) d z ,

where the subscripts E and W denote east (for ηE at 27.875 N and 13.125 W and FE at EB) and west (for ηW at 27.875 N and 74.375 W and FW at West), respectively. The satellite-derived TUMO* is compared with RAPID TUMO time series (r=0.75, statistically significant at the 95 % level) in Fig. 10c. TUMO* follows the weakening apparent in the RAPID TUMO in the year 2006 (positive anomalies indicate weakening in southward transport and vice versa), followed by the sharp intensification of the southward flow, reaching its maximum in 2009–2010. These changes in the transport have been previously documented in Smeed et al. (2014, 2018), and in particular the 2008–2010 intensification in RAPID TUMO is linked to strengthening of southward circulation in the main thermocline (Smeed et al.2014). In general, the satellite-derived TUMO* follows the overall pattern of variability measured by the TUMO, capturing the weakening of southward flow in 2005–2006 and 2011, as well as intensification of southward flow in 2009–2010 and 2016–2017. However, the TUMO* does not capture the full extent of the RAPID TUMO intensification in 2009–2010 or in 2012 by approximately 1 Sv (Fig. 10c, d). Some of the missing variability may be linked to the barotropic component not captured by satellite altimetry, as shown in the previous section (Sect. 4.4). The 90 % confidence interval for TUMO* was estimated using a Monte Carlo technique, whereby 90 % of the Δη time series is randomly sampled 10 000 times, performing the regression with Δϕ to acquire the scale factor as outlined in Sect. 3.3, giving the upper and lower bounds for calculation of TUMO*. Sensitivity of these results to the choice of filtering window was tested. We found that the correlations decreased for TGS* (r=0.6, r=0.68) and TUMO* (r=0.42, r=0.67) when using a 6-month and a 12-month Gaussian filter, respectively (Figs. S5, S6), compared to the 18-month Gaussian filter. Conversely, correlations increased for TGS* (r=0.73) and TUMO* (r=0.79) when using a 24-month Gaussian filter (Fig. S7).

TMOC* is the sum of southward flow TUMO*, the Ekman transport, and northward-flowing TGS. By incorporating the satellite-derived TGS* as per Eq. (17) and the satellite-derived TUMO* as per Eq. (18), TMOC* can be defined by rewriting Eq. (2) such that

(19) T MOC * ( t ) = T GS * ( t ) + T EK ( t ) + T UMO * ( t ) ,

where the Ekman transport (TEK) is derived from ERA5 zonal wind stress as per the RAPID time series method. The satellite-constructed TMOC* and the RAPID TMOC have a correlation of r=0.83 (statistically significant at the 95 % level; Fig. 10e). TMOC* captures the start of the RAPID TMOC in 2004, the weakening in 2005, the subsequent intensification in 2006, and the gradual weakening of RAPID TMOC, reaching a maximum low in 2010 and back to zero in 2011. In 2010 and 2013, the TMOC* noticeably underestimates the magnitude (>1Sv) of the RAPID TMOC fluctuations. The anomalously low event in TUMO during 2009–2010 is attributed to anomalous wind-driven Ekman transport (Roberts et al.2013; Zhao and Johns2014). After the 2009–2010 event, the AMOC recovers and appears to show a positive tendency. The gradual increase in the TMOC apparent from 2010 onwards, however, has not yet been shown to be statistically significant at 26 N (Moat et al.2020b). Using a smaller filtering window, which would allow for higher-frequency variability, we find that the correlation for TMOC* decreases (r=0.71, r=0.81) for 6- and 12-month Gaussian filters (Figs. S5, S6), respectively, but increases (r=0.85) when using a 24-month Gaussian filter (Fig. S7).

The satellite-derived TUMO* and TMOC* were also computed with a time-varying F1(z,t) with almost no improvement: the correlation between TUMO*, using a time-varying F1, and the RAPID TUMO increased slightly from r=0.75 to r=0.76, and for TMOC* and RAPID TMOC there was no change (i.e. r=0.83). This suggests that 69 % of AMOC variability can be estimated from SLA and the time-averaged pressure modes on interannual timescales.

6 Construction of volume transport from historical altimetry data

TOPEX/Poseidon was one of the first satellite altimetry missions to be launched with a major oceanographic focus (Fu et al.1994); since its launch in 1992, it has provided an invaluable source of measurements for studying surface ocean circulation (e.g. Willis2010; Frajka‐Williams2015). In this section, the TUMO* and the TMOC* are constructed for the full available satellite time period of 1993 to 2018 (full annual coverage in the satellite data for the product used here starts in January 1993), yielding 11 years of data that pre-date the RAPID mooring array.

Previously, satellite-derived estimates of the TUMO and the TMOC have been compared with the RAPID 26 N array and constructed for the 1993 to 2013 period in Frajka‐Williams (2015, hereafter EFW15). EFW15 found that SLA could be used to obtain a proxy for mass overturning transport. The statistical relationship between the RAPID TUMO and the SLA on low-frequency (18-month) timescales provided the main method for estimating the EFW15 TUMO. The EFW15 TMOC was estimated by adding the same estimates of TEK (ERA-Interim) and TGS (submarine cable) used by RAPID to the TUMO. Here it is fitting to compare the EFW15 method with the method derived in Sects. 4 and 5, which relies on altimetric data and geostrophic balance and therefore gives a consistent (satellite only) method that incorporates climatological stratification and has the potential to be applied to other latitudes. This is in contrast with EFW15, which was a statistically based method that assumes the existence of in situ observations. Further work is needed to determine the origins of the scale factor required to balance the transports, whether due to sampling issues or uncertainties in satellite sea level anomaly near the boundaries. Additionally, the RAPID TMOC is compared with a third separate TMOC estimated from GloSea5, a global ocean and ice reanalysis product (Blockley et al.2014; Jackson et al.2016, 2019; MacLachlan et al.2015). GloSea5 uses the ocean model NEMO (Nucleus for European Modelling of the Ocean) that has a 1/4 ocean resolution and assimilates observational data including satellite (Megann et al.2014). The GloSea5 TMOC estimates at 26 N are averaged monthly (further details on GloSea5 described in Jackson et al.2019). As in previous sections, annual seasonal climatology is removed from each time series and an 18-month Gaussian filter is applied (with the exception of the EFW15 time series, which uses an 18-month Tukey filter). Since the overlap between the RAPID array data and the EFW15 time series spans 2004 to 2013 only, all correlations presented in this section reflect the r value for the 2004 to 2013 time period.

The TUMO* and the EFW15 TUMO are first compared with the RAPID TUMO time series in Fig. 11a. Over the 2004 to 2013 time period, the TUMO* and EFW15 TUMO have correlations of r=0.79 and r=0.85 (statistically significant at the 90 % level), respectively, with the RAPID TUMO time series. Both time series capture the RAPID TUMO weakening in 2005–2006 (positive anomalies indicate weakening in southward transport and vice versa) and 2010–2011 as well as the strengthening in 2009–2010 and 2013 (Fig. 11a). Both the TUMO* and EFW15 TUMO underestimate the 2009–2010 intensification observed by the RAPID TUMO. In the time period pre-dating the RAPID programme, EFW15 TUMO shows a weakening in 1995 and 2001 and intensification in 1999 and 2004 of the southward flow. In comparison, the TUMO* shows lower-amplitude variability but a higher southward transport in the 1993–2003 time period. The EFW15 TUMO is centred on a single SLA point in the western end of the basin at 70 W and 30 N, while in the satellite estimate derived here, TUMO* is based on the difference in SLA in the east at 13.125 W and 27.875 N and west at 74.375 W and 27.875 N (i.e. variability in the eastern end of the basin is taken into account as well). Thus, holding the eastern component of the TUMO* constant highlights its contribution to total TUMO* variability, which in the pre-RAPID time period acts to strengthen the mean TUMO* southward transport; in other words, holding the eastern component constant leads to a weakening in southward transport that suggests improved agreement with the magnitude of EFW15 TUMO (Fig. 11a).

Figure 11(a) The upper mid-ocean transport, TUMO*, from satellites (η) (blue line) as per Eq. (18), TUMO* with the eastern component held constant (blue dashed line), TUMO from RAPID (orange line), and TUMO from EFW15 (yellow line). (b) MOC transport, TMOC*, from η (blue line) as per Eq. (19); TMOC from RAPID (orange line), EFW15 (yellow line), and GloSea5 (purple line). Correlations indicated in the graphs are over the period of overlap between the RAPID and EFW15 time series, i.e. 2004 to 2013 (units: Sv).


The TMOC*, EFW15 TMOC, and GloSea5 TMOC are compared to the RAPID TMOC over the 2004–2013 time period (Fig. 11b). The TMOC* has the highest correlation with the RAPID TMOC (r=0.87, statistically significant at the 95 % level), while EFW15 TMOC and GloSea5 TMOC show slightly lower correlations with the RAPID TMOC (r=0.84 and r=0.85, respectively, statistically significant for at least the 90 % level). The GloSea5 TMOC proxy underestimates the TMOC mean values in the 2004 to 2010 and does not capture the strengthening apparent in the RAPID TMOC in 2006 (Fig. 11b). In the pre-RAPID time period of 1993–2003, there is surprisingly little agreement between the TMOC proxies. The EFW15 TMOC reconstruction shows a strengthening in 1996–1997 (positive anomaly), while the TMOC* shows weakening in 1997 (negative anomaly). None of the reconstructions quite agree with the general middle to late 90s strengthening in the AMOC, associated with changes in the Atlantic multidecadal variability phase, marked by a positive phase in the late 1990s, as discussed by Zhang et al. (2019). In 2003–2004, GloSea5 TMOC and TMOC* show a strengthening closer to the TMOC values at the start of the RAPID programme, which is not captured by the EFW15 TMOC. In general, it is difficult to determine which of the three TMOC reconstructions might provide a best estimate of the observed TMOC over the 1993–2004 time period. For instance, the correlations between each TMOC reconstruction and RAPID TMOC are very similar; however, there is poor agreement between the three TMOC reconstructions in the time period pre-dating RAPID. The reason for this disagreement may be overfitting. In the case of the satellite-based transport reconstruction, EFW TMOC, the statistical method used to estimate transport relied on RAPID mooring data, and the resulting time series in the period pre-dating RAPID of 1993–2003 may be biased by overfitting of data during the RAPID period of 2004–2013. Similarly, TMOC* may be biased by choice of the scale factor, here based on comparisons between satellites and moorings during the RAPID time period. Further investigation is therefore needed to confidently reconstruct the MOC for the 1993 to 2003 time period.

7 Summary and conclusions

Mooring array programmes such as the RAPID 26 N array have made step-change advancements in our understanding of the AMOC; however, they are limited to measurements at a single latitude and alone cannot be used to infer upstream–downstream changes in the larger-scale structure of the AMOC. Mooring arrays are also a single point of failure should the moorings collapse. In light of the increase in studies linking the AMOC to climate impacts (Zhang et al.2019), it is now more important than ever to find long-term and cost-effective replacements and/or backup systems that can monitor changes in the AMOC.

One previous effort to use satellite altimetry as a proxy for AMOC variability, led by EFW15, relied on a statistical relationship between the RAPID TUMO and the SLA on low-frequency (18-month) timescales to estimate the TUMO. The EFW15 TMOC was estimated by adding the same estimates of TEK (ERA5 reanalysis) and TGS (submarine cable) used by RAPID to the TUMO. Their method successfully recovered 85 % of the TUMO and 92 % of the TMOC on interannual timescales. Earlier studies used more dynamically based methods and combined SLA with the vertical structure of horizontal flow from mooring data and found that the SLA had little skill at capturing the upper mid-ocean transport (e.g. Hirschi et al.2009; Kanzow et al.2009; Szuts et al.2012). The goal of this study was to similarly use satellite altimetry as a proxy for AMOC variability and re-evaluate dynamics-based methods (e.g. Hirschi et al.2009; Kanzow et al.2009; Szuts et al.2012) for using satellite altimetry to estimate the TUMO as well as the TGS components of the AMOC, which could eventually be used at other latitudes. Thus, using principles of geostrophy and normal-mode decomposition a method for constructing the upper mid-ocean transport, TUMO*, and by extension the TMOC* at 26 N from satellite altimetry on low-frequency (18-month) timescales was devised here by combining the first baroclinic mode, derived from time-averaged density profiles from the RAPID moorings, with SLA to reproduce the TUMO* and TMOC* transports. Using this new method, we find that 56 % of RAPID TUMO variability could be captured and 69 % of the RAPID TMOC variability on interannual (18-month) timescales, which is significantly more than at the shorter timescales studied by Kanzow et al. (2009). Sensitivity testing using different filtering windows decreases variance captured by TUMO* and TMOC* to 18 % of TUMO and 50 % of the TMOC, respectively, when using 6-month Gaussian smoothing; the values are 45 % of TUMO and 66 % of the TMOC when using 12-month Gaussian smoothing, with an increase to 62 % of TUMO and 72 % of the TMOC when using 24-month Gaussian smoothing. The TGS component of TMOC was also reproduced separately using the satellite altimetry; we found that it captures 49 % of the GS variability via the Florida Straits as measured from the telephone cable (Meinen et al.2010) on interannual timescales, similar to a recent study by Volkov et al. (2020), who also used satellite data to estimate GS variability in the Florida Straits. Finally, the satellite-derived TUMO* and TMOC* are reconstructed for the full satellite period (1993 to 2018). The 26-year reconstructions of TUMO* and TMOC* were compared with another separate satellite-derived proxy (EFW15) and a TMOC reconstruction from the GloSea5 reanalysis. Though all reconstructions had correlations of at least r=0.79 with RAPID-based transports for the 2004 to 2013 time period, the 1993–2003 time period showed poor agreement between the three TMOC reconstructions, suggesting that further studies are needed to confidently reproduce the fluctuations in the AMOC during times pre-dating the RAPID programme.

Because the altimetry has no knowledge of vertical shear, estimating transport over the upper 1000 m using only satellite altimetry results in an overestimation of the transport's magnitude (Fig. S3). Thus, normal-mode decomposition was investigated using the RAPID moorings to answer two questions: (a) how much does the vertical structure of the flow contribute to the upper mid-ocean transport (TUMO) variability, and (b) can the pressure modes be combined with the satellite to provide an improved way of estimating mass overturning transport? In the first instance, we find that the first baroclinic mode accounts for 83 % of the observed interior geostrophic transport variability, and the combined barotropic and first baroclinic modes account for 98 % of the total variability. In the second instance, we find that combining the satellite altimetry with the vertical structure (from the first baroclinic mode) improves the magnitude of values for the altimetry-based transport. However, a scale factor is still needed to further correct the values to capture the magnitude of the TUMO. We posit that the need for a scale factor is because of the proximity of moorings to land. Close to the coast the SSH signal is not so well resolved by altimetry and is significantly affected by the barotropic mode as well as the baroclinic mode (Kanzow et al.2009). It is also important to note that combining the altimetry with the vertical structure (from the first baroclinic mode) does not improve the amount of variance captured by the satellite-based TUMO* compared to RAPID TUMO, and using time-varying pressure modes (instead of time-averaged) does not necessarily improve correlation between the RAPID TMOC and the satellite-derived TMOC*. Further discrepancies between satellite-based transport and RAPID transport may be due to the barotropic component, which the satellite-based method does not account for.

To summarise, principles of geostrophy can be used with satellite altimetry to effectively capture upper mid-ocean and GS transport and, by extension, AMOC transport. The vertical structure of the flow (from the moorings) does not improve the amount of variability the satellite-based method captures; however, the analysis of the normal-mode decomposition yields insight into the governing modes of variability associated with the TUMO on interannual timescales and the limits this places on any satellite-based transport method. While this method progresses attempts at a dynamically robust method for estimating AMOC transport using satellite altimetry, which can also be used to recreate the AMOC at the start of the satellite time period pre-dating RAPID, the method is still not independent from the mooring array, and without improved understanding of the scale factor, further investigation is needed to produce an altimetry-based method for the AMOC that can be used at other latitudes.

Data availability

The RAPID-MOCHA-WBTS time series (Moat et al.2020a) is available at

The Copernicus Marine Environment Monitoring Service for data access online is available at _047 (last access: 10 December 2020).

ERA5 wind stress is available via!/home (Copernicus Climate Change Service2017). The Gulf Stream cable data are available via The GloSea5 time series is available from Jackson et al. (2019) upon request.


The supplement related to this article is available online at:

Author contributions

ASF wrote the paper with input and contributions from all authors. ASF, EFW, DAS, and BIM contributed to the analysis. EFW, BIM, and DAS contributed to data collection and provided processed data.

Competing interests

The authors declare that they have no conflict of interest.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors thank the reviewers for their helpful comments. The authors also thank the many officers, crew, and technicians who helped to collect these data. Alejandra Sanchez-Franks also thanks Louis Clément for helpful discussions on normal-mode decomposition.

Financial support

This research has been supported by grants from the UK Natural Environment Research Council for the RAPID-AMOC programme and the ACSIS programme (grant no. NE/N018044/1) as well as funding from the European Union Horizon 2020 Research and Innovation programme BLUE-ACTION (grant no. 727852).

Review statement

This paper was edited by Trevor McDougall and reviewed by two anonymous referees.


Baringer, M. O. and Larsen, J. C.: Sixteen years of Florida Current transport at 27 N, Geophys. Res. Lett., 28, 3179–3182, 2001. a, b, c, d

Bingham, R. J. and Hughes, C. W.: Geostrophic dynamics of meridional transport variability in the subpolar North Atlantic, J. Geophys. Res.-Oceans, 114, L02603,, 2009. a

Bingham, R. J., Hughes, C. W., Roussenov, V., and Williams, R. G.: Meridional coherence of the North Atlantic meridional overturning circulation, Geophys. Res. Lett., 34, L23606,, 2007. a

Bisagni, J. J., Gangopadhyay, A., and Sanchez-Franks, A.: Secular change and inter-annual variability of the Gulf Stream position, 1993–2013, 70°–55° W, Deep-Sea Res. Pt. I, 125, 1–10, 2017. a

Blockley, E. W., Martin, M. J., McLaren, A. J., Ryan, A. G., Waters, J., Lea, D. J., Mirouze, I., Peterson, K. A., Sellar, A., and Storkey, D.: Recent development of the Met Office operational ocean forecasting system: an overview and assessment of the new Global FOAM forecasts, Geosci. Model Dev., 7, 2613–2638,, 2014. a

Bryden, H. L., Mujahid, A., Cunningham, S. A., and Kanzow, T.: Adjustment of the basin-scale circulation at 26 N to variations in Gulf Stream, deep western boundary current and Ekman transports as observed by the Rapid array, Ocean Sci., 5, 421–433,, 2009. a

Cabanes, C., Lee, T., and Fu, L.-L.: Mechanisms of interannual variations of the meridional overturning circulation of the North Atlantic Ocean, J. Phys. Oceanogr., 38, 467–480, 2008. a

Caesar, L., Rahmstorf, S., Robinson, A., Feulner, G., and Saba, V.: Observed fingerprint of a weakening Atlantic Ocean overturning circulation, Nature, 556, 191–196, 2018. a

Chelton, D. B., Schlax, M. G., Samelson, R. M., and de Szoeke, R. A.: Global observations of large oceanic eddies, Geophys. Res. Lett., 34, L15606,, 2007. a

Chidichimo, M. P., Kanzow, T., Cunningham, S. A., Johns, W. E., and Marotzke, J.: The contribution of eastern-boundary density variations to the Atlantic meridional overturning circulation at 26.5° N, Ocean Sci., 6, 475–490,, 2010. a, b

Clément, L., Frajka‐Williams, E., Szuts, Z. B., and Cunningham, S. A.: Vertical structure of eddies and R ossby waves, and their effect on the A tlantic meridional overturning circulation at 26.5° N, J. Geophys. Res.-Oceans, 119, 6479–6498, 2014. a, b, c, d, e, f, g

Copernicus Climate Change Service (C3S): ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate, Copernicus Climate Change Service Climate Data Store (CDS), available at:!/home (last access: 7 September 2020), 2017. a

Cunningham, S. A., Kanzow, T., Rayner, D., Baringer, M. O., Johns, W. E., Marotzke, J., Longworth, H. R., Grant, E. M., Hirschi, J. J.-M., Beal, L. M., Meinen, C. S., and Bryden, H. L.: Temporal variability of the Atlantic meridional overturning circulation at 26.5 N, Science, 317, 935–938, 2007. a

Emery, W. J. and Thomson, R. E. (Eds.): Time-series analysis methods, Data Analysis Methods in Physical Oceanography, Elsevier Science: Amsterdam, The Netherlands, chap. 5, 371–567, 2001. a

Evans, D. G., Toole, J., Forget, G., Zika, J. D., Naveira Garabato, A. C., Nurser, A. J. G., and Yu, L.: Recent Wind-Driven Variability in Atlantic Water Mass Distribution and Meridional Overturning Circulation, J. Phys. Oceanogr., 47, 633–647,, 2017. a

Frajka-Williams, E., Ansorge, I. J., Baehr, J., Bryden, H. L., Chidichimo, M. P., Cunningham, S. A., Danabasoglu, G., Dong, S., Donohue, K. A., and Elipot, S.: Atlantic meridional overturning circulation: Observed transport and variability, Front. Mar. Sci., 6, 260,, 2019. a

Frajka‐Williams, E.: Estimating the Atlantic overturning at 26 N using satellite altimetry and cable measurements, Geophys. Res. Lett., 42, 3458–3464, 2015. a, b, c, d, e, f, g

Fratantoni, P. S. and Pickart, R. S.: The western North Atlantic shelfbreak current system in summer, J. Phys. Oceanogr., 37, 2509–2533, 2007. a

Fu, L., Christensen, E. J., Yamarone, C. A., Lefebvre, M., Menard, Y., Dorrer, M., and Escudier, P.: TOPEX/POSEIDON mission overview, J. Geophys. Res.-Oceans, 99, 24369–24381, 1994. a

Gill, A. E.: Atmosphere-Ocean dynamics (International Geophysics Series), Academic Press, 1982. a, b, c

Hall, M. M. and Bryden, H. L.: Direct estimates and mechanisms of ocean heat transport, Deep-Sea Res. Pt. A, 29, 339–359, 1982. a

Herrford, J., Brandt, P., Kanzow, T., Hummels, R., Araujo, M., and Durgadoo, J. V.: Seasonal variability of the Atlantic Meridional Overturning Circulation at 11° S inferred from bottom pressure measurements, Ocean Sci., 17, 265–284,, 2021. a

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz‐Sabater, J., Nicolas, J., Peubey, C., Radu, R., and Schepers, D.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, 2020. a

Hirschi, J. J. M., Killworth, P. D., and Blundell, J. R.: Subannual, seasonal, and interannual variability of the North Atlantic meridional overturning circulation, J. Phys. Oceanogr., 37, 1246–1265, 2007. a, b, c, d, e, f, g, h, i, j, k

Hirschi, J. J. M., Killworth, P. D., Blundell, J. R., and Cromwell, D.: Sea surface height signals as indicators for oceanic meridional mass transports, J. Phys. Oceanogr., 39, 581–601, 2009. a, b, c, d, e, f, g, h, i, j, k

Hummels, R., Brandt, P., Dengler, M., Fischer, J., Araujo, M., Veleda, D., and Durgadoo, J. V.: Interannual to decadal changes in the western boundary circulation in the Atlantic at 11 S, Geophys. Res. Lett., 42, 7615–7622, 2015. a

Jackson, L. C., Peterson, K. A., Roberts, C. D., and Wood, R. A.: Recent slowing of Atlantic overturning circulation as a recovery from earlier strengthening, Nat. Geosci., 9, 518–522, 2016. a

Jackson, L. C., Dubois, C., Forget, G., Haines, K., Harrison, M., Iovino, D., Köhl, A., Mignac, D., Masina, S., Peterson, K. A., Piecuch, C. G., Roberts, C. D., Robson, J., Storto, A., Toyoda, T., Valdivieso, M., Wilson, C., Wang, Y., and Zuo, H.: The Mean State and Variability of the North Atlantic Circulation: A Perspective From Ocean Reanalyses, J. Geophys. Res.-Oceans, 124, 9141–9170,, 2019. a, b, c

Johns, W. E., Kanzow, T., and Zantopp, R.: Estimating ocean transports with dynamic height moorings: An application in the Atlantic Deep Western Boundary Current at 26 N, Deep-Sea Res. Pt. I, 52, 1542–1567, 2005. a

Johns, W. E., Baringer, M. O., Beal, L., Cunningham, S., Kanzow, T., Bryden, H. L., Hirschi, J., Marotzke, J., Meinen, C., Shaw, B., and Curry, R.: Continuous, array-based estimates of Atlantic Ocean heat transport at 26.5 N, J. Climate, 24, 2429–2449, 2011. a, b

Johnson, H. L. and Marshall, D. P.: A theory for the surface Atlantic response to thermohaline variability, J. Phys. Oceanogr., 32, 1121–1132, 2002. a

Kanzow, T., Cunningham, S. A., Rayner, D., Hirschi, J. J.-M., Johns, W. E., Baringer, M. O., Bryden, H. L., Beal, L. M., Meinen, C. S., and Marotzke, J.: Observed flow compensation associated with the MOC at 26.5 N in the Atlantic, Science, 317, 938–941, 2007. a, b

Kanzow, T., Send, U., and McCartney, M.: On the variability of the deep meridional transports in the tropical North Atlantic, Deep-Sea Res. Pt. I, 55, 1601–1623, 2008. a, b, c

Kanzow, T., Johnson, H. L., Marshall, D. P., Cunningham, S. A., Hirschi, J.-M., Mujahid, A., Bryden, H. L., and Johns, W. E.: Basinwide integrated volume transports in an eddy-filled ocean, J. Phys. Oceanogr., 39, 3091–3110, 2009. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Kanzow, T., Cunningham, S. A., Johns, W. E., Hirschi, J. J. M., Marotzke, J., Baringer, M. O., Meinen, C. S., Chidichimo, M. P., Atkinson, C., and Beal, L. M.: Seasonal variability of the Atlantic meridional overturning circulation at 26.5 N, J. Climate, 23, 5678–5698, 2010. a, b, c

Kersalé, M., Meinen, C. S., Perez, R. C., Le Henaff, M., Valla, D., Lamont, T., Sato, O. T., Dong, S., Terre, T., van Caspel, M. Chidichimo, M. P., van den Berg, M., Speich, S., Piola, A. R., Campos, E. J. D., Ansorge, I., Volkov, D. L., Lumpkin, R., and Garzoli, S.: Highly variable upper and abyssal overturning cells in the South Atlantic, Sci. Adv., 6, eaba7573,, 2020. a

Killworth, P. D. and Blundell, J. R.: Long extratropical planetary wave propagation in the presence of slowly varying mean flow and bottom topography. Part II: Ray propagation and comparison with observations, J. Phys. Oceanogr., 33, 802–821, 2003. a

Lozier, M. S., Li, F., Bacon, S., Bahr, F., Bower, A. S., Cunningham, S. A., de Jong, M. F., de Steur, L., deYoung, B., Fischer, J., Gary, S. F., Greenan, B. J. W., Holliday, N. P., Houk, A., Houpert, L., Inall, M. E., Johns, W. E., Johnson, H. L., Johnson, C., Karstensen, J., Koman, G., Le Bras, I. A., Lin, X., Mackay, N., Marshall, D. P., Mercier, H., Oltmanns, M., Pickart, R. S., Ramsey, A. L., Rayner, D., Straneo, F., Thierry, V., Torres, D. J., Williams, R. G., Wilson, C., Yang, J., Yashayaev, I., and Zhao, J.: A sea change in our view of overturning in the subpolar North Atlantic, Science, 363, 516–521, 2019. a

MacLachlan, C., Arribas, A., Peterson, K. A., Maidens, A., Fereday, D., Scaife, A. A., Gordon, M., Vellinga, M., Williams, A., and Comer, R. E.: Global Seasonal forecast system version 5 (GloSea5): A high‐resolution seasonal forecast system, Q. J. Roy. Meteor. Soc., 141, 1072–1084, 2015. a

McCarthy, G. D., Smeed, D. A., Johns, W. E., Frajka-Williams, E., Moat, B. I., Rayner, D., Baringer, M. O., Meinen, C. S., Collins, J., and Bryden, H. L.: Measuring the Atlantic meridional overturning circulation at 26 N, Prog. Oceanogr., 130, 91–111, 2015. a, b, c, d, e, f, g, h

Megann, A., Storkey, D., Aksenov, Y., Alderson, S., Calvert, D., Graham, T., Hyder, P., Siddorn, J., and Sinha, B.: GO5.0: the joint NERC–Met Office NEMO global ocean model for use in coupled and forced applications, Geosci. Model Dev., 7, 1069–1092,, 2014. a

Meinen, C. S., Baringer, M. O., and Garcia, R. F.: Florida Current transport variability: An analysis of annual and longer-period signals, Deep-Sea Res. Pt. I, 57, 835–846, 2010. a, b

Meinen, C. S., Speich, S., Piola, A. R., Ansorge, I., Campos, E., Kersalé, M., Terre, T., Chidichimo, M. P., Lamont, T., Sato, O. T., Perez, R. C., Valla, D., van den Berg, M., Le Hénaff, M., Dong, S., and Garzoli, S. L.: Meridional overturning circulation transport variability at 34.5 S during 2009–2017: Baroclinic and barotropic flows and the dueling influence of the boundaries, Geophys. Res. Lett., 45, 4180–4188, 2018. a

Meinen, C. S., Smith, R. H., and Garcia, R. F.: Evaluating pressure gauges as a potential future replacement for electromagnetic cable observations of the Florida Current transport at 27 N, J. Oper. Oceanogr.,, online first, 2020. a

Moat, B. I., Sinha, B., Josey, S. A., Robson, J., Ortega, P., Sévellec, F., Holliday, N. P., McCarthy, G. D., New, A. L., and Hirschi, J.-M.: Insights into decadal North Atlantic sea surface temperature and ocean heat content variability from an eddy-permitting coupled climate model, J. Climate, 32, 6137–6161, 2019. a

Moat, B. I., Frajka-Williams, E., Smeed, D., Rayner, D., Sanchez-Franks, A., Johns, W. E., Baringer, M. O., Volkov, D., and Collins, J.: Atlantic meridional overturning circulation observed by the RAPID-MOCHA-WBTS (RAPID-Meridional Overturning Circulation and Heatflux Array-Western Boundary Time Series) array at 26N from 2004 to 2018 (v2018.2), British Oceanographic Data Centre, National Oceanography Centre [data set], NERC, UK,, 2020a. a, b, c

Moat, B. I., Smeed, D. A., Frajka-Williams, E., Desbruyères, D. G., Beaulieu, C., Johns, W. E., Rayner, D., Sanchez-Franks, A., Baringer, M. O., Volkov, D., Jackson, L. C., and Bryden, H. L.: Pending recovery in the strength of the meridional overturning circulation at 26 N, Ocean Sci., 16, 863–874,, 2020b. a, b

Peña-Molino, B. and Joyce, T.: Variability in the slope water and its relation to the Gulf Stream path, Geophys. Res. Lett., 35, L03606,, 2008. a

Pérez-Hernández, M. D. and Joyce, T. M.: Two modes of Gulf Stream variability revealed in the last two decades of satellite altimeter data, J. Phys. Oceanogr., 44, 149–163, 2014. a

Petrie, B. and Drinkwater, K.: Temperature and salinity variability on the Scotian Shelf and in the Gulf of Maine 1945–1990, J. Geophys. Res.-Oceans, 98, 20079–20089, 1993. a

Pohlmann, A., Fricke, W. F., Reinecke, F., Kusian, B., Liesegang, H., Cramm, R., Eitinger, T., Ewering, C., Pötter, M., and Schwartz, E.: Genome sequence of the bioplastic-producing “Knallgas” bacterium Ralstonia eutropha H16, Nat. Biotechnol., 24, 1257–1262, 2006. a

Rayner, D., Hirschi, J. J.-M., Kanzow, T., Johns, W. E., Wright, P. G., Frajka-Williams, E., Bryden, H. L., Meinen, C. S., Baringer, M. O., and Marotzke, J.: Monitoring the Atlantic meridional overturning circulation, Deep-Sea Res. Pt. II, 58, 1744–1753, 2011. a, b, c, d

Roberts, C. D., Waters, J., Peterson, K. A., Palmer, M. D., McCarthy, G. D., Frajka-Williams, E., Haines, K., Lea, D. J., Martin, M. J., Storkey, D., Blockley, E. W., and Zuo, H.: Atmosphere drives recent interannual variability of the Atlantic meridional overturning circulation at 26.5°N, Geophys. Res. Lett., 40, 5164–5170,, 2013. a

Sanchez-Franks, A. and Zhang, R.: Impact of the Atlantic meridional overturning circulation on the decadal variability of the Gulf Stream path and regional chlorophyll and nutrient concentrations, Geophys. Res. Lett., 42, 9889–9887,, 2015. a, b

Sanchez-Franks, A., Flagg, C., and Rossby, T.: A comparison of transport and position between the Gulf Stream east of Cape Hatteras and the Florida Current, J. Mar. Res., 72, 291–306, 2014. a

Sanchez-Franks, A., Hameed, S., and Wilson, R. E.: The Icelandic low as a predictor of the Gulf Stream north wall position, J. Phys. Oceanogr., 46, 817–826,, 2016. a

Siegel, D. A., McGillicuddy Jr., D. J., and Fields, E. A.: Mesoscale eddies, satellite altimetry, and new production in the Sargasso Sea, J. Geophys. Res.-Oceans, 104, 13359–13379, 1999.  a

Smeed, D. A., McCarthy, G. D., Cunningham, S. A., Frajka-Williams, E., Rayner, D., Johns, W. E., Meinen, C. S., Baringer, M. O., Moat, B. I., Duchez, A., and Bryden, H. L.: Observed decline of the Atlantic meridional overturning circulation 2004–2012, Ocean Sci., 10, 29–38,, 2014. a, b, c

Smeed, D. A., Josey, S. A., Beaulieu, C., Johns, W. E., Moat, B. I., Frajka‐Williams, E., Rayner, D., Meinen, C. S., Baringer, M. O., and Bryden, H. L.: The North Atlantic Ocean is in a state of reduced overturning, Geophys. Res. Lett., 45, 1527–1533, 2018. a

Srokosz, M., Baringer, M., Bryden, H., Cunningham, S., Delworth, T., Lozier, S., Marotzke, J., and Sutton, R.: Past, present, and future changes in the Atlantic meridional overturning circulation, B. Am. Meteorol. Soc., 93, 1663–1676, 2012. a

Szuts, Z. B., Blundell, J. R., Chidichimo, M. P., and Marotzke, J.: A vertical-mode decomposition to investigate low-frequency internal motion across the Atlantic at 26 N, Ocean Sci., 8, 345–367,, 2012. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Taylor, A. H. and Stephens, J. A.: The North Atlantic oscillation and the latitude of the Gulf Stream, Tellus A, 50, 134–142, 1998. a

Trenberth, K. E. and Caron, J. M.: Estimates of meridional atmosphere and ocean heat transports, J. Climate, 14, 3433–3443, 2001. a

Volkov, D. L., Domingues, R., Meinen, C. S., Garcia, R., Baringer, M., Goni, G., and Smith, R. H.: Inferring Florida Current volume transport from satellite altimetry, J. Geophys. Res.-Oceans, 125, e2020JC016763, e2020JC016763,, 2020. a, b, c

Willis, J. K.: Can in situ floats and satellite altimeters detect long‐term changes in Atlantic Ocean overturning?, Geophys. Res. Lett., 37, L06602,, 2010. a

Wunsch, C.: The vertical partition of oceanic horizontal kinetic energy, J. Phys. Oceanogr., 27, 1770–1794, 1997. a, b, c

Wunsch, C. and Stammer, D.: Atmospheric loading and the oceanic “inverted barometer” effect, Rev. Geophys., 35, 79–107, 1997. a

Zhang, R.: Coherent surface‐subsurface fingerprint of the Atlantic meridional overturning circulation, Geophys. Res. Lett., 35, L20705,, 2008. a

Zhang, R., Sutton, R., Danabasoglu, G., Kwon, Y., Marsh, R., Yeager, S. G., Amrhein, D. E., and Little, C. M.: A review of the role of the Atlantic meridional overturning circulation in Atlantic multidecadal variability and associated climate impacts, Rev. Geophys., 57, 316–375, 2019. a, b, c, d

Zhao, J. and Johns, W.: Wind-forced interannual variability of the Atlantic Meridional Overturning Circulation at 26.5°N, J. Geophys. Res.-Oceans, 119, 2403–2419,, 2014. a, b

Short summary
In the North Atlantic, ocean currents carry warm surface waters northward and return cooler deep waters southward. This type of ocean circulation, known as overturning, is important for the Earth’s climate. This overturning has been measured using a mooring array at 26° N in the North Atlantic since 2004. Here we use these mooring data and global satellite data to produce a new method for monitoring the overturning over longer timescales, which could potentially be applied to different latitudes.