Variability Of Air-sea CO2 Flux Research Articles

Temperature effect on seawater fCO2 revisited: theoretical basis, uncertainty analysis and implications for parameterising carbonic acid equilibrium constants

Abstract. The sensitivity of the fugacity of carbon dioxide in seawater (fCO2) to temperature (denoted υ, reported in % °C−1) is critical for the accurate fCO2 measurements needed to build global carbon budgets and for understanding the drivers of air–sea CO2 flux variability across the ocean. However, understanding and computing υ have been restricted to either using empirical functions fitted to experimental data or determining it as an emergent property of a fully resolved marine carbonate system, and these two approaches are not consistent with each other. The lack of a theoretical basis and an uncertainty estimate for υ has hindered resolving this discrepancy. Here, we develop a new approach for calculating the temperature sensitivity of fCO2 based on the equations governing the marine carbonate system and the van 't Hoff equation. This shows that, to first order, ln (fCO2) should be proportional to 1/tK (where tK is temperature in kelvin), rather than to temperature, as has previously been assumed. This new approach is, to first order, consistent with calculations from a fully resolved marine carbonate system, which we have incorporated into the PyCO2SYS software. Agreement with experimental data is less convincing but remains inconclusive due to the scarcity of direct measurements of υ, particularly above 25 °C. However, the new approach is consistent with field data, performing better than any other approach for adjusting fCO2 by up to 10 °C if spatiotemporal variability in its single fitted coefficient is accounted for. The uncertainty in υ arising from only measurement uncertainty in the main experimental dataset where υ has been directly measured is in the order of 0.04 % °C−1, which corresponds to a 0.04 % uncertainty in fCO2 adjusted by +1 °C. However, spatiotemporal variability in υ is several times greater than this, so the true uncertainty due to the temperature adjustment in fCO2 adjusted by +1 °C using the most widely used constant υ value is around 0.24 %. This can be reduced to around 0.06 % using the new approach proposed here, and this could be further reduced by more measurements. The spatiotemporal variability in υ arises mostly from the equilibrium constants for CO2 solubility and carbonic acid dissociation (K1∗ and K2∗), and its magnitude varies significantly depending on which parameterisation is used for K1∗ and K2∗. Seawater fCO2 can be measured accurately enough that additional experiments should be able to detect spatiotemporal variability in υ and distinguish between different parameterisations for K1∗ and K2∗. Because the most widely used constant υ was coincidentally measured from seawater with roughly global average υ, our results are unlikely to significantly affect global air–sea CO2 flux budgets, but they may have more important implications for regional budgets and studies that adjust by larger temperature differences.

Estimating marine carbon uptake in the northeast Pacific using a neural network approach

Abstract. The global ocean takes up nearly a quarter of anthropogenic CO2 emissions annually, but the variability in this uptake at regional scales remains poorly understood. Here we use a neural network approach to interpolate sparse observations, creating a monthly gridded seawater partial pressure of CO2 (pCO2) data product from January 1998 to December 2019, at 1/12∘ × 1/12∘ spatial resolution, in the northeast Pacific open ocean, a net sink region. The data product (ANN-NEP; NCEI Accession 0277836) was created from pCO2 observations within the 2021 version of the Surface Ocean CO2 Atlas (SOCAT) and a range of predictor variables acting as proxies for processes affecting pCO2 to create nonlinear relationships to interpolate observations at a spatial resolution 4 times greater than leading global products and with better overall performance. In moving to a higher resolution, we show that the internal division of training data is the most important parameter for reducing overfitting. Using our pCO2 product, wind speed, and atmospheric CO2, we evaluate air–sea CO2 flux variability. On sub-decadal to decadal timescales, we find that the upwelling strength of the subpolar Alaskan Gyre, driven by large-scale atmospheric forcing, acts as the primary control on air–sea CO2 flux variability (r2=0.93, p<0.01). In the northern part of our study region, divergence from atmospheric CO2 is enhanced by increased local wind stress curl, enhancing upwelling and entrainment of naturally CO2-rich subsurface waters, leading to decade-long intervals of strong winter outgassing. During recent Pacific marine heat waves from 2013 on, we find enhanced atmospheric CO2 uptake (by as much as 45 %) due to limited wintertime entrainment. Our product estimates long-term surface ocean pCO2 increase at a rate below the atmospheric trend (1.4 ± 0.1 µatm yr−1) with the slowest increase in the center of the subpolar gyre where there is strong interaction with subsurface waters. This mismatch suggests the northeast Pacific Ocean sink for atmospheric CO2 may be increasing.

