In contrast to the open ocean, the sources and sinks for atmospheric carbon dioxide (CO2) in the coastal seas are poorly constrained and understood. Here we address this knowledge gap by analyzing the spatial and temporal variability of the coastal air‐sea flux of CO2 (FCO2) using a recent high‐resolution (0.25°) monthly climatology for coastal sea surface partial pressure in CO2 (pCO2). Coastal regions are characterized by CO2 sinks at temperate and high latitudes and by CO2 sources at low latitude and in the tropics, with annual mean CO2 flux densities comparable in magnitude and pattern to those of the adjacent open ocean with the exception of river‐dominated systems. The seasonal variations in FCO2 are large, often exceeding 2 mol C m−2 year−1, a magnitude similar to the variations exhibited across latitudes. The majority of these seasonal variations stems from the air‐sea pCO2 difference, although changes in wind speed and sea ice cover can also be significant regionally. Globally integrated, the coastal seas act currently as a CO2 sink of −0.20 ± 0.02 Pg C year−1, with a more intense uptake occurring in summer because of the disproportionate influence of high‐latitude shelves in the Northern Hemisphere. Combined with estimates of the carbon sinks in the open ocean and the Arctic, this gives for the global ocean, averaged over the 1998 to 2015 period an annual net CO2 uptake of −1.7 ± 0.3 Pg C year−1.

