Search

We have compared a suite of recent global CO2 atmospheric inversion results to independent airborne observations and to each other, to assess their dependence on differences in northern extratropical (NET) vertical transport and to identify some of the drivers of model spread. We evaluate posterior CO2 concentration profiles against observations from the High-Performance Instrumented Airborne Platform for Environmental Research (HIAPER) Pole-to-Pole Observations (HIPPO) aircraft campaigns over the mid-Pacific in 2009–2011. Although the models differ in inverse approaches, assimilated observations, prior fluxes, and transport models, their broad latitudinal separation of land fluxes has converged significantly since the Atmospheric Carbon Cycle Inversion Intercomparison (TransCom 3) and the REgional Carbon Cycle Assessment and Processes (RECCAP) projects, with model spread reduced by 80 % since TransCom 3 and 70 % since RECCAP. Most modeled CO2 fields agree reasonably well with the HIPPO observations, specifically for the annual mean vertical gradients in the Northern Hemisphere. Northern Hemisphere vertical mixing no longer appears to be a dominant driver of northern versus tropical (T) annual flux differences. Our newer suite of models still gives northern extratropical land uptake that is modest relative to previous estimates (Gurney et al., 2002; Peylin et al., 2013) and near-neutral tropical land uptake for 2009–2011. Given estimates of emissions from deforestation, this implies a continued uptake in intact tropical forests that is strong relative to historical estimates (Gurney et al., 2002; Peylin et al., 2013). The results from these models for other time periods (2004–2014, 2001–2004, 1992–1996) and re-evaluation of the TransCom 3 Level 2 and RECCAP results confirm that tropical land carbon fluxes including deforestation have been near neutral for several decades. However, models still have large disagreements on ocean–land partitioning. The fossil fuel (FF) and the atmospheric growth rate terms have been thought to be the best-known terms in the global carbon budget, but we show that they currently limit our ability to assess regional-scale terrestrial fluxes and ocean–land partitioning from the model ensemble.

Air–sea CO2 fluxes over the Pacific Ocean are known to be characterized by coherent large-scale structures that reflect not only ocean subduction and upwelling patterns, but also the combined effects of wind-driven gas exchange and biology. On the largest scales, a large net CO2 influx into the extratropics is associated with a robust seasonal cycle, and a large net CO2 efflux from the tropics is associated with substantial interannual variability. In this work, we have synthesized estimates of the net air–sea CO2 flux from a variety of products, drawing upon a variety of approaches in three sub-basins of the Pacific Ocean, i.e., the North Pacific extratropics (18–66° N), the tropical Pacific (18° S–18° N), and the South Pacific extratropics (44.5–18° S). These approaches include those based on the measurements of CO2 partial pressure in surface seawater (pCO2sw), inversions of ocean-interior CO2 data, forward ocean biogeochemistry models embedded in the ocean general circulation models (OBGCMs), a model with assimilation of pCO2sw data, and inversions of atmospheric CO2 measurements. Long-term means, interannual variations and mean seasonal variations of the regionally integrated fluxes were compared in each of the sub-basins over the last two decades, spanning the period from 1990 through 2009. A simple average of the long-term mean fluxes obtained with surface water pCO2 diagnostics and those obtained with ocean-interior CO2 inversions are −0.47 ± 0.13 Pg C yr−1 in the North Pacific extratropics, +0.44 ± 0.14 Pg C yr−1 in the tropical Pacific, and −0.37 ± 0.08 Pg C yr−1 in the South Pacific extratropics, where positive fluxes are into the atmosphere. This suggests that approximately half of the CO2 taken up over the North and South Pacific extratropics is released back to the atmosphere from the tropical Pacific. These estimates of the regional fluxes are also supported by the estimates from OBGCMs after adding the riverine CO2 flux, i.e., −0.49 ± 0.02 Pg C yr−1 in the North Pacific extratropics, +0.41 ± 0.05 Pg C yr−1 in the tropical Pacific, and −0.39 ± 0.11 Pg C yr−1 in the South Pacific extratropics. The estimates from the atmospheric CO2 inversions show large variations amongst different inversion systems, but their median fluxes are consistent with the estimates from climatological pCO2sw data and pCO2sw diagnostics. In the South Pacific extratropics, where CO2 variations in the surface and ocean interior are severely undersampled, the difference in the air–sea CO2 flux estimates between the diagnostic models and ocean-interior CO2 inversions is larger (0.18 Pg C yr−1). The range of estimates from forward OBGCMs is also large (−0.19 to −0.72 Pg C yr−1). Regarding interannual variability of air–sea CO2 fluxes, positive and negative anomalies are evident in the tropical Pacific during the cold and warm events of the El Niño–Southern Oscillation in the estimates from pCO2sw diagnostic models and from OBGCMs. They are consistent in phase with the Southern Oscillation Index, but the peak-to-peak amplitudes tend to be higher in OBGCMs (0.40 ± 0.09 Pg C yr−1) than in the diagnostic models (0.27 ± 0.07 Pg C yr−1).