High interannual surface pCO2 variability in the southern Canadian Arctic Archipelago's Kitikmeot Sea

Abstract. Warming of the Arctic due to climate change means the Arctic Ocean is now free from ice for longer, as sea ice melts earlier and refreezes later. Yet, it remains unclear how this extended ice-free period will impact carbon dioxide (CO2) fluxes due to scarcity of surface ocean CO2 measurements. Baseline measurements are urgently needed to understand spatial and temporal air–sea CO2 flux variability in the changing Arctic Ocean. There is also uncertainty as to whether the previous basin-wide surveys are representative of the many smaller bays and inlets that make up the Canadian Arctic Archipelago (CAA). By using a research vessel that is based in the remote Inuit community of Ikaluqtuutiak (Cambridge Bay, Nunavut), we have been able to reliably survey pCO2 shortly after ice melt and access previously unsampled bays and inlets in the nearby region. Here we present 4 years of consecutive summertime pCO2 measurements collected in the Kitikmeot Sea in the southern CAA. Overall, we found that this region is a sink for atmospheric CO2 in August (average of all calculated fluxes over the four cruises was −4.64 mmol m−2 d−1), but the magnitude of this sink varies substantially between years and locations (average calculated fluxes of +3.58, −2.96, −16.79 and −0.57 mmol m−2 d−1 during the 2016, 2017, 2018 and 2019 cruises, respectively). Surface ocean pCO2 varied by up to 156 µatm between years, highlighting the importance of repeat observations in this region, as this high interannual variability would not have been captured by sparse and infrequent measurements. We find that the surface ocean pCO2 value at the time of ice melt is extremely important in constraining the magnitude of the air–sea CO2 flux throughout the ice-free season. However, further constraining the air–sea CO2 flux in the Kitikmeot Sea will require a better understanding of how pCO2 changes outside of the summer season. Surface ocean pCO2 measurements made in small bays and inlets of the Kitikmeot Sea were ∼ 20–40 µatm lower than in the main channels. Surface ocean pCO2 measurements made close in time to ice breakup (i.e. within 2 weeks) were ∼ 50 µatm lower than measurements made > 4 weeks after breakup. As previous basin-wide surveys of the CAA have focused on the deep shipping channels and rarely measure close to the ice breakup date, we hypothesize that there may be an observational bias in previous studies, leading to an underestimate of the CO2 sink in the CAA. These high-resolution measurements constitute an important new baseline for gaining a better understanding of the role this region plays in the uptake of atmospheric CO2.

Tidal mixing of estuarine and coastal waters in the western English Channel is a control on spatial and temporal variability in seawater CO<sub>2</sub>

Abstract. Surface ocean carbon dioxide (CO2) measurements are used to compute the oceanic air–sea CO2 flux. The CO2 flux component from rivers and estuaries is uncertain due to the high spatial and seasonal heterogeneity of CO2 in coastal waters. Existing high-quality CO2 instrumentation predominantly utilises showerhead and percolating style equilibrators optimised for open-ocean observations. The intervals between measurements made with such instrumentation make it difficult to resolve the fine-scale spatial variability of surface water CO2 at timescales relevant to the high frequency variability in estuarine and coastal environments. Here we present a novel dataset with unprecedented frequency and spatial resolution transects made at the Western Channel Observatory in the south-west of the UK from June to September 2016, using a fast-response seawater CO2 system. Novel observations were made along the estuarine–coastal continuum at different stages of the tide and reveal distinct spatial patterns in the surface water CO2 fugacity (fCO2) at different stages of the tidal cycle. Changes in salinity and fCO2 were closely correlated at all stages of the tidal cycle and suggest that the mixing of oceanic and riverine endmembers partially determines the variations in fCO2. The correlation between salinity and fCO2 was different in Cawsand Bay, which could be due to enhanced gas exchange or to enhanced biological activity in the region. The observations demonstrate the complex dynamics determining spatial and temporal patterns of salinity and fCO2 in the region. Spatial variations in observed surface salinity were used to validate the output of a regional high-resolution hydrodynamic model. The model enables a novel estimate of the air–sea CO2 flux in the estuarine–coastal zone. Air–sea CO2 flux variability in the estuarine–coastal boundary region is influenced by the state of the tide because of strong CO2 outgassing from the river plume. The observations and model output demonstrate that undersampling the complex tidal and mixing processes characteristic of estuarine and coastal environment biases quantification of air–sea CO2 fluxes in coastal waters. The results provide a mechanism to support critical national and regional policy implementation by reducing uncertainty in carbon budgets.

