Articles | Volume 19, issue 12
Research article
21 Jun 2022
Research article |  | 21 Jun 2022

Wintertime process study of the North Brazil Current rings reveals the region as a larger sink for CO2 than expected

Léa Olivier, Jacqueline Boutin, Gilles Reverdin, Nathalie Lefèvre, Peter Landschützer, Sabrina Speich, Johannes Karstensen, Matthieu Labaste, Christophe Noisel, Markus Ritschel, Tobias Steinhoff, and Rik Wanninkhof

The key processes driving the air–sea CO2 fluxes in the western tropical Atlantic (WTA) in winter are poorly known. WTA is a highly dynamic oceanic region, expected to have a dominant role in the variability in CO2 air–sea fluxes. In early 2020 (February), this region was the site of a large in situ survey and studied in wider context through satellite measurements. The North Brazil Current (NBC) flows northward along the coast of South America, retroflects close to 8 N and pinches off the world's largest eddies, the NBC rings. The rings are formed to the north of the Amazon River mouth when freshwater discharge is still significant in winter (a time period of relatively low run-off). We show that in February 2020, the region (5–16 N, 50–59 W) is a CO2 sink from the atmosphere to the ocean (−1.7Tg C per month), a factor of 10 greater than previously estimated. The spatial distribution of CO2 fugacity is strongly influenced by eddies south of 12 N. During the campaign, a nutrient-rich freshwater plume from the Amazon River is entrained by a ring from the shelf up to 12 N leading to high phytoplankton concentration and significant carbon drawdown (∼20 % of the total sink). In trapping equatorial waters, NBC rings are a small source of CO2. The less variable North Atlantic subtropical water extends from 12 N northward and represents ∼60 % of the total sink due to the lower temperature associated with winter cooling and strong winds. Our results, in identifying the key processes influencing the air–sea CO2 flux in the WTA, highlight the role of eddy interactions with the Amazon River plume. It sheds light on how a lack of data impeded a correct assessment of the flux in the past, as well as on the necessity of taking into account features at meso- and small scales.

1 Introduction

The North Brazil Current (NBC) is one of the dominant features of the tropical Atlantic circulation. In a region dominated by zonal jets, it flows northward along the coast of South America and separates from the coast around 6–8 N. The NBC seasonally turns back on itself in a tight loop, called the NBC retroflection, and feeds the North Equatorial Counter Current, closing off the equatorial wind-driven gyre (Fig. 1). This retroflection occasionally pinches off some of the world's largest eddies, the North Brazil Current rings (Johns et al., 1990; Richardson et al., 1994).

Figure 1Schematic of the main ocean currents in the western tropical Atlantic superimposed over the SSS field of 7 February 2017 (a) and over the February ΔfCO2 climatology from Landschützer et al. (2020) (b).

After their separation from the NBC retroflection region, the rings travel northwestward toward the Lesser Antilles in a course parallel to the coast of South America. These eddies have been extensively studied using modelling and both in situ (e.g. 1998–2001 NBC Ring experiment; Wilson et al., 2002) and satellite observations (e.g. Goni and Johns, 2001; Fratantoni and Glickson, 2002; Aroucha et al., 2020). They have a mean radius of 200 km, and their diameter can exceed 450 km. Vertically, some of them extend down to more than 1000 m (Fratantoni and Glickson, 2002; Fratantoni and Richardson, 2006; Johns et al., 2003). The NBC rings swirl clockwise and travel with an average northwestward translation speed of 8–15 km d−1 (Johns et al., 1990; Mélice and Arnault, 2017). Most of the anticyclonic eddies detectable by altimetry are rather shallow, extending from the surface to 200–300 m (Garraffo et al., 2003; Wilson et al., 2002). When they reach the Lesser Antilles, they start to coalesce and disintegrate partly due to interactions with the topography (Fratantoni and Richardson, 2006; Jochumsen et al., 2010). There is substantial variability in the number of rings shed per year, ranging from five (Aroucha et al., 2020; Fratantoni and Glickson, 2002; Mélice and Arnault, 2017; Goni and Johns, 2001) to nine (Johns et al., 2003). NBC rings play a crucial role in the inter-hemispheric transport of salt and heat in the Atlantic Ocean and are an important part of the meridional overturning circulation (Johns et al., 2003). The NBC rings disrupt an already complex region located in the vicinity of the Amazon River mouth and at the transition between equatorial and subtropical waters. While most of the studies on rings focused on their physical properties, little is known about their biogeochemical properties and how they affect the air–sea CO2 flux of the western tropical Atlantic.

The global ocean acts as an atmospheric CO2 sink, taking up 23 % of total anthropogenic CO2 emissions (Friedlingstein et al., 2020) and leading to ocean acidification (Masson-Delmotte et al., 2021; Pörtner et al., 2019). The concentration of atmospheric CO2 is increasing due to human activities (IPCC, 2019, 2021), and characterizing the role of the ocean in mitigating climate change through CO2 uptake is thus a key investigation. The equatorial Atlantic Ocean is the second-largest source of CO2 to the atmosphere after the equatorial Pacific (Landschützer et al., 2014; Takahashi et al., 2009). Previous works in this region examined the influence of the equatorial upwelling and of the Amazon plume on the CO2 flux. CO2-rich equatorial waters, originating from the equatorial upwelling (Andrié et al., 1986), strongly contrast with the CO2-undersaturated Amazon River plume waters. The magnitude of the Amazon River discharge is unique in the global ocean. It represents as much freshwater as the next seven largest rivers in the world combined and contributes to almost 20 % of global river freshwater input to the ocean (Dai and Trenberth, 2002). It therefore strongly impacts the physical, biogeochemical and biological properties of the coastal and the open ocean. Often overlooked, the Amazon River plume is an atmospheric CO2 sink of global importance (Ibánhez et al., 2016). The plume carries water rich in silicate, nitrogen and phosphate into the tropical oceanic waters that are strongly depleted in nutrients. As water mixes and turbidity decreases, the primary producer's growth and associated biological drawdown are stimulated (Chen et al., 2012). Nitrogen is rapidly consumed, and nitrogen fixation by diazotrophs becomes the main pathway of carbon sequestration in the plume (Subramaniam et al., 2008). This strong carbon drawdown leads to a significant sink of atmospheric CO2 (Körtzinger, 2003; Lefèvre et al., 2010). Not taking into account the Amazon plume would result in overestimating the tropical Atlantic air–sea CO2 flux by 10 % (Ibánhez et al., 2016).

The Amazon River's discharge reaches a minimum in December and progressively increases from January onwards. The plume extension is minimum from January to March (Fournier et al., 2015), and as a result, it is the period of maximum salinity in the northwestern tropical Atlantic. The Amazon outflow region is particularly hard to reconstruct due to its strong variability and a severe lack of data. Waters located in the southeasternmost part of the domain act as a strong source of CO2 to the atmosphere. The source gradually turns into a sink north of 10 N as waters become colder due to seasonal winter cooling. This situation is typical of a transition zone between equatorial and subtropical waters in winter (Landschützer et al., 2020; Fig. 1b). The northwestern tropical Atlantic is commonly divided into two parts: the northern much less variable part (also called “Trade wind region”) and the southern part, also referred to as Eddy Boulevard (Stevens et al., 2021). The freshwater of the Amazon remains mainly confined to the continental shelf due to winds perpendicular to the coast as it travels northwestward into the Caribbean Sea (Coles et al., 2013). However, it has recently been documented that off-shore freshwater transport is often present in February and significantly alters the physical properties of the region (Reverdin et al., 2021). This is partly due to the interaction of the NBC rings with the Amazon plume (Fig. 1a). The ocean colour signature of the Amazon (Muller-Karger et al., 1988) has been used as a tracer to delineate the rings and better understand their generation, evolution and characteristics (Johns et al., 1990; Fratantoni and Glickson, 2002). The Amazon River also influences the surface temperature and salinity of the rings. For example, Fig. 1a shows a freshwater plume stirred by a large ring. They are considered warm-core rings but have a warm sea surface temperature (SST) anomaly in the first half of the year and a cold one in the second half because the anomaly is relative to the regional SST, with an extensive warm pool in late summer and autumn (Ffield, 2005). Their signature in salinity is therefore plume-dependent as well. Ffield (2005) reported that three out of four rings were surrounded by lower-salinity water. Salinity and chlorophyll a are therefore critical to understand the surface physical and biogeochemical properties of the region, as well as the air–sea fluxes of CO2.

The northwestern tropical Atlantic is a dynamically active region, with eddies several hundred kilometres in diameter and connected to the world's largest river. There are surprisingly few biogeochemical observations available for winter months during low outflow of the Amazon River. Few tropical Atlantic measurements of biogeochemical tracers are available, with one transect in winter in the Surface Ocean CO2 Atlas (SOCAT; Bakker et al., 2016) database south of 10 N crossing the region. This scarcity is a major impediment in understanding the biological and physical processes underlying the oceanic carbon and nutrient cycles in the region. Satellite salinity shows a contrasted spatial structure with eddies and filaments (Fig. 1a). In this study, we take advantage of the physical and biogeochemical data collected during the EUREC4A-OA/ATOMIC (Elucidating the Role of Clouds–Circulation Coupling in Climate Ocean–Atmosphere and Atlantic Tradewind Ocean–Atmosphere Mesoscale Interaction Campaign) experiment in January–February 2020, combined with satellite data, to understand how the NBC rings and their related structures impact the air–sea CO2 flux in winter. The paper is organized as follows. We present the in situ observational data from the EUREC4A-OA/ATOMIC experiment, as well as the satellite data, in Sect. 2. We identify the water masses observed in the region and their physical and biogeochemical properties and estimate the CO2 fluxes at regional scale using empirical relationships in Sect. 3. In Sect. 4, we compare the results with climatologies of air–sea CO2 fluxes to evaluate the added knowledge brought by the intensive surveys of February 2020. We discuss the results and the inter-annual variability in Sect. 4, and we conclude in Sect. 5.

2 Data and methods

2.1 In situ data

The EUREC4A-OA/ATOMIC (Stevens et al., 2021) campaign took place in January and February 2020 and involved research vessels (RVs) from France (RV Atalante; Speich and The Embarked Science Team, 2021), Germany (RV Maria S. Merian, hereby designated as Merian; Karstensen et al., 2020; and RV Meteor, not considered in this study since no CO2 measurements were taken on board) and the United States (RV Ronald H. Brown, hereby designated as Ron Brown; Quinn et al., 2021). These cruises provided numerous in situ measurements, and, in this study, we will focus on the continuous near-surface measurements of temperature, salinity and fCO2.