The Indian Ocean (44° S–30° N) plays an important role in the global carbon cycle, yet it remains one of the most poorly sampled ocean regions. Several approaches have been used to estimate net sea–air CO2 fluxes in this region: interpolated observations, ocean biogeochemical models, atmospheric and ocean inversions. As part of the RECCAP (REgional Carbon Cycle Assessment and Processes) project, we combine these different approaches to quantify and assess the magnitude and variability in Indian Ocean sea–air CO2 fluxes between 1990 and 2009. Using all of the models and inversions, the median annual mean sea–air CO2 uptake of −0.37 ± 0.06 PgC yr−1 is consistent with the −0.24 ± 0.12 PgC yr−1 calculated from observations. The fluxes from the southern Indian Ocean (18–44° S; −0.43 ± 0.07 PgC yr−1 are similar in magnitude to the annual uptake for the entire Indian Ocean. All models capture the observed pattern of fluxes in the Indian Ocean with the following exceptions: underestimation of upwelling fluxes in the northwestern region (off Oman and Somalia), overestimation in the northeastern region (Bay of Bengal) and underestimation of the CO2 sink in the subtropical convergence zone. These differences were mainly driven by lack of atmospheric CO2 data in atmospheric inversions, and poor simulation of monsoonal currents and freshwater discharge in ocean biogeochemical models. Overall, the models and inversions do capture the phase of the observed seasonality for the entire Indian Ocean but overestimate the magnitude. The predicted sea–air CO2 fluxes by ocean biogeochemical models (OBGMs) respond to seasonal variability with strong phase lags with reference to climatological CO2 flux, whereas the atmospheric inversions predicted an order of magnitude higher seasonal flux than OBGMs. The simulated interannual variability by the OBGMs is weaker than that found by atmospheric inversions. Prediction of such weak interannual variability in CO2 fluxes by atmospheric inversions was mainly caused by a lack of atmospheric data in the Indian Ocean. The OBGM models suggest a small strengthening of the sink over the period 1990–2009 of −0.01 PgC decade−1. This is inconsistent with the observations in the southwestern Indian Ocean that shows the growth rate of oceanic pCO2 was faster than the observed atmospheric CO2 growth, a finding attributed to the trend of the Southern Annular Mode (SAM) during the 1990s.

Atmospheric CO2 inversions estimate surface carbon fluxes from an optimal fit to atmospheric CO2 measurements, usually including prior constraints on the flux estimates. Eleven sets of carbon flux estimates are compared, generated by different inversions systems that vary in their inversions methods, choice of atmospheric data, transport model and prior information. The inversions were run for at least 5 yr in the period between 1990 and 2010. Mean fluxes for 2001–2004, seasonal cycles, interannual variability and trends are compared for the tropics and northern and southern extra-tropics, and separately for land and ocean. Some continental/basin-scale subdivisions are also considered where the atmospheric network is denser. Four-year mean fluxes are reasonably consistent across inversions at global/latitudinal scale, with a large total (land plus ocean) carbon uptake in the north (−3.4 Pg C yr−1 (±0.5 Pg C yr−1 standard deviation), with slightly more uptake over land than over ocean), a significant although more variable source over the tropics (1.6 ± 0.9 Pg C yr−1) and a compensatory sink of similar magnitude in the south (−1.4 ± 0.5 Pg C yr−1) corresponding mainly to an ocean sink. Largest differences across inversions occur in the balance between tropical land sources and southern land sinks. Interannual variability (IAV) in carbon fluxes is larger for land than ocean regions (standard deviation around 1.06 versus 0.33 Pg C yr−1 for the 1996–2007 period), with much higher consistency among the inversions for the land. While the tropical land explains most of the IAV (standard deviation ~ 0.65 Pg C yr−1), the northern and southern land also contribute (standard deviation ~ 0.39 Pg C yr−1). Most inversions tend to indicate an increase of the northern land carbon uptake from late 1990s to 2008 (around 0.1 Pg C yr−1, predominantly in North Asia. The mean seasonal cycle appears to be well constrained by the atmospheric data over the northern land (at the continental scale), but still highly dependent on the prior flux seasonality over the ocean. Finally we provide recommendations to interpret the regional fluxes, along with the uncertainty estimates.

The source and sinks of carbon dioxide (CO2) and methane (CH4) due to anthropogenic and natural biospheric activities were estimated for the South Asian region (Bangladesh, Bhutan, India, Nepal, Pakistan and Sri Lanka). Flux estimates were based on top-down methods that use inversions of atmospheric data, and bottom-up methods that use field observations, satellite data, and terrestrial ecosystem models. Based on atmospheric CO2 inversions, the net biospheric CO2 flux in South Asia (equivalent to the Net Biome Productivity, NBP) was a sink, estimated at −104 ± 150 Tg C yr−1 during 2007–2008. Based on the bottom-up approach, the net biospheric CO2 flux is estimated to be −191 ± 193 Tg C yr−1 during the period of 2000–2009. This last net flux results from the following flux components: (1) the Net Ecosystem Productivity, NEP (net primary production minus heterotrophic respiration) of −220 ± 186 Tg C yr−1 (2) the annual net carbon flux from land-use change of −14 ± 50 Tg C yr−1, which resulted from a sink of −16 Tg C yr−1 due to the establishment of tree plantations and wood harvest, and a source of 2 Tg C yr−1 due to the expansion of croplands; (3) the riverine export flux from terrestrial ecosystems to the coastal oceans of +42.9 Tg C yr−1; and (4) the net CO2 emission due to biomass burning of +44.1 ± 13.7 Tg C yr−1. Including the emissions from the combustion of fossil fuels of 444 Tg C yr−1 for the 2000s, we estimate a net CO2 land–atmosphere flux of 297 Tg C yr−1. In addition to CO2, a fraction of the sequestered carbon in terrestrial ecosystems is released to the atmosphere as CH4. Based on bottom-up and top-down estimates, and chemistry-transport modeling, we estimate that 37 ± 3.7 Tg C yr−1 were released to atmosphere from South Asia during the 2000s. Taking all CO2 and CH4 fluxes together, our best estimate of the net land–atmosphere CO2-equivalent flux is a net source of 334 Tg C yr−1 for the South Asian region during the 2000s. If CH4 emissions are weighted by radiative forcing of molecular CH4, the total CO2-equivalent flux increases to 1148 Tg C yr−1 suggesting there is great potential of reducing CH4 emissions for stabilizing greenhouse gases concentrations.

Extensive observations were made during the late Southwest Monsoon of 2004 over the Indian and Omani shelves, and along a transect that extended from the southern coast of Oman to the central west coast of India, tracking the southern leg of the US JGOFS expedition (1994–1995) in the west. The data are used, in conjunction with satellite-derived data, to investigate long-term trends in chlorophyll and sea surface temperature, indicators of upwelling intensity, and to understand factors that control primary production (PP) in the Arabian Sea, focussing on the role of iron. Our results do not support an intensification of upwelling in the western Arabian Sea, reported to have been caused by the decline in the winter/spring Eurasian snow cover since 1997. We also noticed, for the first time, an unexpected development of high-nutrient, low-chlorophyll condition off the southern Omani coast. This feature, coupled with other characteristics of the system, such as a narrow shelf and relatively low iron concentrations in surface waters, suggest a close similarity between the Omani upwelling system and the Peruvian and California upwelling systems, where PP is limited by iron. Iron limitation of PP may complicate simple relationship between upwelling and PP assumed by previous workers, and contribute to the anomalous offshore occurrence of the most severe oxygen (O2) depletion in the region. Over the much wider Indian shelf, which experiences large-scale bottom water O2-depletion in summer, adequate iron supply from reducing bottom-waters and sediments seems to support moderately high PP; however, such production is restricted to the thin, oxygenated surface layer, probably because of the unsuitability of the O2-depleted environment for the growth of oxygenic photosynthesizers.