Importance of El Niño reproducibility for reconstructing historical CO<sub>2</sub> flux variations in the equatorial Pacific

Abstract. Based on a set of climate simulations utilizing two kinds of Earth system models (ESMs) in which observed ocean hydrographic data are assimilated using exactly the same data assimilation procedure, we have clarified that the successful simulation of the observed air–sea CO2 flux variations in the equatorial Pacific is tightly linked to the reproducibility of coupled physical air–sea processes. When an ESM with a weaker ENSO (El Niño–Southern Oscillations) amplitude than that of the observations was used for historical simulations with ocean data assimilation, the observed equatorial anticorrelated relationship between the sea surface temperature (SST) and the air–sea CO2 flux on interannual to decadal timescales could not be represented. The simulated CO2 flux anomalies were upward (downward) during El Niño (La Niña) periods in the equatorial Pacific. The reason for this was that the non-negligible correction term in the governing equation of ocean temperature, which was added via the ocean data assimilation procedure, caused an anomalous, spurious equatorial upwelling (downwelling) during El Niño (La Niña) periods, which brought more (less) subsurface layer water rich in dissolved inorganic carbon (DIC) to the surface layer. On the other hand, in the historical simulations where the observational data were assimilated into the other ESM with a more realistic ENSO representation, the correction term associated with the assimilation procedure remained small enough so as not to disturb an anomalous advection–diffusion balance for the equatorial ocean temperature. Consequently, spurious vertical transport of DIC and the resultant positively correlated SST and air–sea CO2 flux variations did not occur. Thus, the reproducibility of the tropical air–sea CO2 flux variability with data assimilation can be significantly attributed to the reproducibility of ENSO in an ESM. Our results suggest that, when using data assimilation to initialize ESMs for carbon cycle predictions, the reproducibility of the internal climate variations in the model itself is of great importance.

Effects of phytoplankton community composition and productivity on sea surface pCO2 variations in the Southern Ocean

The Southern Ocean is a vast net sink for atmospheric carbon dioxide (CO2), with marine phytoplankton playing a crucial role in CO2 fixation. We assessed how changes in the dominant phytoplankton community and net primary productivity (NPP) affected variations in the partial pressure of CO2 in surface water (pCO2sw) in the Indian sector of the Southern Ocean during austral summer. pCO2sw was negatively correlated with total phytoplankton and diatom abundances, as estimated from pigment signatures, in the zone south of the Antarctic Circumpolar Current; however, pCO2sw was not correlated with haptophyte abundance. Additionally, a stronger correlation was found between pCO2sw and total phytoplankton NPP than between chlorophyll a concentration and pCO2sw. We reconstructed pCO2sw at inter-annual scale using satellite data and assessed the inter-annual variability of air-sea CO2 flux. Over the period from 1997 to 2007, the integrated CO2 fluxes over the study region showed very large variations from a small source to a strong sink. Variations in the integrated CO2 fluxes were also correlated with changes in satellite-derived phytoplankton community in the Indian sector of the Southern Ocean and changes in the dominant phytoplankton community may control CO2 dynamics in the marginal ice zone.

On the role of climate modes in modulating the air–sea CO<sub>2</sub> fluxes in eastern boundary upwelling systems

Abstract. The air–sea CO2 fluxes in eastern boundary upwelling systems (EBUSs) vary strongly in time and space, with some of the highest flux densities globally. The processes controlling this variability have not yet been investigated consistently across all four major EBUSs, i.e., the California (CalCS), Humboldt (HumCS), Canary (CanCS), and Benguela (BenCS) Current systems. In this study, we diagnose the climatic modes of the air–sea CO2 flux variability in these regions between 1920 and 2015, using simulation results from the Community Earth System Model Large Ensemble (CESM-LENS), a global coupled climate model ensemble that is forced by historical and RCP8.5 radiative forcing. Differences between simulations can be attributed entirely to internal (unforced) climate variability, whose contribution can be diagnosed by subtracting the ensemble mean from each simulation. We find that in the CalCS and CanCS, the resulting anomalous CO2 fluxes are strongly affected by large-scale extratropical modes of variability, i.e., the North Pacific Gyre Oscillation (NPGO) and the North Atlantic Oscillation (NAO), respectively. The CalCS has anomalous uptake of CO2 during the positive phase of the NPGO, while the CanCS has anomalous outgassing of CO2 during the positive phase of the NAO. In contrast, the HumCS is mainly affected by El Niño–Southern Oscillation (ENSO), with anomalous uptake of CO2 during an El Niño event. Variations in dissolved inorganic carbon (DIC) and sea surface temperature (SST) are the major contributors to these anomalous CO2 fluxes and are generally driven by changes to large-scale gyre circulation, upwelling, the mixed layer depth, and biological processes. A better understanding of the sensitivity of EBUS CO2 fluxes to modes of climate variability is key in improving our ability to predict the future evolution of the atmospheric CO2 source and sink characteristics of the four EBUSs.