Temperature and salinity from thermosalinographs (TSGs), as well as fCO2, were measured from water pumped ∼5 m below the surface. For each ship, the resulting CO2 data are corrected (Lefèvre et al., 2010; Pierrot et al., 2009) from the temperature difference between the water at the ship's water intake and the one analysed by the instrument. RV Atalante fCO2 measurements started on 30 January and ended on 18 February 2020 (Olivier et al., 2020). The underway oceanic and atmospheric fCO2 were detected by infrared detection using a LI-COR LI-7000 (Takahashi et al., 1993). The fCO2 system was the same as in Lefèvre et al. (2010). It uses a shower air–sea equilibrator described by Poisson et al. (1993). Seawater from the TSG pumping circuit circulates in the equilibrator at a rate of 2 L min−1. A closed loop of about 100 mL of air flows through the equilibrator designed to avoid bubbles at the air–sea interface. To minimize temperature corrections, the equilibrator is thermostated with the same seawater as the one used for CO2 measurements. The temperature difference between the equilibrator and the sea was on the order of 0.5 C.

Furthermore, 138 samples for dissolved inorganic carbon (DIC) and total alkalinity (TA) analysis were collected on board RV Atalante, as well as inorganic nutrients (silicate, phosphate, nitrate and nitrite). DIC and TA were measured at the SNAPO-CO2 (French National Facility for Analysis of Carbonate System Parameter) facility by potentiometric titration using a closed cell, following the method of Edmond (1970). Nutrients were conserved by heat pasteurization and analysed by colourimetry at IRD LAMA service in Brest.

An OceanPack CUBE FerryBox system from SubCtech was installed on RV Merian measuring continuously the oceanic fCO2 from 23 January to 19 February 2020. Water is pumped at a rate of ∼7L min−1 through a debubbler unit subsequently followed by a SeaBird SBE 45 thermosalinograph before it circulates along a membrane through which CO2 diffuses. On the other side of the membrane, the air loop is circulated at a rate of 0.5 L min−1 through a LI-COR LI-840 non-dispersive infrared gas analyser (e.g. Arruda et al., 2020). RV Ron Brown completed two legs, from 10 to 25 January and from 29 January to 13 February 2020. The General Oceanic Inc. 8500 pCO2 instrument installed on board follows a similar methodology to the underway fCO2 system deployed on the Atalante and is detailed in Pierrot et al. (2009).

An inter-comparison of the fCO2 measured by RVs Atalante and Merian is attempted when the ships were navigating together at a distance of less than 5 km (Fig. 2). On average, the fCO2 measured by RV Merian is 6.4 µatm higher than the one measured on RV Atalante, with a standard deviation of 4.8 µatm. RVs Merian and Ron Brown crossed the same water mass at 13–14 N, 57 W on 12 February. On average, RV Merian fCO2 is 6 µatm higher than RV Ron Brown fCO2. In part we link these differences to the slower response time of the membrane system; however, differences also lie within the uncertainties of the fCO2 observing systems (∼5µatm for the membrane system installed on RV Merian, see Arruda et al., 2020, and ∼2µatm for the equilibrator systems installed on RVs Atalante and Ron Brown). The region where RV Merian and RV Atalante were close is very variable (standard deviation of 20 µatm), and RV Merian and RV Ron Brown were never in the same place at the same time, so the observed differences could also be due to the natural variability in fCO2 sampled differently by the various ships. Hence, we did not apply any correction, and we checked that the effect of a 6 µatm systematic bias in RV Merian fCO2 had a minor effect on our resulting interpolations. A comparison between the reconstructed flux with and without a correction for the 6 µatm systematic bias suggests that such a systematic bias would lead to less than 2 µatm difference in our mean interpolated fCO2 and less than 0.1 mmolm-2d-1 in the mean derived air–sea CO2 flux.

Figure 2Ship tracks colour-coded by day of year (a) and by ship name (b).

2.2 Satellite and atmospheric reanalysis data

Daily satellite maps of chlorophyll-a (chl a) and SST, as well as absolute dynamic topography (ADT) and sea surface salinity (SSS), are used in this study.

The salinity maps are a blend of the Soil Moisture Ocean Salinity (SMOS, January 2010–present) and Soil Moisture Active Passive (SMAP, April 2015–present) measurements developed by Reverdin et al. (2021). The European SMOS and US SMAP missions observe the sea surface by L-band radiometry from sun-synchronous polar-orbiting satellites (Entekhabi et al., 2010; Font et al., 2009; Kerr et al., 2010; Piepmeier et al., 2017). Combining 06:00 and 18:00 local time measurements of both missions provides an almost complete coverage each day. When the coverage was not complete over our region, the 06:00 track of the following day was also included. This daily field is available for the first 20 d of February and leaves out only 2 d without sufficient coverage to retrieve salinity data. It has a spatial resolution close to 70 km and an uncertainty on the order of 0.5. This product is optimized for the northwestern tropical Atlantic in February 2020 and has an almost daily resolution. It is an experimental daily product built to have the best representation possible of the Amazon plume variability. The product, its uncertainties and the comparison between the TSG salinity and the satellite SSS are detailed in Reverdin et al. (2021).

Daily chl a concentration maps and SST maps are produced by CLS (Stum et al., 2016) on a spatial grid of 0.02. The chl a concentration maps are composites built from VIIRS (Visible Infrared Imaging Radiometer Suite; on Suomi NPP (Suomi National Polar-orbiting Partnership) and NOAA-20 US platforms) and OLCI (Ocean and Land Colour Instrument; on Sentinel 3A and 3B Copernicus European platforms) satellite sensors. The SST product is a 1 d average of infrared radiometer satellite data. Both datasets are sensitive to the cloud cover, but during our period of interest, they are usually without many gaps except at the end of February. The comparison between the TSG SST and the satellite SST product is detailed in RV Atalante's cruise report (Speich and The Embarked Science Team, 2021).

Daily ADT maps at a 1/4 resolution combine data from all satellites available for the period 1993 to present. From these ADT fields, the TOEddies algorithm, developed by Laxenaire et al. (2018), identifies eddies and their trajectories. The eddy detection is based on the closed contours of ADT, as well as the maximum geostrophic velocity associated with the eddy.

In order to compute the air–sea CO2 fluxes, the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) hourly wind speed and mean sea level pressure, Patm, are used. ERA5 covers the period from January 1950 to present and provides hourly data on a 30 km grid. In addition, the monthly wind speed and SST fields over the period 1998–2015 are used, and a climatology over this period is computed. The wind speed in the region in winter is on average between 6 and 8 m s−1, and its variability is low.

We compare the EUREC4A-OA/ATOMIC observations with the observation-based CO2 partial pressure (pCO2) climatology developed by Landschützer et al. (2020), created using a two-step neural network method (Landschützer et al., 2016) and combining open and coastal ocean datasets. The associated ΔpCO2 and air–sea CO2 flux monthly field climatologies over the 1998–2015 period are computed using the ERA5 climatological wind, SST and Patm fields, as well as the atmospheric CO2 from the Ragged Point, Barbados, station.

2.3 Methods

2.3.1 Air–sea CO2 flux

We compute the air–sea flux (F; mmolm-2d-1) as follows:

(1) F = k K 0 ( f CO 2 - f CO 2atm ) ,

where K0 is the solubility of CO2 in seawater, expressed as a function of SSS and SST by Weiss (1974); fCO2atm is the atmospheric CO2 fugacity; and k is the gas transfer velocity. k is calculated following the relation from Wanninkhof (2014):

(2) k = 0.251 U 2 ( Sc / 660 ) - 0.5 ,

where Sc is the Schmidt number, and U is the wind speed at 10 m above sea level derived from ERA5 wind speed. The ERA5 wind speed is used for the satellite-based analysis and the air–sea CO2 flux climatology. The measured winds from the ships are adjusted to 10 m following a logarithmic profile (Tennekes, 1973) and used to compute the along-track flux for visualization purposes.

In order to compute fCO2atm over that period, we first derived the saturation vapour pressure (PH2O) from SSS and SST and then the atmospheric pCO2 using the monthly averaged CO2 mole fraction (xCO2atm) measured at the NOAA/Earth System Research Laboratory (ESRL) station in Ragged Point, Barbados (13.17 N, 59.43 W):

(3) f CO 2atm = x CO 2atm ( P atm - P H 2 O ) C f ,

where Cf is the fugacity coefficient, a function of the atmospheric pressure and SST (Weiss, 1974).

Figure 3(a) Chlorophyll a, (b) SSS and (c) SST on 6 February 2020 with the contours of NBC rings A1 and A2, their centre, and their trajectory. The NBC retroflection is identified from the 0.51 m contour of the satellite-derived ADT.

2.3.2 Reconstruction of fCO2 from satellite maps

