Articles | Volume 19, issue 2
Biogeosciences, 19, 437–454, 2022
Biogeosciences, 19, 437–454, 2022

Research article 28 Jan 2022

Research article | 28 Jan 2022

Oxygen export to the deep ocean following Labrador Sea Water formation

Oxygen export to the deep ocean following Labrador Sea Water formation
Jannes Koelling1, Dariia Atamanchuk1, Johannes Karstensen2, Patricia Handmann2, and Douglas W. R. Wallace1 Jannes Koelling et al.
  • 1Department of Oceanography, Dalhousie University, Halifax, Nova Scotia, Canada
  • 2GEOMAR Helmholtz Centre for Ocean Research Kiel, Kiel, Germany

Correspondence: Jannes Koelling (


The Labrador Sea in the North Atlantic Ocean is one of the few regions globally where oxygen from the atmosphere can reach the deep ocean directly. This is the result of wintertime deep convection, which homogenizes the water column to a depth of up to 2000 m and brings deep water undersaturated in oxygen into contact with the atmosphere. In this study, we analyze how the intense oxygen uptake during Labrador Sea Water (LSW) formation affects the properties of the outflowing deep western boundary current, which ultimately feeds the upper part of the North Atlantic Deep Water layer in much of the Atlantic Ocean.

Seasonal cycles of oxygen concentration, temperature, and salinity from a 2-year time series collected by sensors moored at 600 m nominal depth in the outflowing boundary current at 53 N show a cooling, freshening, and increase in oxygen content of the water flowing out of the basin between March and August. Analysis of Argo float data suggests that this is preceded by an increased input of LSW into the boundary current about 1 month earlier. This input is the result of newly ventilated LSW entering from the interior, as well as LSW formed directly within the boundary current. Together, these results imply that the southward export of newly formed LSW primarily occurs in the months following the onset of deep convection, from March to August, and that this direct LSW export route controls the seasonal oxygen increase in the outflow at 600 m depth. During the rest of the year, properties of the boundary current measured at 53 N resemble those of Irminger Water, which enters the basin with the boundary current from the Irminger Sea.

The input of newly ventilated LSW increases the oxygen concentration from 298 µmol L−1 in January to a maximum of 306 µmol L−1 in April. As a result of this LSW input, an estimated (1.60 ± 0.42) ×1012mol yr−1 of oxygen are added to the outflowing boundary current, mostly during spring and summer, equivalent to 50 % of the wintertime uptake from the atmosphere in the interior of the basin. The export of oxygen from the subpolar gyre associated with this direct southward pathway of LSW is estimated to supply 42 %–71 % of the oxygen consumed annually in the upper North Atlantic Deep Water layer in the Atlantic Ocean between the Equator and 50 N. Our results show that the formation of LSW is important for replenishing oxygen to the deep oceans, meaning that possible changes in its formation rate and ventilation due to climate change could have wide-reaching impacts on marine life.

1 Introduction

Much of the global supply of oxygen to the deep ocean is concentrated in a few key regions where near-surface water sinks to great depth and spreads away from its source region (Talley2008; Gebbie and Huybers2011; Khatiwala et al.2012). This process ventilates the deep ocean, supplying oxygen to a vast volume of water that would otherwise be barren and making it capable of sustaining life (Rogers2015; Isozaki1997). One of the regions where such deep ocean ventilation occurs is the Labrador Sea, a semi-enclosed marginal sea nestled between eastern Canada and western Greenland. Here, strong wintertime atmospheric cooling leads to buoyancy loss of surface waters that lay above a weakly stratified water column, causing convective overturning that homogenizes the upper 1000–2000 m of the ocean (Marshall and Schott1999; Yashayaev2007). Due to the great depths reached by convection in this region, it is generally referred to as deep convection in order to differentiate it from the much shallower-reaching convection occurring in mixed layers throughout the global ocean, but we will use the terms convection and deep convection interchangeably hereinafter in the interest of brevity. During the convection season, which typically lasts from January to April, deep water masses that are low in oxygen are continuously incorporated into the progressively deepening mixed layer, which leads to a decrease in near-surface oxygen. This results in severe air–sea gradients in oxygen concentration, which, together with extreme atmospheric conditions, drive intense uptake of oxygen during winter (Koelling et al.2017; Wolf et al.2018; Atamanchuk et al.2020).

The climatological flow field in the Labrador Sea features a cyclonic boundary current entering at the southern tip of Greenland and exiting at the southwestern end of the Labrador Sea, with only a weak mean flow in the interior (Fig. 1). At mid-depth, the boundary current entering the basin carries Irminger Water (IW), a water mass originating in the Atlantic ocean and modified in the subpolar gyre, particularly the Irminger Sea (Cuny et al.2002; Pacini et al.2020). The cyclonic circulation in the Labrador Sea is linked to a doming of isopycnals in the center of the basin, which supports deep convection in the interior (Marshall and Schott1999). Different types of eddies, such as Irminger rings and convective eddies, compensate the annual mean surface heat loss in the interior and help restratify the water column above the water mass formed during convection, known as Labrador Sea Water (LSW) (Eden and Böning2002; Straneo2006; de Jong et al.2014; Rieck et al.2019). Variability in atmospheric forcing, lateral fluxes of heat and freshwater, and the evolution of the mixed patch are main drivers for interannual variability in the depth of convection and in properties of LSW (Lazier1973; Yashayaev and Loder2016). Much of the research has been focused on the deep convection region in the interior of the basin, but it has been shown that convection can also occur within the boundary current itself (Pickart et al.1997; Palter et al.2008). The densest water masses formed in the boundary current are similar to the “classical” LSW formed in the interior, while convection over the continental slope forms a lighter water mass known as upper Labrador Sea Water or uLSW (Pickart et al.2002; Cuny et al.2005). While the relative importance of boundary current convection is still uncertain, some model studies have suggested that it could account for a substantial fraction of the LSW that is eventually exported out of the region (Brandt et al.2007; MacGilchrist et al.2020).

The LSW formed each winter in the center of the basin is distinguishable in summertime hydrographic sections as a cold, fresh water mass with high oxygen content, compared to the warmer, saltier IW found in the boundary current near Greenland (see Yashayaev and Loder2016, and Zou et al.2020, for sections across the Labrador Sea). LSW is one of the main water masses making up North Atlantic Deep Water (NADW) (Lazier1973), and it spreads throughout the Atlantic Ocean, both as part of the Deep Western Boundary Current (Molinari et al.1998; Toole et al.2017) and in the interior of the basin (Bower et al.2009; Lozier2010). The signature of LSW is distinguishable in NADW properties far from the source region (Talley and McCartney1982; Rhein et al.2004; Le Bras et al.2017), including elevated oxygen concentrations compared to adjacent water masses (Atkinson et al.2012). Although the importance of ventilation in the Labrador Sea on NADW properties is readily apparent from hydrographic data, the exact timing and mechanisms of the southward export of LSW from the subpolar gyre are not well established. Due to the weak mean currents in the center of the basin (see Fig. 1), some studies have suggested that much of the newly formed LSW remains in the interior for several years, with only a small fraction exported in the months following convection (Rhein et al.2002; Straneo et al.2003; MacGilchrist et al.2021). On the other hand, there is evidence of a rapid export pathway of newly formed LSW which may be associated with convection either within or close to the boundary current (Pickart et al.1997; Brandt et al.2007). An analysis from float data by Palter et al. (2008) showed that both interior and boundary convection can affect the properties of the outflowing boundary current and suggested that eddies may play an important role in the input of newly formed LSW into the boundary current.

Figure 1Mean salinity between 500 and 1000 m in the Labrador Sea, calculated from all available Argo float profiles between 2000 and 2020 (see Sect. 2.2 for details), averaged within overlapping 0.25×0.25 bins. Gray lines show the 1000, 2000, and 3000 m isobaths, vectors show mean currents at the depth of LSW from Fischer et al. (2018), and symbols correspond to positions of the SeaCycler mooring (red) and the 53 N array (magenta).

In this study, we use a novel dataset including moored oxygen concentration measurements from the outflowing boundary current at the southern exit of the Labrador Sea. The data are recorded at the 53 N array (Zantopp et al.2017), which is part of the Overturning in the Subpolar North Atlantic Program (OSNAP) (Lozier et al.2017). We analyze annual cycles of oxygen, temperature, and salinity, along with data from Argo floats, in order to investigate seasonal variations in the input of LSW into the boundary current and its export out of the basin. Our results highlight how changes in properties of the outflowing boundary current over a seasonal cycle, and resulting changes in the export of oxygen, are linked to the formation and export of newly ventilated LSW.

2 Data and methods

2.1 Moored sensor data

Data used in this study were collected from May 2016 to May 2018 on the 53 N array in the boundary current at the exit of the Labrador Sea on moorings K7 (52.86 N, 51.31 W), K8 (52.96 N, 51.31 W), K9 (53.14 N, 50.87 W), and K10 (53.39 N, 50.25 W) (west to east in Fig. 1). Moorings have been deployed at the location since 1997 but in various configurations (Zantopp et al.2017). In the present setup, each mooring is equipped with Aquadopp and RCM current meters and SeaBird SBE 37 instruments measuring conductivity, temperature (T), and pressure, which were also used to derive salinity (S). The depths and positions of the instruments are optimized to measure the strength and properties of the boundary current exiting the Labrador Sea, including transports of different layers (Zantopp et al.2017). Since 2016, Aanderaa 4330 optodes (Aanderaa Data Instruments, AS) have been deployed at select depths alongside TS sensors to measure dissolved oxygen (O2) concentrations (Fig. 2a). We also use data from an Aanderaa 4330 oxygen optode mounted at 500 m depth on the SeaCycler mooring in the central Labrador Sea (see Fig. 1 for location). For the SeaCycler optode, no co-located salinity data are available, and we estimate salinity used in Sect. 3.2 from a climatological TS relationship at 500 m calculated from Argo float data.

Figure 2Oxygen and salinity sections along the array from shipboard measurements collected during four mooring deployment and recovery cruises (see Sect. 2.2). Gray symbols show stations from one representative cruise. Green symbols in (a) show deployment locations of the oxygen sensors, which were co-located with TS sensors. Only data from the instruments at 600 m depth are used in this study. The full coverage of current meters and TS sensors is described in Zantopp et al. (2017). Black contours show the potential density range of 27.68 kg m−3σθ 27.8 kg m−3 for LSW and the σθ=27.74kg m−3 isopycnal which delineates the upper Labrador Sea Water (uLSW) above and classical Labrador Sea Water (cLSW) below.


The deployment depths of the oxygen sensors on the 53 N array are shown in Fig. 2 along with mean oxygen and salinity sections along the array recorded with a conductivity–temperature–depth (CTD) rosette. The majority of the O2 sensors were deployed at nominal depths of approximately 600 m or 1900 m, chosen to represent the LSW core, and a depth near the bottom of the LSW layer. Only the sensors at 600 m depth are used in this study, and the deployment and calibration information for these sensors is summarized in Table 1. LSW is commonly defined in density space by a range of potential density σθ, with a lower boundary of 27.8 kg m−3 and an upper boundary of either 27.68 kg m−3 (Pickart et al.1997; Zantopp et al.2017; Lozier et al.2019) or 27.7 kg m−3 (Zou et al.2020), the former of which is used as the definition of the LSW layer in this study. Within the LSW layer, oxygen concentrations are elevated above 1200 m, and the mean salinity is 34.86. The low-oxygen core below LSW centered at about 2000 m depth coincides with a salinity maximum, which is indicative of Northeast Atlantic Deep Water (NEADW) (Yashayaev2007). Some authors differentiate between “upper” and “classical” LSW (Pickart et al.1997), with the boundary between the two water masses typically at a potential density of σθ= 27.74 kg m−3 (Kieke et al.2006). The mean potential density at 600 m depth is 27.72 kg m−3 for K7 and 27.74 kg m−3 for K8–K10 (Table 1), suggesting that sensors on these three moorings are near the interface between uLSW and cLSW (classical Labrador Sea Water; see also Fig. 2).

The mean salinity of the boundary current in the LSW depth range decreases as it flows cyclonically around the basin (Fig. 1). At the 53 N array, the salinity is closer to values found in the interior Labrador Sea compared to the Greenland side, consistent with the view that there is significant input of LSW from the interior. However, the water here is still warmer and saltier than that found in the convective interior, suggesting that it is a mixture of LSW and IW, which enters the Labrador Sea on the northern side at a mean core depth of about 500 m and can be seen as a local salinity maximum near this depth (Pacini et al.2020). In this study, we will focus on the variability in properties in this LSW layer using data from the sensors mounted at about 600 m on the moorings, representing the full width of the boundary current from the shelf break area (K7, K8), via the core (K9), into the outer edge and recirculation regime (K10).

Table 1Deployment locations, depths, and drifts of the oxygen optodes during 2016–2018.

Download Print Version | Download XLSX

2.1.1 Sensor calibration

Temperature and salinity sensors were calibrated following the procedure outlined in Karstensen (2005). Briefly, the method involves attaching the sensors to a CTD rosette prior to and after each deployment and comparing the data from the different instruments at certain stop depths with the CTD recording, which is calibrated with discrete bottle samples taken throughout the cruise.

Oxygen optodes were supplied with individual multipoint factory calibration, with an absolute accuracy of 1.5 % or 2 µmol L−1 (Tengberg et al.2006; Tengberg and Hovdenes2014), and powered by single-channel loggers (RBR Ltd.). The optodes were further calibrated against Winkler samples collected during CTD casts at the mooring locations during deployment and recovery cruises in 2016 and 2018 to correct for sensor drift. The drift of the sensors during deployment ranged from 2.5 to 6.2 µmol L−1 or 0.4 % yr−1–1 % yr−1 (see Table 1).

Figure 3Example of Argo data used to track LSW input into the boundary current after convection from float ID 6902589. (a) Trajectory showing latest deep mixed layer measurement (blue star) and location of entry into boundary current (red star). Thick green line highlights the part of the trajectory from LSW formation to crossing the 3000 m isobath. (b) Time series of mixed layer depth and maximum layer thickness, with green shading highlighting the time from convection in the interior to entering the boundary current. (c) Consecutive profiles of layer thickness in density space, going from convection (dark blue) to input into boundary current (light green).

2.2 Additional datasets

In addition to the mooring data that are the focus of this study, we use ancillary data from a number of different sources, which will be briefly described in this section

CTD data used in Fig. 2 were collected during four mooring deployment cruises in 2014 (cruise Thalassa MSM40), 2016 (MSM54), 2018 (MSM74), and 2020 (MSM94). The station spacing near the moorings was typically around 10 km, and the stations occupied during the 2016 cruise are indicated in the figure. The sections shown in the figure were produced by interpolating all measurements from the four cruises onto a regular grid with 2.5 km horizontal and 1 dbar vertical spacing using a Gaussian interpolation.

We also use data from the Argo program (Roemmich et al.2009), which are freely available online, including both individual float data and data products. The data from individual floats were downloaded from, using all profiles flagged as “good” on the data selection service that were found in the area of 50–65 N , 65–40 W between 1 January 2000 and 8 September 2020 (Argo2020,, resulting in a total of 41 165 temperature and salinity profiles from 568 floats. These data were used to produce the mean salinity map in Fig. 1 and are used to track the formation and boundary current input of LSW as described in Sect. 2.3. Data from the Argo mixed layer depth database by Holte et al. (2017) are used in Fig. 7 in order to highlight the mean convection area.

Bottom topography shown in Figs. 1, 3a, and 7 comes from the SRTM15+ product (Tozer et al.2019), which is the most recent product based on Smith and Sandwell (1997). These data are also used to differentiate between the interior and boundary current regions for the analysis of Argo data. Current vectors used in Fig. 1 are taken from the dataset described in Fischer et al. (2018).

2.3 Tracking Labrador Sea Water formation and input into the boundary current

To evaluate the input of LSW into the outflowing boundary current, we use data from Argo floats to investigate convection regions and LSW input in Sect. 3.1 and 3.2. The profile data are used to calculate mixed layer depths based on a density difference of Δσθ=0.01kg m−3 from the shallowest measurement, which has been shown to be a suitable threshold in the subpolar North Atlantic (Piron et al.2016; Zunino et al.2020). After the last convective profile for each season, defined here as a mixed layer depth exceeding 600 m, we track the trajectory until the float enters the boundary current. We define LSW input as a float crossing the 3000 m isobath and subsequently staying in the boundary current for at least two out of the next three profiles thereafter, similar to definitions used in previous studies (Georgiou et al.2020). A stricter criterion was also tested, requiring floats to subsequently be exported south of 53 N within the boundary current, reducing the number of floats used for the analysis by about half. However, results are not sensitive to this choice, with the curve shown in Fig. 9b being almost identical if this export criterion is used instead.

During a typical Argo cycle, floats park at a depth of either 1000 m or 1500 m for 10 d before descending to 2000 m, measuring a profile to the surface and staying there to transmit the data, which typically last about 1 h but can be as much as 12 h for older float models (Lebedev et al.2007). As a result, the floats are not strictly following a water parcel, as vertically sheared currents can carry them away from the water mass they had been advected with at depth. To ensure that floats used for the calculation are following LSW, we use layer thickness measurements to track them. Layer thickness is inversely proportional to the vertical density gradient, and newly formed LSW is readily detected as a maximum (Yashayaev and Loder2016). An example of layer thickness measurements from a float measuring convection in the interior of the basin and subsequently entering the boundary current on the southwestern side is shown in Fig. 3. Initially, the time series show mixed layer depth and maximum layer thickness varying in concert, as the thickest layer is the actively convecting mixed layer (Fig. 3b). This is also evident in the last profile measuring convection, as the maximum layer thickness is found at the lowest density measured, suggesting that no lighter water mass is present above the maximum associated with LSW (Fig. 3c, dark blue line). Subsequently, the mixed layer depth decreases to almost zero (Fig. 3b), indicating surface restratification, but a maximum in layer thickness remains. In the profiles (Fig. 3c), this can be seen as lower-density water accumulating above the maximum around σθ= 27.74 kg m−3. The magnitude of the maximum decreases through mixing with surrounding water, but it is still detectable in the profile as the float enters the boundary current.

The evolution of layer thickness and mixed layer depth shown in Fig. 3 is representative of a typical float measuring LSW formation in the interior that later enters the boundary current. Generally, after the initial profile measuring a deep mixed layer, surface restratification and mixing lead to a decrease in maximum layer thickness, but the signature of LSW remains detectable. To ensure that only floats moving with newly formed LSW were used, we discarded data from floats with changes in maximum layer thickness of more than 50 % between consecutive profiles. As a result, data from about 12 % of all floats measuring convection were omitted from the analysis, but the findings discussed in Sect. 3.1 do not change qualitatively if all float data are used instead. The layer thickness criterion is used only for floats measuring convection in the interior. For convection within the boundary current, LSW input occurs at the time of the initial profile, and we only use the constraint that two of the next three profiles have to be in the boundary current.

In order to investigate the seasonal timing of LSW input into the boundary current, we define a LSW input rate which is simply given by the number of floats measuring LSW input at a given time divided by the total number of floats entering the boundary current:

(1) Input LSW ( t ) = n ( t ) N ,

where t is the time given as day of year in 5 d bins, n(t) is the number of floats inputted during each time step, and N is the total number of floats measuring LSW input. This is used in Sect. 3.2 to compare the LSW input to the seasonal cycle of oxygen at K9. One possible caveat of this method is that it does not take into account the volume of LSW associated with each instance of input, which would require some information about the horizontal extent of each patch of LSW measured by the floats. Furthermore, variations in the number of floats present in the region at different times could lead to biases in the resulting estimate of the LSW input rate.

3 Results

Time series of oxygen at 600 m nominal depth (Fig. 4) show an overall similar picture at all sites: Starting in February–March, high-frequency oxygen variability intensifies, and mean oxygen concentrations increase by around 10 µmol L−1, remaining elevated during spring. Subsequently, the short-term variability subsides, and beginning in July, oxygen concentrations decrease slowly throughout the year until the next increase the following winter. The shift in oxygen concentrations appears to occur in concert with changes in temperature and salinity (Fig. 5). TS changes occur largely along isopycnals, meaning that the effect of T and S variations on density compensate each other. Correlation between T and O2 is high for all of the sites, between 0.80 and 0.94. Although some of this can be explained by solubility changes, which are dependent on temperature (Garcia and Gordon1992), correlations between saturation percentage and temperature are still high, between 0.58 and 0.84. If changes in oxygen were simply due to temperature-driven solubility differences at a constant saturation percentage, the correlation of saturation and temperature would be zero. The fact that both saturation percentage and concentration are correlated with temperature suggest that the changes are associated with the presence of different water masses. In fall and early winter, measurements show warmer, saltier, and lower-oxygen water with a lower saturation percentage, while colder, fresher water with higher O2 concentration and saturation is predominantly observed in late winter and spring.

Figure 4Time series of oxygen concentration at 610 m for K7, K8, K9, and K10 (a–d). Gray lines show original 15 min data, and black lines show time series filtered with a 10 d running mean. Note that axis limits are 5 µmol L−1 lower in (a), but the range is the same for all panels.


The majority of the TS measurements lie roughly along a mixing line between typical endmembers for LSW and IW (see Fig. 5 for reference). The clustering of measurements between these two endmembers suggests that variability in properties at 53 N, including oxygen, is chiefly controlled by changes in the fraction of each source water mass found in the boundary current.

In the following sections, we will analyze the main features of the time series in more detail, discussing first the initial oxygen increase in February and March and associated higher-frequency variability using data from the K9 mooring (Sect. 3.1). This is followed by an analysis of a complete seasonal cycle of oxygen concentration at K9, covering convection and the following restratification period in 2016–2017 (Sect. 3.2), and a brief analysis of the differences between the moorings (Sect. 3.3).

3.1 High-frequency oxygen variability in February–April

As oxygen concentrations at the mooring sites begin to increase in late February, there is a concurrent increase in the spread of the measurements (Fig. 4). In 2017, during the months of February, March, and April, oxygen concentrations at K9 vary between values typical of the months prior, between 297 and 301 µmol L−1, and higher values ranging from 305 to 315 µmol L−1, which subsequently persist until late June. During these months, the spread between the minimum and maximum values measured in a single day can be as much as 15 µmol L−1, and the standard deviation increases from 1.1 µmol L−1 for the 3-month period from 1 November 2016 to 31 January 2017 to 4.2 µmol L−1 for 1 February to 30 April 2017.

This wide spread in the measured oxygen concentrations is also reflected in differing temperature and salinity properties for a 20 d period in February (Fig. 6a). As for the overall dataset, the highest oxygen values occur at the lowest temperatures and salinities. Starting from a cluster of points with higher temperatures and salinities for 7 to 12 February, water with properties of LSW begins to emerge in the following days, and by 27 February, some of the measured TS values are similar to those found in the center of the basin during convection (see also Fig. 5). During this time, much of the annual range of temperature, salinity, and oxygen properties found in the 2-year record is observed within just 20 d (see Figs. 5, 6a). High-frequency variability in properties at the beginning of the convection season was also observed in temperature measurements further upstream in the boundary current by Cuny et al. (2005). They concluded that LSW was being formed within the boundary current itself and attributed the high-frequency variability to spatially varying atmospheric forcing producing LSW with differing properties. However, the fact that there is also a large spread in oxygen in Fig. 6a suggests that some of the water observed during this time has not yet been ventilated. Instead, the variability could be a sign that the earliest convective activity is episodic and patchy, resulting in newly ventilated LSW being transported with the boundary current alongside the lower-oxygen IW found prior to the onset of convection.

Figure 5Temperature and salinity diagrams at about 600 m nominal depth from all four boundary current moorings, showing daily averages. Colors correspond to oxygen anomalies relative to the mean at each instrument, ranging from 10 µmol L−1 (blue) to +10 µmol L−1 (red). Purple stars and shading show values from the wintertime mixed layer at the SeaCycler mooring in the interior Labrador Sea (Atamanchuk et al.2020), and green symbols and shading show typical wintertime conditions of IW near Greenland, taken from Pacini et al. (2020). Potential density (σθ) contours are drawn with a spacing of 0.025 kg m−3, and solid black contours show densities of σθ=27.68 and σθ=27.8kg m−3.


Figure 6(a) TS properties at K9 at 15 min resolution for four consecutive 5 d periods during the initial oxygen increase in February. Colors correspond to oxygen concentration. Potential density (σθ) contours are drawn with a spacing of 0.025 kg m−3, and solid black contours show densities of σθ=27.68 and σθ=27.8kg m−3. Note changed axis ranges compared to Fig. 5. (b) Density measurements at K8 at 50 and 608 m nominal depth. Shading for 608 m values shows density ranges of ± 0.01 and ± 0.03 kg m−3 from the measured values.


Although Cuny et al. (2005) showed that LSW formation can occur in the boundary current, there is no evidence of local convection near 53 N in our data. While no near-surface instrument was deployed on the K9 mooring, data from about 50 m depth at K8 suggest that the boundary current at 53 N remains stably stratified (Fig. 6b). Vertical density gradients generally exceed the Δσθ=0.01 kg m−3 threshold often used to determine mixed layer depth in the region (Piron et al.2016), as well as the less stringent threshold of Δσθ=0.03 kg m−3 used by de Boyer Montégut et al. (2004). This suggests that there is no active convection near the mooring sites during the study period. Instead, we conclude that the water must be advected to the 53 N array from upstream. The observations are consistent with LSW being formed in the interior of the basin and entering the boundary current along isopycnals (Georgiou et al.2020). Alternatively, the observed changes could also be explained by convective activity taking place within the boundary current further upstream.

In order to identify the source of the LSW arriving at 53 N in February, we use data from Argo floats to investigate the timing of convection and input into the boundary current. Figure 7 shows the fraction of profiles within 50 km of each 0.5×0.5 bin that measured mixed layers deeper than 600 m based on the Holte et al. (2017) Argo mixed layer depth database. The most salient feature is a large area of deep mixed layers in the interior of the basin, near the mixed patch often defined as the LSW formation region (Lazier1973; Yashayaev2007; Atamanchuk et al.2020). The area extends eastward and indicates a major interior ocean spreading pathway that connects the Labrador Sea and Irminger Sea (Talley and McCartney1982; Sy et al.1997; Fischer et al.2018; Zunino et al.2020). Another region of increased occurrence of deep mixed layers is found inshore of the 3000 m isobath near 55 W, close to where Pickart et al. (2002) first reported evidence of convection in the boundary current. No mixed layers deeper than 600 m are found in the boundary current south of 55 N, consistent with our findings from the mooring data that convection did not occur locally at 53 N.

Also shown in the figure are convection locations for all floats from the 2000–2020 period from the Argo data described in Sect. 2.2, using the last profile measuring a mixed layer deeper than 600 m. The points generally line up with the mean picture from the Holte et al. (2017) data, with differences between the two likely occurring due to continued modification during convection. Overall, the majority of the profiles measuring convection in both datasets are in the interior. Most of those floats that measure deep mixed layers inshore of the 3000 m isobath are located immediately adjacent to the interior patch.

Figure 7Fraction of total profiles from Holte et al. (2017) database with mixed layer depth (MLD) over 600 m. Grey contours show the 2000 and 3000 m isobaths, and the magenta line shows the 53 N array. Overlaid are the convection locations for each float measuring convection as defined in Sect. 2.3, differentiating between floats that are not found in the boundary current within the same year (red dots), those that do enter the boundary current (yellow squares), and those that are found in the boundary current in either January or February (green diamonds).

A quantitative analysis of the relative importance of the different formation regions is beyond the scope of this study as the Lagrangian nature of Argo floats results in undersampling of fast-flowing boundary currents relative to regions with weak currents like the interior Labrador Sea. Nonetheless, our dataset reveals some features about the differences between convection in the interior and within the boundary current. The profile locations in Fig. 7 are color coded by whether or not the float enters the boundary current during the same year using the definition of LSW input described in Sect. 2.3. About half of the floats measuring convection in the interior remain there until the following convective season, consistent with model studies showing that newly formed LSW can remain in the interior for several seasons (Straneo et al.2003; Georgiou et al.2020). Conversely, most of the floats measuring deep mixed layers inshore of the 3000 m isobath stay within the boundary current afterwards, suggesting that a larger fraction of LSW formed in the boundary current is exported out of the basin compared to LSW formed in the interior. Moreover, out of the more than 100 floats measuring convection to depths exceeding 600 m, only six are subsequently found inshore of the 3000 m isobath in the months of January and February. The convection location of these floats is highlighted in Fig. 7 by green diamonds. The floats found in the boundary current during this time measured either convection very close to the 3000 m isobath, which allows them to enter the boundary current quickly, or convection within the boundary current region itself. This suggests that the initial arrival of the elevated oxygen signal at the moorings in February is likely to be a result of convection within, or close to, the boundary current.

3.2 Seasonal cycle

To investigate the broader-scale variability, we look at the relation between oxygen and temperature and salinity for the period of December 2016 to December 2017, covering a full annual cycle including convection and restratification. Figure 8 shows the evolution of the oxygen concentration throughout the year in the form of monthly histograms with a bin size of 1 µmol L−1. The color for each bin shows the mean spiciness (π), which is a tracer of T and S variations that are orthogonal to density, making it useful for differentiating between water masses with different TS properties but similar density (Flament2002). Spiciness is used here as a water mass index to quantify the relative contribution of the two endmember water masses, LSW and IW. Figure 8 uses data from the K9 mooring, and results from all other moorings are discussed in the following section.

Figure 8Monthly histogram of oxygen from December 2016 to December 2017 at about 600 m depth at K9. Colors show the mean spiciness of each O2 bin, with lower values corresponding to LSW and higher values to IW (see text for definition).


Overall, the seasonal cycle is marked by a shift from water with lower oxygen and TS properties closer to those of IW towards high-oxygen LSW, which is most prevalent in the months of March–July. In January, the bulk of the measurements at K9 show similar properties. The most abundant class of oxygen values is at 298 µmol L−1, and the water mass index shows higher temperatures and salinities, indicative of IW most likely advected from upstream within the boundary current. There is little spread around this mean value, with almost all measurements falling in the range of 297–301 µmol L−1. In the month of February, measurements showing water with oxygen concentrations around 297 µmol L−1 still make up the majority of values seen in the record, but the number of observations in the central bin decreases. Moreover, a smaller secondary peak emerges at a concentration of 306 µmol L−1, with lower mean temperature and salinity values, and there is a wide spread of oxygen concentrations measured, ranging from 297 to 315 µmol L−1. This is consistent with the view that, at this time, the boundary current waters comprise non-ventilated water, as well as patches of recently convected and ventilated LSW (see Sect. 3.1).

In the following months, the LSW peak becomes more pronounced and replaces the class of higher temperature and salinity water with low O2 as the most abundant water mass. In April, the most commonly measured O2 value is at its annual maximum of 306 µmol L−1. By May, the boundary current has become more homogenized, as evidenced by the higher number of observations in the central O2 bin. With the increased input of newly formed LSW, the TS properties are at the opposite end of the water mass index (π) compared to December–February. After June, oxygen concentrations are beginning to decrease again, with a concurrent increase in T and S, indicative of a larger fraction of IW. Only isolated instances of elevated oxygen concentrations are observed at K9 by August and September as the water at 600 m depth becomes warmer, saltier, and less oxygenated. In the autumn months, the signature of LSW has all but disappeared, with water found at the mooring once more resembling the properties of IW, before the cycle begins anew the following winter.

Figure 9(a) Monthly histogram of oxygen from December 2016 to December 2017 at 610 m depth at K9, as in Fig. 8. The size of each circle corresponds to the number of observations, the black line shows the bin with the highest number of observations for each month, and colors show the mean spiciness of each O2 bin. (b) Climatological seasonal cycle of LSW input into the boundary current from float data in 5 d bins as a fraction of the total input of LSW over the year (Eq. 1). Gray bars show separately the input from boundary and interior convection in each 5 d bin, and the black line shows the combined input smoothed with a five-point running mean. (c) Monthly histogram of oxygen from December 2016 to July 2017 at 500 m depth for the SeaCycler mooring in the deep convection area in the interior of the basin. The size of each circle corresponds to the number of observations, the black line shows the bin with the highest number of observations for each month, and colors show the mean spiciness of each O2 bin.


Figure 9 compares the seasonal oxygen cycle at K9 from Fig. 8 to the seasonal cycle of LSW input into the boundary current, as well as a partial seasonal cycle at 500 m depth from the SeaCycler mooring in the interior of the basin (see Fig. 1 for location). The input of LSW into the boundary current is an average over all years from the float data used in Sect. 3.1 and is calculated using Eq. (1). The gray bars show the LSW input during each 5 d period resulting from convection within the boundary current and LSW entering the boundary current from the interior, and the black line shows the combined input. Overall LSW input begins to increase in late January, peaks in April, and vanishes in July except for an isolated instance later in the year. Only 47 floats from the dataset measure LSW input into the boundary current, which could result in some bias of our estimate relative to the true underlying variability. However, the timing is consistent with related measures such as the amount of LSW leaving the interior derived from measurements of layer thickness (Yashayaev and Loder2016) and rates of boundary–interior exchange inferred from heat content changes (Straneo2006), suggesting that the overall variability is captured despite the very limited number of data points. The curve has a similar shape to the seasonal cycle of the most commonly measured O2 concentration at K9 for each month, shown as a black line in the top panel, but the two are shifted relative to one another. The O2 concentration maximum increases in March, reaches its peak in April, and then stays at a similar level until June, before starting to decrease in July and August, lagging the increase in LSW input by about 1–2 months. Typical current speeds in the boundary current at the depth of the sensors used here are on the order of 15 cm s−1 (Zantopp et al.2017), and the distance to the region where the interior convection area is closest to the 3000 m isobath and where boundary current convection is most frequent is about 450 km (see Fig. 7). The timescale for LSW to arrive at the 53 N moorings after entering the boundary current would therefore be about 35 d. This is consistent with the time lag of approximately 1 month between the seasonal cycle of LSW input and the O2 maximum time series at the K9 mooring. The data therefore support the hypothesis that the increase in oxygen concentration at the exit of the basin is largely controlled by the rate of formation and subsequent export of LSW.

LSW input from convection that occurs directly within the boundary current precedes the outflow of LSW from the interior by about a month, supporting the qualitative assessment made in Sect. 3.1 that the initial arrival of newly ventilated LSW at 53 N is likely a result of boundary current convection. Overall, 45 % of the input (21 of 47 floats) is associated with boundary current convection and 55 % (26 of 47 floats) with interior convection, with both contributing about equally in March and April when overall LSW input is highest. However, due to the caveats of our method discussed in Sect. 2.3, these numbers should not be interpreted as a measure of the LSW volume flux associated with the two sources.

At the SeaCycler mooring, there are two distinct peaks in the oxygen histogram in January as the mixed layer reaches the sensor depth and ventilates water that had not been in contact with the atmosphere since at least the previous winter. In February, newly ventilated, high-oxygen LSW becomes the most prominent water mass, and O2 remains elevated until the end of the available record in June 2017. The maximum oxygen concentration found in the interior before January is 304 µmol L−1, much lower than the maximum values at 53 N during March–July, confirming that the increased oxygen levels of the outflowing boundary current must be associated with recently ventilated water rather than LSW formed in previous winters.

If the bulk of the LSW entering the boundary current from the interior originates in the center of the basin near the SeaCycler mooring, the shortest distance to the boundary current would be about 230 km. With a time difference of about a month between the oxygen increase at SeaCycler and the increased input of LSW from the interior into the boundary current, this implies a mean advection speed of 9 cm s−1. This is much larger than the climatological advection speed in the interior (see Fig. 1), suggesting that short-lived features such as eddies and current meanders may play a crucial role in the export of newly formed LSW.

Figure 10Monthly oxygen histograms, as in Fig. 8a, for the K7, K8, and K10 moorings. Circle size corresponds to number of observations, the black line shows the bin with the highest number of observations for each month, and colors show the mean spiciness (π) in each bin, highlighting the relative influence of IW and LSW.


3.3 Differences at the 53 N array

Although the seasonal oxygen cycle at K9 discussed in the previous section is representative of the general picture along the 53 N array, there are also some horizontal differences. The four moorings that carry oxygen sensors are deployed in different parts of the boundary current system at the exit of the Labrador Sea (Zantopp et al.2017, their Fig. 1). The K9 mooring, discussed in the previous section, is situated within the core of the Deep Western Boundary Current (DWBC), which has only a weak velocity shear between 400 and 2000 m. K7 sits over the continental shelf beneath the core of the Labrador Current, which is surface intensified and carries colder, fresher water of Arctic origin. The K8 sensors at 600 m depth are within the DWBC but closer to the boundary with the Labrador Current and can be within it when the currents are meandering. At all three of these moorings, the flow is southward, with mean velocities over the deployment period between 13 and 15 cm s−1. In contrast, the K10 mooring is at the boundary between the offshore edge of the DWBC and its northward recirculation. Over the deployment period, there is a weak mean northward flow at about 600 m, with a speed of 3.5 cm s−1, but daily mean values are distributed evenly around this mean, ranging from 20 to +20 cm s−1. The seasonal cycle of oxygen concentrations at all four moorings (Figs. 9a and 10) shows elevated values and water with properties closer to LSW during spring and summer, as well as water with lower oxygen values resembling IW in autumn and winter.

At the border of the Labrador Current and DWBC, the oxygen concentrations measured at K8 show a similar seasonal cycle as those at K9. Properties measured by the water mass index (π) are also comparable, although the spread in T and S is higher (Fig. 5), suggesting that the DWBC here is less homogenized than at K9.

Over the continental slope at K7, the oxygen concentration is lower by about 5 µmol L−1 throughout the year, and the earliest arrival of elevated oxygen occurs in January compared to February at K9. The water is generally warmer and fresher than at the other moorings, which may be the reason for the lower oxygen values and also results in a slightly lower mean density. There also seems to be a stronger seasonal variation in density compared to the other moorings. Between the warmer, saltier, oxygen-poor water in winter and the high-oxygen, low-TS water in spring, density increases by more than 0.05 kg m−3 (Fig. 5). This is similar to variability observed upstream in the boundary current near 55 N by Cuny et al. (2005), who showed that this is a result of isopycnals steepening in winter over the continental slope. This weakened stratification may allow for the formation of uLSW (Pickart et al.1997). At K7, the highest oxygen concentrations occur at lower densities and salinities, seen as a “tail” in the bottom left corner of the TS plot. Moreover, the first high-oxygen measurements occur almost 3 weeks earlier than at K9, which may suggest that these are a product of the formation of uLSW within the Labrador Current and that this formation takes place closer to the mooring section than the production of classical LSW found at the other moorings. Since Argo floats generally stay offshore of the 2000 m isobath, formation of uLSW within the Labrador Current would be difficult to detect with the methods used in Sect. 3.1.

At the offshore edge of the southward boundary current, at mooring K10, the oxygen concentrations are more spread out, and elevated values are not found until the beginning of March. The spread of oxygen concentrations measured during any given month is much larger than for the other moorings, and the TS signature of the two different water masses is less pronounced. This is likely a result of current meanders, with the mooring measuring the DWBC some of the time and its recirculation at other times, resulting in a wider range of ventilation histories. High oxygen concentrations are found during times of both southward and northward flow, suggesting that LSW is carried southward with the offshore edge of the DWBC at K10, but some of it may also enter the northward recirculation back towards the Labrador Sea.

4 Discussion

4.1 LSW formation and export pathways

The analysis in Sect. 3.1 shows that deep mixed layers associated with convection are more commonly found in the interior. In Fig. 7, the area where more than 5 % of profiles measure mixed layers deeper than 600 m is about 15 times larger in the interior than inshore of the 3000 m isobath. However, this does not necessarily translate to 15 times larger oxygen supply as water parcels have a longer residence time in the interior, whereas newly formed LSW in the boundary current is rapidly exported (MacGilchrist et al.2020, 2021). This is also evident in our data analysis, in which a significant number of Argo floats measuring deep mixed layers in the interior do not enter the boundary current during the same season and therefore do not directly contribute to the oxygen export by the outflowing boundary current. Moreover, LSW formation within the boundary current can also occur as slantwise convection (Cuny et al.2005) or in parts of the boundary current that are shallower than 2000 m, neither of which would be detectable with Argo floats. As discussed in Sect. 3.3, some evidence of convection over the continental slope upstream can be seen in data from the K7 mooring at a water depth of 1400 m. However, neither the mooring data, which do not measure local convection, nor the Argo data, biased towards sampling regions with mild currents at depth, provide sufficient information to quantitatively determine the portion of the LSW exported at 53 N that is formed in the boundary current. The importance of boundary versus interior convection for the outflowing LSW remains an open question, with one previous estimate from float data concluding that both could be important for the properties of the outflowing boundary current (Palter et al.2008).

Another crucial aspect concerning the export of LSW from the formation region is how the water formed in the interior enters the boundary current. Palter et al. (2008) showed that both eddies and a mean cross-isobath flow can contribute to LSW export, with floats entering the boundary current on all sides of the basin, which is also the case for the float data used here (not shown). The timescale for a water parcel to round the basin from Cape Farewell to the Labrador side is about 147 d (Cuny et al.2002). This implies that LSW entering the boundary current, or formed within it, on the northern side of the basin during the height of convection in March would not reach 53 N until August, when the influence of LSW at K9 is already decreasing. Thus, the outflow of LSW during summertime observed at 53 N is likely fed mostly by LSW entering the boundary current on the western or southwestern side, allowing for faster export out of the basin. The timing of the seasonal cycles of oxygen concentration in the interior, LSW input from the interior into the boundary current, and oxygen concentration at K9 (Fig. 9) implies that eddies may play an important role in exchanging water between the boundary current and interior, consistent with previous studies (Eden and Böning2002; Straneo2006; de Jong et al.2014).

In a recent study, Georgiou et al. (2020) used model data to analyze the source of water parcels reaching the 53 N section, and they found that 41 % did not leave the boundary current after passing Cape Farewell at the southern tip of Greenland, while a further 44 % are modified in the interior of the basin before re-entering the boundary current, and the remainder follow several less common pathways. These numbers are broadly consistent with the results from Fig. 8, with about equal amounts of LSW and IW found at the moorings over an annual cycle. Our study suggests that there may be a seasonal dependence to the split observed by Georgiou et al. (2020), with the LSW export (interior route) occurring chiefly in the months after convection and IW (water from Cape Farewell) being more ubiquitous during autumn and winter. The input of newly ventilated LSW into the boundary current in the months following convection observed here is consistent with the more direct southward export of lighter LSW found in their study, as well as results from an earlier modeling study (Brandt et al.2007). Georgiou et al. (2020) suggested that pathways are different for the densest LSW formed in the interior, with water entering the boundary current on the northern side of the basin and thus taking longer to reach 53 N. This density dependence may explain why most of the LSW found in the mooring record seems to follow the faster southward export route as the lower part of the LSW layer is not sampled in the current sensor setup. Another model study by Handmann (2019) highlighted different routes that exist for LSW export, including several pathways for LSW entering the Irminger Sea. Their study also proposes that some of the water exported through the 53 N section subsequently leaves the DWBC at Orphan Knoll near 50 N and recirculates back towards the north, reaching both the Labrador and Irminger seas. In our data from the mooring furthest offshore, K10, there is evidence of high-oxygen water with northward velocities (Sect. 3.3), suggesting that newly formed LSW may also be entrained into this northward recirculation.

These studies show that many more LSW export pathways exist in addition to the direct southward route that seems to be responsible for the pronounced seasonal signal at the 53 N moorings. Although the large summertime increase in oxygen (Fig. 8) is associated with LSW formed in the same year, the oxygen concentrations can also be impacted by export of LSW along these alternative routes, as well as further water mass transformation along the pathway. For example, some of the LSW formed in the interior one year may be advected into the Irminger Sea, enter the boundary current on the eastern side of Greenland, and stay within it until it reaches 53 N several years later. Deep convection also takes place in the Irminger Sea (de Jong et al.2018), leading to a somewhat weaker but still significant uptake of oxygen (Maze et al.2012; Palevsky and Nicholson2018). The water mass formed by convection in the Irminger Sea, known there as Irminger Sea Intermediate Water, or ISIW (Le Bras et al.2020), is mixed into the boundary current in the Irminger basin before it enters the Labrador Sea, affecting the properties of what we have referred to here as Irminger Water (IW). Therefore, both the export of LSW into the Irminger Sea and convection occurring there will have affected the properties of the ”low-oxygen” IW seen in Figs. 810, and the annual minimum concentration of 298 µmol L−1 at K9 is likely already elevated compared to what it would have been without ventilation in the Labrador and Irminger seas.

4.2 Global impact of Labrador Sea ventilation

Due to the global importance of LSW, it has long been assumed that changes in its formation would play a crucial role in setting the variability in the Atlantic Meridional Overturning Circulation (Koltermann et al.1999; Yashayaev2007). However, results from the Overturning in the Subpolar North Atlantic Program (OSNAP) have called into question this long-held belief, showing that the bulk of the water mass transformation across isopycnals occurs east of Greenland in the Irminger Sea and Iceland Basin. This water mass transformation in the eastern section appears to be the main source of overturning variability (Lozier et al.2019), although recent modeling results suggest that LSW formation may still be important for variability on multidecadal timescales that are not yet resolved by the OSNAP time series (Yeager et al.2021).

In a recent follow-up study to the OSNAP results, Zou et al. (2020) showed that the disconnect between overturning and LSW formation is in large part due to density-compensated changes in temperature and salinity. Although over 12 Sv of water mass transformation occurs in the Labrador Sea from warm, salty waters to colder and fresher ones, the net transformation into the LSW density range is only about 4 Sv. Our results are consistent with these findings, showing that while the water mass properties of the outflow are altered substantially by the input of newly formed LSW, its density remains largely unchanged. This suggests that LSW from the interior enters the boundary current primarily along isopycnals, consistent with recent results from models (Brüggemann and Katsman2019; Georgiou et al.2020), with some along-isopycnal mixing occurring before it reaches 53 N. Our time series also show for the first time direct evidence that significant changes occur in the oxygen content of the outflow, driven by the rapid export of newly ventilated LSW. While we only use data from a fixed depth, oxygen concentrations in the summer are elevated in much of the LSW layer above 1200 m (Fig. 2). Assuming that the oxygen input is comparable throughout this layer, we can calculate a first estimate of the oxygen transported out of the basin due to the input of newly formed LSW: using a mean transport for the LSW layer at 53 N of 14.5 ± 3.8 Sv (Zantopp et al.2017), the oxygen increase of 6 µmol L−1 at K9 during March–August 2017 relative to the January baseline of 298 µmol L−1 (Fig. 8) corresponds to an oxygen export of 1.37 ± 0.37 Tmol (1 Tmol=1012mol) for these 6 months. Over the whole year, the mean oxygen increase is 3.5 µmol L−1, corresponding to an export of 1.60 ± 0.42 Tmol relative to the baseline value. With an integrated wintertime uptake across the air–sea interface due to gas exchange of 21.2 mol m−2 (Atamanchuk et al.2020) and a convection area of approximately 150 000 km2, the total uptake in the interior is 3.18 Tmol during the same year. Thus, about 43 % of the oxygen taken up during convection is exported out of the basin via this fast export route of LSW along the DWBC during summer and another 7 % through slower export routes during the rest of the year. The good agreement with the air–sea flux estimates from Atamanchuk et al. (2020) provides independent support for their notion that bubble injection plays an important role in oxygen uptake during deep convection. As discussed in Sect. 4.1, LSW is also exported along different pathways than we considered here, and some remains in the basin until the next year, which may account for the remaining 50 % of the oxygen uptake. Conversely, the recirculation of the DWBC towards the formation region after it reaches Orphan Knoll (Fischer and Schott2002) may diminish the net oxygen export, although the recirculation strength is typically only about 10 % of the DWBC transport at 53 N (Zantopp et al.2017).

The deep North Atlantic is one of the few ocean regions that has not experienced significant deoxygenation in recent decades (Schmidtko et al.2017), with oxygen increasing slightly by about 1.7 Tmol yr−1 over their multi-decadal study period compared to a globally integrated oxygen loss of 70 Tmol yr−1 below 1200 m depth. The continued formation and export of LSW and NADW may have contributed to staving off deoxygenation in this basin, and the 1.60 Tmol yr−1 of export calculated for 53 N could be sufficient to supply much of the oxygen consumed in the upper NADW (uNADW) layer in the North Atlantic. Using a global mean apparent oxygen utilization rate of 0.1 µmolkg-1yr-1 at 1500 m depth (Karstensen et al.2008), an area of 26 159 000 km2 for the deep Atlantic Ocean between the Equator and 50 N, and a mean layer thickness of 800 m for the isopycnal range 27.68 kg m−3σθ 27.8 kg m−3, we estimate the annual oxygen consumption for this volume of water to be 2.2 Tmol yr−1. If the depth-varying profile from Karstensen et al. (2008) is used instead, with the depth range occupied by uNADW at each point estimated from density data, the estimate increases to 3.8 Tmol yr−1. Using these two estimates as a range of possible values for the annual oxygen demand in this volume, the rapid southward export of LSW discussed here supplies 42 %–73 % of it, with the rest likely provided by slower export of LSW through other routes, as well as uptake in other deep water formation regions such as the Irminger Sea (Palevsky and Nicholson2018), Iceland Basin (Maze et al.2012), and the Gulf of Lion in the Mediterranean Sea (Ulses et al.2021). Thus, despite playing a less significant role than previously thought in setting the strength of the overturning (Lozier et al.2019), convection in the Labrador Sea does appear to crucially contribute to the supply of oxygen from the subpolar North Atlantic to the deep subtropical and tropical Atlantic Ocean. This highlights the important distinction between the overturning circulation and the ventilation of deep waters, as discussed by Naveira Garabato et al. (2017) and MacGilchrist et al. (2020, 2021). Although the density transformation is weaker in the Labrador Sea, there is a more direct connection between the atmosphere and deep ocean during the formation of LSW compared to the water mass formation and transformation processes in the eastern subpolar gyre, allowing for more intense oxygen uptake.

4.3 Outlook

The focus of the present study was to analyze the seasonal cycle in 2017, the one full calendar year that was resolved by our 2-year dataset. Although only data for the first half of the next year are shown in Fig. 4, some interannual changes in the arrival of LSW in 2018 are already apparent, such as a larger oxygen increase at K7 and a large increase followed by stronger short-term variability at K10. The formation and properties of LSW are well established to be highly variable in time (Lazier1973; Yashayaev2007), and our data were collected during a period which featured some of the deepest-reaching convection since the 1990s (Yashayaev and Loder2017). Although current meter data from the 53 N array show no clear relation between volume transport and the strength of deep convection (Dengler et al.2006; Fischer et al.2010), it might intuitively be expected that weaker LSW formation would lead to less LSW input into the boundary current and consequently lower oxygen concentrations. Accordingly, the oxygen export during a year like 2010 with some of the shallowest convection in the interior (Yashayaev and Loder2016) may be weaker than the values reported here for 2017. The extent to which the oxygen export estimates presented here are subject to interannual variability will be addressed in a future study using data from continued deployments of the oxygen sensors at 53 N that will extend the time series until 2022.

A key topic of research regarding LSW formation is its susceptibility to future climate change. Model studies have suggested that increased freshwater input resulting from the melting of Greenland's ice sheets may weaken or even shut down convection in the Labrador Sea (Manabe and Stouffer1997; Rahmstorf1999; Clark et al.2002). While the re-emergence of LSW formation to near record depths in 2014–2018 shows that the effect of ice melt has not yet critically affected convection (Yashayaev and Loder2016, 2017), recent climate projections and observations continue to point towards a possible shutdown in the future (Böning et al.2016; Oltmanns et al.2018). If deep convection in the Labrador Sea were to cease, the “trap door” of oxygen uptake to the deep ocean would be shut (Atamanchuk et al.2020), cutting off a large part of the oxygen supply discussed above. Moreover, since ventilation of LSW can occur without significant density-space overturning, large changes in oxygen supply are possible even if the overturning circulation remains stable.

In a recent model study on future deep ocean deoxygenation, Oschlies (2021) posited that their work “calls for more research efforts to explore the baseline of these systems before the unavoidable change will hit”. The present study represents a first step towards quantifying the impact of LSW formation and export on the supply of oxygen to the open Atlantic Ocean. However, more measurements will be needed to better understand the complex interplay between different deep water formation regions, as well as the resulting supply of oxygen associated with ventilation in the subpolar North Atlantic as a whole. To that end, oxygen sensors were added to several other moorings in the OSNAP array in 2020, and a total of 70 O2 sensors are now deployed along the path of the boundary current in the Irminger Sea, the Labrador Sea near the southern tip of Greenland, and on the 53 N array. With the resulting data, it will become possible to gain a much more complete understanding of the processes controlling oxygen supply from the subpolar North Atlantic than has been presented here.

5 Conclusions

The results put forth in this study show the effect of Labrador Sea Water formation and export on the oxygen concentration of the outflowing deep western boundary current. During autumn and winter, properties measured at the 53 N array at about 600 m nominal depth are similar to those of Irminger Water, which enters the basin from the Irminger Sea. Following the onset of deep convection, LSW first arrives at the offshore moorings in the second half of February, with a wide range of temperature, salinity, and oxygen properties reflecting a sporadic input of heterogeneous LSW formed in or near the boundary current. As convection becomes more vigorous and sustained, and more newly formed LSW enters the boundary current from the interior, it becomes the dominant water mass found at 53 N in April–May. As a result, the most commonly measured oxygen concentration for these months is 7–8 µmol L−1 higher than in December and January. By September, the signature of LSW has largely disappeared, and properties once again resemble those of IW, with accordingly lower oxygen concentrations. This variability in properties seems to be controlled mostly by a direct southward export of newly ventilated LSW with the outflowing boundary current, which occurs chiefly in the 6 months following convection. However, additional measurements throughout the LSW density range would be needed to evaluate whether the seasonal cycle is the same for the entire layer, or just for the part of it measured with the current sensor setup.

The increase in oxygen content due to the input of LSW into the boundary current was estimated to lead to an added oxygen export out of the basin of 1.60 ± 0.42Tmol yr−1, which is about half of the oxygen taken up during deep convection in the interior. The annual oxygen supply from this fast, direct export of LSW to the Atlantic Ocean is estimated to be equivalent to 42 %–73 % of the oxygen demand by organisms in the upper NADW layer of the North Atlantic Ocean, highlighting the crucial role that the Labrador Sea plays in the supply of oxygen to the deep ocean. When taking into account the complex recirculation and export pathways of LSW through the subpolar gyre, as well as uptake in other deep water formation regions, the total added oxygen export associated with North Atlantic Deep Water formation is likely even higher.

Although the important role of LSW formation for oxygen uptake has long been established, this is the first study that shows a continuously measured time series of the oxygen content of the boundary current over a full year. Likewise, while variations in LSW input into the boundary current over an annual cycle had been reported before, we present the first evidence that these seasonal variations control the properties of the boundary current at the exit of the basin at 600 m depth. Together with the estimates of oxygen export presented above, these results address the question of how quickly and efficiently newly ventilated LSW is exported, which is particularly crucial in the context of estimating storage rates of gases such as oxygen and anthropogenic carbon. However, our study is only a first step towards answering these questions, and more research will be needed to fully quantify the importance of the subpolar North Atlantic for supplying oxygen to the global oceans via the export of NADW and to determine how susceptible this process could be to climate change.

Data availability

The mooring and hydrographic data used in this study will be made freely available on A snapshot of Argo profile data used in this study can be found at (Argo2020).

Author contributions

The manuscript was prepared by JaK, with contributions from DA, JoK, PH, and DWRW. DA was responsible for designing the experiment and calibrating the mooring data, and data analysis was carried out by JaK. DWRW and JoK acquired funding for the project.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


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


The 53 N array is part of the OceanSITES global ocean time series network. Argo float data were collected and made freely available by the International Argo program and the national programs that contribute to it (,, last access: 1 October 2020). The Argo program is part of the Global Ocean Observing System. The authors would like to thank Marilena Oltmanns for help with processing the mooring data.

Financial support

This research has been supported by the Canada Excellence Research Chairs, Government of Canada (Ocean Science and Technology grant), the Canada First Research Excellence Fund (Ocean Frontier Institute grant), the Horizon 2020 (Blue-Action (grant no. 727852) and EuroSea (grant no. 862626)), and the Bundesministerium für Forschung und Technologie (RACE grant). Jannes Koelling was funded by the Ocean Frontier Institute through the International Postdoctoral Fellowship.

Review statement

This paper was edited by Tina Treude and reviewed by Ilaria Stendardo, Graeme MacGilchrist, and three anonymous referees.


Argo: Argo float data and metadata from Global Data Assembly Centre (Argo GDAC) – Snapshot of Argo GDAC of 1 October 2020, SEANOE [data set],, 2020. a, b

Atamanchuk, D., Koelling, J., Send, U., and Wallace, D.: Rapid transfer of oxygen to the deep ocean mediated by bubbles, Nat. Geosci., 13, 232–237,, 2020. a, b, c, d, e, f

Atkinson, C. P., Bryden, H. L., Cunningham, S. A., and King, B. A.: Atlantic transport variability at 25 N in six hydrographic sections, Ocean Sci., 8, 497–523,, 2012. a

Böning, C. W., Behrens, E., Biastoch, A., Getzlaff, K., and Bamber, J. L.: Emerging impact of Greenland meltwater on deepwater formation in the North Atlantic Ocean, Nat. Geosci., 9, 523–527,, 2016. a

Bower, A. S., Lozier, M. S., Gary, S. F., and Böning, C. W.: Interior pathways of the North Atlantic meridional overturning circulation, Nature, 459, 243–247,, 2009. a

Brandt, P., Funk, A., Czeschel, L., Eden, C., and Böning, C. W.: Ventilation and Transformation of Labrador Sea Water and Its Rapid Export in the Deep Labrador Current, J. Phys. Oceanogr., 37, 946–961,, 2007. a, b, c

Brüggemann, N. and Katsman, C. A.: Dynamics of Downwelling in an Eddying Marginal Sea: Contrasting the Eulerian and the Isopycnal Perspective, J. Phys. Oceanogr., 49, 3017–3035,, 2019. a

Clark, P. U., Pisias, N. G., Stocker, T. F., and Weaver, A. J.: The role of the thermohaline circulation in abrupt climate change, Nature, 415, 863–869,, 2002. a

Cuny, J., Rhines, P. B., Niiler, P. P., and Bacon, S.: Labrador Sea Boundary Currents and the Fate of the Irminger Sea Water, J. Phys. Oceanogr., 32, 627–647,<0627:LSBCAT>2.0.CO;2, 2002. a, b

Cuny, J., Rhines, P. B., Schott, F., and Lazier, J.: Convection above the Labrador Continental Slope, J. Phys. Oceanogr., 35, 489–511,, 2005. a, b, c, d, e

de Boyer Montégut, C., Madec, G., Fischer, A. S., Lazar, A., and Iudicone, D.: Mixed layer depth over the global ocean: An examination of profile data and a profile-based climatology, J. Geophys. Res.-Ocean., 109, C12003,, 2004. a

de Jong, M. F., Bower, A. S., and Furey, H. H.: Two Years of Observations of Warm-Core Anticyclones in the Labrador Sea and Their Seasonal Cycle in Heat and Salt Stratification, J. Phys. Oceanogr., 44, 427–444,, 2014. a, b

de Jong, M. F., Oltmanns, M., Karstensen, J., and de Steur, L.: Deep Convection in the Irminger Sea Observed with a Dense Mooring Array, Oceanography, 31, 50–59, 2018. a

Dengler, M., Fischer, J., Schott, F. A., and Zantopp, R.: Deep Labrador Current and its variability in 1996–2005, Geophys. Res. Lett., 33, LS21S06,, 2006. a

Eden, C. and Böning, C.: Sources of Eddy Kinetic Energy in the Labrador Sea, J. Phys. Oceanogr., 32, 3346–3363,<3346:SOEKEI>2.0.CO;2, 2002. a, b

Fischer, J. and Schott, F. A.: Labrador Sea Water Tracked by Profiling Floats–From the Boundary Current into the Open North Atlantic, J. Phys. Oceanogr., 32, 573–584,<0573:lswtbp>;2, 2002. a

Fischer, J., Visbeck, M., Zantopp, R., and Nunes, N.: Interannual to decadal variability of outflow from the Labrador Sea, Geophys. Res. Lett., 37, L24610,, 2010. a

Fischer, J., Karstensen, J., Oltmanns, M., and Schmidtko, S.: Mean circulation and EKE distribution in the Labrador Sea Water level of the subpolar North Atlantic, Ocean Sci., 14, 1167–1183,, 2018. a, b, c

Flament, P.: A state variable for characterizing water masses and their diffusive stability: spiciness, Prog. Oceanogr., 54, 493–501,, 2002. a

Garcia, H. E. and Gordon, L. I.: Oxygen solubility in seawater: Better fitting equations, Limnol. Oceanogr., 37, 1307–1312,, 1992. a

Gebbie, G. and Huybers, P.: How is the ocean filled?, Geophys. Res. Lett., 38, L06604,, 2011. a

Georgiou, S., Ypma, S. L., Brüggemann, N., Sayol, J.-M., van der Boog, C. G., Spence, P., Pietrzak, J. D., and Katsman, C. A.: Direct and indirect pathways of convected water masses and their impacts on the overturning dynamics of the Labrador Sea, J. Geophys. Res.-Ocean., 126, e2020JC016654,, 2020. a, b, c, d, e, f, g

Handmann, P.: Deep Water Formation and Spreading Dynamics in the subpolar North Atlantic from Observations and high-resolution Ocean Models, Ph.D. thesis, Christian-Albrechts-Universität Kiel, Kiel, Germany, available at: (last access: 21 January 2022), 2019. a

Holte, J., Talley, L. D., Gilson, J., and Roemmich, D.: An Argo mixed layer climatology and database, Geophys. Res. Lett., 44, 5618–5626,, 2017. a, b, c, d

Isozaki, Y.: Permo-Triassic Boundary Superanoxia and Stratified Superocean: Records from Lost Deep Sea, Science, 276, 235–238,, 1997. a

Karstensen, J.: Calibration of physical data, Tech. Rep., Internal Report: Animate Report, available at: (last access: 21 January 2022), 2005. a

Karstensen, J., Stramma, L., and Visbeck, M.: Oxygen minimum zones in the eastern tropical Atlantic and Pacific oceans, Prog. Oceanogr., 77, 331–350,, 2008. a, b

Khatiwala, S., Primeau, F., and Holzer, M.: Ventilation of the deep ocean constrained with tracer observations and implications for radiocarbon estimates of ideal mean age, Earth Planet. Sc. Lett., 325–326, 116–125,, 2012. a

Kieke, D., Rhein, M., Stramma, L., Smethie, W. M., LeBel, D. A., and Zenk, W.: Changes in the CFC Inventories and Formation Rates of Upper Labrador Sea Water, 1997–2001, J. Phys. Oceanogr., 36, 64–86,, 2006. a

Koelling, J., Wallace, D. W. R., Send, U., and Karstensen, J.: Intense oceanic uptake of oxygen during 2014–2015 winter convection in the Labrador Sea, Geophys. Res. Lett., 44, 7855–7864,, 2017. a

Koltermann, K., Sokov, A., Tereschenkov, V., Dobroliubov, S., Lorbacher, K., and Sy, A.: Decadal changes in the thermohaline circulation of the North Atlantic, Deep-Sea Res. Pt. II, 46, 109–138,, 1999. a

Lazier, J.: The renewal of Labrador sea water, Deep-Sea Res. Oceanogr., 20, 341–353,, 1973. a, b, c, d

Le Bras, I. A., Yashayaev, I., and Toole, J. M.: Tracking Labrador Sea Water property signals along the Deep Western Boundary Current, J. Geophys. Res.-Ocean., 122, 5348–5366,, 2017. a

Le Bras, I. A.-A., Straneo, F., Holte, J., de Jong, M. F., and Holliday, N. P.: Rapid Export of Waters Formed by Convection Near the Irminger Sea's Western Boundary, Geophys. Res. Lett., 47, e2019GL085989,, 2020. a

Lebedev, K. V., Yoshinari, H., Maximenko, N. A., and Hacker, P. W.: Velocity data assessed from trajectories of Argo floats at parking level and at the sea surface, IPRC Technical Note, 4, 1–16,, 2007. a

Lozier, M. S.: Deconstructing the conveyor belt, Science, 328, 1507–1511,, 2010. a

Lozier, M. S., Bacon, S., Bower, A. S., Cunningham, S. A., de Jong, M. F., de Steur, L., deYoung, B., Fischer, J., Gary, S. F., Greenan, B. J. W., Heimbach, P., Holliday, N. P., Houpert, L., Inall, M. E., Johns, W. E., Johnson, H. L., Karstensen, J., Li, F., Lin, X., Mackay, N., Marshall, D. P., Mercier, H., Myers, P. G., Pickart, R. S., Pillar, H. R., Straneo, F., Thierry, V., Weller, R. A., Williams, R. G., Wilson, C., Yang, J., Zhao, J., and Zika, J. D.: Overturning in the Subpolar North Atlantic Program: A New International Ocean Observing System, Bull. Am. Meteorol. Soc., 98, 737–752, 2017. 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, b, c

MacGilchrist, G. A., Johnson, H. L., Marshall, D. P., Lique, C., Thomas, M., Jackson, L. C., and Wood, R. A.: Locations and Mechanisms of Ocean Ventilation in the High-Latitude North Atlantic in an Eddy-Permitting Ocean Model, J. Clim., 33, 10113–10131,, 2020. a, b, c

MacGilchrist, G. A., Johnson, H. L., Lique, C., and Marshall, D. P.: Demons in the North Atlantic: Variability of Deep Ocean Ventilation, Geophys. Res. Lett., 48, e2020GL092340,, 2021. a, b, c

Manabe, S. and Stouffer, R. J.: Coupled ocean-atmosphere model response to freshwater input: Comparison to Younger Dryas Event, Paleoceanography, 12, 321–336,, 1997. a

Marshall, J. and Schott, F.: Open-ocean convection: Observations, theory, and models, Rev. Geophys., 37, 1–64,, 1999. a, b

Maze, G., Mercier, H., Thierry, V., Memery, L., Morin, P., and Perez, F. F.: Mass, nutrient and oxygen budgets for the northeastern Atlantic Ocean, Biogeosciences, 9, 4099–4113,, 2012. a, b

Molinari, R. L., Fine, R. A., Wilson, W. D., Curry, R. G., Abell, J., and McCartney, M. S.: The arrival of recently formed Labrador sea water in the Deep Western Boundary Current at 26.5 N, Geophys. Res. Lett., 25, 2249–2252,, 1998. a

Naveira Garabato, A. C., MacGilchrist, G. A., Brown, P. J., Evans, D. G., Meijers, A. J., and Zika, J. D.: High-latitude ocean ventilation and its role in Earth's climate transitions, Philos. T. R. Soc. A, 375, 20160324,, 2017. a

Oltmanns, M., Karstensen, J., and Fischer, J.: Increased risk of a shutdown of ocean convection posed by warm North Atlantic summers, Nat. Clim. Change, 8, 300–304,, 2018. a

Oschlies, A.: A committed fourfold increase in ocean oxygen loss, Nat. Commun., 12, 1–8,, 2021. a

Pacini, A., Pickart, R. S., Bahr, F., Torres, D. J., Ramsey, A. L., Holte, J., Karstensen, J., Oltmanns, M., Straneo, F., Bras, I. A. L., Moore, G. W. K., and de Jong, M. F.: Mean Conditions and Seasonality of the West Greenland Boundary Current System near Cape Farewell, J. Phys. Oceanogr., 50, 2849–2871,, 2020. a, b, c

Palevsky, H. I. and Nicholson, D. P.: The North Atlantic Biological Pump: Insights from the ocean observatories initiative Irminger Sea Array, Oceanography, 31, 42–49, 2018. a, b

Palter, J. B., Lozier, M. S., and Lavender, K. L.: How Does Labrador Sea Water Enter the Deep Western Boundary Current?, J. Phys. Oceanogr., 38, 968–983,, 2008. a, b, c, d

Pickart, R. S., Spall, M. A., and Lazier, J. R.: Mid-depth ventilation in the western boundary current system of the sub-polar gyre, Deep-Sea Res. Pt. I, 44, 1025–1054,, 1997. a, b, c, d, e

Pickart, R. S., Torres, D. J., and Clarke, R. A.: Hydrography of the Labrador Sea during Active Convection, J. Phys. Oceanogr., 32, 428–457,<0428:HOTLSD>2.0.CO;2, 2002. a, b

Piron, A., Thierry, V., Mercier, H., and Caniaux, G.: Argo float observations of basin-scale deep convection in the Irminger sea during winter 2011–2012, Deep-Sea Res. Pt. I, 109, 76–90,, 2016. a, b

Rahmstorf, S.: Shifting seas in the greenhouse?, Nature, 399, 523–524,, 1999. a

Rhein, M., Fischer, J., Smethie, W. M., Smythe-Wright, D., Weiss, R. F., Mertens, C., Min, D.-H., Fleischmann, U., and Putzka, A.: Labrador Sea Water: Pathways, CFC Inventory, and Formation Rates, J. Phys. Oceanogr., 32, 648–665,<0648:LSWPCI>2.0.CO;2, 2002. a

Rhein, M., Walter, M., Mertens, C., Steinfeldt, R., and Kieke, D.: The circulation of North Atlantic Deep Water at 16 N, 2000–2003, Geophys. Res. Lett., 31, L14305,, 2004. a

Rieck, J. K., Böning, C. W., and Getzlaff, K.: The Nature of Eddy Kinetic Energy in the Labrador Sea: Different Types of Mesoscale Eddies, Their Temporal Variability, and Impact on Deep Convection, J. Phys. Oceanogr., 49, 2075–2094,, 2019. a

Roemmich, D., Johnson, G. C., Riser, S., Davis, R., Gilson, J., Owens, W. B., Garzoli, S. L., Schmid, C., and Ignaszewski, M.: The Argo Program: Observing the Global Ocean with Profiling Floats, Oceanography, 22, 34–43, 2009. a

Rogers, A. D.: Environmental Change in the Deep Ocean, Annu. Rev. Env. Resour., 40, 1–38,, 2015. a

Schmidtko, S., Stramma, L., and Visbeck, M.: Decline in global oceanic oxygen content during the past five decades, Nature, 542, 335–339,, 2017. a

Smith, W. H. F. and Sandwell, D. T.: Global Sea Floor Topography from Satellite Altimetry and Ship Depth Soundings, Science, 277, 1956–1962,, 1997. a

Straneo, F.: Heat and Freshwater Transport through the Central Labrador Sea, J. Phys. Oceanogr., 36, 606–628,, 2006. a, b, c

Straneo, F., Pickart, R. S., and Lavender, K.: Spreading of Labrador sea water: an advective-diffusive study based on Lagrangian data, Deep-Sea Res. Pt. I, 50, 701–719,, 2003. a, b

Sy, A., Rhein, M., Lazier, J. R., Koltermann, K. P., Meincke, J., Putzka, A., and Bersch, M.: Surprisingly rapid spreading of newly formed intermediate waters across the North Atlantic Ocean, Nature, 386, 675–679,, 1997. a

Talley, L. and McCartney, M.: Distribution and Circulation of Labrador Sea Water, J. Phys. Oceanogr., 12, 1189–1205,<1189:DACOLS>2.0.CO;2, 1982. a, b

Talley, L. D.: Freshwater transport estimates and the global overturning circulation: Shallow, deep and throughflow components, Prog. Oceanogr., 78, 257–303,, 2008. a

Tengberg, A. and Hovdenes, J.: Information on long-term stability and accuracy of Aanderaa oxygen optodes; information about multipoint calibration system and sensor option overview, Aanderaa Data Instruments AS Tech. Note, availlable at: (last access: 21 January 2022), 2014. a

Tengberg, A., Hovdenes, J., Andersson, H. J., Brocandel, O., Diaz, R., Hebert, D., Arnerich, T., Huber, C., Körtzinger, A., Khripounoff, A., Rey, F., Rönning, C., Schimanski, J., Sommer, S., and Stangelmayer, A.: Evaluation of a lifetime-based optode to measure oxygen in aquatic systems, Limnol. Oceanogr.-Method., 4, 7–17,, 2006. a

Toole, J. M., Andres, M., Le Bras, I. A., Joyce, T. M., and McCartney, M. S.: Moored observations of the Deep Western Boundary Current in the NW Atlantic: 2004–2014, J. Geophys. Res.-Ocean., 122, 7488–7505,, 2017. a

Tozer, B., Sandwell, D. T., Smith, W. H. F., Olson, C., Beale, J. R., and Wessel, P.: Global Bathymetry and Topography at 15 arcsec: SRTM15+, Earth Space Sci., 6, 1847–1864,, 2019. a

Ulses, C., Estournel, C., Fourrier, M., Coppola, L., Kessouri, F., Lefèvre, D., and Marsaleix, P.: Oxygen budget of the north-western Mediterranean deep-convection region, Biogeosciences, 18, 937–960,, 2021.  a

Wolf, M. K., Hamme, R. C., Gilbert, D., Yashayaev, I., and Thierry, V.: Oxygen Saturation Surrounding Deep Water Formation Events in the Labrador Sea From Argo-O2 Data, Global Biogeochem. Cy., 32, 635–653,, 2018. a

Yashayaev, I.: Hydrographic changes in the Labrador Sea, 1960–2005, Prog. Oceanogr., 73, 242–276,, 2007. a, b, c, d, e

Yashayaev, I. and Loder, J. W.: Recurrent replenishment of Labrador Sea Water and associated decadal-scale variability, J. Geophys. Res.-Ocean., 121, 8095–8114,, 2016. a, b, c, d, e, f

Yashayaev, I. and Loder, J. W.: Further intensification of deep convection in the Labrador Sea in 2016, Geophys. Res. Lett., 44, 1429–1438,, 2017. a, b

Yeager, S., Castruccio, F., Chang, P., Danabasoglu, G., Maroon, E., Small, J., Wang, H., Wu, L., and Zhang, S.: An outsized role for the Labrador Sea in the multidecadal variability of the Atlantic overturning circulation, Sci. Adv., 7, eabh3592,, 2021. a

Zantopp, R., Fischer, J., Visbeck, M., and Karstensen, J.: From interannual to decadal: 17 years of boundary current transports at the exit of the Labrador Sea, J. Geophys. Res.-Ocean., 122, 1724–1748,, 2017. a, b, c, d, e, f, g, h, i

Zou, S., Lozier, M. S., Li, F., Abernathey, R., and Jackson, L.: Density-compensated overturning in the Labrador Sea, Nat. Geosci., 13, 121–126,, 2020. a, b, c

Zunino, P., Mercier, H., and Thierry, V.: Why did deep convection persist over four consecutive winters (2015–2018) southeast of Cape Farewell?, Ocean Sci., 16, 99–113,, 2020. a, b

Short summary
In this study, we investigate oxygen variability in the deep western boundary current in the Labrador Sea from multiyear moored records. We estimate that about half of the oxygen taken up in the interior Labrador Sea by air–sea gas exchange during deep water formation is exported southward the same year. Our results underline the complexity of the oxygen uptake and export in the Labrador Sea and highlight the important role this region plays in supplying oxygen to the deep ocean.
Final-revised paper