The Spatial and Temporal Variability of Air-Sea CO2 Fluxes and the Effect of Net Coral Reef Calcification in the Indonesian Seas: A Numerical Sensitivity Study

A numerical model system was developed and applied to simulate air-sea fluxes of CO2 and coral reef calcification in the Indonesian Seas and adjacent ocean basin for the period 1960 to 2014 on a fine resolution grid (ca. 11 km) in order to study their response to rising sea water temperatures and CO2 concentrations in the atmosphere. Results were analyzed for different sub-regions on the Sunda Shelf (Gulf of Thailand, Malacca Strait, Java Sea) and show realistic and different levels, signs and pronounced temporal variability in air-sea CO2 flux. The Gulf of Thailand changes from an atmospheric CO2 sink during the boreal winter to a CO2 source in summer due to higher water temperatures, while other sub-regions as well as the entire averaged Sunda Shelf act as a continuous source of CO2 for the atmosphere. However, increasing atmospheric CO2 concentrations weakened this source function during the simulation period. In 2007, the model simulations showed even a first flux inversion, in course of which the Java Sea took up CO2 . The simulated trends suggest that the entire Sunda Shelf will turn into a permanent sink for atmospheric CO2 within the next 30 to 35 years if current trends remain constant. Considering the period between 2010 and 2014, coral reef calcification enhanced the average CO2 emission of the Sunda Shelf by more than 10 % from 15 to 17 Tg C yr−1 due to lowering the pH and increasing the partial pressure of CO2 in surface water. During the entire period of simulation, net reef calcification decreased although increasing seawater temperature mitigated effects of reduced CO2 emission and the resulting decrease of the pH values on reef calcification. Our realistic simulation results already without consideration of any biological processes suggest that biological processes taking up and releasing CO 2 are currently well balanced in these tropical regions. However, the counteracting effects of climate change on the reef calcification, on other biological processes and the carbonate system need to be investigated in more detail. SST ◦increased by about 0.6 C during the last 55 years, while SSS decreased by about 0.7 psu.

Temporal and spatial variability of the carbon cycle in the east of China’s seas: a three-dimensional physical-biogeochemical modeling study

In the east of China’s seas, there is a wide range of the continental shelf. The nutrient cycle and the carbon cycle in the east of China’s seas exhibit a strong variability on seasonal to decadal time scales. On the basis of a regional ocean modeling system (ROMS), a three dimensional physical-biogeochemical model including the carbon cycle with the resolution (1/12)°×(1/12)° is established to investigate the physical variations, ecosystem responses and carbon cycle consequences in the east of China’s seas. The ROMS-Nutrient Phytoplankton Zooplankton Detritus (NPZD) model is driven by daily air-sea fluxes (wind stress, long wave radiation, short wave radiation, sensible heat and latent heat, freshwater fluxes) that derived from the National Centers for Environmental Prediction (NCEP) reanalysis2 from 1982 to 2005. The coupled model is capable of reproducing the observed seasonal variation characteristics over the same period in the East China Sea. The integrated air-sea CO2 flux over the entire east of China’s seas reveals a strong seasonal cycle, functioning as a source of CO2 to the atmosphere from June to October, while serving as a sink of CO2 to the atmosphere in the other months. The 24 a mean value of airsea CO2 flux over the entire east of China’s seas is about 1.06 mol/(m2·a), which is equivalent to a regional total of 3.22 Mt/a, indicating that in the east of China’s seas there is a sink of CO2 to the atmosphere. The partial pressure of carbon dioxide in sea water in the east of China’s seas has an increasing rate of 1.15 μatm/a (1μtm/a=0.101 325 Pa), but pH in sea water has an opposite tendency, which decreases with a rate of 0.001 3 a–1 from 1982 to 2005. Biological activity is a dominant factor that controls the $${p_{C{O_2}air}}$$ in the east of China’s seas, and followed by a temperature. The inverse relationship between the interannual variability of air-sea CO2 flux averaged from the domain area and Nino3 SST Index indicates that the carbon cycle in the east of China’s seas has a high correlation with El Nino-Southern Oscillation (ENSO).