Our approach is to derive from the large EUREC4A-OA dataset a relationship linking fCO2 to SST, SSS and chl a in order to provide maps of fCO2 based on the satellite maps of SST, SSS and chl a. Chl a was not measured on board; thus, we use satellite surface chl a co-located along the ship track. This set of parameters is used as proxy to describe the influence of ocean dynamics, chemistry and also marine biology on fCO2. The biological carbon pump is one of the major components of the oceanic and global carbon cycles as the photosynthetic production of organic carbon by marine phytoplankton accounts for about half of the carbon fixation associated with global primary production (Arrigo, 2007; Behrenfeld et al., 2006; Field et al., 1998). However, while chl a is an indicator of biological activity, it is also a very good tracer of ocean circulation, so, depending on the water masses origin, fCO2 and chl a are not expected to be systematically negatively correlated. Waters rich in detrital material tend to limit the phytoplankton growth, and microbial respiration of riverine material on the continental shelf likely dominates (Aller and Blair, 2006; Medeiros et al., 2015; Mu et al., 2021). Even considering the extent of the EUREC4A-OA atomic cruise, the dataset is still sparse and cannot fully represent the small-scale variability it highlights. In order to understand the fluxes at regional scale the need for a good spatial resolution arises. For that, the surface T (temperature), S (salinity) and chl a diagram computed from the ship measurements (and co-located satellite chl a) is interpolated using a linear 3D interpolation on a grid of SST, SSS and chl a. The grid has a resolution of 0.01 C in SST, 0.1 in SSS and 0.01 in log(chl a (mg m−3)). Using a 3D linear interpolation to map the fCO2 data over a grid is a simple yet effective solution for a dataset that is still relatively sparse. The method is presented in more details in the Supplement (Text S1 in the Supplement). Using a linear fit prevents oscillations between two data points and yields good results. Along the ship track, the standard deviation between the measured and reconstructed fCO2 is ∼4µatm. Each triplet of surface T, S and log(chl a) in the range of the values measured by the ship is therefore associated with a value of fCO2 based on the 3D linear interpolation of the in situ values. In order to cover the whole range of TS–log(chl a) present in the region, we extrapolate to lower temperatures and lower salinities than the ones measured by the ship. In order to do so, we add four points to the TS–chl a diagram at lower salinities and lower temperatures based on previous knowledge of the region. For the low-salinity domain (SSS<30), fCO2 is strongly dominated by salinity, and the influence of temperature is weak (Lefèvre et al., 2010). The SSS–fCO2 relation developed by Lefèvre et al. (2010) is in good agreement with the SSS–fCO2 relationship computed from this study data (Fig. S1 in the Supplement) in the common range, and we therefore use it to compute fCO2 at a salinity of 26 (fCO2(S=26)=251.4µatm). The lower temperature is mostly located in the northern part of the domain, which is the least variable and where the variations in fCO2 are dominated by the ones in temperature. From this dataset we compute a variation of 15 µatmC-1, which is consistent with the 4.23 % C-1 expected variation in fCO2 with temperature due to the temperature sensitivity of the carbonate dissociation constants and CO2 solubility (Takahashi et al., 2002; Wanninkhof et al., 1999). We use this dependency to compute the fCO2 at a temperature of 24 C to cover the whole range of temperature in the region.

We combine the interpolated fCO2 with satellite maps of SST, SSS and chl a to obtain daily high-resolution maps of fCO2. Some days, either the presence of clouds altering the chl a and SST or the lack of salinity coverage prevents the retrieval of fCO2. In order to limit the error in fCO2, we only keep 9 d out of the first 20 days of February (2, 4, 6, 7, 9, 11, 12, 17 and 19 February) when the coverage is sufficient. Then, daily mean sea level pressure maps and wind fields are used to compute the air–sea CO2 flux over the region in a similar way as described in Sect. 2.3.1. The salinity maps present major errors near islands because no correction of the island effect was applied on the SMAP maps (Grodsky et al., 2018). Therefore, the reconstructed flux will be studied over a region excluding the close vicinity of the islands (5–16 N, 59–50 W).

3 Results

3.1 A transition region presenting strong mesoscale activity

Figure 3 shows how in February 2020 the ocean currents of the western tropical Atlantic (WTA) strongly influenced the variability in SSS, SST and chl a at many scales, with an NBC ring stirring a plume of fresher water rich in chlorophyll a toward the open ocean. This is also supported by the measurements done during the EUREC4A-OA/ATOMIC campaign in January–February 2020 (Fig. 4). They show a complex environment – for example ΔfCO2 presents similar large-scale features as the climatology – but it also reveals numerous smaller-scale structures (Fig. 4a). Among the latter, two stand out. These are the very low fCO2 in the southeastern part of the domain and the high fCO2 around 11 N. In early 2020, the NBC retroflection was very variable and shed two large anticyclonic rings (Fig. 3). They are long-lived 250 km large eddies travelling northwestward toward the Lesser Antilles in the Eddy Boulevard region. The ring detection algorithm TOEddies based on ADT indicates that NBC ring A2 separated from the retroflection in late December, was fairly stationary during the cruise period (February 2020) and was located around 11 N, 58 W. NBC ring A1 separated from the retroflection in early February 2020 and then stayed around 10 N, 54 W for 10 d before translating northwestward toward the Lesser Antilles after 20 February. These eddies contribute to the variability in the region in two ways. As they travel, they transport the water trapped in their core during their formation (eddy trapping), but they also stir the surrounding waters (eddy stirring).

Figure 4In situ measurements of (a) ΔfCO2, (b) salinity and (c) temperature.

3.2 Surface water mass identification

In an effort to understand how biogeochemistry is forced by physical processes in the ocean, we used surface chl a to complement SST and SSS data in defining surface water masses. We observed that in the northwestern tropical Atlantic in winter in situ fCO2 strongly depends on these three variables (Figs. 4 and 5). There is a strong positive dependence of fCO2 on SSS, with low fCO2 for low SSS (Fig. 5c and d). Across the whole EUREC4A-OA/ATOMIC region, SST did not vary much (mean SST of 27 C and standard deviation of 0.5 C), but warmer waters present higher fCO2. The dependence on chl a allows for the discrimination of water masses with the same surface T and S properties but not the same fCO2 (Fig. 5). Satellite-based chl a is hard to discriminate from detrital material using ocean colour where both are present as they have close spectral characteristics. Figure 5 shows that waters with an SST of 26.5 C and SSS between 35–36 can either be rich in chl a and have a high fCO2 or be low in chl a and have a low fCO2. By combining SST and SSS with chl a and using information from the dynamical structures of the region, we identified six upper-ocean water masses (Fig. 5a and b). The along-track ΔfCO2 for each ship is presented in Fig. 6, colour-coded with the identified water mass, highlighting the link between the surface TS–chl a relation and ΔfCO2. The way we defined water masses, considering time-varying boundaries, is relatively similar to the one used by Longhurst (2010) and some of the surface water masses compare well with Longhurst's (2010) biogeochemical provinces. He identified three provinces in the northwestern tropical Atlantic: the North Atlantic Tropical Gyre province (NATR), the Western Tropical Atlantic province (WTRA) and the Guiana Coastal province. In this first part, we will present two surface water masses that are usually identified in the region (e.g. Longhurst et al., 1995, 2010) and their physical properties.

Figure 5(a) Map representing Atalante, Merian and Ron Brown ship tracks colour-coded with the identified water masses. (b) TS diagram colour-coded with the water masses; the grey colour corresponds to points that do not fit into the definition of the identified water masses.

North of Barbados, the domain is mostly dominated by North Atlantic Subtropical Waters (NASW), similar to Longhurst's NATR. They have a SSS in the range of 35 to 36 and are relatively cold (Table 1). Their SST diminishes over time towards the end of February. These waters are less influenced by coastal dynamics and therefore are not very productive at the surface (chl a levels inferior to 0.14 mg m−3) due to low nutrient levels. They are mainly located north of 13 N and get progressively colder toward the northeast. RV Ron Brown stayed in that Trade Wind region for almost a month (Fig. 5). The observations collected from this ship show lower fCO2 in NASW with respect to the atmosphere (ΔfCO2-40µatm). Similar results are found on one of RV Merian's transects, with ΔfCO2<-30µatm.

Table 1Thresholds in SSS, SST and chl a used to define the six water masses identified in the winter WTA.

Download Print Version | Download XLSX

The NBC is surface-intensified and fed by the central branch of the South Equatorial Current (SEC; Schott et al., 1998). As the cold and saline water from the upwelling region is transported westward by the SEC, it warms up (SST>27C) but retains its saline characteristic (SSS>36, Table 1) as it reaches the NBC retroflection region (Fig. 4). This NBC water mass is oligotrophic (Fig. 3b) and therefore in our area of interest distinguished itself by its low level of surface chl a (chl a<0.14mg m−3). These waters are found in the retroflection area and were sampled by both RVs Merian and Atalante at the beginning of February (Fig. 5). The NBC water mass is in some way similar to the WTRA but does not extend beyond the retroflection, therefore representing a more limited part of the WTRA. We will introduce four new water masses in the following parts, as well as their associated dynamical structures. They can be considered as subsets of the WTRA and Guiana Coastal provinces, these two provinces being too large and not well suited to represent small-scale processes.

3.3 North Brazil Current rings

The extension of the NBC retroflection varies depending on the state of eddy formation. It moves northwestward up until 9 N as an eddy is forming and then retracts to the southeastern part of the region. During our period of interest, the retroflection shed anticyclone A1 at the beginning of February. It is difficult to estimate the date of shedding as the area is highly dynamic, and detecting the first closed contour of ADT is complicated and may be inaccurate. It is however interesting that the two ships sampled the retroflection when it was expanding to generate A1, and this northwestward expansion is well observed on several physical and biogeochemical parameters (Fig. 7a). RV Merian crossed the retroflection on 27 January and stayed in the area until 2 February (Fig. 6). Chl a present on the shelf is advected by the strong currents on the periphery of the retroflection and delineates well its southwestern side (Fig. 7a). The NBC waters stand out on the surface TS diagram as they are the most saline waters observed in the region (Fig. 7b). They are also high in fCO2 which reflects their equatorial origin. Their SST is relatively warm, varying from 27.8 C at the crossing of the first retroflection front to 27.2 C. The region is rather homogeneous, with an almost constant SSS of 36.3 and ΔfCO2 along the multiple crossings, as observed on the Merian and Atalante transects (Fig. 6). On average along those transects, the NBC fCO2 is higher than fCO2atm by 20 µatm.

Figure 6RVs Atalante (a), Merian (b) and Ron Brown (c) ΔfCO2 time series. The background colour indicates the crossed water mass domains (see definition in legend of Fig. 5).


Anticyclone A1 is further crossed by RV Atalante on 6 February, just a few days after its separation from the retroflection (Fig. 7c). The surface signal is almost lost, both in SST and in fCO2 (Fig. 7c and d). It is mainly composed of modified NBC water, whose properties are close to the NBC water (high SSS, high SST, low chl a) but not as pronounced. This water mass covers a larger area, which mainly encompasses the Eddy Boulevard region. It is defined here as SSS>35.6, 27.16C<SST<27.6C and 0.11mgm-3<chla<0.25mg m−3 (Table 1). While the high chl a water delimits well the retroflection area, it partly covers eddy A1.

Figure 7(a) RVs Atalante and Merian ship tracks in the NBC retroflection (Merian: 27 January to 2 February; Atalante: 2 February), (c) in NBC ring A1 (Atalante: 6 February) and (e) in NBC ring A2 (Atalante: 12–13 February; Merian: 13–14 February) colour-coded with fCO2. The background represents the chl a on 2 (a), 6 (c) and 12 February (e), and the contours of NBC rings A1 and A2 are indicated. (b, d and f) Corresponding TS diagrams colour-coded with fCO2.

