Articles | Volume 17, issue 9
Research article
14 May 2020
Research article |  | 14 May 2020

Ocean carbonate system variability in the North Atlantic Subpolar surface water (1993–2017)

Coraline Leseurre, Claire Lo Monaco, Gilles Reverdin, Nicolas Metzl, Jonathan Fin, Solveig Olafsdottir, and Virginie Racapé

The North Atlantic is one of the major ocean sinks for natural and anthropogenic atmospheric CO2. Given the variability of the circulation, convective processes or warming–cooling recognized in the high latitudes in this region, a better understanding of the CO2 sink temporal variability and associated acidification needs a close inspection of seasonal, interannual to multidecadal observations. In this study, we investigate the evolution of CO2 uptake and ocean acidification in the North Atlantic Subpolar Gyre (50–64 N) using repeated observations collected over the last 3 decades in the framework of the long-term monitoring program SURATLANT (SURveillance de l'ATLANTique). Over the full period (1993–2017) pH decreases (−0.0017 yr−1) and fugacity of CO2 (fCO2) increases (+1.70 µatm yr−1). The trend of fCO2 in surface water is slightly less than the atmospheric rate (+1.96 µatm yr−1). This is mainly due to dissolved inorganic carbon (DIC) increase associated with the anthropogenic signal. However, over shorter periods (4–10 years) and depending on the season, we detect significant variability investigated in more detail in this study. Data obtained between 1993 and 1997 suggest a rapid increase in fCO2 in summer (up to +14µatm yr−1) that was driven by a significant warming and an increase in DIC for a short period. Similar fCO2 trends are observed between 2001 and 2007 during both summer and winter, but, without significant warming detected, these trends are mainly explained by an increase in DIC and a decrease in alkalinity. This also leads to a pH decrease but with contrasting trends depending on the region and season (between −0.006 and −0.013 yr−1). Conversely, data obtained during the last decade (2008–2017) in summer show a cooling of surface waters and an increase in alkalinity, leading to a strong decrease in surface fCO2 (between −4.4 and −2.3µatm yr−1; i.e., the ocean CO2 sink increases). Surprisingly, during summer, pH increases up to +0.0052 yr−1 in the southern subpolar gyre. Overall, our results show that, in addition to the accumulation of anthropogenic CO2, the temporal changes in the uptake of CO2 and ocean acidification in the North Atlantic Subpolar Gyre present significant multiannual variability, not clearly directly associated with the North Atlantic Oscillation (NAO). With such variability it is uncertain to predict the near-future evolution of air–sea CO2 fluxes and pH in this region. Thus, it is highly recommended to maintain long-term observations to monitor these properties in the next decade.

1 Introduction

The ocean plays an important role in climate regulation by absorbing between one-quarter and one-third of anthropogenic carbon dioxide (CO2) emitted to the atmosphere (Le Quéré et al., 2018; Gruber et al., 2019). The North Atlantic (NA) is one of the strongest ocean sinks for natural (Takahashi et al., 2009) and anthropogenic atmospheric CO2 (Khatiwala et al., 2013; Sabine et al., 2004; Perez et al., 2008). Although many observational studies were conducted in the NA to evaluate how the CO2 sink varies from seasonal to multi-decadal scales and modeling studies attempted to reproduce the observed changes, there are still open questions regarding the processes that control the contemporary CO2 sink variability in this region (e.g., Schuster et al., 2013, for a synthesis). The uptake of CO2 in the NA is mainly due to extensive biological activity during spring–summer and considerable heat loss during winter, both processes being subject to significant spatiotemporal variability related to climate mode such as the North Atlantic Oscillation (NAO) and the Atlantic Multidecadal Variability (AMV), also called the Atlantic Multidecadal Oscillation (AMO). The recent decades have indeed witnessed large variations in the NAO index (Hurrell index, 2013), which displayed strong positive values in the early 1990s, tended to be positive from 1996 to the early 2000, and was then more neutral, albeit with extreme anomalies, such as the negative 2010 event, until again witnessing large positive anomalies in the mid-2010s. The NA has also been at the core of large decadal to multi-decadal variability of the AMV with a positive peak by the mid-2000s and a decrease since then (Robson et al., 2018). In the North Atlantic Subpolar Gyre (NASPG, around 50–60 N) where the oceanic CO2 sink is particularly strong (e.g., Watson et al., 2009), recent years in the mid-2010s have also witnessed particularly large negative surface temperature anomalies (Josey et al., 2018) and an extreme freshening (Holliday et al., 2020). Associated changes in the ocean circulation and water masses in the NASPG are for example presented in Chafik et al. (2014, 2019), Desbruyères et al. (2015) and Nigam et al. (2018). These physical processes linked to the NAO and/or AMV directly impact the sea surface fCO2 and air–sea CO2 fluxes through warming–cooling or deeper convection as this has been observed in the NASPG for specific periods (Corbière et al., 2007; Fröb et al., 2019). However, a direct link of the CO2 uptake variability with the NAO depends on the period investigated, and this is not always clearly revealed from observations (Metzl et al., 2010; Schuster et al., 2013).

Ocean models and Earth system models (ESMs) can help to understand the link between the NAO and biogeochemical cycles, but results from models in the NA are still controversial (Thomas et al., 2008; Ullman et al., 2009; Keller et al., 2012; Tjiputra et al., 2012). Keller et al. (2012) investigated simulations from six ESMs and found that the interannual variability of the CO2 sink in the North Atlantic is the largest in the subpolar gyre. They also conclude that for winter (i.e., not the productive season) on-site entrainment in the NASPG (mixing and upwelling) is the main driver of carbon sink variability as opposed to advection (Thomas et al., 2008), but the magnitude and responses of the carbon uptake to the NAO significantly differ between the ESMs (Keller et al., 2012). Thomas et al. (2008) suggest that negative or neutral NAO conditions result in a substantial decline of the CO2 uptake for the years 1997–2004 along the North Atlantic Current (NAC) region and in the eastern subpolar gyre. Conversely, Ullman et al. (2009) conclude that the CO2 uptake increased over 1992–2006. During the transition of the NAO (from positive to neutral), their model simulates a decline of the convection and vertical DIC supply to the surface in the subpolar region, counteracting the increase in pCO2 due to warming. This leads to a relatively small pCO2 increase compared with the atmospheric CO2 trend and an increasing CO2 sink, while observations in 1993–2003 suggest a reduced CO2 sink in the NASPG after the NAO shift in the mid-1990s (Corbière et al., 2007).

For longer timescales, 2 decades or more, observations in the NASPG show a gradual fCO2 increase slightly smaller or indistinguishable from trends in atmospheric CO2 and that the long-term oceanic fCO2 trend is not correlated with the NAO (Takahashi et al., 2009; McKinley et al., 2011; Fay and McKinley, 2013). On a decadal scale the link between the AMO and fCO2 variability in the NASPG appears more robust (Breeden and McKinley, 2016; Landschützer et al., 2019). When the AMO enters a positive phase, these studies indicate a reduction of DIC in the NASPG due to reduced mixing, a process that dominates the effect of warming on fCO2 (Breeden and McKinley, 2016). Consequently, since the mid-1990s the fCO2 long-term trend lagged the one in the atmosphere, and the CO2 uptake is increasing in this region. Such a signal is not captured by current ESM CMIP5 models (Tjiputra et al., 2014; Lebehot et al., 2019) likely due to inadequate representation of biogeochemical cycles (here DIC and/or TA) and in part to bad representation of mixing leading to pCO2 seasonality that is opposed to observations (Goris et al., 2018). This leads to uncertainties on the evolution of the oceanic CO2 uptake in the North Atlantic in the future (Lebehot et al., 2019). A better knowledge of DIC and TA trends (not only fCO2) is needed to correct and validate biogeochemical representation in 35 ESMs.

The responses of the marine biological processes to the NAO and AMO are much less studied than their physical counterpart. However, and somehow, as for fCO2, the biology link seems to also be better identified with the AMO than with the NAO. Based on the biomass variability detected from SeaWiFS over 1998–2007 in the North Atlantic and supported by biogeochemical simulations, McKinley et al. (2018) conclude that “nowhere is the NAO correlated with biomass variability”. The same is true, i.e., no correlation with the NAO, when analyzing long-term in situ observations of marine species in the NA since the 1960s (Beaugrand et al., 2013; Rivero-Calle et al., 2015). However, these authors detect a significant increase in calcifying species (coccolithophores, foraminifera) associated with the warming in the NA since the mid-1990s and correlated with the AMO. This increase in calcification would impact the DIC and TA concentrations in summer and thus the fCO2 trend and carbon uptake in the NA. On the other hand, this also suggests that these phytoplanktonic groups are not yet altered by acidification, a phenomenon known as the other CO2 problem (Doney et al., 2009). Indeed, the accumulation of anthropogenic CO2 in the ocean has led to a decrease in pH in surface waters by 0.1 unit since the industrial revolution (IPCC, 2013 – Hartmann et al., 2013). It is projected that surface pH would be reduced by 0.2 to 0.4 units by the end of the century depending on the anthropogenic emission scenario (Orr et al., 2005; Doney et al., 2014; Jiang et al., 2019). This would be associated with a dramatic reduction of the aragonite saturation horizon at high latitudes, from 2000 m at present up to 100 m by 2100 in the North Atlantic (north of 50 N). Indeed, these northern surface waters (north of 50 N) have initially low CaCO3 saturation because of high natural (preindustrial) DIC concentrations due to low temperatures giving high CO2 solubility. Such fast lysocline shoaling based on a simulation (Orr et al., 2005) is starting to be observed in this region (e.g., Vázquez-Rodríguez et al., 2012). Due to the potential threat to marine life, the decrease in pH is now recognized as a true indicator of global change, similarly to warming and sea level rise (World Meteorological Organization, 2018). The urgent need to document and understand changes in oceanic pH and its impact on marine life has motivated more studies in recent years. Among them, the data syntheses from Lauvset et al. (2015) and Bates et al. (2014) indicated that pH is decreasing in most of the surface ocean as a result of the increase in oceanic CO2. In the North Atlantic they reported a decrease in pH ranging between −0.001 yr−1 in the subtropical region and −0.0026 yr−1 in the Irminger Sea. These results have been confirmed at a regional scale. Based on 12 cruises conducted in summer between 1991 and 2015 in the Irminger and Iceland basins, García-Ibáñez et al. (2016) identified a pH decrease in all water masses, with the largest change in surface and intermediate waters of between −0.0016 and −0.0018yr−1. They conclude that it is mainly due to anthropogenic CO2 and modulated by change in alkalinity. Recently Omar et al. (2019) analyzed regular surface observations during 2004–2015 in the North Sea; when selecting winter data in the North Atlantic surface waters they found significant long-term trends of fCO2 and pH of +2.4µatm yr−1 and −0.0024 yr−1, consistent with results obtained in the subpolar North Atlantic (Bates et al., 2014; Lauvset and Gruber, 2014). A driver analysis indicates that fCO2 and pH trends are almost entirely explained by increasing DIC due to anthropogenic CO2 uptake (Omar et al., 2019).

As anthropogenic CO2 accumulates each year in the surface and interior ocean with a direct impact on both air–sea CO2 exchange and pH, it is important to conduct and maintain regular observations to follow the evolution of the ocean carbonate system properties and interpret how they vary from interannual to multi-decadal scales in relation with physical and biological changes. The aim of this study is to document the recent evolution of the surface CO2 system parameters in the NASPG and to evaluate the drivers for the evolution of surface fCO2 and pH for winter and summer. We mostly use total alkalinity (TA) and dissolved inorganic carbon (DIC) discrete observations, while other studies are based on observations of partial pressure of CO2 (pCO2) and TA–salinity relations (Lauvset and Gruber, 2014; Lauvset et al., 2015). The observations were collected in 1993–2017 between Iceland and Newfoundland, and this extends previous analyses based on winter SURATLANT observations in 1993–2008 to both summer and winter (Corbière et al., 2007; Metzl et al., 2010).

Figure 1SURATLANT cruise tracks (in black, TA–DIC data over the period 1993–2017) with the five boxes (box A: 46–50 N (in yellow), box B: 50–54 N (in red), box C: 54–58 N (in green), box D: 58–62 N (in blue), box E: 62–64 N (in black)) and surface currents in the North Atlantic Subpolar Gyre (NASPG) circulation resulting from the Gulf Stream and the North Atlantic Current (in red).

The NASPG as defined here extends from the east coast of Canada (65 W) to the west of northern Europe (5 W) mostly north of 45 N (Fig. 1). It consists mostly of the Labrador, Irminger and Iceland basins, and the Rockall Plateau and Trough. The NASPG is a region of mixing between subtropical, subpolar and polar surface waters. The inflow of these surface waters is dominated by three surface currents: the North Atlantic Current (NAC), also known in this sector as the North Atlantic Drift, mostly of subtropical origin; the East Greenland Current (EGC); and the Labrador Current (LC), which are the conduits that bring waters from further north and the Arctic (Fig. 1). The North Atlantic Drift delimits the NASPG in the south and transports warm (8 C to 15 C) and salty (35 to 36) subtropical surface waters to the north. These waters are then cooled as they circulate in the NASPG due to a loss of heat to the atmosphere, while they become less salty due to excess local precipitation and freshwater inputs from the Arctic and the ice sheet. The East Greenland Current (EGC) and Labrador Current (LC) are the western branches of the NASPG and carry cold (<4C) and fresh (<34.6) waters of polar origin. The study area is mostly in the central part of the NASPG between Iceland and Newfoundland.

2 Materials and method

2.1 Data collection and measurements

This study is based on observations collected in the framework of the long-term monitoring program SURATLANT (SURveillance de l'ATLANTique) initiated in 1993. The main objective of this program is to monitor hydrological and biogeochemical properties in surface waters of the NASPG, notably sea surface salinity (SSS), to improve the understanding of the role of salinity in the variability and predictability of climate as well as the water cycle (Reverdin et al., 2018a; Holliday et al., 2020). To this end, two to four cruises per year are conducted between Reykjavík (Iceland) and Newfoundland (Canada) on board merchant ships. Seawater samples are collected every 3 to 4 h from a pumping system (at approximately 5 m deep) in order to measure salinity, as well as alkalinity (measured since 2001), DIC, silicate, nitrate and phosphate concentrations.

Underway measurements of sea surface temperature (SST) are obtained using a Sea-Bird thermosalinograph, with a precision of ±0.1C. SSS values in this paper are from sample measurements of conductivity using a Guildline AUTOSAL salinometer since 1997 at the Marine and Freshwater Research Institute (MFRI) in Reykjavík (the associated error is estimated at ±0.005). For TA and DIC, the 500 mL glass bottles are rinsed three times before introducing seawater (avoiding introducing air bubbles) with an overflow in order to remove the water in contact with air during filling. The samples are then poisoned with mercuric chloride and stored in a cool, dark place. Since 2001, these samples are measured within 3 months of collection at the SNAPO-CO2 (Service National d'Analyse des Paramètres Océaniques du CO2) located at the laboratory LOCEAN in Paris. The accuracy is on the order of ±3µmol kg−1 for both DIC and TA based on Certified Reference Material (CRM) analysis and occasional intercomparisons as explained by Reverdin et al. (2018b). Samples for the analysis of nutrient concentrations are frozen just after sampling. Spring and summer samples are filtered before analysis. The measurements are carried out by the MFRI team in Reykjavík according to the standard colorimetric method described by Ólafsson et al. (2010) with an accuracy of ±0.2µmol L−1 for nitrate and silicate and ±0.03µmol L−1 for phosphate. Details of the sampling and analyses for all properties are provided in Reverdin et al. (2018b). The SURATLANT dataset is freely available and is accessible at (last access: 11 March 2019). In the SURATLANT data file, each sample and data are associated with a quality flag (based on WOCE flag criteria, i.e., the same as used in GLODAP database; Olsen et al., 2016). In this study we use only data qualified with flag 2 (acceptable data).

2.2 Calculations of the carbonate system parameters and contributions

Carbonate system parameters such as pH, the fugacity of CO2 (fCO2) and the saturation states for the calcium carbonate minerals calcite (ΩCa) and aragonite (ΩAr) are calculated using SST, SSS, TA, DIC, silicate and phosphate data described by Reverdin et al. (2018b). For the last two, we used monthly climatological values derived from the SURATLANT data over the period 1993–2017 because nutrient interannual variability does not have a significant impact on the carbonate system parameter calculations (and because nutrients were missing for some cruises). The calculation program used is CO2SYS originally developed by Lewis et al. (1998) available in MATLAB (van Heuven et al., 2011) that now includes error propagation (Orr et al., 2018). The constants of the thermodynamic equilibrium of CO2 in seawater used are K1 (for the dissociation of carbonic acid) and K2 (for the bicarbonate ion) defined by Mehrbach et al. (1973), refitted by Dickson and Millero (1987). The total boron value is calculated according to Uppström (1974), and the KHSO4 dissociation constant is from Dickson (1990). The adopted pH scale is total scale.

When TA was not measured (notably before 2001), it was calculated from salinity data. The correlation between sea surface alkalinity and salinity in the open ocean can be described with an empirical linear relationship (Millero et al., 1998; Friis et al., 2003). For the NASPG, the original TA–SSS relation used to analyze fCO2 variability and trend in 1993–2007 (e.g., Corbière et al., 2007; Metzl et al., 2010; McKinley et al., 2011) was revisited by Reverdin et al. (2018b) based on SURATLANT data over the period 2001–2016 (Eq. 1). This equation was obtained for samples collected north of 50 N with SSS >34.

(1) TA = 45.5337 × SSS + 713.58

To evaluate the air–sea CO2 difference (delta-fCO2) for each sample, we used atmospheric fCO2 values for the period 1993–2017 calculated from the molar fraction (xCO2) data at Mace Head provided by the Cooperative Global Atmospheric Data Integration Project (Dlugokencky et al., 2019). The xCO2 data are available at (last access: 6 March 2019). xCO2 data were converted to fCO2 at 100 % humidity following Weiss and Price (1980), with standard atmospheric pressure (i.e., 1013.25 hPa).

The trends in sea surface pH, fCO2 and Ω are driven by changes in SST, SSS, TA and DIC (Keeling et al., 2004; Fröb et al., 2019). For DIC, this includes anthropogenic CO2 concentration (Cant) and the DIC natural component (DIC-nat = DIC-Cant). The contribution of each term for a specific period and season is evaluated by allowing a change in only one parameter according to their observed trend, while setting the other parameters to their climatological seasonal values (Eq. 2). The uncertainty on the contributions was evaluated by performing 1000 random perturbations within the range of the standard deviation of the observed trends in SST, SSS, TA and DIC. The same method was applied when separating the effect of DIC-nat and anthropogenic CO2 (Cant).

(2) d X d t = X SST dSST d t + X SSS dSSS d t + X DIC sDIC SSS 0 dSSS d t + SSS SSS 0 dsDIC d t + X TA sTA SSS 0 dSSS d t + SSS SSS 0 dsTA d t

Here, X corresponds to pH, fCO2 and Ω. SSS0 is the reference salinity, which is set to 35 (Normal Standard Seawater; Millero et al., 2008) and very close to the mean salinity observed in the NASPG (SSS = 34.84). sDIC and sTA were computed by normalizing DIC and TA to a salinity of 35 and assuming a nonzero freshwater end-member DIC and TA concentration (Friis et al., 2003; Reverdin et al., 2018b). The trends (dX∕dt) were estimated for each season, each box and each period, and SST,SSSTA and DIC correspond to their climatological values calculated over the same period.

The uncertainties of measured parameters (SST, SSS, DIC and TA; Fig. 4a–f) used in random calculations were based on the standard deviation associated with the averaging of the data in each box and month. For the calculated parameters (pH, fCO2, ΩAr; Fig. 4d–j), the uncertainties consider the errors associated with each measurement of the calculation parameters (Orr et al., 2018). The uncertainty on the trends is linked to the discrete sampling of spatiotemporal variability within the season that we wish to analyze. The uncertainty on the trends is linked either to the discrete sampling of a spatial–temporal variability within the season analyzed or to the interannual variability (real or linked to insufficient sampling). If the uncertainty associated with the annual averages is low, then the trend uncertainty will be low and will correspond to the first part. On the other hand, in the case where the interannual (real) variability is high, the uncertainty will correspond to the second part, which will decrease the trend significance.

2.3 Data selection: regions and seasons

The sampled region was separated into five 4 latitude boxes, from 46 to 64 N (Fig. 1) according to Reverdin et al. (2018b). The southern box (box A: 46–50 N) covers the shelf and continental margin of North America and is excluded in the present long-term trend analysis because of insufficient sampling and large interannual variations due to freshwater inputs from continental runoff. Box B (50–54 N) incorporates only samples where SSS is between 34 and 35 (to avoid including continental shelf or NAC waters). Thus, whereas over 470 samples with DIC data are found in box B for 1993–2017, only 399 data in the salinity range 34–35 are considered (we excluded 38 data for S<34 and 33 data for S>35). For boxes C (54–58 N), D (58–62 N) and E (62–64 N), we used all available data.

Table 1SURATLANT sampling since 1993 for each box in Fig. 1.

Download Print Version | Download XLSX

The trend analysis and contribution calculations (Eq. 2) were first performed for each box but results of trends and drivers for boxes C, D and E were found to be quasi-identical. Therefore, we present and discuss the results by combining the data collected north of 54 N (i.e., boxes C, D and E). The trends in box B are evaluated separately because properties and their trends (SST, SSS, DIC, TA) differ significantly from the other three boxes. We separate seasons, and the trends for winter (January to March) and summer (June to August) are analyzed separately as they correspond to seasonal extrema of surface carbonate properties (described below in Sect. 3.1) and these seasons were more regularly sampled along the SURATLANT line (Table 1). Finally, we choose February and July as a reference for winter and summer, respectively. When February or July data were missing for a specific year, the data collected in January or March (for winter) and in June or August (for summer) were used after correcting the observed seasonal anomalies from the climatological cycles described in Sect. 3.1. We tested this reconstruction for February and July data using SST and SSS data based on the binned products from Reverdin et al. (2018a) and averaging the bins in the same boxes (B and CDE) used for discrete sampling (Fig. S1 and Table S1 in the Supplement). As expected, the differences for SST and SSS are most pronounced for summer and box B (Fig. S1a, c) due to higher variability observed during this season, both in time and at the small scale. For winter the reconstruction based on discrete sampling appears more robust when compared to the binned products, especially for box CDE (Fig. S1b, c). To interpret our results, we thus are more confident with the trends based on winter data than summer and for box CDE. No attempt to reconstruct data was made when the SURATLANT line was not sampled for at least 3 months (e.g., winter 2016 and 2017 not included in this analysis and trends limited to 1993–2015 for this season).

3 Results

3.1 Seasonal cycle

As mentioned above, the SURATLANT sampling is not available each year for the reference months of February and July selected for the summer and winter trend analysis (Table 1). When no data are available in February or July, we adjust the representative winter and summer data based on the deviations from the seasonal cycle observed in January or March (for winter) and June or August (for summer) according to Eq. (3).

(3) SST month ( A ) year = SST month ( A ) clim + SST month ( B ) year - SST month ( B ) clim

Here, month(A) corresponds to February or July, month(B) to January or March (for winter) and June or August (for summer), and clim to climatology.

Figure 2Mean seasonal cycle of (a) SST, (b) SSS, (c) TA, (d) DIC, (e) fCO2, (f) pH and (g) Ω(Ar) over the period 1993–2017 in boxes B, C, D and E identified with corresponding colors in Fig. 1. Mean seasonal cycle of atmospheric fCO2 from Mace Head station is represented by grey dots in (e). Errors bars result from both interannual and spatial variability within a box.


The mean seasonal cycles constructed from data collected over the period 1993–2017 (only 2001–2017 for TA) are portrayed in Fig. 2. For nutrients (not shown here) the seasonal cycle was described by Reverdin et al. (2018b). The seasonal cycles deduced here from the SURATLANT data represent classical variability observed in the North Atlantic north of 50 N and are coherent with other climatology of the oceanic carbonate system constructed from different data and methods (Lee et al., 2006; Takahashi et al., 2009, 2014; Becker et al., 2018). Except for salinity and TA, all properties show a marked seasonality in all regions with maxima and minima identified at the same period for boxes B, C, D and E. For the northern box E (62–64 N, black lines in Fig. 2), we note some short time variability during summer for salinity and TA (minimum in July) that impacts the fCO2, pH and Ω cycles in this region.

The mean monthly DIC concentration (Fig. 2d) is at a maximum in winter due to deep vertical mixing, presents a steep decline from April to May when phytoplankton blooms start to occur and shows a minimum at the end of summer. The mean DIC concentrations are very similar each month north of 54 N (boxes C, D, E) and significantly lower in the southern region (box B, red line in Fig. 2) because salinity is also lower and the temperature is higher from June to November in the south. The seasonal DIC amplitude of around 70–80 µmol kg−1 in the north is slightly more pronounced in the south at around 90 µmol kg−1. Given the variability from year to year, the mean seasonal cycles of fCO2, pH and ΩAr (Fig. 2e, f, g) show fairly similar variations in all boxes. As the alkalinity seasonality is low, the seasonal changes in fCO2 and pH are anticorrelated and the seasonal variability of Ω is relatively similar to that of pH.

Although the temperature is +3 to +6C warmer in summer compared to winter (Fig. 2a), the seasonal fCO2 variability is mainly driven by DIC, i.e., the biological production in spring–summer dominates the effect of warming and deep mixing in winter dominates the cooling effect. The seasonal amplitude of fCO2 that varies between 50 and 70 µatm (Fig. 4e) is much larger than in the atmosphere (around 14 µatm at these latitudes). Consequently, for climatology, the NASPG is near equilibrium in winter and a strong CO2 sink in summer. How and why this CO2 sink changes over 1993–2017 will be investigated in the next sections.

Table 2Trends (per year) evaluated from data presented in Fig. 4. In bold are represented the significant trends (Student test).

Download Print Version | Download XLSX

3.2 Long-term trend and anthropogenic CO2

When describing the SURATLANT dataset for 1993–2017, Reverdin et al. (2018b) evaluated the long-term trend of properties over a broad region (50–63 N). Using data over 24 years (and all seasons) they estimate an increase in DIC of +0.77 (±0.11) µmol kg−1 yr−1, a fCO2 trend of +1.95 (±0.12) µatm yr−1 and a pH decrease of −0.0021 (±0.0001) yr−1. Here we re-evaluate these trends in each box (B, CDE merged) and only for summer (July) or winter (February) using observations and reconstructed data as described in Sect. 2.3. The trends for 1993–2017 in each box and season are listed in Table 2 (first lines). The DIC increasing rate ranges between +0.5 and +0.9µmol kg−1 yr−1 depending on the season and region. For fCO2 we evaluate trends ranging between +1.5 and +1.7µatm yr−1 and for pH between −0.0016 and −0.0019 yr−1. Not surprisingly, because the same original data are used, the summer and winter trends for each box (B, CDE) are consistent with preliminary results from Reverdin et al. (2018b). From 1993 to 2017, we also observed a warming, most pronounced in summer (up to +0.05C yr−1, Table 2), but the fCO2 increase and pH decrease are mainly explained by DIC increase. The long-term fCO2 and pH changes are thus mainly attributed to anthropogenic CO2 uptake. This is also supported by a decrease in δ13CDIC (Suess effect) along the SURATLANT track but for a shorter period, 2005–2017 (Reverdin et al., 2018b), a signal observed at depth in the Irminger Basin between 1981 and 2006 and linked to anthropogenic CO2 (Racapé et al., 2013).

To gain insight into the DIC trends and separate natural versus anthropogenic contributions to the fCO2 and pH trends, we evaluate the anthropogenic DIC (hereafter denoted C-ant) in the NASPG region. Here we use the TrOCA method (Tracer Combining Oxygen, Inorganic Carbon and Total Alkalinity; Touratier et al., 2007) applied to data available in the GLODAPv2 database (Key et al., 2015; Olsen et al., 2016). In the NASPG, repeated observations since 1997 were mostly conducted during summer (13 cruises in June–September). Because indirect methods such as TrOCA are not suitable to evaluate C-ant concentrations in surface waters (due to biological activity and gas exchange), we calculate C-ant in the layer 150–200 m only. For the period 1997–2010, we estimate a C-ant trend of +0.6µmol kg−1 yr−1 in the NASPG. The anthropogenic signal would explain 65 % of the DIC trend of +0.9µmol kg−1 yr−1 observed in GLODAPv2 subsurface data for the same period in the region 54–64 N, 40–20 W. Note that we obtain similar results when selecting different periods (C-ant =+0.53µmol kg−1 yr−1 for 1997–2007 whereas it is +0.58µmol kg−1 yr−1 when restricted to 2001–2007).

Recently, a new data-based method, eMLR(C), developed by Clement and Gruber, (2018) was used to evaluate the accumulation of C-ant from 1994 to 2007 in the global ocean (Gruber et al., 2019a). To compare with our estimates based on TrOCA, we explore the C-ant concentrations from eMLR(C) in the NASPG (data extracted from Gruber et al., 2019b). Along the SURATLANT line, the accumulated C-ant between 1994 and 2007 is +8.5 (±1.7) µmol kg−1 in the layer 150–200 m. This would correspond to a C-ant trend of +0.65µmol kg−1 yr−1 close to our C-ant estimate based on GLODAPv2 and the TrOCA method in the same layer. This is not surprising as Gruber et al. (2019a, b) used GLODAPv2 data as well. In the surface layer, the eMLR(C) method leads to a C-ant accumulation in the layer 0–50 m along the SURATLANT line (55–64 N) of +10.1 (±0.8) µmol kg−1 yr−1 between 1994 and 2007, i.e., about +0.8µmol kg−1 yr−1. Again, this is in the range of the long-term (1993–2017) DIC surface trends that we report (Table 2) between +0.7 and +0.9µmol kg−1 yr−1 in boxes CDE (depending on the season). Interestingly, in the region 55–64 N winter DIC from SURATLANT averaged 2134.1 (±2.8) µmol kg−1 in 1994 and 2145.1 (±2.1) µmol kg−1 in 2007, i.e., an increase of +11µmol kg−1 almost equal to the C-ant accumulation deduced from eMLR(C) over the same period. These independent estimates confirm that in the NASPG the observed long-term DIC trend and derived trends for fCO2 and pH in surface waters (Table 2, 1993–2017) are mostly linked with anthropogenic CO2 uptake. The same conclusion was drawn for the subpolar mode waters found in subsurface layers in the NASPG (Vázquez-Rodríguez et al., 2012; García-Ibáñez et al., 2016; Fröb et al., 2018). For example, based on data collected in 1981–2008, Vázquez-Rodríguez et al. (2012) estimated that 75 % of the observed pH decrease (−0.0019 yr−1) in the subarctic intermediate waters (SAIW) in the Irminger Basin is due to anthropogenic DIC. García-Ibáñez et al. (2016), incorporating four additional cruises (2010–2015), estimated a pH decrease of around −0.0018 yr−1 over 1991–2015 in the subpolar mode waters (SPMW) found at around 200–300 m in the Irminger Basin. They also conclude that in the SPMW the pH decrease is dominated by DIC increase (+0.82µmol kg−1 yr−1) through anthropogenic CO2 (C-ant trend of +0.9µmol kg−1 yr−1) but this is modulated by a pH increase over time due to an increase in alkalinity (about +0.11µmol kg−1 yr−1). From the SURATLANT time series (1993–2017) we also observe a small increase in TA during winter (+0.1µmol kg−1 yr−1, Table 2) coherent with the TA signal observed in the SPMW layer from data collected in summer (García-Ibáñez et al., 2016). The TA increase in the SPMW layer in the Irminger Basin is directly linked with salinity increase (García-Ibáñez et al., 2016) attributed to advection of higher salinity of subtropical origin and also associated with the contraction of the subpolar gyre since the mid-1990s (Häkkinen and Rhines, 2004). However, we observed a much stronger TA increase in surface water during summer, up to +0.4µmol kg−1 yr−1 in the northern region (Table 2, box CDE). Such a TA increase would lead to a pH change of around +0.0009 yr−1, counteracting the impact of DIC increase. In addition, when TA is normalized in salinity, we evaluate summer sTA trends of +0.2µmol kg−1 yr−1 in the north and +0.3 µmol kg−1 yr−1 in the south (Table 2), suggesting that the TA changes are not solely linked to salinity and advective process.

At high latitudes in the open ocean, such a long-term positive TA trend was only observed (to our knowledge) in the Pacific Western Subarctic Gyre (Wakita et al., 2013, 2017). There, in the winter mixed layer, Wakita et al. (2017) estimate positive TA and sTA trends of +0.4 and +0.34µmol kg−1 yr−1 over 1999–2015. This leads to a long-term pH decrease in this region (−0.0008 yr−1) much slower than observed in other open-ocean times series stations (Bates et al., 2014) and in the NASPG for 1993–2017 (Table 2). The TA increase also leads to slow fCO2 increase +0.9µatm yr−1 (Wakita et al., 2017). Wakita et al. (2017) suggested the TA increase might be due to weakened calcification (but this process was not quantified). Therefore, in addition to advection or convective processes, biological process such as calcification cannot be ruled out to explain TA and sTA long-term positive trends observed during summer in the NASPG (Table 2). At the high latitudes of the North Atlantic, blooms of coccolithophorid (Emiliania huxleyi) are well captured from space (Brown and Yoder, 1994), and there is indication that these blooms were more pronounced in the 1990s compared to the last 2 decades, 2000–2017 (Loveday and Smyth, 2018). Thus, in addition to DIC and anthropogenic CO2, a better knowledge of temporal TA dynamics is relevant for an understanding of fCO2 and pH changes. It is worth noting that ocean biogeochemical models (OBGMs) used to quantify the carbon cycle at a large scale are not able to fully capture TA variability when calcifying phytoplankton is not explicitly included (Ullman et al., 2009, who used 1993–2005 SURATLANT data for model comparison). When multiple phytoplankton functional groups including coccolithophores are explicitly parameterized in ecosystem carbon models, the seasonal cycle of both DIC and TA are better represented in the NASPG (Signorini et al., 2012). The same is likely true to interpret interannual to decadal TA variability and explain the drivers of observed fCO2 and pH trends along the SURATLANT line.

Because interannual variability is more pronounced in summer than in winter because of the added influence of biological activity, most studies focused on winter to analyze and interpret the decadal trends of the carbonate system in surface waters (Olafsson et al., 2009; Metzl et al., 2010; Fröb et al., 2019; Omar et al., 2019). Only a few observational studies conducted at high latitudes (Olafsson et al., 2009; Munro et al., 2015; Wakita et al 2017) showed that the trends of DIC, fCO2 and pH are seasonally different and thus driven by different processes not yet fully explained. In a global context, taking into account a seasonal view for the trend analysis is also relevant to better understand the strengthening of fCO2 seasonality observed over the last 3 decades (Landschützer et al., 2018) and how this would change the degree of acidification in the future (Kwiatkowski and Orr, 2018). Note that the observed trend of fCO2 seasonality (winter minus summer) in the North Atlantic appears less pronounced than in other regions with significant sub-decadal variability (see Figs. 1b and 4a in Landschützer et al., 2018). In this context we now investigate the observed trends for winter and summer in more detail and over specific periods to quantify the effect of temperature, salinity, DIC and TA that drive the fCO2 and pH variations in the NASPG.

Figure 3SST (a) and SSS (b) anomalies for winter (January–February) along the SURATLANT track (north of 50 N), based on the binned products constructed by Reverdin et al. (2018a) and updated data to 2018, available at (last access: 25 March 2020). The mean trends for each period selected for the carbonate data analysis (when available) are represented as black lines (1: 1993–1997; 2: 2001–2007; 3: 2008–2017). The NAO index (in blue, Hurrell data available at, last access: 30 September 2019) and AMO index (in red, Friedman ERSST V5 data available at, last access: 20 November 2019) are represented (c).


The periods 1993–1997, 2001–2007 and 2008–2017 are selected for several reasons. First, as the trends and quantification of drivers are sensitive to the data selection, we only select periods when time series for carbonate data were available (e.g., no DIC and TA data available in 1998–2000; see Table 1). This is in contrast with previous work where trends based on SURATLANT data were evaluated over 1993–2003 (Corbière et al., 2007; Metzl et al., 2010). Second, we use new binned products based on more regular SST and SSS observations (Reverdin et al., 2018a) to separate periods that present significant interannual variability and different trends in temperature and salinity (Fig. 3a, b). The three selected periods also present contrasting SST trends well identified at a regional scale in the North Atlantic, at monthly or seasonal scales (Fig. S2). Finally, in the NASPG it is well recognized that decadal and multi-decadal variability in temperature, salinity, winter convection and large-scale gyre circulation is associated with the phase of the NAO or AMO (Fig. 3c).

During the first observational period 1993–1997 (rather short), the NAO shifted from a positive to negative phase while the AMO was in a transitional stage and regularly increased during the 1990s (Fig. 3c). During the positive NAO in the early 1990s, the NASPG experienced a deep convection regime in the western subpolar gyre (Pickart et al., 2003). During the second period 2001–2007, the NAO was negative or neutral and the AMO reached a high index. The convection in the western NASPG was relatively shallow during this period and the ocean surface warmed (Fig. 3a). In the last decade, 2008–2017, the NAO was highly variable compared to previous periods, with the lowest phase in 2010 and highest phase in 2015 (Fig. 3c). Conversely, the AMO has strongly decreased since the late 2000s. During this period, both SST and SSS anomalies present clear negative trends: the NASPG becomes colder and fresher (Fig. 3a, b) during this period, also associated with several deep convection events occurring in the western subpolar gyre during the winters 2008, 2012 and 2015 (Våge et al., 2008; de Jong and de Steur, 2016; Fröb et al., 2016; Piron et al., 2017).

Figure 4Evolution of (a) SST, (b) SSS, (c) TA, (d) sTA, (e) DIC, (f) sDIC, (g) fCO2, (h) ΔfCO2, (i) pH and (j) ΩAr between 1993 and 2017 is obtained in box B (open symbols) and by combining all data in boxes C, D and E (filled symbols) during summer (in red) and winter (in blue). Data from February and July are indicated with circles, and the reconstructed data are depicted with triangles.


3.3 Winter trends for different periods

The trends and drivers in winter are identified as blue symbols in Figs. 4 and 5. During the first period, 1993–1997, the results are contrasted between the southern and northern regions. The strongest change occurs in the south (box B), where we observe a rapid cooling over 4 years (−0.19C yr−1), an increase in salinity (+0.045 yr−1) and an increase in TA (+2µmol kg−1 yr−1) that is directly linked to salinity (recall that for 1993–1997 TA was not measured and the TA values are based on salinity, Eq. 1). Over this short period, this leads to a strong decrease in fCO2 (−7µatm yr−1) and an increase in pH, both attributed to the cooling and TA changes (Fig. 5a, d), the effect of DIC being minor. In the north (boxes CDE), we do not observe any significant trend of winter properties in 1993–1997 for SST, sTA, sDIC, pH and fCO2. As a result, in the NASPG (all region), the oceanic fCO2 trend does not follow the atmospheric CO2 increase, and the air–sea disequilibrium increases; i.e., the region was a CO2 source in winter 1994 and reached near equilibrium afterwards (Fig. 4h).

Figure 5Decomposition of the trends in surface pH (a–c), fCO2 (d–f) and ΩAr (g–i) according to Eq. (2). The effect of the changes in SSS, SST, DIC, TA and the normalized part (sDIC and sTA) is shown for the three periods: 1993–1997 (left column), 2001–2007 (middle column) and 2008–2017 (right column). Color coding is the same as in Fig. 4. Here, DIC (or TA) represents the total carbon (or total alkalinity) part, whereas sDIC (or sTA) is the part of the DIC (or TA) driver not related to salinity (see explanations for Eq. 2).


As opposed to the 1990s, data collected in winter during the second period (2002–2007) show a gradual DIC increase in the two regions (+1.7 and +2.1µmol kg−1 yr−1, Table 2) higher than expected from the anthropogenic CO2 increase. A decrease in TA is also observed in the northern region but not directly related to salinity. This leads to a significant pH decrease (−0.006 to −0.008 yr−1) and fast fCO2 increase (+6 to +7µatm yr−1). For this period, the DIC increase explains the temporal changes in pH and fCO2 in box B, whereas both DIC and TA drive these changes in box CDE (Fig. 5b, e). In contrast to 1993–1997, the winter oceanic CO2 source increases during 2002–2007. It is near equilibrium in 2002 and a source in 2004–2007 (Fig. 4h), a result previously confirmed with independent sea surface fCO2 observations (Metzl et al., 2010). At the end of this period, in winter 2007, oceanic fCO2 was higher than 400 µatm and pH was low, near 8.04.

The last period, 2008–2015, was not investigated in previous studies (Corbière et al., 2007; Metzl et al., 2010). During 2008–2015, the DIC concentrations in winter continue to rise (+1.0 to +1.3µmol kg−1 yr−1, Table 2) but at a slower pace compared to previous years and still higher than the anthropogenic DIC trend. Contrary to 2001–2007, we do not observe large changes in TA. In the northern region (box CDE) the sea surface cools (−0.05C yr−1), which is not found in the southern region, inducing different fCO2 and pH trends in the southern and northern sectors. In the north, the low fCO2 and pH trends are due to a DIC increase compensated for by cooling and a small TA increase, whereas in the south the fCO2 trend (+2.5µatm yr−1) and pH trend (−0.0024 yr−1) are mainly driven by DIC increase (Fig. 5c, f). These trends deviate from the ones in 2001–2007 with different processes possibly related to the NAO presenting much more variability during the last period (Fig. 3c). At the end of our time series and in both regions fCO2 in winter 2015 reaches the highest value (410 µatm) recorded since 1993 and the lowest pH values (8.03), with the same for ΩAr (1.7). This is mainly explained by DIC concentrations approaching 2160 µmol kg−1 in February 2015 (Fig. 4e). These high DIC concentrations observed in winter 2015 correspond to a high positive NAO phase (Fig. 3c) and to negative SST anomalies (Fig. 3a). It was also associated with strengthening of fresh water coming from the western NASPG (Holliday et al., 2020) and likely resulted in increased vertical mixing over the Reykjanes Ridge (De Boisséson et al., 2012). The high DIC and sDIC concentrations observed in winter 2015 are not unique to the Irminger Sea. Recently, based on sea surface fCO2 measurements conducted in winter 2004–2017 and DIC calculated from fCO2 data and reconstructed TA (using an adapted TA–SSS relation), Fröb et al. (2019) report average winter sDIC concentrations of around 2160 µmol kg−1 in 2012 and 2014. These authors also observed high sDIC concentrations in winter 2017, above 2165 µmol kg−1. Unfortunately, we cannot compare these values as we have no winter data after 2015 along the SURATLANT line (Table 1). Although the periods are not exactly the same and methods to evaluate trends are different, the sDIC trends reported by Fröb et al. (2019) are in the same range as here: between +1.3 and +1.75µmol kg−1 yr−1 for the period 2004–2017, when we find (i.e., boxes CDE including the eastern Irminger Sea) a sDIC trend of +1.8µmol kg−1 yr−1 for 2001–2007 and +1.3µmol kg−1 yr−1 for 2008–2015 (Table 2). Both results show a clear increase in sDIC in the NASPG or Irminger Sea since 2001 (or 2004), and the differences in trends might also be modulated by interannual variability.

During 2008–2015, we also observe significant interannual variability in winter, with marked DIC minima in 2010 and 2013 also found in sDIC concentrations as well as in TA and sTA (Fig. 4c–f). Interestingly, Fröb et al. (2019) also report lower sDIC concentrations in 2010 and 2013 in the Irminger Sea. In 2010, the NAO shifted to a negative phase (Fig. 3c) supporting relatively shallow convection in the NASPG and inducing positive SST anomalies (Fig. 3a); the absence of deep mixed layers during these years would explain the relatively low DIC and TA concentrations observed in winter (i.e., less input of high DIC and TA from subsurface layers). In 2013, although the NAO was neutral, we also observed lower DIC and TA in winter. It is suggested that mixing was relatively shallow in the Irminger Sea and near the Reykjanes Ridge in 2013, which could drive less renewal of the surface layers by the enriched sDIC deeper waters and thus negative winter DIC and TA anomalies for this particular year. However, as both DIC and TA decrease, fCO2 and pH do not strongly change for these particular years (2010 and 2013), and associated air–sea CO2 disequilibrium remains stable (Fig. 4h). The strong negative DIC anomalies observed in 2010 and 2013 nonetheless influence the DIC trend for the period 2008–2015: not including those 2 years increases the trend in boxes CDE from +1.0 to +1.3µmol kg−1 yr−1. However, this is still lower compared to 2002–2007 (+1.7µmol kg−1 yr−1), and thus a much slower trend for fCO2 and pH is estimated in recent years. As a result, in 2008–2015 we interpret the fCO2 and pH winter changes equally due to the contribution of anthropogenic CO2 (C-ant) and the natural component (DIC-nat), whereas in 2001–2007 the natural component dominates due to the dynamics in the NASPG (Fig. 6c, d).

3.4 Summer trends for different periods

We examine the trends and drivers in summer over the same periods selected for winter (except for the last period, with 2001–2017 in summer and 2002–2015 in winter). Results for summer are identified by the red symbols in Figs. 4 and 5. As mentioned earlier, a trend analysis in summer was not performed in previous work based on SURATLANT data (Metzl et al., 2010). Here we attempt for the first time to detect trends and processes occurring during this season and how they impact fCO2, air–sea CO2 equilibrium and pH changes. Not surprisingly, the temporal variability of the carbonate properties is much more pronounced during summer when biological processes generally starting in spring imprint large DIC variations at the seasonal scale (Fig. 2). The biological bloom and its timing may also lead to significant interannual variability, but, as we use July as a reference for summer, we expect to record the low DIC concentrations each year due to accumulated carbon uptake though production occurring in spring–summer.

During the first period, 1993–1996, the most remarkable feature is a very rapid increase in DIC in the northern and southern sectors (+4 and +9µmol kg−1 yr−1). This contrasts with the DIC trends in winter over 1994–1997 (Table 2). Because the northern box (CDE) also experienced a sharp warming of +0.6C yr−1, there the fCO2 increase (+12 to +14µatm yr−1) and pH decrease (−0.015 to −0.017 yr−1) are particularly fast. However, the period is too short to clearly interpret which are the main processes at play and derive any conclusion on the trends. We notice that in 1993–1996, the NAO changed from a positive to negative phase and the AMO progressively increased in the 1990s (Fig. 3c), but no particular anomalies were revealed in the winter observations. Thus, we have no direct explanation for the fast DIC, fCO2 and pH changes observed during the summer in 1993–1996, except that relatively higher DIC in 1996 might have resulted from a decrease in primary production compared with previous years. This is rather speculative as we have no other information on nutrients or on the biological activity (e.g., no nutrient data for summer 1996 and no remote sensing data at the beginning of the 1990s). However, using continuous plankton recorder (CPR) data over 1960–2010 in both the western and eastern NASPG regions, Martinez et al. (2016) identified negative anomalies of Chl a around 1996–1997 after a sharp increase in Chl a during 1986–1995. The fCO2 data show that the region was a large CO2 sink in summer in 1993–1995 but abruptly increased to near-equilibrium values in 1996 (Fig. 4h). In the northern region (box CDE), the mean delta-fCO2 value in summer 1996 was close to 0 µatm, i.e., about the same as observed during the winter in 1996–1997.

The SURATLANT regular sampling for DIC, TA and nutrients was restarted in June 2001 (Table 1). During 2001–2007, and as in summer in 1993–1996, we observe a rapid increase in fCO2 (+8 to +12µatm yr−1) and pH decrease (−0.010 to −0.013 yr−1). Despite having reconstructed July from the data obtained in August (for boxes B and CDE) for both the starting (2001) and ending (2007) years, this result is fairly robust. Indeed, if we do not normalize summer to July, but directly include the June–July–August data, we also observe a rapid increase in fCO2 of +6.2µatm yr−1. The summer DIC and TA trends in the northern and southern areas have the same sign (similar to during winter), positive for DIC and negative for TA but not directly related to salinity (Table 2). However, the DIC increasing and TA decreasing rates are regionally significantly different. The DIC increase is most pronounced in the north (+5µmol kg−1 yr−1), and the TA decrease is more pronounced in the south (−3.9µmol kg−1 yr−1). Although the fCO2 and pH trends are similar in the two areas, this suggests different drivers for boxes B and CDE (Fig. 5b, e). In the north, DIC explains most of the fCO2 and pH change, whereas in the south the TA contribution dominates. In both regions, as in winter, temperature and salinity changes have a small effect for this period (Fig. 5b, e). The rapid oceanic fCO2 increase strongly impacts the variations in air–sea CO2 disequilibrium (Fig. 4h). In the north, the region was a sink in 2001–2004 but reaches equilibrium in 2005 and 2007 (no data in July 2006 in the north). In July 2007, temperature was relatively low (around 10 C) and the high fCO2 at that period was mainly due to high DIC concentrations linked to low productivity that year.

For the last period, 2008–2017, results in summer are very different compared with previous decades as we now observe a decrease in fCO2 and an increase in pH in both regions (Fig. 4g, i; Table 2). This also contrasts with the winter trends in 2008–2015. The variation in fCO2 leads to a strong ocean CO2 sink with low delta-fCO2 in 2010–2017 (−40 to −80µatm, Fig. 4h). During the summer in 2008–2017 sea surface properties are very variable and this might be related to the NAO variability in this period (Fig. 3c) but also to changes in the ship's route in some years or to productivity occurring at the mesoscale. This was the case in July 2013 where the sampling took place further northwest (Fig. 1), impacting SST (colder, Fig. 4a) and DIC (stronger, Fig. 4e, f). However, these anomalies are less noticeable for fCO2 and pH (Fig. 4g, i) because the effects of SST and DIC partially compensate. Thus, we chose to not consider July 2013 for calculating trends and drivers.

Although the interannual variability is large in summer, the trends are evaluated over the period 2008–2017 using all available data (except July 2013). For this period all terms (SST, DIC and TA) contribute to the fCO2 and pH changes (Fig. 5c, f). The DIC increase is only revealed in the north (box CDE, +1.2µmol kg−1 yr−1); this is slightly higher than in winter (+1.0µmol kg−1 yr−1) and thus higher than anthropogenic signal (around +0.6µmol kg−1 yr−1). However, as the surface ocean cools (−0.11C yr−1) and TA increases (+1.7µmol kg−1 yr−1), the net effect is a decrease in fCO2 (−2.3µatm yr−1) and increase in pH (+0.0026 yr−1). During 2008–2017, the progressive cooling in the NASPG, found here in both winter and summer data (Fig. 4a) is a large-scale signal (Fig. 3a, Robson et al., 2016; Reverdin et al., 2018a). An intriguing signal in 2008–2017 is the increase in TA (Fig. 4c) opposed to the decease in salinity (Fig. 3e, Fig. 4b); this leads to a positive sTA trend opposite to summer in 2001–2007. Therefore, we cannot interpret the observed TA increase directly linked to salinity. The shift of the sTA trend, negative in 2001–2007 and positive in 2008–2017, is an important signal that drives opposite trends for fCO2 and pH between the two periods (Fig. 5b, c, e, f).

Figure 6Decomposition of the trends in surface pH according to the effect of the changes in DIC-tot (filled), DIC-nat (dotted) and C-ant (open) during summer (a, b) and winter (c, d) for the long-term period and the three periods.


This contrasting TA signal in summer cannot be attributed to error or drift in laboratory analyses as the measurements were performed using the same methods for the whole time series (Reverdin et al., 2018b), and no such signal is identified for winter cruises in 2008–2015. In addition, during OVIDE cruises conducted in summer in 2006–2014 in the North Atlantic (e.g., García-Ibáñez et al., 2016), we performed regular TA and DIC intercomparisons with the ICM/CSIC group in Vigo (F. Pérez) and certified our results to within around ±4µmol kg−1 for both TA and DIC. Results based on the surface TA samples measured at LOCEAN for the OVIDE cruises in summer in 2006–2018 (data in Metzl et al., 2018) also suggest that sTA increases in the NASPG at a rate of around +1.5µmol kg−1 yr−1.

To support these results, we have evaluated sTA trends in the NASPG from independent data available in the most recent GLODAPv2.2019 version (Olsen et al., 2019). For this we selected the data in the layer 0–20 m for all cruises conducted in June–August in 1997–2014. We found a significant difference of sTA trends during this period: −0.52µmol kg−1 yr−1 in 1997–2006 against +0.54µmol kg−1 yr−1 in 2006–2014. Although the periods are not exactly the same and sTA trends from GLODAPv2.2019 are smaller than those deduced from the SURATLANT time series, the changing sTA trends from negative to positive in recent years are also features in that dataset. We are thus confident with TA data over time and need to find a process that explains why the sTA trend was negative in summer in 2001–2007 (−1.9µmol kg−1 yr−1) and positive in summer in 2008–2017 (+2.3µmol kg−1 yr−1).

The variability of calcification through the production of calcifying species (e.g., coccolithophores, foraminifera) is a possible mechanism that would impact TA trends as suggested by Wakita et al. (2017) in the Pacific Western Subarctic Gyre. In the northeast Atlantic, long-term in situ continuous plankton recorder (CPR) observations showed a significant increase in the calcifying species starting in the mid-1990s (Beaugrand et al., 2013). These authors suggest that the temperature (a warming) was the main driver of the positive trend of calcifying plankton. In addition they show that no correlation was identified between the NAO and species variability but that a positive trend of the calcifying species was correlated with the AMO, a result also confirmed by Rivero-Calle et al. (2015). The TA decrease in summer we observe during the period 2001–2007 after the AMO moved to its positive state (Fig. 3c) might be explained by the increase in calcifying species as identified from in situ CPR observations in the North Atlantic (Beaugrand et al., 2013; Rivero-Calle et al., 2015). Unfortunately, we do not yet have direct in situ observational evidence of a reduced calcifying species after 2010. However, if one follows the proposed scenario for 2001–2017, the TA increase in summer in 2008–2017 during the cooling phase in the NASPG could be linked to a weakening of calcification. In this region, this seems supported by the absence of coccolithophore blooms in recent years (2010–2017) as identified from remote-sensing reflectance records (Loveday and Smyth, 2018).

4 Discussion

Our observations collected over the last 3 decades show an abrupt change in the evolution of hydrological and biogeochemical properties in the NASPG around the year 2007 (Fig. 4) when the AMO reached a maximum and the NAO was around neutral (Fig. 3c). Following a warming since the mid-1990s, the region experienced a large cooling and freshening after 2005 (Fig 3a, b; Robson et al., 2016; Holliday et al., 2020). This change is associated with an increase in the size of the gyre and increased currents along the gyre's southern rim (Chafik et al., 2014, 2019; Desbruyères et al., 2015), as well as by a very large heat loss during positive NAO years, particularly in 2015 (Fig. 3c, Josey et al., 2018). Indeed, the NAO was previously recognized as a possible cause of the rapid fCO2 increase when the NAO shifted from the positive to negative phase in 1995–1996 (Fig. 3c; Corbière et al., 2007). However, this was not confirmed for the period 2001–2007 when the NAO did not vary so much around a neutral value (Fig. 3c; Metzl et al., 2010). For the whole NA, Schuster et al. (2009) synthesized fCO2 observations over 1990–2006, i.e., before 2007 when we observe abrupt changes in property trends. North of 45 N, Schuster et al. (2009) evaluate a trend of fCO2 exceeding +3µatm yr−1 and a significant decrease in the CO2 sink in the NA by over 50 % between 1990 and 2006. They also predict an increasing sink in the subpolar regions following the increasing NAO index in 2007. Although the impact of the climatic modes (NAO and/or AMO) on the oceanic physical properties and circulation has been well established, their effects, if any, on the fCO2 and pH trends need to be clarified.

Whatever the NAO variability, the long-term fCO2 and pH trends we evaluate over 1993–2017 in the NASPG for summer or winter are mainly explained by the increase in DIC associated with the uptake of anthropogenic CO2 (Fig. 6); that is, the DIC or sDIC trends (Table 2) are not significantly different from the anthropogenic DIC trend estimated between 1994 and 2007 in this region (Gruber et al., 2019a). For fCO2, here calculated from DIC–TA pairs, the long-term trends between +1.5 and +1.7µatm yr−1 in the NASPG (Table 2) are slightly higher than the mean trend of +1.47 (±0.06) µatm yr−1 evaluated for 1992–2014 in the whole NA (Lebehot et al., 2019) based on monthly reconstructed fCO2 using a multiple linear regression (MLR) approach and SOCAT v4 fCO2 data (Bakker et al., 2016). Lebehot et al. (2019) show that the fCO2 trends based on observations are much lower than those derived from 19 ESM CMIP5 models, +1.90 (±0.09) µatm yr−1. By performing several sensitivity test analyses on the ocean ESMs, Lebehot et al. (2019) conclude that the discrepancy between observed and simulated fCO2 trends originates mainly in the model's biogeochemistry, e.g., biases in simulated TA and also related to the way ESMs do or do not represent winter mixing (Goris et al., 2018). This might be especially relevant for the NA subpolar region where some ESMs project faster change in ocean fCO2 in the future (2061–2100) with 60 % due to the DIC increase and up to 29 % due to the TA changes (Tjiputra et al., 2014). We also suspect that ESMs used to predict future change of the oceanic CO2 sink and ocean acidification would produce faster pH trends than what we observed.

If the long-term fCO2 and pH trends could be mainly explained by anthropogenic CO2 uptake in the NASPG, at shorter timescales (4–10 years) the trends are very different (Table 2). Indeed, the gradual changes in fCO2 and pH caused each year by the uptake of anthropogenic CO2 can be significantly masked by the natural variability of DIC (DIC-nat), temperature and/or TA. To better identify when and why the DIC-nat dominates, we separately compute the impact of anthropogenic CO2 (C-ant) and DIC-nat on pH and fCO2 trends for each season and period. The results for pH are presented in Fig. 6. The same results are obtained for fCO2 (not shown), with an opposite sign for each bar plotted in Fig. 6. For C-ant we adopt a value of +0.6µmol kg−1 yr−1 as described in Sect. 3.2.

For all sub-periods and both seasons, the effect of DIC-nat on pH trend is significant and with similar or higher magnitude than the effect of C-ant (Fig. 6). Occasionally it opposes C-ant, i.e., DIC-nat decreases with a positive effect on pH trend (box B, 2008–2017 in summer, Fig. 6a, and box CDE, 1994–1997 in winter, Fig. 6d). Over 1993–1997 the DIC-nat effect on pH trends for this short period is opposed in summer and winter, suggesting that this is not linked to changes in regional circulation, e.g., less or more input of DIC from different water masses. For 1993–1997 when the largest DIC-nat effect is observed in summer (in the northern and southern boxes), one might suggest that it is linked to primary productivity; unfortunately, prior to SeaWiFS in 1998, we have no direct or indirect information on biological changes to explain why DIC-nat increased in summer in 1993–1997. However, the length of the period is short and trend results are very sensitive to interannual anomalies, especially in the first and last years.

Apart from the first short period (1993–1997), the largest effect of DIC-nat is observed in summer in 2001–2007 in the northern region (Fig. 6b). This period includes a decade (1995–2005) where the strength of the subpolar gyre circulation decreased (Häkkinen and Rhines, 2004, 2009; Häkkinen et al., 2011, 2013), suggesting that more water of subtropical origin penetrates the NASPG. This would have decreased the DIC surface concentrations (and increase TA as well), but for 2001–2007 we observed the opposite, including in winter (Table 2). We thus eliminate the effect of advection to interpret the DIC-nat increase in the NASPG that impacts pH trend significantly. As mixed-layer depths are not deeper than 30–40 m in summer and present no significant interannual changes (according to the dataset Armor which presents the same variations as the data of the product GLOBAL-REANALYSIS-PHY-001 030 distributed on, last access: 28 October 2019), the variation in the vertical mixing during this season is not a likely candidate to explain the changes in DIC-nat. A possible explanation for the large contribution of DIC-nat to pH (and fCO2) in 2001–2007 could be a decline in biomass identified in SeaWiFS data in the eastern subpolar region for the years 1998–2007, although the trends for net primary productivity (NPP) appear to not be statistically significant over 10 years (McKinley et al., 2018). For the NPP subject to high-frequency variability, the evaluation of trends over nearly 10 years is challenging. Annual biomass anomalies (based on SeaWiFS and MODIS sensors) changed from positive values in 1998–2004 to negative ones in 2005–2009 (McKinley et al., 2018). To explain the biomass decline, McKinley et al. (2018) used a coupled physical–biogeochemical ocean model (OBGM) to reproduce these changes and found that nutrient concentrations decline significantly in the region north of 50 N due to reduced physical supply through horizontal and vertical fluxes; their model suggests that enhanced phosphate and silicate limitation over time dominates the light limitation in this region.

Due to limited sampling months and the high variability of nutrients in spring–summer, we have not been able to detect such a mechanism from the SURATLANT nutrient data. However, the reduced productivity due to nutrient limitations (McKinley et al., 2018) is supported by independent observations of silicate concentrations in the subpolar NA over 25 years (McKinley et al., 2018). They showed a decline in pre-bloom silicate concentrations in the winter mixed layer since the 1990s until 2010 and attributed it to the decrease in winter convection depth and the weakening of the NASPG. This decline of silicates observed over 25 years would negatively impact diatom blooms in spring and favor coccolithophore blooms occurring in summer. Such a scenario is coherent with the increase in calcifying species since the mid-1990s as observed from the CPR data in the NA during a warm period associated with the AMO (Beaugrand et al., 2013; Rivero-Calle et al., 2015). The shift of phytoplankton species and local or regional intensified calcification might also explain the low sTA concentrations occasionally observed in 2003 and 2005, i.e., lower than for summer 2001 and also for winter (Fig. 4d). The interannual variability of TA of around ±8µmol kg−1 in 2001–2007 is twice the TA change of ±4µmol kg−1 due to coccolithophores in a biological model applied at 60 N–30 W (Signorini et al., 2012). The TA and sTA negative trends observed in 2001–2007 (Table 2) contribute less to pH and fCO2 trends than with DIC (Fig. 5b, e) but nearly the same magnitude compared to the contribution of C-ant. This is not the case after 2007.

During the last decade, 2008–2017, observations were obtained during a strong negative NAO in 2010 and a positive NAO phase in 2015 (Fig. 3c). The NAO presents large interannual variation compared to 1995–2007 (Fig. 3c) and a decline in AMO index after 2010 (Fig. 3c). The 2010 event was associated with a warming and freshening (Fig. 4a, b) found in both SURATLANT discrete sampling (in August 2010) and monthly reconstructed binned products (Reverdin et al., 2018a). In August 2010 the DIC concentrations were very low in the north (<2070µmol kg−1) compared with other years since 1993 (Fig. 4e). This was not associated with particular signals in TA but with high δ13CDIC as reported in Racapé et al. (2014) and also with higher Chl a concentration identified from MODIS data (McKinley et al., 2018). Indeed, Henson et al. (2013) showed that the physical forcing caused by the very negative NAO recorded in winter 2009–2010 stimulated spring blooms and not the eruption of the volcano Eyjafjallajökull in Iceland, which erupted in spring 2010, depositing large amounts of iron in the North Atlantic Subpolar region. The effect of the higher productivity in summer 2010 is due to lower fCO2 and increased pH compared to previous summers. Thus, there was a rapid drop in delta-fCO2 in 2009–2010, such that summer 2010 was a strong CO2 sink. In 2015, when the NAO was in a positive phase, the SST anomaly was on the order of −1C (Fig. 3a): for that year, DIC was high in winter, close to the maximum observed in winter in our time series (Fig. 4e, f). As the temperature also lowers fCO2, the fCO2 values (and pH) were not so different from other winters, illustrating the competing effects of cooling and deeper vertical mixing on fCO2. Although significant changes are observed for DIC variability during NAO events in 2010 and 2015, they seem to have a small impact on the trends. After 2007 the positive DIC trends are less important than in 2001–2007 (Table 2) and the contributions of the natural and anthropogenic parts of DIC have a similar magnitude in summer and winter (Fig. 6). This result alone does not explain the decrease in fCO2 and the increase in pH observed during the last period in summer. Indeed, there is also a significant impact of cooling and, as opposed to 2001–2007, an intriguing increase in TA, both leading to decreasing fCO2 and increasing pH trends (Figs. 4, 5).

5 Conclusion and perspectives

Based on sea surface observations of DIC and TA collected in the NASPG over 1993–2017, we have analyzed the variability and trends of the carbonate system properties including fCO2, pH, ΩAr and ΩCa calculated from the DIC–TA pairs. This study extends to summer, and for pH and ΩAr there is previous work based on winter data (Corbière et al., 2007; Metzl et al., 2010). It also extends the analysis for the last decade, 2008–2017, after the AMO reached a maximum and then decline and when the NAO was highly variable.

In the last decade we observed a continuous surface DIC increase in winter in the northern and southern NASPG. In February 2015 when the NAO was in a positive state, DIC reached a maximum concentration (DIC >2150µmol kg−1, i.e., more than +20µmol kg−1 higher than in the 1990s). In 2015 pH was at a minimum (8.03), and fCO2 exceeded 400 µatm and was close to the atmospheric fCO2. Such high fCO2 was also observed in the NASPG in January–February 2015 from direct underway measurements (range 405–415 µatm for cruises AGFO20150115 and AGFO20150212 (measured by Wanninkhof NOAA/AOML on the MV Skogafoss) in SOCAT v5; Bakker et al., 2016). As opposed to the wintertime DIC continuous increase observed since the 1990s, the TA decadal variability is not uniform. In 2001–2007 the decreasing TA, added to the DIC changes, reinforced a rapid fCO2 increase of up to +7µatm yr−1 and a strong pH decline of around −0.007 yr−1 confirming previous studies (Metzl et al., 2010). Conversely, in 2008–2015 the increasing TA and the cooling in the northern NASPG compensated for the effect of DIC increase, leading to much smaller winter trends for fCO2 and pH, on the order of +1µatm yr−1 and −0.001 yr−1, respectively.

During summer, the interannual variability of all properties is much more pronounced due to active primary productivity in the NASPG in spring–summer and shallow mixed layers. Consequently, the decadal trends of the carbonate system properties in summer are more difficult to detect compared with winter or other oceanic regions such as the subtropics (Bates et al., 2014; Ono et al., 2019). In addition, primary production often occurs at small spatiotemporal scales, and the DIC–TA heterogeneous sampling may have occasionally missed planktonic blooms. In 2001–2007 the summer trends of SST, DIC and TA have the same sign as in winter. This confirms the fast increase in fCO2 and the strong decline of pH observed during this period. The natural variability of DIC dominates the effect of the anthropogenic uptake on fCO2 and pH trends, especially in the northern region (Fig. 6). The DIC increase in summer is likely due to a reduced productivity during this period (Tilstone et al., 2014; McKinley et al., 2018). After 2007, the fCO2 and pH trends for summer are drastically different from the previous decade. In both southern and northern parts of the transect, the fCO2 trend becomes negative while the pH trend is positive (up to +0.0052 yr−1 in the south). This is driven by a complex interplay of cooling (−0.1C yr−1), DIC increase and significant increase in TA more pronounced in summer than in winter.

Before 2007, we evaluate a rapid fCO2 increase faster than in the atmosphere. As a result, the ocean CO2 uptake decreased in the NASPG from 1993 to 2007, in agreement with other studies (Schuster et al., 2009; Landschützer et al., 2013). However, at larger scale, here for the NA-SPSS regional biome (Fay and McKinley, 2014), the decreasing CO2 sink in the NASPG for this period is not always resolved from fCO2 data-based methods that evaluate an increasing CO2 sink in the North Atlantic Subpolar region (Rödenbeck et al., 2015; Denvil-Sommer et al., 2019). This disagreement might be in part due to missing fCO2 data for the period 1994–2003 in the NASPG (Bakker et al., 2016) and need to be investigated in further studies. After 2007, our results show that winter fCO2 increased at a lower rate than in the atmosphere, whereas in summer we observe a decrease in fCO2. Thus, the ocean CO2 sink increases over 2008–2017 (Fig. 4h). For this period this is coherent with the results derived from data-based methods in the NA-SPSS biome but only for 2007–2013. Indeed after 2013 the indirect methods produce either an increasing CO2 sink or a decreasing CO2 sink and with a very high variability noticeable in 2015 when the NAO was positive (four indirect methods present in Fig. S5h in Denvil-Sommer et al., 2019).

The observed change in all sea surface properties around 2007 in the NASPG, as well as a decreasing CO2 sink before 2007 and an increasing CO2 sink after 2008 (Fig. 4h), seems linked to the AMO rather than directly to the NAO (Fig. 3c) as also suggested by data-based reconstructed fCO2 fields in the North Atlantic (Landschützer et al., 2019).

The temporal change in fCO2 and CO2 uptake around 2007 that we deduce from SST, DIC and TA is also clearly observed in pH variability. In summer, the rapid pH decline in 2001–2007 is followed by a significant pH increase in 2008–2017. This is due to a changing trend in TA that was negative before 2007 and positive after 2008, a signal observed in both seasons. This unexpected decadal change in TA, also observed in the northwestern subpolar Pacific, might be linked to changes in calcification processes (Wakita et al., 2017). Indeed, it has been recognized that calcifying species in the North Atlantic present significant variations since the mid-1990s (Beaugrand et al., 2013). Quantification of how this process impacts fCO2 and pH variability deserves further studies to investigate the coupling of chemical measurements with species observations such as those obtained from CPR in the NASPG.

Although we identified significant interannual to decadal variability in surface carbonate system properties over the full time series (1993–2017), the long-term trends of fCO2 and pH in winter and summer are almost entirely explained by the DIC increase and anthropogenic CO2 uptake (Fig. 6). The long-term trend of fCO2 in the NASPG (between +1.5 and +1.7µatm yr−1) is slightly higher than the mean trend for 1992–2014 (+1.47µatm yr−1) evaluated at a large scale in the NA (Lebehot et al., 2019). Our results now extended to 2017 confirm the fCO2 trends evaluated by McKinley et al. (2011) for the period 1981–2009, and our new estimate for 1993–2017 in summer (+1.7µatm yr−1) is not significantly different from the trend of +1.8µatm yr−1 estimated for 1981–2002 based on a few observations in August (Corbière et al., 2007). However, it is worth noting that in other NA sectors and periods, such a fCO2 increase is not always observed. In a recent study focused on the midlatitude NA (40–50 N), Macovei et al. (2019) showed that fCO2 is highly variable with a small trend of +0.37 (±0.22) µatm yr−1 in 2002–2016, implying a significant increase in the CO2 uptake in this region, especially after 2010 when we also observed a sudden drop of fCO2 in summer in the NASPG (Fig. 4g, h). Our results added to those from Macovei et al. (2019) suggest the ocean CO2 sink increased in the North Atlantic (north of 40 N) at least since 2007 and until 2017, contributing to the increase in the global ocean CO2 sink (Friedlingstein et al., 2019). At the NA basin scale this is consistent with results deduced from fCO2 reconstructed data-based methods (Denvil-Sommer et al., 2019; Gregor et al., 2019) but opposed to the ocean CO2 sink variability generated by current ESMs, suggesting uncertainties for predicting the evolution of the NA CO2 sink in the future (Lebehot et al., 2019).

For pH the long-term trend in the NASPG of −0.0017 yr−1 is in the range of what is observed in other oceanic regions but we note some differences. This trend is equal to the trend estimated at BATS station but is lower than in the Irminger Sea and higher than in the Iceland Sea (Bates et al., 2014). It is also lower than the pH trend of −0.0020 yr−1 estimated in the North Atlantic Subpolar region in 1991–2011 (Lauvset et al., 2015) and −0.0024 yr−1 recently observed in the NA waters in the North Sea in 2004–2015 (Omar et al., 2019).

Given the large differences of fCO2 and pH trends at the regional scale listed above, our results among many others highlight the need for acquiring sustained time series of the ocean carbonate system in different regions as was strongly recommended at the recent Ocean-Obs19 conference (Tilbrook et al., 2019; Wanninkhof et al., 2019).

An understanding of these differences also calls for a comprehensive analysis including a synthesis of all DIC and TA sea surface observations collected in different regions. This should be achieved at an international level as is done for sea surface fCO2 in the frame of SOCAT (Bakker et al., 2016) or for CLIVAR GO-SHIP cruises assembled in GLODAP (Olsen et al., 2019). In addition to ship-based observations, the analysis of data from BGC Argo floats equipped with pH sensors (together with temperature and salinity sensors, from which TA, DIC and fCO2 can be estimated (e.g., Williams et al., 2017) will help to better constrain spatial, seasonal and interannual variability. Such data synthesis would also help to validate ocean and Earth system models that at present do not correctly represent the temporal change in marine biogeochemistry as demonstrated by Lebehot et al. (2019) for the North Atlantic.

Data availability

The dataset is freely available and is accessible at ( (last access: 10 March 2020), Reverdin et al. (2018b).


The supplement related to this article is available online at:

Author contributions

CL produced the data analyses and wrote the manuscript with input from NM, GR and CLM. GR, NM and VR produced the data synthesis. SO provided the nutrient data. JF provided the DIC and TA data.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “The 10th International Carbon Dioxide Conference (ICDC10) and the 19th WMO/IAEA Meeting on Carbon Dioxide, other Greenhouse Gases and Related Measurement Techniques (GGMT-2017) (AMT/ACP/BG/CP/ESD inter-journal SI)”. It is a result of the 10th International Carbon Dioxide Conference, Interlaken, Switzerland, 21–25 August 2017.


The SURATLANT project is supported by the French institute INSU (Institut National des Sciences de l'Univers), the Lamont-Doherty Earth Observatory (LDEO), the National Oceanic and Atmospheric Administration (NOAA) – Atlantic Oceanographic and Meteorological Laboratory (AOML) and the Climate Program Office (CPO). We thank the EIMSKIP Company and the MFRI team, both based in Reykjavík (Iceland), for their cooperation in seawater sampling and analysis. We also thank the numerous scientific volunteers who worked at sea as well as the crew and captain of the vessels for their help. Support from the European Integrated Projects CARBOOCEAN (511176) and CARBOCHANGE (264879) is also acknowledged. We thank the reviewers Are Olsen and the anonymous reviewer for their helpful comments on this work. We would like to pay tribute to our late colleague Taro Takahashi, who contributed to initiate this sampling in 1993 and was a strong source of motivation for maintaining this long-time monitoring.

Financial support

This research has been supported by the European Integrated Projects CARBOOCEAN (511176) and CARBOCHANGE (264879).

Review statement

This paper was edited by Christoph Heinze and reviewed by Are Olsen and one anonymous referee.


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. 

Bates, N., Astor, Y., Church, M., Currie, K., Dore, J., Gonaález-Dávila, M., Lorenzoni, L., Muller-Karger, F., Olafsson, J. and Santa-Casiano, M.: A Time-Series View of Changing Ocean Chemistry Due to Ocean Uptake of Anthropogenic CO2 and Ocean Acidification, Oceanography, 27, 126–141,, 2014. 

Beaugrand, G., McQuatters-Gollop, A., Edwards, M., and Goberville, E.: Long-term responses of North Atlantic calcifying plankton to climate change, Nat. Clim. Change, 3, 263–267,, 2013. 

Becker, M., Steinhoff, T., and Körtzinger, A.: A Detailed View on the Seasonality of Stable Carbon Isotopes Across the North Atlantic, Global Biogeochem. Cy., 32, 1406–1419,, 2018. 

Breeden, M. L. and McKinley, G. A.: Climate impacts on multidecadal pCO2 variability in the North Atlantic: 1948–2009, Biogeosciences, 13, 3387–3396,, 2016. 

Brown, C. W. and Yoder, J. A.: Distribution pattern of coccolithophorid blooms in the western North Atlantic Ocean, Cont. Shelf Res., 14, 175–197,, 1994. 

Chafik, L., Rossby, T., and Schrum, C.: On the spatial structure and temporal variability of poleward transport between Scotland and Greenland, J. Geophys. Res.-Oceans, 119, 824–841,, 2014. 

Chafik, L., Nilsen, J. E. Ø., Dangendorf, S., Reverdin, G., and Frederikse, T.: North Atlantic Ocean Circulation and Decadal Sea Level Change During the Altimetry Era, Sci. Rep., 9, 1–9,, 2019. 

Clement, D. and Gruber, N.: The eMLR(C*) Method to Determine Decadal Changes in the Global Ocean Storage of Anthropogenic CO2, Global Biogeochem. Cy., 32, 654–679,, 2018. 

Corbière, A., Metzl, N., Reverdin, G., Brunet, C., and Takahashi, T.: Interannual and decadal variability of the oceanic carbon sink in the North Atlantic subpolar gyre, Tellus B, 59, 168–178,, 2007. 

De Boisséson, E., Thierry, V., Mercier, H., Caniaux, G., and Desbruyères, D.: Origin, formation and variability of the Subpolar Mode Water located over the Reykjanes Ridge, J. Geophys. Res.-Oceans, 117, C12005,, 2012. 

Denvil-Sommer, A., Gehlen, M., Vrac, M., and Mejia, C.: LSCE-FFNN-v1: a two-step neural network model for the reconstruction of surface ocean pCO2 over the global ocean, Geosci. Model Dev., 12, 2091–2105,, 2019. 

Desbruyères, D., Mercier, H., and Thierry, V.: On the mechanisms behind decadal heat content changes in the eastern subpolar gyre, Prog. Oceanogr., 132, 262–272,, 2015. 

Dickson, A. G.: Standard potential of the reaction: AgCl(s) + 12H2(g) = Ag(s) + HCl(aq), and the standard acidity constant of the ion HSO4 – in synthetic sea water from 273.15 to 318.15 K, J. Chem. Thermodyn., 22, 113–127,, 1990. 

Dickson, A. G. and Millero, F. J.: A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media, Deep-Sea Res. Pt. I, 34, 1733–1743,, 1987. 

Dlugokencky, E. J., Lang, P. M., Mund, J. W., Crotwell, M. J., and Thoning, K. W.: Atmospheric Carbon Dioxide Dry Air Mole Fractions from the NOAA ESRL Carbon Cycle Cooperative Global Air Sampling Network, 1968–2017, Version: 2018-07-31, available at:, last access: 6 March 2019. 

Doney, S., Bopp, L., and Long, M.: Historical and Future Trends in Ocean Climate and Biogeochemistry, Oceanog, 27, 108–119,, 2014. 

Doney, S. C., Fabry, V. J., Feely, R. A., and Kleypas, J. A.: Ocean Acidification: The Other CO2 Problem, Ann. Rev. Mar. Sci., 1, 169–192,, 2009. 

Fay, A. R. and McKinley, G. A.: Global trends in surface ocean pCO2 from in situ data, Global Biogeochem. Cy., 27, 541–557,, 2013. 

Fay, A. R. and McKinley, G. A.: Global open-ocean biomes: mean and temporal variability, Earth Syst. Sci. Data, 6, 273–284,, 2014. 

Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Hauck, J., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Bakker, D. C. E., Canadell, J. G., Ciais, P., Jackson, R. B., Anthoni, P., Barbero, L., Bastos, A., Bastrikov, V., Becker, M., Bopp, L., Buitenhuis, E., Chandra, N., Chevallier, F., Chini, L. P., Currie, K. I., Feely, R. A., Gehlen, M., Gilfillan, D., Gkritzalis, T., Goll, D. S., Gruber, N., Gutekunst, S., Harris, I., Haverd, V., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Joetzjer, E., Kaplan, J. O., Kato, E., Klein Goldewijk, K., Korsbakken, J. I., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Marland, G., McGuire, P. C., Melton, J. R., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Neill, C., Omar, A. M., Ono, T., Peregon, A., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Séférian, R., Schwinger, J., Smith, N., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Werf, G. R., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2019, Earth Syst. Sci. Data, 11, 1783–1838,, 2019. 

Friis, K., Körtzinger, A., and Wallace, D. W. R.: The salinity normalization of marine inorganic carbon chemistry data, Geophys. Res. Lett., 30, 1085,, 2003. 

Fröb, F., Olsen, A., Våge, K., Moore, G. W. K., Yashayaev, I., Jeansson, E., and Rajasakaren, B.: Irminger Sea deep convection injects oxygen and anthropogenic carbon to the ocean interior, Nat. Commun., 7, 1–8,, 2016. 

Fröb, F., Olsen, A., Pérez, F. F., García-Ibáñez, M. I., Jeansson, E., Omar, A., and Lauvset, S. K.: Inorganic carbon and water masses in the Irminger Sea since 1991, Biogeosciences, 15, 51–72,, 2018. 

Fröb, F., Olsen, A., Becker, M., Chafik, L., Johannessen, T., Reverdin, G. and Omar, A.: Wintertime fCO2 Variability in the Subpolar North Atlantic Since 2004, Geophys. Res. Lett., 46, 1580–1590,, 2019. 

García-Ibáñez, M. I., Zunino, P., Fröb, F., Carracedo, L. I., Ríos, A. F., Mercier, H., Olsen, A., and Pérez, F. F.: Ocean acidification in the subpolar North Atlantic: rates and mechanisms controlling pH changes, Biogeosciences, 13, 3701–3715,, 2016. 

Goris, N., Tjiputra, J. F., Olsen, A., Schwinger, J., Lauvset, S. K. and Jeansson, E.: Constraining projection-based estimates of the future North Atlantic carbon uptake, J. Climate, 31, 3959–3978,,2018. 

Gregor, L., Lebehot, A. D., Kok, S., and Scheel Monteiro, P. M.: A comparative assessment of the uncertainties of global surface ocean CO2 estimates using a machine-learning ensemble (CSIR-ML6 version 2019a) – have we hit the wall?, Geosci. Model Dev., 12, 5113–5136,, 2019. 

Gruber, N., Clement, D., Carter, B. R., Feely, R. A., Heuven, S. van, Hoppema, M., Ishii, M., Key, R. M., Kozyr, A., Lauvset, S. K., Lo Monaco, C., Mathis, J. T., Murata, A., Olsen, A., Perez, F. F., Sabine, C. L., Tanhua, T. and Wanninkhof, R.: The oceanic sink for anthropogenic CO2 from 1994 to 2007, Science, 363, 1193–1199,, 2019a. 

Gruber, N., Clement, D., Carter, B. R., Feely, R. A., Heuven, S. van, Hoppema, M., Ishii, M., Key, R. M., Kozyr, A., Lauvset, S. K., Lo Monaco, C., Mathis, J. T., Murata, A., Olsen, A., Perez, F. F., Sabine, C. L., Tanhua, T., and Wanninkhof, R.: The oceanic sink for anthropogenic CO2 from 1994 to 2007 – the data (NCEI Accession 0186034), NOAA National Centers for Environmental Information. Dataset,, 2019b. 

Häkkinen, S. and Rhines, P. B.: Decline of Subpolar North Atlantic Circulation During the 1990s, Science, 304, 555–559,, 2004. 

Häkkinen, S. and Rhines, P. B.: Shifting surface currents in the northern North Atlantic Ocean, J. Geophys. Res.-Oceans, 114, C04005,, 2009. 

Häkkinen, S., Rhines, P. B., and Worthen, D. L.: Atmospheric Blocking and Atlantic Multidecadal Ocean Variability, Science, 334, 655–659,, 2011. 

Häkkinen, S., Rhines, P. B., and Worthen, D. L.: Northern North Atlantic sea surface height and ocean heat content variability, J. Geophys. Res.-Oceans, 118, 3670–3678,, 2013. 

Hartmann, D. L., Klein Tank, A. M. G., Rusticucci, M., Alexander, L. V., Brönnimann, S., Charabi, Y., Dentener, F. J., Dlugokencky, E. J., Easterling, D. R., Kaplan, A., Soden, B. J., Thorne, P. W., Wild, M., and Zhai, P. M.: Observations: Atmosphere and Surface. In: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 159–254,, 2013. 

Hátún, H., Azetsu-Scott, K., Somavilla, R., Rey, F., Johnson, C., Mathis, M., Mikolajewicz, U., Coupel, P., Tremblay, J.-É., Hartman, S., Pacariz, S. V., Salter, I., and Ólafsson, J.: The subpolar gyre regulates silicate concentrations in the North Atlantic, Sci. Rep., 7, 1–9,, 2017. 

Henson, S. A., Painter, S. C., Holliday, N. P., Stinchcombe, M. C., and Giering, S. L. C.: Unusual subpolar NorthAtlantic phytoplankton bloom in 2010: Volcanic fertilization or North Atlantic Oscillation?, J. Geophys. Res.-Oceans, 118, 4771–4780,, 2013. 

van Heuven, S., Pierrot, D., Rae, J. W. B., Lewis, E., and Wallace, D. W. R.: MATLAB Program Developed for CO2 System Calculations, ORNL/CDIAC-105b, 2011. 

Holliday, N. P., Bersch, M., Berx, B., Chafik, L., Cunningham, S.,Florindo-Lopez, C., Hátún, H., Johns, W., Josey, S., A., Larsen, K., M., H., Mulet, S., Oltmanns, M., Reverdin, G., Rossby, T., Thierry, V., Valdimarsson, H.: Ocean circulation causes the largest freshening event for 120 years in eastern subpolar North Atlantic, Nat. Commun., 11, 585,, 2020. 

Jiang, L.-Q., Carter, B. R., Feely, R. A., Lauvset, S. K., and Olsen, A.: Surface ocean pH and buffer capacity: past, present and future, Sci. Rep., 9, 1–11,, 2019. 

de Jong, M. F. and de Steur, L.: Strong winter cooling over the Irminger Sea in winter 2014–2015, exceptional deep convection, and the emergence of anomalously low SST, Geophys. Res. Lett., 43, 7106–7113,, 2016. 

Josey, S. A., Hirschi, J. J.-M., Sinha, B., Duchez, A., Grist, J. P., and Marsh, R.: The Recent Atlantic Cold Anomaly: Causes, Consequences, and Related Phenomena, Ann. Rev. Mar. Sci., 10, 475–501,, 2018. 

Keeling, C. D., Brix, H., and Gruber, N.: Seasonal and long-term dynamics of the upper ocean carbon cycle at Station ALOHA near Hawaii, Global Biogeochem. Cy., 18, GB4006,, 2004. 

Keller, K., Joos, F., Raible, C., Cocco, V., Frolicher, T., Dunne, J., Gehlen, M., Bopp, L., Orr, J., Tjiputra, J., Heinze, C., Segscheider, J., Roy, T., and Metzl, N.: Variability of the Ocean Carbon Cycle in Response to the North Atlantic Oscillation, Tellus B, 64, 18738,, 2012. 

Key, R. M., Olsen, A., van Heuven, S., Lauvset, S. K., Velo, A., Lin, X., Schirnick, C., Kozyr, A., Tanhua, T., Hoppema, M., Jutterström, S., Steinfeldt, R., Jeansson, E., Ishii, M., Perez, F. F., and Suzuki, T.: Global Ocean Data Analysis Project, Version 2 (GLODAPv2), ORNL/CDIAC-162, NDP-093,, NDP093_GLODAPv2, 2015. 

Khatiwala, S., Tanhua, T., Mikaloff Fletcher, S., Gerber, M., Doney, S. C., Graven, H. D., Gruber, N., McKinley, G. A., Murata, A., Ríos, A. F., and Sabine, C. L.: Global ocean storage of anthropogenic carbon, Biogeosciences, 10, 2169–2191,, 2013. 

Kwiatkowski, L. and Orr, J. C.: Diverging seasonal extremes for ocean acidification during the twenty-first century, Nat. Clim. Change, 8, 141–145,, 2018. 

Landschützer, P., Gruber, N., Bakker, D. C. E., Schuster, U., Nakaoka, S., Payne, M. R., Sasse, T. P., and Zeng, J.: A neural network-based estimate of the seasonal to inter-annual variability of the Atlantic Ocean carbon sink, Biogeosciences, 10, 7793–7815,, 2013. 

Landschützer, P., Gruber, N., Bakker, D. C. E., Stemmler, I. and Six, K. D.: Strengthening seasonal marine CO2 variations due to increasing atmospheric CO2, Nat. Clim. Change, 8, 146–150,, 2018. 

Landschützer, P., Ilyina, T., and Lovenduski, N. S.: Detecting Regional Modes of Variability in Observation-Based Surface Ocean pCO2, Geophys. Res. Lett., 46, 2670–2679,, 2019. 

Lauvset, S. K. and Gruber, N.: Long-term trends in surface ocean pH in the North Atlantic, 0304-4203,, 2014. 

Lauvset, S. K., Gruber, N., Landschützer, P., Olsen, A., and Tjiputra, J.: Trends and drivers in global surface ocean pH over the past 3 decades, Biogeosciences, 12, 1285–1298,, 2015. 

Le Quéré, C., Andrew, R. M., Friedlingstein, P., Sitch, S., Hauck, J., Pongratz, J., Pickers, P. A., Korsbakken, J. I., Peters, G. P., Canadell, J. G., Arneth, A., Arora, V. K., Barbero, L., Bastos, A., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Doney, S. C., Gkritzalis, T., Goll, D. S., Harris, I., Haverd, V., Hoffman, F. M., Hoppema, M., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Johannessen, T., Jones, C. D., Kato, E., Keeling, R. F., Goldewijk, K. K., Landschützer, P., Lefèvre, N., Lienert, S., Liu, Z., Lombardozzi, D., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., Neill, C., Olsen, A., Ono, T., Patra, P., Peregon, A., Peters, W., Peylin, P., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rocher, M., Rödenbeck, C., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Steinhoff, T., Sutton, A., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., Wright, R., Zaehle, S., and Zheng, B.: Global Carbon Budget 2018, Earth Syst. Sci. Data, 10, 2141–2194,, 2018. 

Lebehot, A. D., Halloran, P. R., Watson, A. J., McNeall, D., Ford, D. A., Landschützer, P., Lauvset, S. K., and Schuster, U.: Reconciling Observation and Model Trends in North Atlantic Surface CO2, Global Biogeochem. Cy., 33, 1204–1222,, 2019. 

Lee, K., Tong, L. T., Millero, F. J., Sabine, C. L., Dickson, A. G., Goyet, C., Park, G.-H., Wanninkhof, R., Feely, R. A., and Key, R. M.: Global relationships of total alkalinity with salinity and temperature in surface waters of the world's oceans, Geophys. Res. Lett., 33, L19605,, 2006. 

Lewis, E., Wallace, D., and Allison, L. J.: Program developed for CO2 system calculations, Brookhaven National Lab., Dept. of Applied Science, Upton, NY (United States); Oak Ridge National Lab., Carbon Dioxide Information Analysis Center, TN (United States), 1998. 

Loveday, B. R. and Smyth, T.: A 40-year global data set of visible-channel remote-sensing reflectances and coccolithophore bloom occurrence derived from the Advanced Very High Resolution Radiometer catalogue, Earth Syst. Sci. Data, 10, 2043–2054,, 2018. 

Macovei, V. A., Hartman, S. E., Schuster, U., Torres-Valdés, S., Moore, C. M., and Sanders, R. J.: Impact of physical and biological processes on temporal variations of the ocean carbon sink in the mid-latitude North Atlantic (2002–2016), Prog. Oceanogr., 180, 102223,, 2019. 

Martinez, E., Raitsos, D. E., and Antoine, D.: Warmer, deeper, and greener mixed layers in the North Atlantic subpolar gyre over the last 50 years, Glob. Chang. Biol. 22, 604–612,, 2016. 

McKinley, G. A., Fay, A. R., Takahashi, T., and Metzl, N.: Convergence of atmospheric and North Atlantic carbon dioxide trends on multidecadal timescales, Nat. Geosci., 4, 606–610,, 2011. 

McKinley, G. A., Ritzer, A. L., and Lovenduski, N. S.: Mechanisms of northern North Atlantic biomass variability, Biogeosciences, 15, 6049–6066,, 2018. 

Mehrbach, C., Culberson, C. H., Hawley, J. E., and Pytkowicx, R. M.: Measurement of the Apparent Dissociation Constants of Carbonic Acid in Seawater at Atmospheric Pressure1, Limnol. Oceanogr., 18, 897–907,, 1973. 

Metzl, N., Corbière, A., Reverdin, G., Lenton, A., Takahashi, T., Olsen, A., Johannessen, T., Pierrot, D., Wanninkhof, R., Ólafsdóttir, S. R., Olafsson, J., and Ramonet, M.: Recent acceleration of the sea surface fCO2 growth rate in the North Atlantic subpolar gyre (1993–2008) revealed by winter observations, Global Biogeochem. Cy., 24, GB4004,, 2010. 

Metzl, N., Ferron, B., Lherminier, P., Sarthou, G., and Thierry, V.: Discrete profile measurements of dissolved inorganic carbon (DIC), total alkalinity (TALK), temperature and salinity during the multiple ships Observatoire de la variabilité interannuelle et décennale en Atlantique Nord (OVIDE) project, OVIDE-2006, OVIDE-2008, OVIDE-2010, OVIDE-2012, OVIDE-2014 cruises in the North Atlantic Ocean from 2006-05-23 to 2014-06-30 (NCEI Accession 0177219), Version 1.1, NOAA National Centers for Environmental Information, Dataset,, 2018. 

Millero, F. J., Lee, K., and Roche, M.: Distribution of alkalinity in the surface waters of the major oceans, Mar. Chem., 60, 111–130,, 1998. 

Millero, F. J., Feistel, R., Wright, D. G., and McDougall, T. J.: The composition of standar seawater and the definition of the reference-composition salinity scale, Deep-Sea Res. Pt. I, 55, 50–72,, 2008. 

Munro, D. R., Lovenduski, N. S., Takahashi, T., Stephens, B. B., Newberger, T., and Sweeney C.: Recent evidence for a strengthening CO2 sink in the Southern Ocean from carbonate system measurements in the Drake Passage (2002–2015), Geophys. Res. Lett., 42, 7623–7630,, 2015. 

Nigam, S., Ruiz-Barradas, A., and Chafik, L.: Gulf Stream Excursions and Sectional Detachments Generate the Decadal Pulses in the Atlantic Multidecadal Oscillation, J. Climate, 31, 2853–2870,, 2018. 

Olafsson, J., Olafsdottir, S. R., Benoit-Cattin, A., Danielsen, M., and Takahashi, T.: North Atlantic Ocean acidification from time series measurements, IOP Conf. Ser.: Earth Environ. Sci., 6, 462005,, 2009. 

Olafsson, J., Olafsdottir, S. R., Benoit-Cattin, A., and Takahashi, T.: The Irminger Sea and the Iceland Sea time series measurements of sea water carbon and nutrient chemistry 1983–2008, Earth Syst. Sci. Data, 2, 99–104,, 2010. 

Olsen, A., Key, R. M., van Heuven, S., Lauvset, S. K., Velo, A., Lin, X., Schirnick, C., Kozyr, A., Tanhua, T., Hoppema, M., Jutterström, S., Steinfeldt, R., Jeansson, E., Ishii, M., Pérez, F. F., and Suzuki, T.: The Global Ocean Data Analysis Project version 2 (GLODAPv2) – an internally consistent data product for the world ocean, Earth Syst. Sci. Data, 8, 297–323,, 2016. 

Olsen, A., Lange, N., Key, R. M., Tanhua, T., Álvarez, M., Becker, S., Bittig, H. C., Carter, B. R., Cotrim da Cunha, L., Feely, R. A., van Heuven, S., Hoppema, M., Ishii, M., Jeansson, E., Jones, S. D., Jutterström, S., Karlsen, M. K., Kozyr, A., Lauvset, S. K., Lo Monaco, C., Murata, A., Pérez, F. F., Pfeil, B., Schirnick, C., Steinfeldt, R., Suzuki, T., Telszewski, M., Tilbrook, B., Velo, A., and Wanninkhof, R.: GLODAPv2.2019 – an update of GLODAPv2, Earth Syst. Sci. Data, 11, 1437–1461,, 2019. 

Omar, A. M., Thomas, H., Olsen, A., Becker, M., Skjelvan, I., and Reverdin, G.: Trends of Ocean Acidification and pCO2 in the Northern North Sea, 2003–2015, J. Geophys. Res.-Biogeo., 124, 3088–3103,, 2019. 

Ono, H., Kosugi, N., Toyama, K., Tsujino, H., Kojima, A., Enyo, K., Iida, Y., Nakano, T., and Ishii, M.: Acceleration of Ocean Acidification in the Western North Pacific, Geophys. Res. Lett., 46, 13161–13169,, 2019. 

Organization (WMO), W. M., (UNESCO) United Nations Educational, S. and C. O., Commission, (IOC) Intergovernmental Oceanographic, Programme, (UNEP) United Nations Environment, Science, (ICSU) International Council for, Climate (OOPC-21), 21st Session of the Ocean Observations Panel for and World Meteorological Organization (WMO): GCOS, 217, 21st Session of the Ocean Observations Panel for Climate (the GOOS Physics and Climate Panel) (OOPC-21), WMO, 2018. 

Orr, J. C., Fabry, V. J., Aumont, O., Bopp, L., Doney, S. C., Feely, R. A., Gnanadesikan, A., Gruber, N., Ishida, A., Joos, F., Key, R. M., Lindsay, K., Maier-Reimer, E., Matear, R., Monfray, P., Mouchet, A., Najjar, R. G., Plattner, G.-K., Rodgers, K. B., Sabine, C. L., Sarmiento, J. L., Schlitzer, R., Slater, R. D., Totterdell, I. J., Weirig, M.-F., Yamanaka, Y., and Yool, A.: Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms, Nature, 437, 681–686,, 2005. 

Orr, J. C., Epitalon, J.-M., Dickson, A. G., and Gattuso, J.-P.: Routine uncertainty propagation for the marine carbon dioxide system, Mar. Chem., 207, 84–107,, 2018. 

Pickart, R. S., Spall, M. A., Ribergaard, M. H., Moore, G. W. K., and Milliff, R. F.: Deep convection in the Irminger Sea forced by the Greenland tip jet, Nature, 424, 152–156,, 2003. 

Piron, A., Thierry, V., Mercier, H., and Caniaux, G.: Gyre-scale deep convection in the subpolar North Atlantic Ocean during winter 2014–2015, Geophys. Res. Lett., 44, 1439–1447,, 2017. 

Racapé, V., Pierre, C., Metzl, N., Lo Monaco, C., Reverdin, G., Olsen, A., Morin, P., Vázquez-Rodríguez, M., Ríos, A. F., and Pérez, F. F.: Anthropogenic carbon changes in the Irminger Basin (1981–2006): Coupling δ13CDIC and DIC observations, J. Mar. Syst., 126, 24–32,, 2013. 

Racapé, V., Metzl, N., Pierre, C., Reverdin, G., Quay, P. D., and Olafsdottir, S. R.: The seasonal cycle of δ3CDIC in the North Atlantic subpolar gyre, Biogeosciences, 11, 1683–1692,, 2014. 

Reverdin, G., Valdimarsson, H., Alory, G., Diverres, D., Bringas, F., Goni, G., Heilmann, L., Chafik, L., Szekely, T., and Friedman, A. R.: North Atlantic subpolar gyre along predetermined ship tracks since 1993: a monthly data set of surface temperature, salinity, and density, Earth Syst. Sci. Data, 10, 1403–1415,, 2018a. 

Reverdin, G., Metzl, N., Olafsdottir, S., Racapé, V., Takahashi, T., Benetti, M., Valdimarsson, H., Benoit-Cattin, A., Danielsen, M., Fin, J., Naamar, A., Pierrot, D., Sullivan, K., Bringas, F., and Goni, G.: SURATLANT: a 1993–2017 surface sampling in the central part of the North Atlantic subpolar gyre, Earth Syst. Sci. Data, 10, 1901–1924,, 2018b. 

Rivero-Calle, S., Gnanadesikan, A., Del Castillo, C. E., Balch, W. M., and Guikema, S. D.: Multidecadal increase in North Atlantic coccolithophores and the potential role of rising CO2, Science, 350, 1533–1537,, 2015. 

Robson, J., Ortega, P., and Sutton, R.: A reversal of climatic trends in the North Atlantic since 2005, Nat. Geosci., 9, 513–517,, 2016. 

Robson, J., Polo, I., Hodson, D. L. R., Stevens, D. P., and Shaffrey, L. C.: Decadal prediction of the North Atlantic subpolar gyre in the HiGEM high-resolution climate model, Clim. Dynam., 50, 921–937,, 2018. 

Rödenbeck, C., Bakker, D. C. E., Gruber, N., Iida, Y., Jacobson, A. R., Jones, S., Landschützer, P., Metzl, N., Nakaoka, S., Olsen, A., Park, G.-H., Peylin, P., Rodgers, K. B., Sasse, T. P., Schuster, U., Shutler, J. D., Valsala, V., Wanninkhof, R., and Zeng, J.: Data-based estimates of the ocean carbon sink variability – first results of the Surface Ocean pCO2 Mapping intercomparison (SOCOM), Biogeosciences, 12, 7251–7278,, 2015. 

Sabine, C. L., Feely, R. A., Gruber, N., Key, R. M., Lee, K., Bullister, J. L., Wanninkhof, R., Wong, C. S., Wallace, D. W. R., Tilbrook, B., Millero, F. J., Peng, T.-H., Kozyr, A., Ono, T. and Rios, A. F.: The Oceanic Sink for Anthropogenic CO2, Science, 305, 367–371,, 2004. 

Schuster, U., Watson, A. J., Bates, N. R., Corbiere, A., Gonzalez-Davila, M., Metzl, N., Pierrot, D., and Santana-Casiano, M.: Trends in North Atlantic sea-surface fCO2 from 1990 to 2006, Deep-Sea Res. Pt. II, 56, 620–629,, 2009. 

Schuster, U., McKinley, G. A., Bates, N., Chevallier, F., Doney, S. C., Fay, A. R., González-Dávila, M., Gruber, N., Jones, S., Krijnen, J., Landschützer, P., Lefèvre, N., Manizza, M., Mathis, J., Metzl, N., Olsen, A., Rios, A. F., Rödenbeck, C., Santana-Casiano, J. M., Takahashi, T., Wanninkhof, R., and Watson, A. J.: An assessment of the Atlantic and Arctic sea-air CO2 fluxes, 1990–2009, Biogeosciences, 10, 607–627,, 2013. 

Signorini, S. R., Häkkinen, S., Gudmundsson, K., Olsen, A., Omar, A. M., Olafsson, J., Reverdin, G., Henson, S. A., McClain, C. R., and Worthen, D. L.: The role of phytoplankton dynamics in the seasonal and interannual variability of carbon in the subpolar North Atlantic – a modeling study, Geosci. Model Dev., 5, 683–707,, 2012. 

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. 

Takahashi, T., Sutherland, S. C., Chipman, D. W., Goddard, J. G., Ho, C., Newberger, T., Sweeney, C., and Munro, D. R.: Climatological distributions of pH, pCO2, total CO2, alkalinity, and CaCO3 saturation in the global surface ocean, and temporal changes at selected locations, Mar. Chem., 164, 95–125,, 2014. 

Thomas, H., Prowe, A. E. F., Lima, I. D., Doney, S. C., Wanninkhof, R., Greatbatch, R. J., Schuster, U., and Corbière, A.: Changes in the North Atlantic Oscillation influence CO2 uptake in the North Atlantic over the past 2 decades, Global Biogeochem. Cy., 22, GB4027,, 2008. 

Tilbrook, B., Jewett, E. B., DeGrandpre, M. D., Hernandez-Ayon, J. M., Feely, R. A., Gledhill, D. K., Hansson, L., Isensee, K., Kurz, M. L., Newton, J. A., Siedlecki, S. A., Chai, F., Dupont, S., Graco, M., Calvo, E., Greeley, D., Kapsenberg, L., Lebrec, M., Pelejero, C., Schoo, K. L., and Telszewski, M.: An Enhanced Ocean Acidification Observing Network: From People to Technology to Data Synthesis and Information Exchange, Front. Mar. Sci., 6, 337,, 2019. 

Tilstone, G. H., Miller, P. I., Brewin, R. J. W., and Priede, I. G.: Enhancement of primary production in the North Atlantic outside of the spring bloom, identified by remote sensing of ocean colour and temperature, Remote Sens. Environ., 146, 77–86,, 2014. 

Tjiputra, J. F., Olsen, A., Assmann, K., Pfeil, B., and Heinze, C.: A model study of the seasonal and long-term North Atlantic surface pCO2 variability, Biogeosciences, 9, 907–923,, 2012. 

Tjiputra, J. F., Olsen, A., Bopp, L., Lenton, A., Pfeil, B., Roy, T., Segschneider, J., Totterdell, I., and Heinze, C.: Long-term surface pCO2 trends from observations and models, Tellus B, 66, 23083,, 2014. 

Touratier, F., Azouzi, L., and Goyet, C.: CFC-11, Δ14C and 3H tracers as a means to assess anthropogenic CO2 concentrations in the ocean, Tellus B, 59, 318–325,, 2007. 

Ullman, D. J., McKinley, G. A., Bennington, V., and Dutkiewicz, S.: Trends in the North Atlantic carbon sink: 1992–2006, Global Biogeochem. Cy., 23, GB4011,, 2009. 

Uppström, L. R.: The boron/chlorinity ratio of deep-sea water from the Pacific Ocean, Deep Sea Research and Oceanographic Abstracts, 21, 161–162,, 1974. 

Våge, K., Pickart, R. S., Moore, G. W. K., and Ribergaard, M. H.: Winter Mixed Layer Development in the Central Irminger Sea: The Effect of Strong, Intermittent Wind Events, J. Phys. Oceanogr., 38, 541–565,, 2008. 

Vázquez-Rodríguez, M., Pérez, F. F., Velo, A., Ríos, A. F., and Mercier, H.: Observed acidification trends in North Atlantic water masses, Biogeosciences, 9, 5217–5230,, 2012. 

Wakita, M., Watanabe, S., Honda, M., Nagano, A., Kimoto, K., Matsumoto, K., Kitamura, M., Sasaki, K., Kawakami, H., Fujiki, T., Sasaoka, K., Nakano, Y., and Murata, A.: Ocean acidification from 1997 to 2011 in the subarctic western North Pacific Ocean, Biogeosciences, 10, 7817–7827,, 2013. 

Wakita, M., Nagano, A., Fujiki, T., and Watanabe, S.: Slow acidification of the winter mixed layer in the subarctic western North Pacific, J. Geophys. Res.-Oceans, 122, 6923–6935,, 2017. 

Wanninkhof, R., Pickers, P. A., Omar, A. M., Sutton, A., Murata, A., Olsen, A., Stephens, B. B., Tilbrook, B., Munro, D., Pierrot, D., Rehder, G., Santana-Casiano, J. M., Müller, J. D., Trinanes, J., Tedesco, K., O'Brien, K., Currie, K., Barbero, L., Telszewski, M., Hoppema, M., Ishii, M., González-Dávila, M., Bates, N. R., Metzl, N., Suntharalingam, P., Feely, R. A., Nakaoka, S., Lauvset, S. K., Takahashi, T., Steinhoff, T., and Schuster, U.: A Surface Ocean CO2 Reference Network, SOCONET and Associated Marine Boundary Layer CO2 Measurements, Front. Mar. Sci., 6, 400,, 2019. 

Watson, A. J., Schuster, U., Bakker, D. C. E., Bates, N. R., Corbière, A., González-Dávila, M., Friedrich, T., Hauck, J., Heinze, C., Johannessen, T., Körtzinger, A., Metzl, N., Olafsson, J., Olsen, A., Oschlies, A., Padin, X. A., Pfeil, B., Santana-Casiano, J. M., Steinhoff, T., Telszewski, M., Rios, A. F., Wallace, D. W. R., and Wanninkhof, R.: Tracking the Variable North Atlantic Sink for Atmospheric CO2, Science, 326, 1391–1393,, 2009.  

Weiss, R. F. and Price, B. A.: Nitrous oxide solubility in water and seawater, Mar. Chem., 8, 347–359,, 1980. 

Williams, N. L., Juranek, L. W., Feely, R. A., Johnson, K. S., Sarmiento, J. L., Talley, L. D., Dickson, A. G., Gray, A. R., Wanninkhof, R., Russell, J. L., Riser, S. C., and Takeshita, Y.: Calculating surface ocean pCO2 from biogeochemical Argo floats equipped with pH: An uncertainty analysis, Global Biogeochem. Cy., 31, 591–604,, 2017. 

Short summary
In this study, we investigate the evolution of CO2 uptake and ocean acidification in the North Atlantic Subpolar surface water. Our results show an important reduction in the capacity of the ocean to absorb CO2 from the atmosphere (1993–2007), due to a rapid increase in the fCO2 and associated with a rapid decrease in pH. Conversely, data obtained during the last decade (2008–2017) show a stagnation of fCO2 (increasing the ocean sink for CO2) and pH.
Final-revised paper