Global Air–Sea CO2 Flux in 22 CMIP5 Models: Multiyear Mean and Interannual Variability*

Abstract To assess the capability of the latest Earth system models (ESMs) in representing historical global air–sea CO2 flux, 22 models from phase 5 of the Coupled Model Intercomparision Project (CMIP5) are analyzed, with a focus on the spatial distribution of multiyear mean and interannual variability. Results show that the global distribution of air–sea CO2 flux is reasonable in most of the models and that the main differences between models and observationally based results exist in regions with strong vertical movement. The annual mean flux in the 18-member multimodel ensemble (MME; four models were excluded because of their poor performances) mean during 1996–2004 is 1.95 Pg C yr−1 (1 Pg = 1015 g; positive values mean into the ocean), and all but one model describe the rapid increasing trend of air–sea CO2 flux observed during 1960–2000. The first mode of the global air–sea CO2 flux variability during 1870–2000 in six of the models represents the El Niño–Southern Oscillation (ENSO) mode. The remaining 12 models fail to represent this important character for the following reasons: in five models, the tropical Pacific does not play a dominant role in the interannual variability of global air–sea CO2 flux because of stronger interannual variability in the Southern Ocean; two models poorly represent the interannual fluctuation of dissolved inorganic carbon (DIC) in the surface ocean of the tropical Pacific; and four models have shorter periods of the air–sea CO2 flux, which are out of the period range of ENSO events.

Modeling Spatial and Temporal pCO2 Variability at 49°N/16.5°W and 56.5°N/52.6°W in the North Atlantic Ocean

Numerical models have been employed in understanding and capturing spatiotemporal CO2 trends, inter seasonal to decadal variability, and the characterization of thermal (pCO2-T) and non-thermal (pCO2-nonT) components of surface ocean pCO2 and air-sea CO2 flux variability. We employed the MIT Ocean General Circulation Model and available data from two sub polar North Atlantic observatories at 49°N, 16.5°W and 56.5°N, 52.6°W to capture in situ pCO2 observations and deconvolute bio-physical controlling processes. The model suggests that the pCCO2 cycle is marked by a summertime minimum and wintertime maximum. The physical-chemical and biological response pattern of the model is in good accordance with the observed pCO2, pCO2-T and pCO2-nonT trends. Model outputs suggest that the CO2 cycle is governed by contrasting effects of seasonal cooling and warming on one side and spring and summer net biology activities on the other side. It also predicts year-round under saturation, indicating that the region is a weak to moderate sink for CO2.

Climate impacts on the structures of the North Pacific air-sea CO<sub>2</sub> flux variability

Abstract. Some dominant spatial and temporal structures of the North Pacific air-sea CO2 fluxes in response to the Pacific Decadal Oscillation (PDO) are identified in three data products from three independent sources: an assimilated CO2 flux product and two forward model solutions. The interannual variability of CO2 flux is found to be an order of magnitude weaker compared to the seasonal cycle of CO2 flux in the North Pacific. A statistical approach is employed to quantify the signal-to-noise ratio in the reconstructed dataset to delineate the representativity errors. The dominant variability with a signal-to-noise ratio above one is identified and its correlations with PDO are examined. A tentative four-pole pattern in the North Pacific air-sea CO2 flux variability linked to PDO emerges in which two positively correlated poles are oriented in the northwest and southeast directions and contrarily, the negatively correlated poles are oriented in the northeast and southwest directions. This pattern is identified in three products, providing CO2 and pCO2. Its relations to the interannual El Niño-Southern Oscillation (ENSO) and lower-frequency PDO are separately identified. A combined EOF analysis between air-sea CO2 flux and key variables representing ocean-atmosphere interactions is carried out to elicit robust oscillations in the North Pacific CO2 flux in response to the PDO. The proposed spatial and temporal structures of the North Pacific CO2 fluxes are insightful since they separate the secular trends of the surface ocean carbon from the interannual variability. The regional characterization of the North Pacific in terms of PDO and CO2 flux variability is also instructive for determining the homogeneous oceanic domains for the Regional Carbon Cycle and Assessment Processes (RECCAP).