NBC ring A2 presents a different situation. Detached from the retroflection in early December 2019 (as defined from altimetry detection), it travelled northwestward while retaining an intense coherent core. Coastal waters identified by their high chl a content were less present and mostly entrained at the northwestward edge of the eddy. After 2 months, A2 almost reached Trinidad and Tobago and was located around 11 N, 58 W when it was sampled by the two ships (Fig. 7e). The SST signal is eroded, and most of the eddy is mostly made of modified NBC water with relatively low fCO2. However, high SSS (36.5) and fCO2 (415 µatm) are still visible near the eddy centre on the two crossings made by RV Atalante on 12 and 13 February. This is confirmed by the two sections of the Merian that crossed the altimetric eddy centre and measured SSS at 0.5 higher in the 50 km radius around the centre and above 36. The NBC water mass is therefore found close to the centre of eddy A2, as well as its associated high fCO2.

From the collected observations, it appears that the surface signature of NBC rings is relatively variable and complex. It is well marked in their formation area in the NBC retroflection, where waters brought north by the NBC are warmer, saltier and higher in fCO2 than the water of the northwestern tropical Atlantic. As the eddies travel northwestward, further away from the retroflection, they may be subject to various processes that modify the surface signal. Unfortunately, the data collected are not sufficient to shed light on which processes are involved in this situation. South of Barbados, away from the retroflection, the modified NBC is therefore the most common water mass. Nevertheless, the NBC water is still sometimes observed months after the separation from the retroflection in the eddy centre.

3.4 Freshwater plume

The NBC rings form and evolve in an area highly influenced by the Amazon River plume. Even if February is a period of low Amazon River outflow (Dai and Trenberth, 2002), freshwater events are relatively common. In February 2020, a freshwater plume detached from the Guiana plateau and spread out into the northwestern tropical Atlantic. The off-shelf plume was steered northward by the retroflection and NBC ring A1 up to 12 N and then extended westward toward the Caribbean Sea. Waters carried by the plume strongly contrast with the saline waters of the retroflection. They include water from the Amazon and present low SSS (SSS<34.5), low fCO2 (fCO2<380µatm) and high chl a (Table 1). The plume was crossed three times: twice on 2 February and once on 5 February (Fig. 6). Freshwater from the Amazon arrived on the plateau on 1–2 February and was then entrained northwestward by Ekman transport and geostrophic currents (Reverdin et al., 2021). On 2 February RVs Atalante and Merian left the retroflection area to cross the adjacent nascent plume. SSS rapidly decreased, reaching 33, which is associated with a strong decrease in fCO2 (Fig. 8). From 2 to 5 February, the plume formed, and on 5 February the plume was approximately 100 km wide, with the lowest salinities around 30. Based on satellite SSS data of the following days, the plume appears to have reached even lower salinities and then spread out over the northwestern tropical Atlantic. It first spread northward, steered by A1, and then northwestward, channelled between A1 and A2 and reaching all the way up to 12 N and extending over more than 100 000 km2 (Reverdin et al., 2021). The plume can be followed by satellite SSS and chl a maps (Fig. 7a, c and e). Indeed, the low SSS is also accompanied by high chl a as water from the Amazon is considered highly productive. The northwestern tropical Atlantic is in general nutrient-limited, but the nutrients brought by the Amazon can support the occurrence of a bloom. The plume is also characterized by high silicate (between 4 and 10 µmol kg−1 in the plume, almost 0 elsewhere), while nitrate and phosphate are rapidly consumed. Traces of inorganic phosphorus were observed in the plume, while nitrates were absent from surface waters (Fig. S2 in the Supplement). Low salinity combined with high biological productivity led to low fCO2 and a strong carbon drawdown in the plume as the ΔfCO2 reached −73µatm on 5 February (Fig. 6).

Figure 8(a) RVs Atalante and Merian ship tracks in the freshwater plume (Atalante: 2 and 5 February; Merian: 2 February) colour-coded with fCO2. The background represents the chl a on 2 February. (b) Corresponding TS diagram colour-coded with fCO2.

In an area highly influenced by the NBC waters, through rings or the retroflection, the plume stands out and modifies the biogeochemical dynamics of the region.

3.5 Shelf water and filaments

The freshwater plume is not the only water stirred by the NBC rings travelling from the NBC retroflection towards the Lesser Antilles. The shelf water is very different from the plume water and was only sampled sparsely on the way in and out of the plume (Figs. 4 and 5). On the Guiana plateau water is very rich in chl a and detrital material, rather saline (SSS∼35.5), and relatively cold (SST∼26.5C) (Fig. 5 and Table 1). Since the water sampled on the edge of the plume was cold due to a local upwelling event (and/or vertical mixing event) detailed in the Supplement (Fig. S3 in the Supplement), temperature is not homogenous on the shelf.

Further north, a filament is stirred on the western side of NBC ring A2 (Fig. 9a). It is a small-scale structure, approximately 10 km wide and easily identifiable due to its high chl a. The filament is continuously stirred by A2 and so is already visible on chl a maps on 2 February (Fig. 7). It followed A2's westward translation and was crossed on 6 and 17 February by RV Merian and on 14 February by RV Atalante (Fig. 6). It has a SSS close to 36 and an SST between 27 and 27.5 C; thus, it is slightly colder and more saline than its surrounding waters (Fig. 9b and Table 1). It stands out because of its high chl a content (chl a>0.25mg m−3) even if this is lower than close to the coast or in the freshwater plume. The strongest signal is observed on the ocean carbon parameters. In contrast to the freshwater plume, this filament presents very high fCO2 (>430µatm), highlighting different origins. It stands out from the ship track time series by also having a larger positive ΔfCO2 (50 µatm). Whereas the freshwater plume observed more southeastward carries water recently arrived on the plateau from the Amazon, the northwestward filament contains shelf waters.

3.6 Air–sea CO2 flux

In order to better characterize the impact of each structure on the regional flux, we computed air–sea CO2 flux maps from satellite data, at a resolution of 2.5 km (Fig. 10), averaged over the period of the cruise (2–19 February). The along-track flux represented in Fig. 10a and the reconstructed regional field (Fig. 10b) show the importance of the small-scale dynamical structures and highlight two strong regimes that are found on the reconstructed map. The air–sea CO2 flux in the northeastern part of the domain, characterized by the NASW, is mainly dominated by temperature effects, while further south the presence of NBC rings and their interactions with coastal waters create a strong dependence of the CO2 flux on SSS and on the biological and biogeochemical processes highlighted by the chl a.

Figure 9(a) RV Merian ship track in the shelf water filament (6 February) colour-coded with fCO2. The background represents the chl a on 6 February. (b) Corresponding TS diagram colour-coded with fCO2.

We evaluate the integrated air–sea CO2 flux over the region. In February, waters are the coldest, and the region is a strong CO2 sink of −1.7Tg C pr month (Fig. 11). Three biogeochemical domains mainly contribute to the air–sea CO2 flux: the NASW, the freshwater plume and the NBC retroflection. The impact on the flux of the small-scale coastal filament is evident along the ship tracks (Fig. 10b). However, its contribution to the total flux is weak as the signal is smoothed when averaging over 2–19 February as the filament moves following the A2 ring's northwestward translation. Each of the main three regions is identified based on its averaged SST, SSS and chl a properties in February, and the region-specific flux is determined (Fig. 11).

NASW contributes about 60 % of the total sink due to the relatively cold temperature and to strong winds that enhance the air–sea exchanges. These waters extend from Barbados northward and eastward, cover more than one third of the domain, and show a weak variability over 2–19 February.

The NBC retroflection is a source of CO2 to the atmosphere. In February, the strongest signal is observed in the southeastern part of the domain up to 8 N, 53 W. The retroflection nevertheless impacts the region as far as 10 N, 54 W as it is spatially variable, reaching up to 10 N when shedding an eddy. The two NBC rings present a small positive February air–sea CO2 flux average. Eddy A1 is almost stationary from its formation date (around 6 February) until 20 February, and its neutral to slightly positive CO2 flux is centred around 10 N, 54.5 W (Fig. 10). Eddy A2 translates rapidly westward at the beginning of February and then northward from 15 February. Its signal is therefore not as visible as on the ship tracks as it is averaged over 19 d. The retroflection is the main region with a positive air–sea CO2 flux; even if the region is too small to have a big global impact, understanding small-scale features may be significant for the total flux. The NBC rings carry part of the signal, which is heavily modified as they travel northwestward. As a result, on average in early to late February only the retroflection maintains a positive flux, while a large part of the domain dominated by modified NBC waters (not influenced by the plume) behaves as a small sink.

The freshwater plume with Amazon water is nascent when crossed by the ships (Fig. 10b) but is already the strongest signal of the time series. As the plume develops, it is entrained by NBC ring A1, then A2, and spreads out into the open ocean, as observed on SSS and chl a maps. The plume generates a strong CO2 sink that is amplified by strong winds and reaches up to 12 N (Fig. 10). The freshwater plume covers only 10 % of the total area but contributes to almost 20 % of the sink. In winter, this region either is not characterized in previous studies or is considered as dominated by high fCO2 waters brought by the NBC to the climatology. We observe here that the trapping of CO2-rich water by the NBC ring effect is relatively weak in winter, and the main signal is associated with the filaments they stir.

The northwestern tropical Atlantic therefore behaves as a sink of CO2 in early to mid-February, driven by the cold North Atlantic subtropical waters and the Amazon freshwater plume stirred by NBC rings.

Figure 10(a) Air–sea CO2 flux measured in January–February 2020 during the EUREC4A-OA/ATOMIC cruise. (b) Air–sea CO2 flux reconstructed over February 2020. (c) February climatology of the air–sea CO2 flux over 1998–2015 (Landschützer et al., 2020).

4 Discussion

4.1 Integrated air–sea CO2 flux

The northwestern tropical Atlantic presents a strong seasonal variability in air–sea CO2 fluxes (Landschützer et al., 2016). In February, waters are the coldest, and we estimate the 5–16 N, 59–50 W domain to be a strong CO2 sink of −1.7Tg C per month (Fig. 11). This region, located at tropical latitudes but combining characteristics of subtropical waters and river outflow, is difficult to represent in large-scale climatologies. Indeed, the sink for the month of February is smaller by a factor of 10 in Landschützer et al. (2020) and is also considerably smaller in Takahashi et al. (2009), but the low spatial resolution of this last product does not allow for a good quantitative comparison. This region has been rarely observed, and the inter-annual variability described in Landschützer et al.'s (2020) climatology is therefore rather uncertain. The compensating effect of different years cannot explain entirely the difference of signal observed in February 2020 with respect to the two climatologies.