Abstract. Air–sea flux of carbon dioxide (CO2) is a critical component of the global carbon cycle and the climate system with the ocean removing about a quarter of the CO2 emitted into the atmosphere by human activities over the last decade. A common approach to estimate this net flux of CO2 across the air–sea interface is the use of surface ocean CO2 observations and the computation of the flux through a bulk parameterization approach. Yet, the details for how this is done in order to arrive at a global ocean CO2 uptake estimate vary greatly, enhancing the spread of estimates. Here we introduce the ensemble data product, SeaFlux (Gregor and Fay, 2021, https://doi.org/10.5281/zenodo.5482547​​​​​​​, https://github.com/luke-gregor/pySeaFlux, last access: 9 September 2021​​​​​​​); this resource enables users to harmonize an ensemble of products that interpolate surface ocean CO2 observations to near-global coverage with a common methodology to fill in missing areas in the products. Further, the dataset provides the inputs to calculate fluxes in a consistent manner. Utilizing six global observation-based mapping products (CMEMS-FFNN, CSIR-ML6, JENA-MLS, JMA-MLR, MPI-SOMFFN, NIES-FNN), the SeaFlux ensemble approach adjusts for methodological inconsistencies in flux calculations. We address differences in spatial coverage of the surface ocean CO2 between the mapping products, which ultimately yields an increase in CO2 uptake of up to 17 % for some products. Fluxes are calculated using three wind products (CCMPv2, ERA5, and JRA55). Application of a scaled gas exchange coefficient has a greater impact on the resulting flux than solely the choice of wind product. With these adjustments, we present an ensemble of global surface ocean pCO2 and air–sea carbon flux estimates. This work aims to support the community effort to perform model–data intercomparisons which will help to identify missing fluxes as we strive to close the global carbon budget.

Abstract. The calculation of the air-water CO 2 exchange (F CO 2 ) in the ocean not only depends on the gradient in CO 2 partial pressure at the air-water interface but also on the parameterization of the gas exchange transfer velocity (k) and the choice of wind product. Here, we present regional and global-scale quantifications of the uncertainty in F CO 2 induced by several widely used k formulations and four wind speed data products (CCMP, ERA, NCEP1 and NCEP2). The analysis is performed at a 1 • × 1 • resolution using the sea surface pCO 2 climatology generated by Landschützer et al. (2015a) for the 1991-2011 period, while the regional assessment relies on the segmentation proposed by the Regional Carbon Cycle Assessment and Processes (REC-CAP) project. First, we use k formulations derived from the global 14 C inventory relying on a quadratic relationship between k and wind speed (k = c · U 10 2 ; Sweeney et al., 2007;Takahashi et al., 2009;Wanninkhof, 2014), where c is a calibration coefficient and U 10 is the wind speed measured 10 m above the surface. Our results show that the range of global F CO 2 , calculated with these k relationships, diverge by 12 % when using CCMP, ERA or NCEP1. Due to differences in the regional wind patterns, regional discrepancies in F CO 2 are more pronounced than global. These global and regional differences significantly increase when using NCEP2 or other k formulations which include earlier relationships (i.e., Wanninkhof, 1992;Wanninkhof et al., 2009) as well as numerous local and regional parameterizations derived experimentally. To minimize uncertainties associated with the choice of wind product, it is possible to recalculate the coefficient c globally (hereafter called c * ) for a given wind product and its spatiotemporal resolution, in order to match the last evaluation of the global k value. We thus performed these recalculations for each wind product at the resolution and time period of our study but the resulting global F CO 2 estimates still diverge by 10 %. These results also reveal that the Equatorial Pacific, the North Atlantic and the Southern Ocean are the regions in which the choice of wind product will most strongly affect the estimation of the F CO 2 , even when using c * .

Abstract. In this study, we present the first combined open- and coastal-ocean pCO2 mapped monthly climatology (Landschützer et al., 2020b, https://doi.org/10.25921/qb25-f418, https://www.nodc.noaa.gov/ocads/oceans/MPI-ULB-SOM_FFN_clim.html, last access: 8 April 2020) constructed from observations collected between 1998 and 2015 extracted from the Surface Ocean CO2 Atlas (SOCAT) database. We combine two neural network-based pCO2 products, one from the open ocean and the other from the coastal ocean, and investigate their consistency along their common overlap areas. While the difference between open- and coastal-ocean estimates along the overlap area increases with latitude, it remains close to 0 µatm globally. Stronger discrepancies, however, exist on the regional level resulting in differences that exceed 10 % of the climatological mean pCO2, or an order of magnitude larger than the uncertainty from state-of-the-art measurements. This also illustrates the potential of such an analysis to highlight where we lack a good representation of the aquatic continuum and future research should be dedicated. A regional analysis further shows that the seasonal carbon dynamics at the coast–open interface are well represented in our climatology. While our combined product is only a first step towards a true representation of both the open-ocean and the coastal-ocean air–sea CO2 flux in marine carbon budgets, we show it is a feasible task and the present data product already constitutes a valuable tool to investigate and quantify the dynamics of the air–sea CO2 exchange consistently for oceanic regions regardless of its distance to the coast.

<p>The spatio-temporal variability and the underlying drivers of the carbon dioxide (CO<sub>2</sub>) exchange at the air-water interface (FCO<sub>2</sub>) of the global coastal ocean are still poorly understood and their quantification remains highly uncertain. Here, we present an analysis of the spatial and seasonal variability of FCO<sub>2</sub> using a high-resolution (0.25 degree) monthly climatology (1998-2015 period) for coastal sea surface partial pressure in CO<sub>2</sub> (pCO<sub>2</sub>), globally.</p><p>Overall, a clear latitudinal pattern emerges from our analysis regarding sources/sinks distribution of atmospheric CO<sub>2</sub> and we find that in most regions, annual mean CO<sub>2</sub> flux densities are comparable in sign and magnitude to those of the adjacent open ocean except for river dominated systems. Globally, coastal regions act as a CO<sub>2</sub> sink with a more intense uptake occurring in summer because of the disproportionate influence of high latitude coastal seas in the Northern Hemisphere. The majority of the coastal seasonal FCO<sub>2</sub> variations stems from the air-sea pCO<sub>2</sub> gradient, although changes in wind speed and sea-ice cover can also be significant regionally. To investigate further the drivers of the spatio-seasonal variability, our observation-based pCO<sub>2</sub> climatology is used in conjunction with global ocean biogeochemistry model MOM6-COBALT. The model outputs allow us to quantify the respective contributions of thermal effects, biology, and non-thermal physical processes (circulation and freshwater inputs) to seasonal variations in coastal pCO<sub>2</sub>. Generally, biological activity is the dominant driver of the pCO<sub>2</sub> seasonal variability in temperate and high latitudes while thermal and non-thermal physical processes dominate in low latitudes.</p>

Abstract. The temporal variability of the sea surface partial pressure of CO2 (pCO2) and the underlying processes driving this variability are poorly understood in the coastal ocean. In this study, we tailor an existing method that quantifies the effects of thermal changes, biological activity, ocean circulation and freshwater fluxes to examine seasonal pCO2 changes in highly variable coastal environments. We first use the Modular Ocean Model version 6 (MOM6) and biogeochemical module Carbon Ocean Biogeochemistry And Lower Trophics version 2 (COBALTv2) at a half-degree resolution to simulate coastal CO2 dynamics and evaluate them against pCO2 from the Surface Ocean CO2 Atlas database (SOCAT) and from the continuous coastal pCO2 product generated from SOCAT by a two-step neuronal network interpolation method (coastal Self-Organizing Map Feed-Forward neural Network SOM-FFN, Laruelle et al., 2017). The MOM6-COBALT model reproduces the observed spatiotemporal variability not only in pCO2 but also in sea surface temperature, salinity and nutrients in most coastal environments, except in a few specific regions such as marginal seas. Based on this evaluation, we identify coastal regions of “high” and “medium” agreement between model and coastal SOM-FFN where the drivers of coastal pCO2 seasonal changes can be examined with reasonable confidence. Second, we apply our decomposition method in three contrasted coastal regions: an eastern (US East Coast) and a western (the Californian Current) boundary current and a polar coastal region (the Norwegian Basin). Results show that differences in pCO2 seasonality in the three regions are controlled by the balance between ocean circulation and biological and thermal changes. Circulation controls the pCO2 seasonality in the Californian Current; biological activity controls pCO2 in the Norwegian Basin; and the interplay between biological processes and thermal and circulation changes is key on the US East Coast. The refined approach presented here allows the attribution of pCO2 changes with small residual biases in the coastal ocean, allowing for future work on the mechanisms controlling coastal air–sea CO2 exchanges and how they are likely to be affected by future changes in sea surface temperature, hydrodynamics and biological dynamics.

Abstract. The temporal variability of the sea surface partial pressure of CO2 (pCO2) and the underlying processes driving this variability are poorly understood in the coastal ocean. In this study, we tailor an existing method that quantifies the effects of thermal changes, biological activity, ocean circulation and fresh water fluxes to examine seasonal pCO2 changes in highly-variable coastal environments. We first use the Modular Ocean Model version 6 (MOM6) and biogeochemical module Carbon Ocean Biogeochemistry And Lower Trophics version 2 (COBALTv2) at a half degree resolution to simulate the coastal CO2 dynamics and evaluate it against pCO2 from the Surface Ocean CO2 Atlas database (SOCAT) and from the continuous coastal pCO2 product generated from SOCAT by a two-step neuronal network interpolation method (coastal-SOM-FFN, Laruelle et al., 2017). The MOM6-COBALT model not only reproduces the observed spatio-temporal variability in pCO2 but also in sea surface temperature, salinity, nutrients, in most coastal environments except in a few specific regions such as marginal seas. Based on this evaluation, we identify coastal regions of ‘high’ and ‘medium’ model skill where the drivers of coastal pCO2 seasonal changes can be examined with reasonable confidence. Second, we apply our decomposition method in three contrasted coastal regions: an Eastern (East coast of the U.S) and a Western (the Californian Current) boundary current and a polar coastal region (the Norwegian Basin). Results show that differences in pCO2 seasonality in the three regions are controlled by the balance between ocean circulation, biological and thermal changes. Circulation controls the pCO2 seasonality in the Californian Current, biological activity controls pCO2 in the Norwegian Basin, while the interplay between biology, thermal and circulation changes is key in the East coast of the U.S. The refined approach presented here allows the attribution of pCO2 changes with small residual biases in the coastal ocean, allowing future work on the mechanisms controlling coastal air-sea CO2 exchanges and how they are likely to be affected by future changes in sea surface temperature, hydrodynamics and biological dynamics.

The coastal ocean contributes to regulating atmospheric greenhouse gas concentrations by taking up carbon dioxide (CO2) and releasing nitrous oxide (N2O) and methane (CH4). Major advances have improved our understanding of the coastal air-sea exchanges of these three gasses since the first phase of the Regional Carbon Cycle Assessment and Processes (RECCAP in 2013), but a comprehensive view that integrates the three gasses at the global scale is still lacking. In this second phase (RECCAP2), we quantify global coastal ocean fluxes of CO2, N2O and CH4 using an ensemble of global gap-filled observation-based products and ocean biogeochemical models. The global coastal ocean is a net sink of CO2 in both observational products and models, but the magnitude of the median net global coastal uptake is ~60% larger in models (-0.72 vs. -0.44 PgC/yr, 1998-2018, coastal ocean area of 77 million km2). We attribute most of this model-product difference to the seasonality in sea surface CO2 partial pressure at mid- and high-latitudes, where models simulate stronger winter CO2 uptake. The global coastal ocean is a major source of N2O (+0.70 PgCO2-e /yr in observational product and +0.54 PgCO2-e /yr in model median) and of CH4 (+0.21 PgCO2-e /yr in observational product), which offsets a substantial proportion of the net radiative effect of coastal \co uptake (35-58% in CO2-equivalents). Data products and models need improvement to better resolve the spatio-temporal variability and long term trends in CO2, N2O and CH4 in the global coastal ocean.