Interannual variability of net community production and air-sea CO2 flux in a naturally iron fertilized region of the Southern Ocean (Kerguelen Plateau)

AbstractThe interannual variability of net community production (NCP) and air-sea CO2 flux in a naturally iron fertilized and productive area of the Southern Ocean (Kerguelen plateau) was investigated using a 1D biogeochemical model driven by satellite chlorophyll, sea surface temperature and wind speed data for the 1997–2007 period. The model simulates the low fCO2 and dissolved inorganic carbon (DIC) measured during summers 2004–05, 2005–06, 2006–07 and the high NCP derived from a seasonal carbon budget in the surface waters of these blooms. Although satellite data show high interannual variability in the dynamics and magnitude of the bloom during the 1997–2007 decade, the simulated interannual variability of the NCP was only ± 14%. This unexpected result could be due to the combined effect of both the duration and the start date of the bloom, the latter determining the depth of the mixed layer used to compute the NCP. In the productive area, the interannual variability of air-sea CO2 flux (± 13%) was not only driven by the biological effect but also by the solubility effect. Our results contrast with previous studies in the high nutrient, low chlorophyll regions of the Southern Ocean.

Seasonal and interannual variability of oceanic carbon cycling in the western and central tropical-subtropical pacific: A physical-biogeochemical modeling study

The oceanic carbon cycle in the tropical-subtropical Pacific is strongly affected by various physical processes with different temporal and spatial scales, yet the mechanisms that regulate air-sea CO2 flux are not fully understood due to the paucity of both measurement and modeling. Using a 3-D physical-biogeochemical model, we simulate the partial pressure of CO2 in surface water (pCO2sea) and air-sea CO2 flux in the tropical and subtropical regions from 1990 to 2004. The model reproduces well the observed spatial differences in physical and biogeochemical processes, such as: (1) relatively higher sea surface temperature (SST), and lower dissolved inorganic carbon (DIC) and pCO2sea in the western than in the central tropical-subtropical Pacific, and (2) predominantly seasonal and interannual variations in the subtropical and tropical Pacific, respectively. Our model results suggest a non-negligible contribution of the wind variability to that of the air-sea CO2 flux in the central tropical Pacific, but the modeled contribution of 7% is much less than that from a previous modeling study (30%; McKinley et al., 2004). While DIC increases in the entire region SST increases in the subtropical and western tropical Pacific but decreases in the central tropical Pacific from 1990 to 2004. As a result, the interannual pCO2sea variability is different in different regions. The pCO2sea temporal variation is found to be primarily controlled by SST and DIC, although the role of salinity and total alkalinity, both of which also control pCO2sea, need to be elucidated by long-term observations and eddy-permitting models for better estimation of the interannual variability of air-sea CO2 flux.

Role of wind stress and heat fluxes in interannual‐to‐decadal variability of air‐sea CO2 and O2 fluxes in the North Atlantic

A coupled ecosystem‐circulation model of the North Atlantic is used to examine the individual contributions by wind stress and surface heat fluxes to naturally driven interannual‐to‐decadal variability of air‐sea fluxes of CO2 and O2 during 1948–2002. The model results indicate that variations in O2 fluxes are mainly driven by variations in surface heat fluxes in the extratropics (15°N to 70°N), and by wind stress in the tropics (10°S to 15°N). Conversely, variations in simulated CO2 fluxes are predominantly wind‐stress driven over the entire model domain (18°S to 70°N); while variability in piston velocity and surface heat fluxes is less important. The simulated uptake of O2 by the North Atlantic amounts to 70 ± 11 Tmol yr−1 to which the subpolar region (45°N to 70°N) contributes by 62 ± 10 Tmol yr−1. Whereas the subpolar North Atlantic takes up more than 2/3 of the total carbon absorbed by the North Atlantic in our model (about 0.3 Pg C yr−1), interannual variability of air‐sea CO2 fluxes reaches similar values (about 0.01 Pg C yr−1 each) in the subpolar (45°N to 70°N), the subtropical (15°N to 45°N) and the equatorial (10°S to 15°N) Atlantic.

Air-sea CO2 flux variability in the equatorial Pacific Ocean near 100oW

Air-sea CO2 flux variability in the equatorial Pacific Ocean near 100oW