Three water masses mainly contribute to the air–sea CO2 flux: the NASW, the fresh plume and the NBC retroflection. The NASW contributes about 60 % of the total sink and is not well captured in climatologies, with noticeable differences of more than 20 µatm between the measured ΔfCO2 in 2020 and the one computed from Landschützer et al. (2020) and Takahashi et al. (2009) (the closest grid point is considered for this comparison). The retroflection is the main region with a positive air–sea CO2 flux. Its influence is observed up to 10 N, 55 W, but its area is small, so its impact on the regional flux is weak. The positive flux of the retroflection is slightly overestimated in the climatologies, but it could also be due to the difficulty to detect the retroflection at the beginning of February. The main difference is that the NBC waters rich in CO2 are localized in the retroflection area and are heavily modified when spreading into the Eddy Boulevard.

The freshwater plume is a feature previously not well described for this region in winter, and we found a contribution of almost 20 % to the sink. The impact of the Amazon River has been overlooked so far in winter, but it accounts for a large part of the salinity and biogeochemical variability. Freshwater from the Amazon is not just located on the shelf, but it can spread northward, advected by the strong current variability associated with the NBC rings (Reverdin et al., 2021). These rings are the largest, faster rotating and the most energetic during boreal winter compared to other seasons (Aroucha et al., 2020). Combined with a seasonal increase in the Amazon's outflow, it induces a large variability in SSS, chl a and fCO2. The occurrence of freshwater export from the shelf to the open ocean has a strong influence on the salinity and therefore on the mixed layer depth and air–sea heat exchanges (Reverdin et al., 2021). It also strongly impacts the biogeochemistry of the region as the low fCO2 is due both to the low salinity of the plume waters and to the biological activity. The plume stirred into the open ocean by the NBC rings brings nutrients in a region strongly nutrient-limited and generates a local winter bloom. This in turn plays an important role in the air–sea CO2 flux and is a crucial feature of the southern part of the northwestern tropical Atlantic.

4.2 Extension to other years and inter-annual variability

Few tropical Atlantic measurements of biogeochemical tracers are available, in particular in the northwestern tropical Atlantic. The EUREC4A-OA/ATOMIC campaign provides the first in situ comprehensive measurements of fCO2 in this region for the boreal winter season. The reconstruction of fCO2 maps likely provides a good understanding of the spatial evolution of fCO2 and air–sea CO2 fluxes and is fitted for the months of January–February 2020. Although the processes described here are specific to winter and thus cannot be extended to other seasons, they will be useful to understand the winter variability in other years.

Only a few cruises crossed the region according to the SOCAT database between 2010 and 2019 (period with satellite SSS data), and investigating inter-annual variability is not possible. However, we can test the relation developed for 2020 for other years by using selected cruises from the SOCAT database. We thus first use the relationship to reconstruct fCO2 along the ship tracks (using in situ SSS and SST and colocalized chl a) and then over the whole region based on satellite products (OSTIA (Operational Sea Surface Temperature and Sea Ice Analysis) SST, GlobColour chl a, CCI+SSS (Climate Change Initiative Sea Surface Salinity), detailed in Appendix A). A comparison between the measured and reconstructed fCO2 for the water masses sampled by the SOCAT cruises (NASW, fresh plume, NBC retroflection, modified NBC) is presented in the appendix (Table A1). Good agreement is found between the fCO2 from the SOCAT database and the one reconstructed from the in situ temperature and salinity, as well as co-localized chl a for the four water masses (averaged difference of 5.5 µatm). When comparing reconstructed fCO2 maps with fCO2 on ship tracks (Fig. A1), the agreement between fCO2 in various water masses is very clear even though the spatial structures are sometimes a bit misplaced. This is attributable to the slightly coarser resolution of satellite products not designed specifically for each campaign, to the high spatio-temporal variability in fCO2, and to missing chl a and SST observations in cloudy areas. February 2020 was mainly cloud free, so we were able to use high-resolution daily SST and chl a. The SSS product used in 2020 is also a daily product. However, for the other years, the satellite chl a (if clouds) and SSS products have a weekly temporal resolution which smears the fast-moving structures. The gradients between water masses are therefore not always well represented, but we find a good agreement between the fCO2 of each structure, which is encouraging for future studies on inter-annual variability in winter.

Figure 11Integrated flux for the (5–16 N, 59–50 W) domain and for three water masses. For each bar duet, the one on the left in faded colours represents the integrated flux from Landschützer et al.'s (2020) February climatology, while the one on the right is computed from the reconstructed flux. Same colour code as in Fig. 5.


By identifying the main processes responsible for the variability in the air–sea CO2 flux in 2020, we can better understand the inter-annual variability in the region. Indeed, each of the main water masses has its own inter-annual variability that shapes the CO2 variability. The northern part of the domain is dominated by the variation in temperature, and therefore its inter-annual variability is mainly linked to the one of SST. From 32 years of monthly mean SST data, the SST standard deviation in the area is relatively weak and does not exceed 0.5 C. The northern sink of CO2 is therefore rather similar from year to year, coherent with the low standard deviation of the air–sea CO2 flux computed from Landschützer et al. (2020). Some variability is still observed in the snapshot of the reconstructed fCO2 (Fig. 12) but to a much smaller extent than south of Barbados. For example, the strong sink observed in March 2014 is caused by cold SST anomalies over the whole domain. Some small-scale variability in the northern part of the domain is sometimes correlated to SSS anomalies, as in 2019.

The freshwater plume sampled during EUREC4A-OA is a common feature in February. During the 2010–2019 period, events of freshwater reaching the open ocean were observed each year, and freshwater plumes similar to the one described in this paper were observed during 7 out of 10 years of satellite salinity data (Reverdin et al., 2021). Two of the main mechanisms driving the occurrence of the plume are the winds near the Amazon estuary that can induce along-shelf transport to the Guyana plateau and the presence of NBC rings. Most of the plume events similar to the one in this study suggest the presence of an anticyclone to its east. This region is commonly crossed by several NBC rings during winter (Jochumsen et al., 2010; Johns et al., 2003; Mélice and Arnault, 2017), but it also is subject to a strong year-to-year variability that has linkages with the variability in the Amazon River outflow (Aroucha et al., 2020). Therefore, identifying and understanding the processes happening in 2020 should contribute to the better assessment of the inter-annual variability in fCO2, as well as air–sea CO2 fluxes, in the northwestern tropical Atlantic during winter. Using a combination of SSS, SST and chl a brings information on the biogeochemistry of the area in winter and represents well the mesoscale structure.

Figure 12Snapshot of reconstructed fCO2 for all occurrences of fresh plumes extending at least to 10 N and east of 56 W in January–March 2010–2019 (2010, 2011 and 2013 do not present this type of event).

5 Conclusions

The EUREC4A-OA/ATOMIC campaign provides for the first time synoptic measurements related to the air–sea fluxes of CO2 in the northwestern tropical Atlantic in winter. Six main surface water masses are identified, one of them north of Barbados (North Atlantic Subtropical Water) and the other five (the NBC retroflection, modified NBC waters, the freshwater plume, the shelf water and the shelf filament) south of Barbados. The investigation highlights the two different regimes of the region. In the northern part, the variability in the CO2 flux is low, and the area is covered by relatively cold, saline and low-chlorophyll NASW. The southern part is highly variable due to the presence of large mesoscale anticyclonic eddies. In January and February 2020, two NBC rings influence the physical and biogeochemical properties of the region. The NBC retroflection is characterized by waters with equatorial origins that are relatively warm, saline and high in fCO2. As the rings separate from the retroflection, they interact with the surrounding waters, and the initial signal in fCO2 is dampened. The main impact of the rings is therefore not necessarily on the surface water they transport in their core (eddy trapping) but rather on the filaments they stir off the coast (eddy stirring). A fresh plume from the Amazon River is transported by the coastal current up to the French Guiana shelf in the beginning of February. The NBC rings entrain the plume of freshwater up to 12 N. This plume is fresh, rich in chl a, and low in fCO2, strongly contrasts with the surrounding waters, and spreads over ∼100 000km2. On the shelf not influenced by the plume, water is relatively saline and high in fCO2 and chl a probably due to a high concentration of detrital material. As ring A2 propagates westward, it continuously stirs a thin (10 km wide) filament of high fCO2 shelf water up to 12 N.

Based on the ship observations we identify distinct regimes in fCO2 linked to certain combinations of SST, SSS and chl a properties. We use this information to construct high-resolution maps of fCO2 and air–sea CO2 flux using satellite maps of SSS, SST and chl a. On average over early to mid-February, the region acts as a strong sink of CO2 (−1.7Tg C per month), the sink being 10 times smaller in air–sea CO2 flux climatologies. The NASW is responsible for most of the flux (60 %) due to low temperatures associated with winter cooling and strong winds. South of Barbados, the region also acts as a sink of CO2. The influence of equatorial water is localized in the retroflection region that acts as a small source of CO2. The main feature in this part of the domain is the fresh plume that contributes almost 20 % of the total sink.

The processes described here highlight the high variability in air–sea CO2 fluxes in winter that are quite different from the ones in summer. These features are relatively common in winter and can be used to better understand the inter-annual variability in air–sea CO2 fluxes. The northern part of the domain is driven by the variability in SST, while the southern one is a combination of the inter-annual variability in temperature, salinity and chlorophyll. It is therefore linked to the year-to-year variability in the NBC rings and the Amazon outflow.

This study is limited by the paucity of data in the region and for this time period. More fCO2 data closer to the coast would help us to better quantify the influence of shelf water on the flux. The signature of the NBC rings has been described for only two rings that had different signatures. In order to reach more robust conclusions on the transport of surface NBC water by the rings, more eddies should be observed. The variability in fCO2 occurs at large and small scales. Salinity is one of the most valuable predictors of fCO2 south of 10 N, but the satellite salinity resolution is much lower than those of temperature and chlorophyll. To have a more accurate prediction of the fCO2, a high-resolution SSS product would also be very useful.

Appendix A

Due their long time series, the following SST, chl a and SSS products are used to reconstruct fCO2 maps in winter in the northwestern tropical Atlantic for years other than 2020. Results are shown in Figs. A1 and 12. They are different than the satellite products used in the main study that were only available for a short period.

Figure A1fCO2 reconstructed from OSTIA SST, CCI+SSS and GlobColour chl a for 23 December 2015 superimposed with the fCO2 from cruise 642B20151209.

Table A1Comparison for the four main water masses between the fCO2 from SOCAT transect and the fCO2 reconstructed from in situ SSS and SST and colocalized chl a.

Download Print Version | Download XLSX

The OSTIA SST product, distributed by the Copernicus Marine Service (CMEMS) is used here. Daily maps of SST are produced at a resolution of 1/20 and are available from 1981 to present. OSTIA SST uses most SST data available for a day from both infrared and microwave inferred SST.

Surface chl a from GlobColour dataset derived from ocean colour at a 1/24 resolution is used. It is a merged product from multiple satellite missions' observations (SeaWiFS, MERIS, MODIS, VIIRS NPP, OLCI-A, VIIRS JPSS-1 and OLCI-B). GlobColour data are developed and validated by ACRI-ST and distributed by the CMEMS.

We also use SMOS and SMAP combined weekly SSS generated by the Climate Change Initiative Sea Surface Salinity (CCI+SSS) project (Boutin et al., 2021, 10.5285/5920a2c77e3c45339477acd31ce62c3c). It provides weekly level-3 SSS data from 2010 to 2019 at a spatial resolution of 50 km, sampled daily on a 25 km×25 km grid, by combining data from the SMOS, Aquarius and SMAP missions.

Code availability

Code used in this study can be made available upon reasonable request to the corresponding author.

Data availability

We benefited from numerous data sets made freely available and listed here: the ADT produced by Ssalto/Duacs distributed by CMEMS (, the chl a and SST maps produced by CLS (
url and, the SMOS L2Q field produced by CATDS (CATDS, 2019) (, the SMAP maps produced by Remote Sensing System (RSS v4 40 km), the CCI+SSS maps produced in the frame of the ESA CCI+SSS project (, and the OSTIA SST and Copernicus GlobColour chl a distributed by the CMEMS (SST_GLO_SST_L4_REP_OBSERVATIONS_010_011 and OCEANCOLOUR_GLO_CHL_L4_REP_OBSERVATIONS_009_082).

RV Atalante fCO2 data are available on the SEANOE website: RV Ron Brown and RV Merian fCO2 data can be found on the SOCAT database (expocodes 33RO20200106 and 06M220200117 respectively). The Surface Ocean CO2 Atlas (SOCAT) is an international effort, endorsed by the International Ocean Carbon Coordination Project (IOCCP), the Surface Ocean Lower Atmosphere Study (SOLAS) and the Integrated Marine Biosphere Research programme, to deliver a uniformly quality-controlled surface ocean CO2 database. The many researchers and funding agencies responsible for the collection of data and quality control are thanked for their contributions to SOCAT.


The supplement related to this article is available online at:

Author contributions

LO, JB, GR and NL conceptualized the project. LO carried out the measurements and data analysis. LO, JB, GR, NL and PL contributed to result interpretation. PL, MR and RW provided the crucial datasets. LO, MR, SS, JK, ML, TS and CN conducted field work. LO wrote the manuscript with input from all co-authors.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Biogeosciences. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.


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


This research has been supported by the European Research Council (ERC) advanced grant EUREC4A (grant agreement no. 694768) under the European Union's Horizon 2020 research and innovation programme (H2020), with additional support from CNES (the French National Centre for Space Studies) through the TOSCA SMOS-Ocean, TOEddies, and EUREC4A-OA proposals, the French national programme LEFE INSU, IFREMER, the French research fleet, the French research infrastructures AERIS and ODATIS, IPSL, the Chaire Chanel programme of the Geosciences Department at ENS, and the EUREC4A-OA JPI Ocean and Climate programme. Léa Olivier was supported by a scholarship from ENS and Sorbonne Université. We thank Jonathan Fin at the Service National d'Analyse des paramètres Océaniques du CO2 (SNAPO-CO2) at LOCEAN for the analysis of DIC and TA samples, as well as François Baurand at the US IMAGO for the nutrient analysis. Kevin Sullivan performed the data reduction and quality control of data on the Ronald H. Brown. We also warmly thank the captain and crew of RVs Atalante, Maria S. Merian and Ronald H. Brown. The measurements on the Ronald H. Brown were supported by the Global Ocean Monitoring and Observation (GOMO) programme (fund Ref. 100007298).

Financial support

This research has been supported by the Centre National d'Etudes Spatiales (grant no. 6146) and the European Research Council, H2020 European Research Council (grant no. EUREC4A (694768)), and Global Ocean Monitoring and Observation (GOMO) programme (fund ref. 100007298).

Review statement

This paper was edited by Manmohan Sarin and reviewed by Peter Land and three anonymous referees.


Aller, R. C. and Blair, N. E.: Carbon remineralization in the Amazon–Guianas tropical mobile mudbelt: A sedimentary incinerator, Cont. Shelf Res., 26, 2241–2259,, 2006. 

Andrié, C., Oudot, C., Genthon, C., and Merlivat, L.: CO2 fluxes in the tropical Atlantic during FOCAL cruises, J. Geophys. Res.-Oceans, 91, 11741–11755,, 1986. 

Aroucha, L. C., Veleda, D., Lopes, F. S., Tyaquiçã, P., Lefèvre, N., and Araujo, M.: Intra- and Inter-Annual Variability of North Brazil Current Rings Using Angular Momentum Eddy Detection and Tracking Algorithm: Observations From 1993 to 2016, J. Geophys. Res.-Oceans, 125, e2019JC015921,, 2020. 

Arrigo, K. R.: Marine manipulations, Nature, 450, 491–492,, 2007. 

Arruda, R., Atamanchuk, D., Cronin, M., Steinhoff, T., and Wallace, D. W. R.: At-sea intercomparison of three underway pCO2 systems, Limnol. Oceanogr.-Meth., 18, 63–76,, 2020. 

Bakker, D. C. E., Pfeil, B., Landa, C. S., Metzl, N., O'Brien, K. M., Olsen, A., Smith, K., Cosca, C., Harasawa, S., Jones, S. D., Nakaoka, S., Nojiri, Y., Schuster, U., Steinhoff, T., Sweeney, C., Takahashi, T., Tilbrook, B., Wada, C., Wanninkhof, R., Alin, S. R., Balestrini, C. F., Barbero, L., Bates, N. R., Bianchi, A. A., Bonou, F., Boutin, J., Bozec, Y., Burger, E. F., Cai, W.-J., Castle, R. D., Chen, L., Chierici, M., Currie, K., Evans, W., Featherstone, C., Feely, R. A., Fransson, A., Goyet, C., Greenwood, N., Gregor, L., Hankin, S., Hardman-Mountford, N. J., Harlay, J., Hauck, J., Hoppema, M., Humphreys, M. P., Hunt, C. W., Huss, B., Ibánhez, J. S. P., Johannessen, T., Keeling, R., Kitidis, V., Körtzinger, A., Kozyr, A., Krasakopoulou, E., Kuwata, A., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lo Monaco, C., Manke, A., Mathis, J. T., Merlivat, L., Millero, F. J., Monteiro, P. M. S., Munro, D. R., Murata, A., Newberger, T., Omar, A. M., Ono, T., Paterson, K., Pearce, D., Pierrot, D., Robbins, L. L., Saito, S., Salisbury, J., Schlitzer, R., Schneider, B., Schweitzer, R., Sieger, R., Skjelvan, I., Sullivan, K. F., Sutherland, S. C., Sutton, A. J., Tadokoro, K., Telszewski, M., Tuma, M., van Heuven, S. M. A. C., Vandemark, D., Ward, B., Watson, A. J., and Xu, S.: A multi-decade record of high-quality fCO2 data in version 3 of the Surface Ocean CO2 Atlas (SOCAT), Earth Syst. Sci. Data, 8, 383–413,, 2016. 

Behrenfeld, M. J., O'Malley, R. T., Siegel, D. A., McClain, C. R., Sarmiento, J. L., Feldman, G. C., Milligan, A. J., Falkowski, P. G., Letelier, R. M., and Boss, E. S.: Climate-driven trends in contemporary ocean productivity, Nature, 444, 752–755,, 2006. 

Boutin, J., Vergely, J.-L., Reul, N., Catany, R., Koehler, J., Martin, A., Rouffi, F., Arias, M., Chakroun, M., Corato, G., Estella-Perez, V., Guimbard, S., Hasson, A., Josey, S., Khvorostyanov, D., Kolodziejczyk, N., Mignot, J., Olivier, L., Reverdin, G., Stammer, D., Supply, A., Thouvenin-Masson, C., Turiel, A., Vialard, J., Cipollini, P., and Donlon, C.: ESA Sea Surface Salinity Climate Change Initiative (Sea_Surface_Salinity_cci): weekly and monthly sea surface salinity products, v03.21, for 2010 to 2020, 2021. 

CATDS: CATDS-PDC L3OS 2Q – Debiased daily valid ocean salinity values product from SMOS satellite, CATDS (CNES, IFREMER,LOCEAN, ACRI),, 2019. 

Chen, C.-T. A., Huang, T.-H., Fu, Y.-H., Bai, Y., and He, X.: Strong sources of CO2 in upper estuaries become sinks of CO2 in large river plumes, Curr. Opin. Env. Sust., 4, 179–185,, 2012. 

Coles, V. J., Brooks, M. T., Hopkins, J., Stukel, M. R., Yager, P. L., and Hood, R. R.: The pathways and properties of the Amazon River Plume in the tropical North Atlantic Ocean, J. Geophys. Res.-Oceans, 118, 6894–6913,, 2013. 

Dai, A. and Trenberth, K. E.: Estimates of Freshwater Discharge from Continents: Latitudinal and Seasonal Variations, J. Hydrometeorol., 3, 660–687,<0660:EOFDFC>2.0.CO;2, 2002. 

Edmond, J. M.: High precision determination of titration alkalinity and total carbon dioxide content of sea water by potentiometric titration, Deep-Sea Res. Oceanogr. Abstr., 17, 737–750,, 1970. 

Entekhabi, D., Njoku, E. G., O'Neill, P. E., Kellogg, K. H., Crow, W. T., Edelstein, W. N., Entin, J. K., Goodman, S. D., Jackson, T. J., and Johnson, J.: The soil moisture active passive (SMAP) mission, Proc. IEEE, 98, 704–716, 2010. 

Ffield, A.: North Brazil current rings viewed by TRMM Microwave Imager SST and the influence of the Amazon Plume, Deep-Sea Res. Pt. I, 52, 137–160,, 2005. 

Field, C. B., Behrenfeld, M. J., Randerson, J. T., and Falkowski, P.: Primary Production of the Biosphere: Integrating Terrestrial and Oceanic Components, Science, 281, 237–240,, 1998. 

Font, J., Camps, A., Borges, A., Martín-Neira, M., Boutin, J., Reul, N., Kerr, Y. H., Hahne, A., and Mecklenburg, S.: SMOS: The challenging sea surface salinity measurement from space, Proc. IEEE, 98, 649–665, 2009. 

Fournier, S., Chapron, B., Salisbury, J., Vandemark, D., and Reul, N.: Comparison of spaceborne measurements of sea surface salinity and colored detrital matter in the Amazon plume, J. Geophys. Res.-Oceans, 120, 3177–3192,, 2015. 

Fratantoni, D. M. and Glickson, D. A.: North Brazil Current Ring Generation and Evolution Observed with SeaWiFS, J. Phys. Oceanogr., 32, 1058–1074,<1058:NBCRGA>2.0.CO;2, 2002. 

Fratantoni, D. M. and Richardson, P. L.: The Evolution and Demise of North Brazil Current Rings, J. Phys. Oceanogr., 36, 1241–1264,, 2006. 

Friedlingstein, P., O'Sullivan, M., Jones, M. W., Andrew, R. M., Hauck, J., Olsen, A., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S., Aragão, L. E. O. C., Arneth, A., Arora, V., Bates, N. R., Becker, M., Benoit-Cattin, A., Bittig, H. C., Bopp, L., Bultan, S., Chandra, N., Chevallier, F., Chini, L. P., Evans, W., Florentie, L., Forster, P. M., Gasser, T., Gehlen, M., Gilfillan, D., Gkritzalis, T., Gregor, L., Gruber, N., Harris, I., Hartung, K., Haverd, V., Houghton, R. A., Ilyina, T., Jain, A. K., Joetzjer, E., Kadono, K., Kato, E., Kitidis, V., Korsbakken, J. I., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Liu, Z., Lombardozzi, D., Marland, G., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Niwa, Y., O'Brien, K., Ono, T., Palmer, P. I., Pierrot, D., Poulter, B., Resplandy, L., Robertson, E., Rödenbeck, C., Schwinger, J., Séférian, R., Skjelvan, I., Smith, A. J. P., Sutton, A. J., Tanhua, T., Tans, P. P., Tian, H., Tilbrook, B., van der Werf, G., Vuichard, N., Walker, A. P., Wanninkhof, R., Watson, A. J., Willis, D., Wiltshire, A. J., Yuan, W., Yue, X., and Zaehle, S.: Global Carbon Budget 2020, Earth Syst. Sci. Data, 12, 3269–3340,, 2020. 

Garraffo, Z. D., Johns, W. E., Chassignet, E. P., and Goni, G. J.: North Brazil Current rings and transport of southern waters in a high resolution numerical simulation of the North Atlantic, in: Elsevier Oceanography Series, vol. 68, edited by: Goni, G. J. and Malanotte-Rizzoli, P., Elsevier, 375–409,, 2003. 

Goni, G. J. and Johns, W. E.: A census of North Brazil Current Rings observed from TOPEX/POSEIDON altimetry: 1992–1998, Geophys. Res. Lett., 28, 1–4,, 2001. 

Grodsky, S. A., Vandemark, D., and Feng, H.: Assessing Coastal SMAP Surface Salinity Accuracy and Its Application to Monitoring Gulf of Maine Circulation Dynamics, J. Geophys. Res.-Oceans, 10, 1232,, 2018. 

Ibánhez, J. S. P., Araujo, M., and Lefèvre, N.: The overlooked tropical oceanic CO2 sink, J. Geophys. Res.-Oceans, 43, 3804–3812,, 2016. 

Jochumsen, K., Rhein, M., Hüttl-Kabus, S., and Böning, C. W.: On the propagation and decay of North Brazil Current rings, Proc. IEEE, 115,, 2010. 

Johns, W. E., Lee, T. N., Schott, F. A., Zantopp, R. J., and Evans, R. H.: The North Brazil Current retroflection: Seasonal structure and eddy variability, Geophys. Res. Lett., 95, 22103–22120,, 1990. 

Johns, W. E., Zantopp, R. J., and Goni, G. J.: Cross-gyre transport by North Brazil Current rings, in: Elsevier Oceanography Series, vol. 68, edited by: Goni, G. J. and Malanotte-Rizzoli, P., Elsevier, 411–441,, 2003. 

Karstensen, J., Lavik, G., Acquistapace, C., Bagheri, G., Begler, C., Bendinger, A., Bodenschatz, E., Böck, T., Güttler, J., Hall, K., Körner, M., Kopp, A., Lange, D., Mehlmann, M., Nordsiek, F., Reus, K., Ribbe, J., Philippi, M., Piosek, S., Ritschel, M., Tschitschko, B., and Wiskandt, J.: EUREC4A Campaign, Cruise No. MSM89, 17 January–20 February 2020, Bridgetown (Barbados) – Bridgetown (Barbados), The ocean mesoscale component in the EUREC4A field study, Gutachterpanel Forschungsschiffe, Bonn, 70 pp.,, 2020. 

Kerr, Y. H., Waldteufel, P., Wigneron, J.-P., Delwart, S., Cabot, F., Boutin, J., Escorihuela, M.-J., Font, J., Reul, N., and Gruhier, C.: The SMOS mission: New tool for monitoring key elements ofthe global water cycle, Global Biogeochem. Cy., 98, 666–687, 2010. 

Körtzinger, A.: A significant CO2 sink in the tropical Atlantic Ocean associated with the Amazon River plume, Global Biogeochem. Cy., Geophys. Res. Lett., 30,, 2003. 

Landschützer, P., Gruber, N., Bakker, D. C. E., and Schuster, U.: Recent variability of the global ocean carbon sink, Global Biogeochem. Cy., 28, 927–949,, 2014. 

Landschützer, P., Gruber, N., and Bakker, D. C. E.: Decadal variations and trends of the global ocean carbon sink, Global Biogeochem. Cy., 30, 1396–1417,, 2016. 

Landschützer, P., Laruelle, G. G., Roobaert, A., and Regnier, P.: A uniform pCO2 climatology combining open and coastal oceans, Earth Syst. Sci. Data, 12, 2537–2553,, 2020. 

Laxenaire, R., Speich, S., Blanke, B., Chaigneau, A., Pegliasco, C., and Stegner, A.: Anticyclonic Eddies Connecting the Western Boundaries of Indian and Atlantic Oceans, J. Geophys. Res.-Oceans, 123, 7651–7677,, 2018. 

Lefèvre, N., Diverrés, D., and Gallois, F.: Origin of CO2 undersaturation in the western tropical Atlantic, Tellus B, 62, 595–607,, 2010. 

Longhurst, A. R.: Ecological geography of the sea, Elsevier, 2010. 

Longhurst, A., Sathyendranath, S., Platt, T., and Caverhill, C.: An estimate of global primary production in the ocean from satellite radiometer data, J. Plankton Res., 17, 1245–1271, 1995. 

Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, Ö., Yu, R., and Zhou, B. (Eds.): Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, 2021. 

Medeiros, P. M., Seidel, M., Ward, N. D., Carpenter, E. J., Gomes, H. R., Niggemann, J., Krusche, A. V., Richey, J. E., Yager, P. L., and Dittmar, T.: Fate of the Amazon River dissolved organic matter in the tropical Atlantic Ocean, Global Biogeochem. Cy., 29, 677–690,, 2015. 

Mélice, J.-L. and Arnault, S.: Investigation of the Intra-Annual Variability of the North Equatorial Countercurrent/North Brazil Current Eddies and of the Instability Waves of the North Tropical Atlantic Ocean Using Satellite Altimetry and Empirical Mode Decomposition, J. Atmos. Ocean. Tech., 34, 2295–2310,, 2017. 

Mu, L., Gomes, H. D. R., Burns, S. M., Goes, J. I., Coles, V. J., Rezende, C. E., Thompson, F. L., Moura, R. L., Page, B., and Yager, P. L.: Temporal Variability of Air-Sea CO2 flux in the Western Tropical North Atlantic Influenced by the Amazon River Plume, Global Biogeochem. Cy., 35, e2020GB006798,, 2021. 

Muller-Karger, F. E., McClain, C. R., and Richardson, P. L.: The dispersal of the Amazon's water, Nature, 333, 56–59,, 1988. 

Olivier, L., Labaste, M., Noisel, C., and Lefèvre, N.: Underway fCO2 distribution during the EUREC4A-OA experiment,, 2020. 

Piepmeier, J. R., Focardi, P., Horgan, K. A., Knuble, J., Ehsan, N., Lucey, J., Brambora, C., Brown, P. R., Hoffman, P. J., and French, R. T.: SMAP L-band microwave radiometer: Instrument design and first year on orbit, IEEE T. Geosci. Remote Sens., 55, 1954–1966, 2017. 

Pierrot, D., Neill, C., Sullivan, K., Castle, R., Wanninkhof, R., Lüger, H., Johannessen, T., Olsen, A., Feely, R. A., and Cosca, C. E.: Recommendations for autonomous underway pCO2 measuring systems and data-reduction routines, Deep-Sea Res. Pt. II, 56, 512–522,, 2009. 

Poisson, A., Metzl, N., Brunet, C., Schauer, B., Bres, B., Ruiz-Pino, D., and Louanchi, F.: Variability of sources and sinks of CO2 in the western Indian and southern oceans during the year 1991, J. Geophys. Res.-Oceans, 98, 22759–22778,, 1993. 

Pörtner, H.-O., Roberts, D. C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., and Weyer, N. M.: The ocean and cryosphere in a changing climate, 2019. 

Quinn, P. K., Thompson, E. J., Coffman, D. J., Baidar, S., Bariteau, L., Bates, T. S., Bigorre, S., Brewer, A., de Boer, G., de Szoeke, S. P., Drushka, K., Foltz, G. R., Intrieri, J., Iyer, S., Fairall, C. W., Gaston, C. J., Jansen, F., Johnson, J. E., Krüger, O. O., Marchbanks, R. D., Moran, K. P., Noone, D., Pezoa, S., Pincus, R., Plueddemann, A. J., Pöhlker, M. L., Pöschl, U., Quinones Melendez, E., Royer, H. M., Szczodrak, M., Thomson, J., Upchurch, L. M., Zhang, C., Zhang, D., and Zuidema, P.: Measurements from the RV Ronald H. Brown and related platforms as part of the Atlantic Tradewind Ocean-Atmosphere Mesoscale Interaction Campaign (ATOMIC), Earth Syst. Sci. Data, 13, 1759–1790,, 2021. 

Reverdin, G., Olivier, L., Foltz, G. R., Speich, S., Karstensen, J., Horstmann, J., Zhang, D., Laxenaire, R., Carton, X., Branger, H., Carrasco, R., and Boutin, J.: Formation and Evolution of a Freshwater Plume in the Northwestern Tropical Atlantic in February 2020, J. Geophys. Res.-Oceans, 126, e2020JC016981,, 2021. 

Richardson, P. L., Hufford, G. E., Limeburner, R., and Brown, W. S.: North Brazil Current retroflection eddies, J. Geophys. Res.-Oceans, 99, 5081–5093,, 1994. 

Schott, F. A., Fischer, J., and Stramma, L.: Transports and Pathways of the Upper-Layer Circulation in the Western Tropical Atlantic, J. Phys. Oceanogr., 28, 1904–1928,<1904:TAPOTU>2.0.CO;2, 1998. 

Speich, S. and The Embarked Science Team: EUREC4A-OA. Cruise Report. 19 January–19 February 2020. Vessel: L'ATALANTE,, 2021. 

Stevens, B., Bony, S., Farrell, D., Ament, F., Blyth, A., Fairall, C., Karstensen, J., Quinn, P. K., Speich, S., Acquistapace, C., Aemisegger, F., Albright, A. L., Bellenger, H., Bodenschatz, E., Caesar, K.-A., Chewitt-Lucas, R., de Boer, G., Delanoë, J., Denby, L., Ewald, F., Fildier, B., Forde, M., George, G., Gross, S., Hagen, M., Hausold, A., Heywood, K. J., Hirsch, L., Jacob, M., Jansen, F., Kinne, S., Klocke, D., Kölling, T., Konow, H., Lothon, M., Mohr, W., Naumann, A. K., Nuijens, L., Olivier, L., Pincus, R., Pöhlker, M., Reverdin, G., Roberts, G., Schnitt, S., Schulz, H., Siebesma, A. P., Stephan, C. C., Sullivan, P., Touzé-Peiffer, L., Vial, J., Vogel, R., Zuidema, P., Alexander, N., Alves, L., Arixi, S., Asmath, H., Bagheri, G., Baier, K., Bailey, A., Baranowski, D., Baron, A., Barrau, S., Barrett, P. A., Batier, F., Behrendt, A., Bendinger, A., Beucher, F., Bigorre, S., Blades, E., Blossey, P., Bock, O., Böing, S., Bosser, P., Bourras, D., Bouruet-Aubertot, P., Bower, K., Branellec, P., Branger, H., Brennek, M., Brewer, A., Brilouet , P.-E., Brügmann, B., Buehler, S. A., Burke, E., Burton, R., Calmer, R., Canonici, J.-C., Carton, X., Cato Jr., G., Charles, J. A., Chazette, P., Chen, Y., Chilinski, M. T., Choularton, T., Chuang, P., Clarke, S., Coe, H., Cornet, C., Coutris, P., Couvreux, F., Crewell, S., Cronin, T., Cui, Z., Cuypers, Y., Daley, A., Damerell, G. M., Dauhut, T., Deneke, H., Desbios, J.-P., Dörner, S., Donner, S., Douet, V., Drushka, K., Dütsch, M., Ehrlich, A., Emanuel, K., Emmanouilidis, A., Etienne, J.-C., Etienne-Leblanc, S., Faure, G., Feingold, G., Ferrero, L., Fix, A., Flamant, C., Flatau, P. J., Foltz, G. R., Forster, L., Furtuna, I., Gadian, A., Galewsky, J., Gallagher, M., Gallimore, P., Gaston, C., Gentemann, C., Geyskens, N., Giez, A., Gollop, J., Gouirand, I., Gourbeyre, C., de Graaf, D., de Groot, G. E., Grosz, R., Güttler, J., Gutleben, M., Hall, K., Harris, G., Helfer, K. C., Henze, D., Herbert, C., Holanda, B., Ibanez-Landeta, A., Intrieri, J., Iyer, S., Julien, F., Kalesse, H., Kazil, J., Kellman, A., Kidane, A. T., Kirchner, U., Klingebiel, M., Körner, M., Kremper, L. A., Kretzschmar, J., Krüger, O., Kumala, W., Kurz, A., L'Hégaret, P., Labaste, M., Lachlan-Cope, T., Laing, A., Landschützer, P., Lang, T., Lange, D., Lange, I., Laplace, C., Lavik, G., Laxenaire, R., Le Bihan, C., Leandro, M., Lefèvre, N., Lena, M., Lenschow, D., Li, Q., Lloyd, G., Los, S., Losi, N., Lovell, O., Luneau, C., Makuch, P., Malinowski, S., Manta, G., Marinou, E., Marsden, N., Masson, S., Maury, N., Mayer, B., Mayers-Als, M., Mazel, C., McGeary, W., McWilliams, J. C., Mech, M., Mehlmann, M., Meroni, A. N., Mieslinger, T., Minikin, A., Minnett, P., Möller, G., Morfa Avalos, Y., Muller, C., Musat, I., Napoli, A., Neuberger, A., Noisel, C., Noone, D., Nordsiek, F., Nowak, J. L., Oswald, L., Parker, D. J., Peck, C., Person, R., Philippi, M., Plueddemann, A., Pöhlker, C., Pörtge, V., Pöschl, U., Pologne, L., Posyniak, M., Prange, M., Quiñones Meléndez, E., Radtke, J., Ramage, K., Reimann, J., Renault, L., Reus, K., Reyes, A., Ribbe, J., Ringel, M., Ritschel, M., Rocha, C. B., Rochetin, N., Röttenbacher, J., Rollo, C., Royer, H., Sadoulet, P., Saffin, L., Sandiford, S., Sandu, I., Schäfer, M., Schemann, V., Schirmacher, I., Schlenczek, O., Schmidt, J., Schröder, M., Schwarzenboeck, A., Sealy, A., Senff, C. J., Serikov, I., Shohan, S., Siddle, E., Smirnov, A., Späth, F., Spooner, B., Stolla, M. K., Szkółka, W., de Szoeke, S. P., Tarot, S., Tetoni, E., Thompson, E., Thomson, J., Tomassini, L., Totems, J., Ubele, A. A., Villiger, L., von Arx, J., Wagner, T., Walther, A., Webber, B., Wendisch, M., Whitehall, S., Wiltshire, A., Wing, A. A., Wirth, M., Wiskandt, J., Wolf, K., Worbes, L., Wright, E., Wulfmeyer, V., Young, S., Zhang, C., Zhang, D., Ziemen, F., Zinner, T., and Zöger, M.: EUREC4A, Earth Syst. Sci. Data, 13, 4067–4119,, 2021. 

Stum, J., Tebri, H., Lehodey, P., Senina, I., Greiner, E., and Lucas, M.: NRT operational chlorophyll maps calculation for marine applications, 1, n.d.,, 2015. 

Subramaniam, A., Yager, P. L., Carpenter, E. J., Mahaffey, C., Björkman, K., Cooley, S., Kustka, A. B., Montoya, J. P., Sañudo-Wilhelmy, S. A., Shipe, R., and Capone, D. G.: Amazon River enhances diazotrophy and carbon sequestration in the tropical North Atlantic Ocean, P. Natl. Acad. Sci. USA, 105, 10460–10465,, 2008. 

Takahashi, T., Olafsson, J., Goddard, J. G., Chipman, D. W., and Sutherland, S. C.: Seasonal variation of CO2 and nutrients in the high-latitude surface oceans: A comparative study, Global Biogeochem. Cy., 7, 843–878,, 1993. 

Takahashi, T., Sutherland, S. C., Sweeney, C., Poisson, A., Metzl, N., Tilbrook, B., Bates, N., Wanninkhof, R., Feely, R. A., Sabine, C., Olafsson, J., and Nojiri, Y.: Global sea–air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects, Deep-Sea Res. Pt. II, 49, 1601–1622,, 2002.  

Takahashi, T., Sutherland, S. C., Wanninkhof, R., Sweeney, C., Feely, R. A., Chipman, D. W., Hales, B., Friederich, G., Chavez, F., Sabine, C., Watson, A., Bakker, D. C. E., Schuster, U., Metzl, N., Yoshikawa-Inoue, H., Ishii, M., Midorikawa, T., Nojiri, Y., Körtzinger, A., Steinhoff, T., Hoppema, M., Olafsson, J., Arnarson, T. S., Tilbrook, B., Johannessen, T., Olsen, A., Bellerby, R., Wong, C. S., Delille, B., Bates, N. R., and de Baar, H. J. W.: Climatological mean and decadal change in surface ocean pCO2, and net sea–air CO2 flux over the global oceans, Deep-Sea Res. Pt. II, 56, 554–577,, 2009. 

Tennekes, H.: The Logarithmic Wind Profile, J. Atmos. Sci., 30, 234–238,<0234:TLWP>2.0.CO;2, 1973. 

Wanninkhof, R.: Relationship between wind speed and gas exchange over the ocean revisited, Limnol. Oceanogr.-Meth., 12, 351–362,, 2014. 

Wanninkhof, R., Lewis, E., Feely, R. A., and Millero, F. J.: The optimal carbonate dissociation constants for determining surface water pCO2 from alkalinity and total inorganic carbon, Mar. Chem., 65, 291–301,, 1999. 

Weiss, R. F.: Carbon dioxide in water and seawater: the solubility of a non-ideal gas, Mar. Chem., 2, 203–215,, 1974. 

Wilson, W. D., Johns, W. E., and Garzoli, S. L.: Velocity structure of North Brazil Current rings, Geophys. Res. Lett., 29, 114-1–114-4,, 2002. 

Short summary
We investigate the impact of the interactions between eddies and the Amazon River plume on the CO2 air–sea fluxes to better characterize the ocean carbon sink in winter 2020. The region is a strong CO2 sink, previously underestimated by a factor of 10 due to a lack of data and understanding of the processes responsible for the variability in ocean carbon parameters. The CO2 absorption is mainly driven by freshwater from the Amazon entrained by eddies and by the winter seasonal cooling.
Final-revised paper