Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
Atmos. Chem. Phys., 19, 9269–9285, 2019
https://doi.org/10.5194/acp-19-9269-2019
Atmos. Chem. Phys., 19, 9269–9285, 2019
https://doi.org/10.5194/acp-19-9269-2019

Research article 19 Jul 2019

Research article | 19 Jul 2019

# Global carbon budgets estimated from atmospheric O2∕N2 and CO2 observations in the western Pacific region over a 15-year period

Global carbon budgets estimated from atmospheric O2∕N2 and CO2 observations in the western Pacific region over a 15-year period
Yasunori Tohjima1, Hitoshi Mukai2, Toshinobu Machida2, Yu Hoshina1, and Shin-Ichiro Nakaoka2 Yasunori Tohjima et al.
• 1Center for Environmental Measurement and Analysis, National Institute for Environmental Studies, Tsukuba, 305-8506, Japan
• 2Center for Global Environmental Research, National Institute for Environmental Studies, Tsukuba, 305-8506, Japan

Correspondence: Yasunori Tohjima (tohjima@nies.go.jp)

Abstract

Time series of the atmospheric O2∕N2 ratio and CO2 mole fraction of flask samples obtained from the National Institute for Environmental Studies' (NIES's) flask sampling network are presented. The network includes two ground sites, Hateruma island (HAT; 24.05 N, 123.81 E) and Cape Ochiishi (COI; 43.17 N, 145.50 E), and cargo ships regularly sailing in the western Pacific. Based on temporal changes in fossil-fuel-derived CO2 emissions, global atmospheric CO2 burden and atmospheric potential oxygen (APO), which were calculated from the observed O2∕N2 ratio and CO2 mole fraction according to APO =O2$+\mathrm{1.1}×{\mathrm{CO}}_{\mathrm{2}}$, we estimated the global carbon sinks of the ocean and land biosphere for a period of more than 15 years. In this carbon budget calculation, we adopted a correction for the time-varying ocean O2 outgassing effect with an average of 0.54 PgC yr−1 for 2000–2016. The outgassing effect, attributed mainly to global ocean warming, was evaluated under the assumption that the net ocean gas flux is proportional to the change in the ocean heat content for the 0–2000 m layer. The resulting oceanic and land biotic carbon sinks were 2.6±0.7 and 1.5±0.9 PgC yr−1, respectively, for a 17-year period (2000–2016) and 2.4±0.7 and 1.9±0.9 PgC yr−1, respectively, for a 14-year period (2003–2016). Despite the independent approaches, the sink values of this study agreed with those estimated by the Global Carbon Project (GCP) within a difference of about ±0.4 PgC yr−1. We examined the carbon sinks for an interval of 5 years to assess the temporal trends. The pentad (5-year) ocean sinks showed an increasing trend at a rate of 0.08±0.02 PgC yr−2 during 2001–2014, while the pentad land sinks showed an increasing trend at a rate of 0.23±0.04 PgC yr−2 for 2001–2009 and a decreasing trend at a rate of $-\mathrm{0.22}±\mathrm{0.04}$ PgC yr−2 during 2009–2014. Although there is good agreement in the trends of the pentad sinks between this study and that of GCP, the increasing rate of the pentad ocean sinks of this study was about 2 times larger than that of GCP.

1 Introduction

In spite of various international efforts to reduce anthropogenic carbon dioxide (CO2) emissions, the atmospheric CO2 levels observed around the world have shown a steady increase and exceeded the benchmark of 400 parts per million (ppm) mole fraction in past years (Betts et al., 2016). The Carbon Dioxide Information Analysis Center (CDIAC) reported that global fossil-fuel-derived CO2 emissions in recent years still increased gradually and rose to 9.9 PgC yr−1 by 2014 (Boden et al., 2017). Under these circumstances, the Paris Agreement adopted at the 21st Conference of Parties (COP) to the United Nations Framework Convention on Climate Change in 2015 aimed to reduce the anthropogenic greenhouse gas emissions to maintain the increase in global mean surface temperatures at well below 2 by 2100 and, if possible, to limit the increase to 1.5. To achieve this goal, it is crucially important to quantitatively understand the natural sink strengths or land biosphere and ocean sinks. A variety of approaches have so far been applied to the quantification of ocean or land sinks or both, including process-based land and ocean models, bottom-up emission estimates based on flux measurements, and top-down estimates based on atmospheric measurements. Developing process-based models to enhance the accuracy of the global carbon budget is crucially important because they are expected to predict the future global carbon cycle in a warmer world. However, carbon budget estimates based on observations are still important in validating and improving the process-based models.

The budget estimation based on atmospheric CO2 and O2 observations is a simple and straightforward approach, and it has historically settled the controversy of whether the land biosphere is a net carbon sink or source (Keeling and Shertz, 1992). Although several techniques based on the interferometer (Keeling, 1988), mass spectrometry (Bender et al., 1994), paramagnetic analyzer (Manning et al., 1999), fuel cell analyzer (Stephens et al., 2007), vacuum ultraviolet absorption photometer (Stephens et al., 2003) and so on have been developed to detect the parts-per-million-level changes in the atmospheric O2 concentration, the accurate quantification of the O2 change is still challenging. The carbon budget is evaluated by simultaneously solving the mass balance equations of the atmospheric CO2 and O2 as follows (Manning and Keeling, 2006):

$\begin{array}{}\text{(1)}& \mathrm{\Delta }{\mathrm{CO}}_{\mathrm{2}}=F-B-O,\text{(2)}& \mathrm{\Delta }{\mathrm{O}}_{\mathrm{2}}=-{\mathit{\alpha }}_{\mathrm{F}}F+{\mathit{\alpha }}_{\mathrm{B}}B+Z,\end{array}$

where ΔCO2 and ΔO2 represent the changes in the atmospheric CO2 and O2 burdens based on atmospheric observations, respectively, F represents the fossil-fuel-derived emissions, and B and O represent the uptake by the land biosphere and the ocean, respectively. αF and αB are the O2∕C exchange ratio for the globally averaged fossil-fuel combustions and land biotic processes, respectively. The estimated value for αB is about 1.1 (Severinghaus, 1995), and that for αF is about 1.4 (Keeling, 1988). These equations mean that the CO2 and O2 fluxes associated with fossil-fuel combustion and land biotic processes are tightly coupled. In contrast, the ocean CO2 uptake, O, and ocean O2 emissions, denoted as Z, are decoupled because the ocean acts as a carbon sink by physicochemically dissolving the CO2. Since the values of F and αF can be evaluated from energy statistics (Keeling, 1988), we can evaluate the ocean and land uptake by solving the above equations if we could evaluate the value of Z.

The global carbon budget can also be related to tracer atmospheric potential oxygen (APO), which is defined by the equation of APO =O2$+{\mathit{\alpha }}_{\mathrm{B}}×$CO2 (Stephens et al., 1998). Since the APO is defined to be invariant with respect to the land biotic exchange, the secular trend in the APO is determined by fossil-fuel combustions which cause a gradually decreasing trend in APO and the air–sea gas exchange. Combining Eq. (1), multiplied by αB, and Eq. (2), in accordance with the APO definition, results in the following equation for the APO budget (Manning and Keeling, 2006):

$\begin{array}{}\text{(3)}& \mathrm{\Delta }\mathrm{APO}=-\left({\mathit{\alpha }}_{\mathrm{F}}-{\mathit{\alpha }}_{\mathrm{B}}\right)F-{\mathit{\alpha }}_{\mathrm{B}}O+Z.\end{array}$

Since observation sites for atmospheric O2 are still limited compared with those for atmospheric CO2, Manning and Keeling (2006) proposed an alternative approach that the global carbon budgets could be obtained by simultaneously solving Eqs. (1) and (3) and using globally averaged CO2 data based on the measurements of the Global Monitoring Division of the National Oceanic and Atmospheric Administration's Earth System Research Laboratory (NOAA/ESRL/GMD). This approach, making maximum use of the available data, is expected to give the most reliable estimation. This APO approach has been adopted for the estimation of the global carbon budget based on atmospheric O2 and CO2 measurements (e.g., Manning and Keeling, 2006; Tohjima et al., 2008; Ishidoya et al., 2012; Goto et al., 2017).

To evaluate the carbon budgets based on O2 and CO2 measurements, we need to quantify the magnitude of Z and its temporal variation, if possible. It is considered that Z has a large interannual variability because observed trends of APO generally show large interannual variations which would result in unrealistic variations in the ocean uptake if the variability in Z were rather small (e.g., Bender et al., 2005). Probably, an imbalance of the air–sea seasonal O2 exchanges, outgassing flux associated with primary production in spring and summer, and ingassing flux associated with ocean ventilation in autumn and winter cause the interannual variations in Z. The results of ocean model simulations also support this mechanism (e.g., McKinley et al., 2003; Nevison et al., 2008). Therefore, it is difficult to estimate the short-term carbon budgets unless the temporal variations in Z are accurately evaluated. Additionally, as for long timescales, it is considered that the present ocean acts as an O2 source because of the global ocean warming (Keeling and Garcia, 2002). The increase in surface ocean temperature not only reduces the solubility of gases in seawater but also strengthens the ocean stratification, which reduces the ventilation of interior water masses. The reduction of ventilation reduces the ingassing flux of O2. In contrast, the reduction of ventilation also causes a reduction of the nutrient supply from deep water, which might decrease the primary production and O2 outgassing in summer. Therefore, the influence of the ocean warming on the net air–sea gas exchange is rather complicated.

Unfortunately, there is little observational evidence to quantify the magnitude of Z and its temporal variations. The long-term average values of Z are inferred under the assumption that Z is proportional to the air–sea heat flux (Keeling and Garcia, 2002). The change in the global ocean heat content has been evaluated based on the large dataset of ocean observations (Levitus et al., 2012). Keeling and Garcia (2002) estimated the ratio of the O2 flux to heat flux from the relationship between the dissolved O2 corrected for the mineralization effect and the potential temperature. This approach was basically adopted by most of the studies to evaluate the long-term global carbon budgets (e.g., Bender et al., 2005; Manning and Keeling, 2006; Tohjima et al., 2008). On the other hand, Ishidoya et al. (2012) evaluated the instantaneous variations in the land and ocean sinks based on the APO data at Ny-Ålesund, Svalbard and Syowa, Antarctica, for the period 2001–2009 by calculating the interannual variation in Z from the temporal variation in the ocean heat content. They concluded that the above-mentioned Z values adequately suppressed artifacts caused by the imbalance of the seasonal air–sea O2 exchange.

We have been conducting air sampling into glass flasks for the measurement of the atmospheric O2∕N2 ratio and CO2 mole fraction at two ground sites in Japan since the late 1990s (Tohjima et al., 2003) and evaluated the global carbon budgets for up to 7 years (1999–2005) based on the APO data from the flask observations (Tohjima et al., 2008). To extend the observation area, we started additional flask sampling aboard commercial cargo ships regularly sailing in the Pacific region in 2002 (Tohjima et al., 2005b, 2012). About a decade has passed since we previously reported the global carbon budgets, and now we have more than a 15-year-long record of atmospheric O2∕N2 and CO2 of the flask samples. In this study, we estimated the ocean and land biotic carbon sinks for over a decade by using the temporal changes in the APO based on these flask data. In addition, we sequentially computed the ocean and land sinks for an interval of 5 years and examined the changing trends of both sinks. In these budget calculations, we estimated the values of Z for the corresponding period by using the temporal changes in the global ocean heat content. Finally, the estimated ocean and land carbon sinks of this study were compared with those of the Global Carbon Project (GCP).

2 Data and analysis

We started air samplings for the measurement of the atmospheric O2∕N2 ratio and CO2 mole fraction at two monitoring stations located on Hateruma island (HAT; 24.05 N, 123.81 E) in July 1997 and at Cape Ochiishi (COI; 43.17 N, 145.50 E) in December 1998 (Tohjima et al., 2003). In addition, we have been collecting air samples from the Pacific regions by using commercial cargo vessels equipped with automated flask sampling systems (Tohjima et al., 2005b; 2012). The shipboard flask samplings were started between Japan and North America in December 2001, between Japan and Australia and New Zealand in December 2001, and between Japan and southeastern Asia in September 2007. The flask sampling sites are depicted in Fig. 1. Unfortunately, the shipboard data in southeastern Asia, the northern North Pacific (north of 30 N) and the eastern North Pacific are spatiotemporally rather sporadic. (See the inserted figure in Fig. 1 showing time–latitude plots of the shipboard flask samples.) Thus, in the following analysis, we only used the dataset obtained at HAT, COI and the western Pacific region between 40 S and 30 N and between 130 and 180 E.

Air samples were collected in glass flasks hermetically sealed by two glass valves with Viton O-rings. The volumes of the flasks were 2 L for the samplings at HAT and COI and 2.5 L for the shipboard samplings. It should be noted that glass flasks with a volume of 1 L were also used in the early period from the start to March 2006, and only 1 L flasks were used from the start to January 1999 at HAT.

Figure 1Map showing the air sampling locations in the Pacific region. The light blue and red squares represent the monitoring stations of COI and HAT, respectively. The orange, blue and green circles correspond to the positions where flask samplings were taken aboard cargo ships in southeastern Asia, western Pacific and eastern Pacific, respectively. The inserted figure shows the time–latitude distribution of onboard flask samples. The flask data from COI, HAT and the regions in the black rectangle (130–180 E, 40 S–30 N) were used in the budget calculations.

## 2.2O2∕N2 and CO2 analytical methods

In this study, we used a gas chromatograph (GC) equipped with a thermal conductivity detector (TCD) for the measurements of atmospheric O2 (Tohjima, 2000). In this GC–TCD method, O2∕N2 ratios of sample air and working reference air were alternately measured, and the atmospheric O2 change was determined as the relative difference in the O2∕N2 ratio from an arbitrary reference. We used the delta notation according to Keeling and Shertz (1992) to express the relatively small difference in the O2∕N2 ratio as follows:

$\begin{array}{}\text{(4)}& \mathit{\delta }\left({\mathrm{O}}_{\mathrm{2}}/{\mathrm{N}}_{\mathrm{2}}\right)=\frac{{\left({\mathrm{O}}_{\mathrm{2}}/{\mathrm{N}}_{\mathrm{2}}\right)}_{\mathrm{sam}}}{{\left({\mathrm{O}}_{\mathrm{2}}/{\mathrm{N}}_{\mathrm{2}}\right)}_{\mathrm{ref}}}-\mathrm{1},\end{array}$

where subscripts “sam” and “ref” refer to sample and reference, respectively, and the δ(O2∕N2) value multiplied by 106 is expressed in per meg units. The change in 1 µmol of O2 per mole of dry air changed the O2∕N2 ratio by 4.77 per meg, which corresponds to a 1 ppm change in the atmospheric trace gas abundance. APO was calculated from the CO2 mole fraction (X${}_{{\mathrm{CO}}_{\mathrm{2}}}$) in parts per million and δ(O2∕N2) in per meg according to

$\begin{array}{}\text{(5)}& \mathit{\delta }\mathrm{APO}=\mathit{\delta }\left({\mathrm{O}}_{\mathrm{2}}/{\mathrm{N}}_{\mathrm{2}}\right)+{\mathit{\alpha }}_{\mathrm{B}}{X}_{{\mathrm{CO}}_{\mathrm{2}}}/{S}_{{\mathrm{O}}_{\mathrm{2}}}-\mathrm{1850},\end{array}$

where ${S}_{{\mathrm{O}}_{\mathrm{2}}}$ is the mole fraction of O2 in the air (${S}_{{\mathrm{O}}_{\mathrm{2}}}=\mathrm{0.2094}$, Tohjima et al., 2005a) and the value of 1850 is an arbitrary reference point of δAPO in per meg. The values of δ(O2∕N2) were determined against the National Institute for Environmental Studies (NIES) O2∕N2 scale (Tohjima et al., 2008). Its temporal stability is examined in the following section.

A nondispersive infrared (NDIR) analyzer (LI-COR, Lincoln, Nebraska, model LI-6252) was used for the CO2 measurement of the flask samples. The CO2 mole fractions were determined against the NIES 09 scale, which is based on a set of gravitationally prepared CO2-in-air standard gases (Machida et al., 2011). The relationship between the NIES 09 scale and the NOAA scale were repeatedly compared through the WMO round-robin intercomparison program. The results showed that the differences of the NIES 09 scale from the NOAA scale were kept within ±0.15 ppm during the period from 1996 to 2014 (https://www.esrl.noaa.gov/gmd/ccgg/wmorr/wmorr_results.php?, last access: 16 July 2019).

## 2.3O2∕N2 scale stability

As details of the NIES O2∕N2 scale are given elsewhere (Tohjima et al., 2008), here we describe only briefly the outline of the scale and add some new information below. The zero point of NIES O2∕N2 scale was related to an ambient dry air stored in a high-pressure cylinder (HDA-1). The O2∕N2 scale was maintained by three cylinders during the first 4 years (1997–2001) of our O2∕N2 measurement program. In 2001, the NIES O2∕N2 scale was transferred to 11 other high-pressure cylinders (five 10 L cylinders and six 48 L cylinders), in which the δ(O2∕N2) values were carefully determined against the original O2∕N2 scale. Another primary reference gas (48 L cylinder, CQB-07080) was added in 2002, and now 12 primary reference gases keep the NIES O2∕N2 scale. The air samples delivered from glass flasks or high-pressure cylinders were measured against the working reference airs stored in 48 L aluminum cylinders, from which the δ(O2∕N2) values were repeatedly determined against the individual primary gas cylinders at intervals of a few months. The working reference gas cylinders were replaced by new ones every 1–2 years.

Figure 2Temporal changes in the O2∕N2 ratio of primary reference gases relative to the NIES O2∕N2 scale. The differences of the O2∕N2 ratio from the average are plotted for HDA-2 along with the differences of the O2∕N2 ratio from the initial values for the individual cylinders except HDA-2. Solid and dashed horizontal bars in the lower part of the figure indicate the periods when the working reference gases were used.

Figure 2 is the extended version of the previously reported figure (Fig. 1 in Tohjima et al., 2008), showing the temporal changes in the O2∕N2 ratio of primary reference gases relative to the NIES O2∕N2 scale. In the figure, the deviations of the O2∕N2 ratio from the average value for another high-pressure cylinder (HDA-2) and from the initially determined values for the second set of 12 cylinders are plotted. The averages and the standard deviations (1σ) of the differences for the individual 12 cylinders range from 4.2 to 3.3 per meg and from 3.1 to 5.2 per meg, respectively. The changing rates of the deviations for the 12 reference gases during 2001–2017, determined by least square linear regression, range from 0.34 to 0.2 per meg yr−1. Solid and dashed horizontal bars in the bottom of the figure indicate the durations of use of the individual working reference gases.

To assess the stability of the NIES O2∕N2 scales, we have continued to measure the reference gases in two 48 L aluminum cylinders (CQB-15645 and CQB-15649) since 2003, which are independent of the reference gases for the NIES scale. The results are shown in Fig. 3. The average changing rates for the whole period, evaluated by linear regression analysis, are $-\mathrm{0.14}±\mathrm{0.06}$ per meg yr−1 for CQB-15645 and $-\mathrm{0.05}±\mathrm{0.06}$ per meg yr−1 for CQB-15649. Therefore, we conclude that the stability of the NIES O2∕N2 scale has been maintained within ±0.2 per meg yr−1, at least during the period of 2003–2016. However, this stability test cannot exclude the possibility that the O2∕N2 ratios of the reference gases drift across all the cylinders rather uniformly. There are several mechanisms that affect the O2∕N2 ratios of the gases within the high-pressure cylinders, including corrosion of the inner surface, leakage, thermal diffusion and gravitational fractionation. Keeling et al. (2007) assessed carefully and comprehensively the influences of those potential mechanisms on the long-term stability of the O2∕N2 ratio of the reference gases and obtained an estimated uncertainty of ±0.4 per meg yr−1. We also treated the reference cylinders, which were kept horizontally in a thermally insulated box, with the greatest care (Tohjima et al., 2008). Therefore, we adopted the value of ±0.4 per meg yr−1 as the long-term drift of the reference gases caused by the above degradation effects. Consequently, we assumed that the total uncertainty of the long-term stability of the O2∕N2 reference scale was ±0.45 per meg yr−1 (or $\left({\mathrm{0.2}}^{\mathrm{2}}+{\mathrm{0.4}}^{\mathrm{2}}{\right)}^{\mathrm{1}/\mathrm{2}}$) in this study.

Figure 3Temporal changes in the O2∕N2 ratio of reference airs in two aluminum cylinders which are independent of the NIES primary reference gases, relative to the NIES O2∕N2 scale. The dashed lines represent the linear regression lines.

## 2.4 Carbon budget calculation

The ocean uptake and land uptake, O and B, respectively, are given by the following equations (Manning and Keeling, 2006; Tohjima et al., 2008):

$\begin{array}{}\text{(6)}& O=\left[-\left({\mathit{\alpha }}_{\mathrm{F}}-{\mathit{\alpha }}_{\mathrm{B}}\right)F-\left(\frac{{S}_{{\mathrm{O}}_{\mathrm{2}}}}{\mathit{\beta }}\right)×\mathrm{\Delta }\mathrm{APO}+{Z}_{\mathrm{eff}}\right]×\frac{\mathrm{1}}{{\mathit{\alpha }}_{\mathrm{B}}},\text{(7)}& B=\left[{\mathit{\alpha }}_{\mathrm{F}}F+\left(\frac{{S}_{{\mathrm{O}}_{\mathrm{2}}}}{\mathit{\beta }}\right)×\mathrm{\Delta }\mathrm{APO}-\left(\frac{{\mathit{\alpha }}_{\mathrm{B}}}{\mathit{\beta }}\right)×\mathrm{\Delta }{X}_{{\mathrm{CO}}_{\mathrm{2}}}-{Z}_{\mathrm{eff}}\right]×\frac{\mathrm{1}}{{\mathit{\alpha }}_{\mathrm{B}}},\end{array}$

where β is the coefficient converting petagrams of carbon to parts per million CO2 in the atmosphere (β=0.470 ppm PgC−1; Tohjima et al., 2008), and Zeff represents the net effect of the air–sea O2 and N2 exchange on the atmospheric O2∕N2 ratio. In these equations, O, B, F and Zeff are given in units of petagrams of carbon per year; ΔAPO in units of per meg per year; and $\mathrm{\Delta }{X}_{{\mathrm{CO}}_{\mathrm{2}}}$ in units of parts per million per year. Note that F and αF include the CO2 emissions associated with cement manufacturing. The values of αF were calculated from the CO2 emission amounts and the O2∕CO2 molar exchange ratios of the individual fuel types (Keeling, 1988). Since αF slightly varies from year to year, the value of αF for the relevant period was used for the carbon budget calculations. The values of Zeff were calculated based on the effects of global ocean warming and anthropogenic nitrogen deposition in accordance with the approach of Keeling and Manning (2014). The details of the Zeff calculation is discussed in the following section.

We used the same dataset of the fossil-fuel-derived CO2 emissions and the global average of the atmospheric CO2 mole fractions as those used by the GCP for the estimation of global carbon budget 2018 (Le Quéré et al., 2018). The fossil CO2 emissions from fossil-fuel combustion and cement production were basically based on the dataset from the CDIAC and other energy statistics (Boden et al., 2017). The change in the atmospheric CO2 burden was calculated based on the global observation by the US National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA/ESRL; Dlugokencky and Tans, 2018). The temporal variations in the fossil CO2 emissions and the atmospheric CO2 accumulation rate are depicted in Fig. 4.

Annual means of APO centered on 1 January were computed by using the same procedure as Tohjima et al. (2008). First, smooth-curve fits to the data were computed in accordance with the methods of Thoning et al. (1989), with a cut-off frequency of 4.6 cycles yr−1. Then the flask APO data were modified to represent the values at the center of the individual months by shifting them in parallel with the smooth-curve fits. This procedure aimed to reduce the influence from biases of the sampling timings within the individual months. The monthly averages were calculated from the modified APO data. When there were no flask data in the monthly time frame, the monthly average of the smooth curve was used. The annual means were calculated from the consecutive 12 monthly averages from July to June of the following year. The standard errors of the differences between the flask data and the smooth-curve fits for the corresponding annual periods were adopted as the uncertainties for the annual averages. The averages and ranges (minimum–maximum) of the errors for the annual means of APO were 0.8 per meg (0.6–1.3 per meg) for HAT and 1.1 per meg (0.9–1.4 per meg) for COI.

In this study, we adopted the value of 1.10 of Severinghaus (1995) for αB, in accordance with a series of previous studies (e.g., Bender et al., 2005; Manning and Keeling, 2006; Ishidoya et al., 2012; Keeling and Manning, 2014; Goto et al., 2017). However, several studies (e.g., Randerson et al., 2006; Worrall et al., 2013), investigating the elemental compositions of organic matters in soil and plants, indicated that the value of 1.1 is rather large for the globally averaged net O2∕CO2 exchange ratio for the terrestrial biosphere. These studies suggest that the value of 1.05 is much more appropriate for αB. Although it is beyond the scope of this study to discuss which value is better for αB, it is useful to mention that the use of 1.05 for αB results in larger decreasing rates of APO by about 5 % and an increase in land sinks and a decrease in ocean sinks by about 0.06 PgC yr−1 on average in the following results. Considering recent reports about the global net O2∕CO2 exchange ratio, Keeling and Manning (2014) revised the uncertainty of αB to be increased, from ±0.05 (Severinghaus, 1995) to ±0.10. Thus, we also adopted ±0.10 for the uncertainty of αB in this study.

Figure 4Temporal changes in fossil-fuel CO2 emissions (red squares), atmospheric CO2 accumulation rate (orange circles) and ocean outgassing effect Zeff divided by land biotic O2∕CO2 exchange rate (1.1; purple diamonds). The 5-year averages of Zeff∕1.1 used for the pentad (5-year) carbon sink calculations are also depicted as purple lines.

## 2.5 Evaluation of outgassing effect (Zeff)

As is discussed in the Introduction, today's ocean is considered to act as a net source of atmospheric O2 because of global ocean warming, which also affects the air–sea N2 exchange. Since the atmospheric O2 change is measured as the change in the atmospheric O2∕N2 ratio, the outgassing effect caused by the global ocean warming, which is denoted as Zgow, should include the influences from not only ocean O2 outgassing but also ocean N2 outgassing. Assuming that the relationship is proportional between the gas fluxes and heat fluxes across the air–sea interface, Manning and Keeling (2006) gave the equation for Zgow as

$\begin{array}{}\text{(8)}& {Z}_{\mathrm{gow}}=\left({\mathit{\gamma }}_{{\mathrm{O}}_{\mathrm{2}}}-\frac{{S}_{{\mathrm{O}}_{\mathrm{2}}}}{{S}_{{\mathrm{N}}_{\mathrm{2}}}}{\mathit{\gamma }}_{{\mathrm{N}}_{\mathrm{2}}}\right)×Q×\phantom{\rule{0.125em}{0ex}}{m}_{\mathrm{c}}×{\mathrm{10}}^{-\mathrm{15}},\end{array}$

where Q represents the changing rate of the global ocean heat storage (in units of J yr−1), ${\mathit{\gamma }}_{{\mathrm{O}}_{\mathrm{2}}}$ and ${\mathit{\gamma }}_{{\mathrm{N}}_{\mathrm{2}}}$ are ratios of the gas flux to heat flux between the air and sea (in units of mol J−1), ${S}_{{\mathrm{N}}_{\mathrm{2}}}$ is the mole fraction of atmospheric nitrogen (${S}_{{\mathrm{N}}_{\mathrm{2}}}=\mathrm{0.7809}$; Tohjima et al., 2005b), and mc is the atomic mass of carbon (mc=12.01). Zgow is given in units of petagrams of carbon per year.

The primary mechanism that affects the air–sea gas exchange is a reduction of gas solubility caused by the increase in the ocean temperature. Therefore, the ratio of the gas flux to heat flux derived from the above thermal effect can be evaluated from the temperature dependence of gas solubility in the seawater and the specific heat of the seawater. Since the air–sea N2 exchange is predominantly driven by the thermal effect, we adopted the estimated ${\mathit{\gamma }}_{{\mathrm{N}}_{\mathrm{2}}}$ of 2.2 nmol J−1 in this study, in accordance with previous studies (Keeling and Garcia, 2002; Manning and Keeling, 2006).

In contrast to the air–sea N2 exchange, the changes in the ocean circulation and ocean primary production also affect the air–sea O2 exchange, as mentioned in the Introduction. Examining the ratio of the seasonal ocean outgassing of O2 to the seasonal ocean heating and the negative linear relationship between the dissolved O2 concentrations corrected for ocean biological processes and the potential temperature in the main thermocline based on archived global observation data, Keeling and Garcia (2002) obtained the estimate of 4.9 nmol J−1 for ${\mathit{\gamma }}_{{\mathrm{O}}_{\mathrm{2}}}$. The value of ${\mathit{\gamma }}_{{\mathrm{O}}_{\mathrm{2}}}$ was also investigated by using ocean biogeochemical models to revise the global carbon budgets based on O2 observations (e.g., Plattner et al., 2002; Bopp et al., 2002). Keeling et al. (2010) summarized the model-based values of ${\mathit{\gamma }}_{{\mathrm{O}}_{\mathrm{2}}}$ ranging from 5.9 to 6.7 nmol J−1. On the other hand, Stendardo and Gruber (2012) examined a huge archived dataset of observations in the northern Antarctic Ocean during the past 5 decades and obtained changing ratios of O2 inventory to heat content of $-\mathrm{4.3}±\mathrm{2.4}$ nmol J−1 in the upper 700 m and $-\mathrm{1.6}±\mathrm{1.9}$ nmol J−1 between 700 and 2750 m. These basin-scale ocean O2 heat changing ratios seem to suggest that the global ocean acts as a net O2 source due to global ocean warming. Therefore, in this study, we used the value of 4.9 nmol J−1 for ${\mathit{\gamma }}_{{\mathrm{O}}_{\mathrm{2}}}$ in accordance with previous studies.

To compute Q, we used the estimates of the world ocean heat content (OHC) based on a variety of oceanographic data (Levitus et al., 2000, 2012). Time series of OHC for the 0–700 and 0–2000 m layers are available from the NOAA's website (https://www.nodc.noaa.gov/OC5/3M_HEAT_CONTENT/, last access: 19 February 2019). In previous carbon budget estimations based on atmospheric O2∕N2 measurements, the values of Q were estimated from the OHC for the 0–700 m layer (e.g., Manning and Keeling, 2006; Tohjima et al., 2008; Ishidoya et al., 2012). Levitus et at. (2012) showed, however, that the ocean heat storage of the 700–2000 m layer contributes to about one-third of the total heat storage of the 0–2000 m layer. Keeling and Manning (2014) estimated the value of Q by considering not only Q for the depths above 700 m but also Q for the depths below 700, which contributed to 30 % of the total value of Q. Therefore, the time series of OHC for the 0–2000 m layer was used in this study. Note that since the annual average of OHC for the 0-2000 m layer is available only after 2005, we used the pentad (5-year) averages before 2005.

In addition to the ocean warming effect, Keeling and Manning (2014) introduced recently another ocean outgassing effect caused by atmospheric deposition of excess anthropogenic nitrogen to the open ocean. The excess nitrogen is considered to enhance the ocean biotic production of organic matter, which is associated with the O2 production. Keeling and Manning (2014) evaluated the anthropogenic nitrogen-induced outgassing as being about 0.1×1014 mol O2 yr−1. Since the outgassing effect caused by the anthropogenic nitrogen deposition is small but rather significant, we adopted the effect as ZanthN, with a magnitude of 0.12 PgC yr−1 ($=\mathrm{0.1}×{\mathrm{10}}^{\mathrm{14}}$ mol O2 yr${}^{-\mathrm{1}}×\mathrm{12.01}$ gC mol−1). Eventually, the total outgassing effect, Zeff, is expressed as the summation of Zgow and ZanthN:

$\begin{array}{}\text{(9)}& {Z}_{\mathrm{eff}}={Z}_{\mathrm{gow}}+{Z}_{\mathrm{anthN}}.\end{array}$

The time series of the annual Zeff divided by αB are depicted in Fig. 4 as purple lines. The value of Zeff/αB ranges from 0.2 to 1.0 Pg yr−1, and the 19-year average for 1998–2016 is 0.52 Pg yr−1. There are not many differences in Zeff between this study and previous studies (e.g., Manning and Keeling, 2006; Keeling and Manning, 2014; Tohjima et al., 2008; Ishidoya et al., 2012; Goto et al., 2017).

3 Results and discussion

## 3.1 Observed O2∕N2, CO2 and APO

The time series of the atmospheric CO2 mole fraction, O2∕N2 ratio and APO of the air samples collected at HAT and COI and aboard cargo ships sailing between 40 S and 30 N in the western Pacific are depicted in Fig. 5 together with the smooth-curve fits. The ship data were binned into 10 latitudinal bands (40–30 S, 30–20 S, 20–10 S, 10–0 S, 0–10 N, 10–20 N and 20–30 N). Note that there are no data gaps with more than 50 days in the time series at HAT and COI, while the time series of the cargo ships have data gaps during the 7-month period from October 2006 to April 2007. The ship data during the 7-month period were significantly contaminated by the inboard air due to the failure of diaphragm of the sampling pump.

Figure 5Time series of the atmospheric CO2 mole fraction (a), O2∕N2 ratio (b) and APO (c) of the flask samples obtained from the NIES flask sampling network shown in Fig. 1. Observed data from COI, HAT and cargo ships operating between 40 S and 30 N were used for the global carbon budget calculation. The time series of CO2, O2∕N2 and APO are offset by 20 ppm, 150 per meg and 100 per meg, respectively, to allow them to be plotted on the same panels. The numbers on the y axis represent the values for the data at HAT.

The temporal variations in the annual means of the atmospheric CO2, O2∕N2 and APO, centered on 1 January, are shown in Fig. 6a, b and c, respectively, where linear trends obtained from least square fitting to the data of HAT were subtracted from the individual time series to emphasize the interannual variations. The standard errors of the annual means for HAT and the 10–0 S bin are depicted as vertical bars for typical examples. Note that the annual means of the atmospheric CO2, O2∕N2 and APO and the corresponding standard errors for HAT, COI and the 10 bins are summarized in Table S1–S9 in the Supplement. The annual means of CO2 and O2∕N2 show latitudinal gradients of northward increase and southward increase, respectively, because the fossil-fuel-derived CO2 emissions and O2 consumptions occur predominantly in the northern midlatitudes. In contrast, the highest values of the annual mean APO were generally observed around the Equator, as previously reported (Battle et al., 2006; Tohjima et al., 2005b, 2012). This equatorial peak is mainly attributed to large-scale air–sea gas exchanges: ingassing in the midlatitudes and high latitudes and outgassing in the equatorial region. Conducting atmospheric inversion analyses based on the APO data from the Scripps observation network, Rödenbeck et al. (2008) suggested anomalous outgassing of APO from the equatorial region during El Niño periods, while Tohjima et al. (2015) found a suppressed equatorial peak during El Niño periods based on the western Pacific observations. Eddebbar et al. (2017) reconciled these conflicting results by predicting the existence of a zonal dipole-like El Niño–Southern Oscillation (ENSO) response in the equatorial Pacific based on several ocean process-based models and an atmospheric transport model. These results suggest that enhanced zonal coverage of the atmospheric observations in the equatorial Pacific is needed to constrain the full basin-scale ENSO response. We can see a considerable suppression of the equatorial peak during the strong 2015–2016 El Niño event in Fig. 6c, which was not reported in Tohjima et al. (2015). Any detailed discussion about the temporal variation of the equatorial peak during the 2015/2016 El Niño event is, however, beyond the scope of this study and will be given elsewhere.

Figure 6Temporal variations of the (a) annual mean CO2, (b) annual mean O2∕N2, (c) annual mean APO and (d) annual changing rate of APO based on the flask samples collected from HAT, COI and cargo ships in the western Pacific (40 S–30 N). The differences in the annual means from the linear trends fitted to the data at HAT are depicted in the figures to emphasize the interannual variations. Vertical bars for the plots of HAT and 10–0 S bin correspond to the standard errors of the annual means.

Figure 6d shows the time series of the annual changes in the annual mean APO, which are the annual changing rates of APO for a 1-year interval (Δt=1 year). As you can see, there are considerable differences in the annual changing rates among the observation sites in the same years; the standard deviations range from 1.6 to 4.4 per meg and the average is 2.8 per meg. We also depict the averages of the annual changing rates of APO of HAT and COI and of all the shipboard data as a thick grey line in Fig. 6d. Note that these average annual changing rates of APO were used for calculation of the global carbon budget in the following sections. The average annual changing rates show also a large interannual variability with a standard deviation of 4.7 per meg yr−1 for the entire observation period.

The differences in the changing rate of APO among the sites in the same years decrease with the increase in the interval for the calculation (Δt) as shown in Fig. 7, where the average (red circles) and the minimum and maximum (red dashed lines) standard deviations of the changing rates are plotted against Δt. The differences among the sites decrease almost inversely with Δt; the average standard deviation for the Δt of 5 years is 0.54 per meg yr−1. The temporal variability in the changing rate also decreases inversely with Δt as depicted in Fig. 7 (blue circles); the standard deviation is reduced to 1.2 per meg yr−1 for the Δt of 5 years. The above results seem to suggest that the temporal variability in the APO fluxes exceeds the spatial variability. As is indicated by Eq. (3), the temporal variability in the APO changing rate should be attributed mostly to those in O and Zeff. Therefore, the above results also indicate that an interval of 5 years could suppress the temporal variability in Zeff to the level of ±1.2 per meg yr−1, which corresponds to a carbon budget of about ±0.5 PgC yr−1.

Table 1Comparison of changing rate of the atmospheric CO2, O2 and APO.

* Average and standard deviation of the changing rates for the three sites (ALT, MLO and SPO) listed in Table 1 of Goto et al. (2017) are given in this table.

Figure 7Relationship between the standard deviation (SD) of the APO changing rate and the time interval to calculate the changing rate. The red circles represent the averages of the standard deviations of the APO changing rates from the different sites for the same year. The dashed lines represent the minimum and maximum of the standard deviations. The blue circles represent the temporal variability in the average APO changing rate of the different sites.

The changing rates of the atmospheric CO2, O2∕N2 and APO for several combinations of time periods and the observed data (HAT, COI and shipboard) are summarized in Table 1. Here, the uncertainties of the changing rates were computed from the uncertainties of the corresponding annual means at both ends of the periods, the estimated uncertainty of the O2∕N2 scale stability (±0.45 per meg yr−1; Sect. 2.3) and the uncertainty in the O2∕N2 span sensitivity (±3 %; see below). The time periods of 2000–2010, 2001–2010 and 2001–2014 were selected to compare the observational results of this study with those of the Scripps Institution of Oceanography (SIO) and Tohoku University (TU; Keeling and Manning, 2014; Ishidoya et al., 2012; Goto et al., 2017). As discussed in the above section, the differences in the long-term changing rates of APO between HAT, COI and shipboard data are less than 0.3 per meg yr−1, while the increasing rates of CO2 and the decreasing rates of O2∕N2 for HAT are slightly larger than those for COI and other sites. The monitoring station of HAT is located at the marginal region of continental East Asia, and the anthropogenic CO2 emissions from China often influence the observations at HAT during winter due to the East Asian monsoon (c.f. Minejima et al., 2012; Tohjima et al., 2010, 2014). Additionally, for the period of 2000–2014, the fossil-fuel-derived CO2 emissions from China show a rapid increase associated with the unprecedented economic growth. These situations may explain the rather large increase in CO2 and decrease in O2∕N2 at HAT. In contrast to CO2 and O2∕N2, the emissions from fossil-fuel combustion and land biotic processes contribute less to the APO variations, resulting in relatively small differences in the long-term APO changing rates among the sites.

It should be noted that the decreasing rates of APO of our study are 0.5–1.1 per meg yr−1 smaller than those of SIO and TU. Except for the differences of the observation sites, we can offer two explanations for the discrepancy. First, the calculation methods of the changing rate adopted by Goto et al. (2017) are different from those adopted in this study, which might partially explain the discrepancy. This is understandable when comparing the changing rates of CO2, O2∕N2 and APO for the individual studies. In this study, the APO changing rates are almost consistent with those calculated from the CO2 and O2∕N2 changing rates according to the APO definition. However, in the study of Goto et al. (2017), the CO2 and O2∕N2 changing rates of Ny-Ålesund give APO decreasing rates of −9.4 per meg yr−1, which is 0.7 per meg yr−1 smaller than the originally reported values. Second, inter-laboratory comparison of flask samples and high-pressure cylinders suggests a possibility that the span sensitivity of the O2∕N2 measurements of NIES is about 3 % lower than that of SIO, which can almost explain the differences in the APO decreasing rates. However, to obtain an accurate conclusion, we need many more studies.

Table 2Comparison of global carbon budgets based on APO with those from GCPa,b.

a Figures are given (in units of per meg yr−1 for ΔAPO, mol mol−1 for αF and PgC yr−1 for the others). b Figures in parentheses represent the uncertainties. c ΔAPO is based on the data from HAT, COI and cargo ships (40 S–30 N). The uncertainty in parentheses includes the uncertainty associated with the observations, stability in the O2∕N2 scale, uncertainty derived from limited sampling and uncertainty in the O2∕N2 span sensitivity (see text). d These figures are computed from the dataset summarized by the GCP. The uncertainties are ±0.2 PgC yr−1 for the atmospheric CO2 and ±5 % for the fossil-fuel emissions, ±0.5 PgC yr−1 for the ocean sinks, and ±0.9 PgC yr−1 for the land sinks (Le Quéré et al., 2018). e The uncertainties for αF are ±0.04 mol mol−1 (Tohjima et al., 2008). f The values of Zeff include both global ocean warming and anthropogenic nitrogen deposition effects, and the uncertainties are assumed to be ±100 % (see text).

## 3.2 Calculation of global carbon budgets

The rates of the global carbon uptake by the ocean and land biosphere were calculated from the average changing rate of APO based on observations at COI, HAT and on cargo ships (40 S–30 N). The results for several time periods are summarized in Table 2 together with the average changing rates of APO, globally averaged atmospheric CO2 accumulation rates, fossil-fuel-derived CO2 emission rates and the ocean outgassing effect divided by αB. The 19-year (1998–2016), 17-year (2000–2016) and 14-year (2003–2016) periods correspond to the individual maximum observation periods for HAT, COI and the western Pacific, respectively. For example, the estimated ocean and land sinks for 2000–2016 were found to be 2.6±0.7 and 1.5±0.9 PgC yr−1, respectively.

The uncertainties in the parameters used for the carbon budget calculation (Eqs. 6 and 7), which are also listed in Table 2, are briefly discussed here. Note that in this study the estimated uncertainties are ±1σ. Since the ocean outgassing effect is rather speculative, we assumed that the values of Zeff for the individual periods had ±100 % uncertainties, in accordance with previous studies (e.g., Manning and Keeling, 2006; Tohjima et al., 2008). We adopted uncertainties of ±5 % for the fossil-fuel-derived CO2 emission rate and ±0.2 PgC yr−1 for the atmospheric CO2 increasing rate from Le Quéré et al. (2018). As for the uncertainties of the observed APO changing rates, we adopted the standard deviations among the sites shown in Fig. 7 (±0.37 per meg yr−1 for longer than 10 years and ±0.54 per meg yr−1 for 5 years). The estimated uncertainty of the O2∕N2 scale stability (±0.45 per meg yr−1) discussed in Sect. 2.3, the uncertainty of the O2∕N2 span sensitivity (±3 %) and the uncertainty in the globally averaged APO associated with the limited atmospheric sampling (±0.2 PgC yr−1) discussed in Nevison et al. (2008) were also included in the calculation of the uncertainties in ΔAPO. The uncertainties of αB and αF were ±0.10 (Keeling and Manning, 2014) and ±0.04 (Tohjma et al., 2008), respectively. Finally, these uncertainties were propagated to the ocean and land sink uncertainties, in accordance with Eqs. (6) and (7).

We compared our global carbon budget estimations with those of the GCP (Global Carbon Budget, 2018) updated by Le Quéré et al. (2018). In the GCP carbon budget assessment, the ocean and land sinks were estimated by combining multiple results from a variety of models, including global ocean biogeochemistry models (GOBMs) and dynamic global vegetation models (DGVMs). Since the sum of the model-based ocean and land sinks was not necessarily balanced with the difference between fossil-fuel emissions and atmospheric accumulation, Le Quéré et al. (2018) listed the discrepancies as budget imbalances. The ocean sinks, the land sinks which are net land sinks computed as the differences between land uptake and emissions associated with land-use change, and the budget imbalances for the corresponding periods are listed in Table 2. Note that the uncertainties of the sinks of the GCP are ±0.5 PgC yr−1 for ocean and ±0.9 PgC yr−1 for land. The carbon sinks of this study and the GCP for the three long periods are consistent with each other; the largest difference in sink strength is 0.35 PgC yr−1, which is smaller than the uncertainties associated with the individual estimations. A 3 % higher span sensitivity of the O2∕N2 measurements, which corresponds to the difference in the span sensitivity between SIO and NIES (Sect. 3.1), would result in an increase and decrease of 0.27 PgC yr−1 in the ocean and land sinks, respectively, for the three long periods. Although these changes would enlarge the differences in sinks between the GCP and this study, the differences are still not significant given the uncertainties of both this study and the GCP.

The carbon budgets for four pentad periods (2000–2004, 2004–2008, 2008–2012 and 2012–2016) are also listed in Table 2. Here, we consider that the 5-year interval effectively reduced the apparent errors caused by the imbalance of the seasonal ocean O2 fluxes, as discussed in Sect. 3.1. Again, the discrepancies of the pentad ocean and land sinks between this study and the GCP are within ±0.5 PgC yr−1, which is also less than the estimated uncertainties. The land sink during 2008–2012 and the ocean sink during 2012–2016 of this study are about 0.5 PgC yr−1 larger than those of the GCP. These discrepancies in the carbon sinks, partly explained by the rather large values of the carbon budget imbalances of the GCP estimation, might give a clue about how to partition the imbalance values between the land and ocean sinks.

Examining the temporal variations in the pentad sink strengths of this study, we found a gradual increase in the ocean sinks for the latter three pentad periods and a rapid increase and decrease in the land sinks for the former and latter two pentad periods, respectively. The pentad averages of the GCP sinks seem to show similar temporal variations: a steady increase in the ocean sinks for the whole period and a rapid increase and decrease in the land sinks for the former and latter two pentad periods, respectively. These results suggest that the carbon sinks for the pentad periods can be used to evaluate the temporal changes. In the following section, we will examine the temporal change in the carbon sinks in more detail.

## 3.3 Temporal change in the carbon sinks

Figure 8 shows the temporal variations in the ocean and land sinks for the annual (dashed red lines) and pentad (red lines) intervals calculated from the average of the APO changing rates based on the observations from HAT, COI and cargo ships in the western Pacific. The uncertainties for the pentad sinks (±1σ), which were calculated as described in the previous section, are shown as grey shading. To clearly understand the effect of Zeff correction, the pentad sinks without corrections are also depicted as purple lines in the figure. Additionally, the annual and pentad sinks of the GCP are also depicted in the figures for comparison. Although the annual sinks show considerable variability, especially for the first several years with peak-to-peak differences of more than 10 PgC yr−1, the variability in the pentad sinks is effectively suppressed. Only the pentad budgets for 2000 show a rather large ocean uptake (3.21 PgC yr−1) and a rather weak land emission (0.56 PgC yr−1), which are depicted as dotted lines. These anomalous values may be explained by the fact that the influences from the considerable drawdown of APO in 2000–2001 cannot be compensated for in the pentad APO changing rate for 2000. Hamme and Keeling (2008) reported that the APO drawdown in 2000–2001, which was also observed in the SIO observations, may be attributed to deep ventilation associated with the unprecedented cooling of the western Pacific, and the variations in the ocean heat content exerted only secondary influence. Therefore, we do not use the anomalous pentad ocean and land sinks for 2000 in the following discussions.

Figure 8Temporal variations in (a) ocean and (b) land biospheric sinks estimated from APO variations of this study (red) and process-based models of GCP (blue). The thin dashed lines represent the annual sinks, and the thick lines represent the pentad sinks. The purple lines represent the pentad sinks based on APO without ocean outgassing correction (Zeff), and the light blue lines represent the sinks of GCP with the imbalance sinks added. The uncertainty associated with the pentad sinks with Zeff corrections is shown as shaded area.

Figure 9Comparison of the temporal variations in the ocean sinks based on the APO data of this study (red), global ocean biogeochemistry models (GOBMs) of GCP (blue), and pCO2 data of Landschützer et al. (2016; light blue) and Rödenbeck et al. (2014; orange). The dashed lines represent the regression lines for the corresponding data during 2001–2014. Note that the pCO2-based ocean sinks are adjusted for the pre-industrial ocean CO2 emissions (±0.78 PgC yr−1) caused by riverine CO2 input to the ocean (Resplandy et al., 2018).

In contrast, the pentad land sinks of both this study and the GCP study show an increasing trend during 2001–2009 followed by a decreasing trend during 2009–2014, although the range of variations of this study is about 2 times larger than that of the GCP. The linear trends for the former period are 0.23±0.04 PgC yr−2 for this study and 0.10±0.03 PgC yr−2 for the GCP, and those for the latter period are $-\mathrm{0.22}±\mathrm{0.04}$ PgC yr−2 for this study and $-\mathrm{0.12}±\mathrm{0.04}$ PgC yr−2 for the GCP. An enhancement of the land uptake during the 2000s has been reported recently by several studies based on atmospheric inversions and biosphere models (Keenan et al., 2016; Kondo et al., 2018; Piao et al., 2018). Although there is an ongoing discussion about the detailed mechanisms of the enhanced net land uptake, the accelerated land uptake may partially explain the stagnation of the growth rate of the atmospheric CO2 in the 2000s despite the increasing anthropogenic CO2 emissions. Examining the atmospheric inversion studies and the previous version of the GCP (Le Quéré et al., 2015), in which the net land uptake was computed as residuals among the other carbon budget components, Piao et al. (2018) found that the linear increasing trend of the net land carbon sink during 1998–2012 was 0.17±0.05 PgC yr−2. The linear trend of this study during 2001–2009 is close to the above value within the uncertainty. Although the corresponding linear trend of the latest GCP estimation is about half that of the present study, the sum of the net land sink and the budget imbalances of the GCP, plotted as light blue lines in Fig. 8, shows a much larger increasing trend, 0.20±0.03 PgC yr−2, which is almost identical to our trend.

The land sinks of both this study and the GCP study exhibit decreasing trends for the period 2009–2014, which are partially compensated for by the steady increase in the ocean uptake. The atmospheric accumulation rate of CO2 significantly increased in 2015 and 2016 (see Fig. 5), when one of the strongest El Niño events occurred. Studies based on atmospheric CO2 observations from stations and by satellite indicated that a reduction in biospheric uptake and an increase in biomass burning contributed to the CO2 increase during the El Niño event (Chatterjee et al., 2017; Patra et al., 2017). The decreasing trend of the pentad land uptake also reflects the change in the global carbon cycle associated with the El Niño event.

Shown as discrepancies between the pentad sinks of the GCP with and without budget imbalances (Fig. 8), the magnitude of the budget imbalances seems to increase after 2007. For the pentad sinks centered on 2007, 2008 and 2009, the ocean and land sinks of this study agree with those of the GCP without and with the budget imbalances, respectively. For the pentad sinks between 2010 and 2014, both the ocean and land sinks of this study are plotted between those of the GCP with and without the budget imbalances. Although we cannot show any definitive partitioning of the budget imbalance between ocean and land sinks because of a rather large uncertainty associated with the sink estimations, the above results seem to suggest that the budget imbalances are allocated to land sinks for the former period and to both sinks for the latter period.

From the above discussions, we feel that a 5-year duration effectively suppresses to some extent the anomalous variations in the carbon budget estimations based on APO, which are considered to be caused by the imbalance of the seasonal air–sea O2 exchange. Probably, the 5-year average suppresses the variability in Zeff to a level of ±0.5 PgC yr−1, as discussed in Sect. 3.1. To reduce uncertainty in the carbon budget estimation, we need more effort to improve the quantification of the net O2 outgassing associated with global ocean warming because the quantification of the Zeff at this state is still very speculative. Applying the approach of Stendardo and Gruber (2012), who examined the long-term changes in dissolved O2 and heat content by using archived oceanographic data of the Atlantic Ocean, to other ocean basins would improve our understanding of the long-term net ocean ratio of the O2 flux to the heat flux.

4 Conclusions

We evaluated the global carbon budgets based on the APO data computed from the O2∕N2 and CO2 of the flask samples collected in the Pacific region since 1997. In the carbon budget calculation, we corrected the ocean and land sinks with the ocean O2 outgassing effect, Zeff, based on the ocean heat increment for the 0–2000 m layer. Eventually, we obtained the following conclusions:

1. The long-term oceanic and land biotic carbon sinks were 2.6±0.7 PgC yr−1 and 1.5±0.9 PgC yr−1, respectively, for a 17-year period (2000–2016), and 2.4±0.7 and 1.9±0.9 PgC yr−1, respectively, for a 14-year period (2003–2016). These long-term carbon sinks agreed well with those of the latest GCP estimation (Le Quéré et al., 2018); the differences of the individual estimations are less than ±0.35 PgC yr−1.

2. The ocean and land sinks for the four pentad (5-year) periods (2000–2004, 2004–2008, 2008–2012 and 2012–2016) of this study also showed good agreement with those of the GCP within a difference of ±0.5 PgC yr−1. The land and ocean sinks of this study showed larger values by about 0.5 PgC yr−1 than those of the GCP for 2008–2012 and 2012–2016, respectively, when rather large carbon budget imbalances (> 0.5 PgC yr−1) were found. Therefore, the discrepancies in the sinks between this study and the GCP might give a clue about how to partition the imbalance values between the land and ocean sinks.

3. Calculating the carbon budgets for the pentad periods consecutively, we examined the changing trend of the ocean and land sinks during a 14-year period (2001–2014). In general, the changing trends of both land and ocean sinks of this study agreed well with those of the GCP, although the range of variations of this study was about 2 times larger than that of the GCP study. The pentad ocean sinks showed an overall increasing trend for the entire period (2001–2014), with a linear increasing rate of 0.08±0.02 PgC yr−2. This increasing rate was about 2 times larger than that for the GCP ocean sinks (0.04±0.01 PgC yr−2) but was consistent with that for the global ocean sinks based on pCO2 observations and interpolation techniques (Landschützer et al., 2016; Rödenbeck et al., 2014). In contrast, the pentad land sinks showed an increasing trend for 2001–2009 and a decreasing trend for 2009–2014. The linear trends of the land sinks for this study and the GCP (in parentheses) were 0.23±0.04 PgC yr−2 (0.10±0.03 PgC yr−2) for the former period and $-\mathrm{0.22}±\mathrm{0.04}$ PgC yr−2 ($-\mathrm{0.12}±\mathrm{0.04}$ PgC yr−2) for the latter period. Enhancement of the land carbon uptake was reported also by previous studies (Keenan et al., 2016; Kondo et al., 2018; Piao et al., 2018). In addition, the recent decreasing trend of the land uptake was found to be partially related to the global carbon cycle variation associated with the strong El Niño event in 2015 and 2016.

Code and data availability
Code and data availability.

All the observed data needed to evaluate the global carbon budget are provided in the tables in the Supplement. Other data are available upon request to the corresponding author (tohjima@nies.go.jp).

Supplement
Supplement.

Author contributions
Author contributions.

YT designed the study and drafted the paper. YT and YH conducted the O2∕N2 measurements and analyzed the observed data. HM and TM organized the flask sampling at the monitoring stations and conducted the CO2 measurements. SN managed the shipboard flask sampling. HM, TM, SN and YH examined the results and provided feedback on the paper. All the authors approved the final paper.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
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.

Acknowledgements
Acknowledgements.

We thank Nobukatsu Oda, Fujio Shimano, Kousei Yumoto, Shigeru Kariya, Tomoyasu Yamada and other members of the Global Environment Forum (GEF) for their continued support in flask sampling of air samples. We would like to acknowledge the assistance of the caretakers and local staff of the Ochiishi and Hateruma monitoring stations as well as the owners, operators and crew of the cargo ships. We are grateful to Hisayo Sandanbata, Eri Matsuura, Yoko Kajita and Motoki Sasakawa of the National Institute for Environmental Studies for their continued support in the CO2 measurements of the flask samples. We also thank Keiichi Katsumata for preparing the working standard gases for the O2∕N2 measurements. We gratefully acknowledge the generous cooperation of Toyofuji Shipping Co. and Kagoshima Senpaku Co. for providing us the opportunity to make the aboard atmospheric observations. We would like to thank the two anonymous referees for helping us to greatly improve the paper.

Financial support
Financial support.

This research has been supported by the Global Environmental Research Coordinate System from the Ministry of the Environment, Japan (grant nos. E0955, E1451).

Review statement
Review statement.

This paper was edited by Rachel Law and reviewed by two anonymous referees.

References

Battle, M., Fletcher, S. M., Bender, M. L., Keeling, R. F., Manning, A. C., Gruber, N., Tnas, P. P., Hendricks, M. B., Ho, D. T., Simonds, C., Mika, R., and Paplawsky, B.: Atmospheric potential oxygen: New observations and their implications for some atmospheric and oceanic models, Global Biogeochem. Cy., 20, GB1010, https://doi.org/10.1029/2005GB002534, 2006.

Bender, M., Tans, P. P., Ellis, J. T., Orchardo, J., and Habfast, K.: A high-precision isotope ratio mass-spectrometry method for measuring the O2∕N2 ratio of air, Geochim. Cosmochim. Ac., 58, 4851–4758, 1994.

Bender, M. L., Ho, D. T., Hendricks, M. B., Mika, R., Battle, M. O., Tans, P. P., Conway, T. J., Sturtevant, B., and Casser, N.: Atmospheric O2∕N2 changes, 1993-2002: Implications for the partitioning of fossil fuel CO2 sequestration, Global Biogeochem. Cy., 19, GB4017, https://doi.org/10.1029/2004GB002410, 2005.

Betts, R. A., Jones, C. D., Knight, J. R., Keeling, R. F., and Kennedy, J. J.: El Niño and a record CO2 rise, Nat. Clim. Change, 6, 806–810, 2016.

Boden, T. A., Marland, G., and Andres, R. J.: Global, regional, and national fossil-fuel CO2 Emissions, Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, U.S. Department of Energy, Oak Ridge, Tenn., USA, available at: https://cdiac.ess-dive.lbl.gov/ftp/ndp030/global.1751_2014.ems (last access: 16 July 2019), 2017.

Bopp, L., Le Quéré, C., Heimann, M., Manning, A. C., and Monfray, P.: Climate-induced oceanic oxygen fluxes: implications for the contemporary carbon budget, Global Biogeochem. Cy., 16, 1022, https://doi.org/10.1029/2001GB001445, 2002.

Chatterjee, A., Cierach, M. M., Sutton, A. J., Feely, R. A., Crisp, D., Eldering, A., Gunson, M. R., O'Dell, C. W., Stephens, B. B., and Schimel, D. S.: Influence of El Niño on atmospheric CO2 over the tropical Pacific Ocean: findings from NASA's OCO-2 mission, Science, 358, eaam5776, https://doi.org/10.1126/science.aam5776, 2017.

DeVries, T., Holzer, M., and Primeau, F.: Recent increase in oceanic carbon uptake driven by weaker upper-ocean overturning, Nature, 542, 215–218, https://doi.org/10.1038/nature21068, 2017.

Dlugokencky, E. and Tans, P.: Trends in atmospheric carbon dioxide, National Oceanic & Atmospheric Administration, Earth System Research Laboratory (NOAA/ESRL), available at: http://www.esrl.noaa.gov/gmd/ccgg/trends/global.html, last access: 1 December 2018.

Eddebbar, Y. A., Long, M. C., Resplandy, L., Rödenbeck, C., Rodgers, K. B., Manizza, M., and Keeling, R. F.: Impacts of ENSO on air-sea oxygen exchange: Observations and mechanisms, Global Biogeochem. Cy., 31, 901–921, https://doi.org/10.1002/2017GB005630, 2017.

Goto, D., Morimoto, S., Ishidoya, S., Aoki, S., and Nakazawa, T.: Terrestrial biospheric and oceanic CO2 uptakes estimated from long-term measurements of atmospheric CO2 mole fraction, δ13C, and δ(O2∕N2) at Ny-Ålesund, Svalbard, J. Geophys. Res.-Biogeo., 122, 1192–1202, https://doi.org/10.1002/2017JG003845, 2017.

Hamme, R. C. and Keeling, R. F.: Ocean ventilation as a driver of interannual variability in atmospheric potential oxygen, Tellus B, 60, 706–717, https://doi.org/10.1111/j.1600-0889.2008.00376.x, 2008.

Ishidoya, S., Morimoto, S, Aoki, S., Taguchi, S., Goto, D., Murayama, S., and Nakazawa, T.: Oceanic and terrestrial biospheric CO2 uptake estimated from atmospheric potential oxygen observed at Ny-Ålesund, Svalbard, and Syow, Antactica, Tellus B, 64, 18924, https://doi.org/10.3402/tellusb.v64i0.18924, 2012.

Keeling, R. F.: Development of an interferometric oxygen analyzer for precise measurement of the atmospheric O2 mole fraction, PhD thesis, Harvard Univ., Cambridge, Mass., USA, 178 pp., 1988.

Keeling, R. F. and Garcia, H. E.: The change in oceanic O2 inventory associated with recent global warming, P. Natl. Acad. Sci. USA, 99, 7848–7853, 2002.

Keeling, R. F. and Manning, A. C.: Studies of recent changes in atmospheric O2 content, in Treatise on Geochemistry, vol. 5, 2nd ed., Elsevier, Amsterdam, 385–404, 2014.

Keeling, R. F. and Shertz, S. R.: Seasonal and interannual variations in atmospheric oxygen and implications for the global carbon cycle, Nature, 358, 723–727, 1992.

Keeling, R. F., Manning, A. C., Paplawsky, W. J., and Cox, A. C.: On the long-term stability of reference gases for atmospheric O2∕N2 and CO2 measurements, Tellus B, 59, 3–14, https://doi.org/10.1111/j.1600-0889.2006.00228.x, 2007.

Keeling, R. R., Kortzinger, A., and Gruber, N.: Ocean deoxygenation in a warming world, Annu. Rev. Mar. Sci., 2, 199–229, https://doi.org/10.1146/annurev.marine.010908.163855, 2010.

Keenan, T. F., Prentice, I. C., Canadell, J. G., Williams, C. A., Wang, H., Raupach, M., and Collatz, G. J.: Recent pause in the growth rate of atmospheric CO2 due to enhance terrestrial carbon uptake, Nat. Commun. 7, 13428, https://doi.org/10.1038/ncomms13428, 2016.

Kondo, M., Ichii, K., Patra, P. K., Poulter, B., Calle, L., Koven, C., Pugh, T. A. M., Kato, E., Harper, A., Zaehle, S., and Wiltshire, A.: Plant regrowth as a driver of recent enhancement of terrestrial CO2 uptake, Geophys. Res. Lett., 45, 4820–4830, https://doi.org/10.1029/2018GL077633, 2018.

Landschützer, P., Gruber, N., Haumann, F. A., Rödenbeck, C., Bakker, D. C. E., van Heuven, S., Hoppema, M., Metzl, N., Sweeney, C., Takahashi, T., Tilbrook, B., and Wanninkhof, R.: The reinvigoration of the Southern Ocean carbon sink, Science, 349, 1221–1224, https://doi.org/10.1126/science.aab2620, 2015.

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

Le Quéré, C., Moriarty, R., Andrew, R. M., Canadell, J. G., Sitch, S., Korsbakken, J. I., Friedlingstein, P., Peters, G. P., Andres, R. J., Boden, T. A., Houghton, R. A., House, J. I., Keeling, R. F., Tans, P., Arneth, A., Bakker, D. C. E., Barbero, L., Bopp, L., Chang, J., Chevallier, F., Chini, L. P., Ciais, P., Fader, M., Feely, R. A., Gkritzalis, T., Harris, I., Hauck, J., Ilyina, T., Jain, A. K., Kato, E., Kitidis, V., Klein Goldewijk, K., Koven, C., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lima, I. D., Metzl, N., Millero, F., Munro, D. R., Murata, A., Nabel, J. E. M. S., Nakaoka, S., Nojiri, Y., O'Brien, K., Olsen, A., Ono, T., Pérez, F. F., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Rödenbeck, C., Saito, S., Schuster, U., Schwinger, J., Séférian, R., Steinhoff, T., Stocker, B. D., Sutton, A. J., Takahashi, T., Tilbrook, B., van der Laan-Luijkx, I. T., van der Werf, G. R., van Heuven, S., Vandemark, D., Viovy, N., Wiltshire, A., Zaehle, S., and Zeng, N.: Global Carbon Budget 2015, Earth Syst. Sci. Data, 7, 349–396, https://doi.org/10.5194/essd-7-349-2015, 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.-I., 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, https://doi.org/10.5194/essd-10-2141-2018, 2018.

Levitus, S., Antonov, J., Boyer, T. P., and Stephens, C.: Warming of the world ocean, Science, 287, 2225–2229, https://doi.org/10.1126/science.287.5461.2225, 2000.

Levitus, S., Antonov, J. I., Boyer, T. P., Baranova, O. K., Garcia, H. E., Locarnini, R. A., Mishonov, A. V., Reagan, J. R., Yarosh, E. S., and Zweng, M. M.: World ocean heat content and thermosteric sea level change (0–2000m), 1955–2010, Gephys. Res. Lett., 39, L10603, https://doi.org/10.1029/2012GL051106, 2012.

Machida, T., Tohjima, Y., Katsumata, K., and Mukai, H.: A new CO2 calibration scale based on gravimetric one-step dilution cylinders in National Institute for Environmental Studies-NIES 09 CO2 scale. Paper presented at: Report of the 15th WMO Meeting of Experts on Carbon Dioxide Concentration and Related Tracer Measurement Techniques; September 2009; Jena, Germany, WMO/GAW Rep. 194, edited by: Brand, W., 165–169, WMO, Geneva, Switzerland, 2011.

Mannning, A. C. and Keeling R. F.: Global oceanic and biotic carbon sinks from the Scripps atmospheric oxygen flask sampling network, Tellus, 58, 95–116, 2006.

Manning, A. C., Keeling, R. F., and Severinghaus, J. P.: Precise atmospheric oxygen measurements with a paramagnetic oxygen analyzer, Glob. Biogeochem. Cy., 13, 1107–1115, 1999.

McKinley, G. A., Follows, M. J., Marshall, J., and Fan, S. M.: Interannual variability of air-sea O2 fluxes and the determination of CO2 sinks using atmospheric O2∕N2, Geophys. Res. Lett., 30, 1101, https://doi.org/10.1029/2002GL016044, 2003.

Minejima, C., Kubo, M., Tohjima, Y., Yamagishi, H., Koyama, Y., Maksyutov, S., Kita, K., and Mukai, H.: Analysis of ΔO2∕ΔCO2 ratios for the pollution events observed at Hateruma Island, Japan, Atmos. Chem. Phys., 12, 2713–2723, https://doi.org/10.5194/acp-12-2713-2012, 2012.

Nevison, C. D., Mahowald, N. M., Doney, S. C., Lima, I. D., and Cassar, N.: Impact of variable air-sea O2 and CO2 fluxes on atmospheric potential oxygen (APO) and land-ocean carbon sink partitioning, Biogeosciences, 5, 875–889, https://doi.org/10.5194/bg-5-875-2008, 2008.

Patra, P. K., Crisp, D., Kaiser, J. W., Wunch, D., Saeki, T., Ichii, K., Sekiya, T., Wennberg, P. O., Feist, D. G., Pollard, D. F., Griffith, D. W., Velazco, V. A., De Maziere, M., Sha, M. K., Roehl, C., Chatterjee, A., and Ishijima, K.: The orbiting carbon observatory (OCO-2) tracks 2–3 peta-gram increase in carbon release to the atmosphere during the 2014–2016 El Niño, Sci. Rep., 7, 13567, https://doi.org/10.1038/s41598-017-13459-0, 2017.

Piao, S., Huang, M., Liu, Z., Wang, X., Ciais, P., Canadell, J. G., Wang, K., Bastos, A., Friedlingstein, P., Houghton, R. A., Le Quéré, C., Liu, Y., Myneni, R., Peng, S., Pongratz, J., Sitch, S., Yan, T., Wang, Y., Zhu, Z., Wu, D., and Wang, T.: Lower land-use emissions responsible for increased net land carbon sink during the slow warming period, Nat. Geosci., 11, 739–743, https://doi.org/10.1038/s41561-018-0204-7, 2018.

Plattner, G. K., Joos, F., and Stocker, T. F.: Revision of the global carbon budget due to changing air-sea oxygen fluxes, Global Biogeochm. Cy., 16, 1096, https://doi.org/10.1029/2001GB001746, 2002.

Randerson, J. T., Masiello, C. A., Still, C. J., Rahn, T., Poorter, H., and Field, C. B.: Is carbon within the global terrestrial biosphere becoming more oxidized? Implications for trends in atmospheric O2, Glob. Change Biol., 12, 260–271, https://doi.org/10.1111/j.1365-2486.2006.01099.x, 2006.

Resplandy, L., Keeling, R. F., Rödenbeck, C., Stephens, B. B., Khatiwala, S., Rodgers, K. B., Long, M. C. Long, Bopp, L., and Tans, P. P.: Revision of global carbon fluxes based on a reassessment of oceanic and riverine carbon transport, Nat. Geosci., 11, 504–509, https://doi.org/10.1038/s41561-018-0151-3, 2018.

Rödenbeck, C., Le Quéré, C., Heimann, M., and Keeling, R. F.: Interannual variability inoceanic biogeochemical processes inferred by inversion of atmospheric O2∕N2 and CO2 data, Tellus B, 60, 685–705, https://doi.org/10.1111/j.1600-0889.2008.00375.x, 2008.

Rödenbeck, C., Bakker, D. C. E., Metzl, N., Olsen, A., Sabine, C., Cassar, N., Reum, F., Keeling, R. F., and Heimann, M.: Interannual sea–air CO2 flux variability from an observation-driven ocean mixed-layer scheme, Biogeosciences, 11, 4599–4613, https://doi.org/10.5194/bg-11-4599-2014, 2014.

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, https://doi.org/10.5194/bg-12-7251-2015, 2015.

Severinghaus, J. P.: Studies of the terrestrial O2 and carbon cycles in sand dunes gases and in Biosphere 2, PhD thesis, Columbia Univ., New York, NY, pp. 148, 1995.

Stendardo, I. and Gruber, N.: Oxygen trends over five decades in the North Atlantic, J. Geophys. Res., 117, C11004, https://doi.org/10.1029/2012JC007909, 2012.

Stephens, B. B., Keeling, R. F., Heimann, M., Six, K. D., Murnane, R., and Caldeira, K.: Testing global ocean carbon cycle models using measurements of atmospheric O2 and CO2 concentration, Glob. Biogeochem. Cy., 12, 213–230, 1998.

Stephens, B. B., Keeling, R. F., and Paplawsky, W.: Shipboard measurements of atmospheric oxygen using a vacuum-ultraviolet absorption technique, Tellus B, 55, 857–878, 2003.

Stephens, B. B., Bakwin, P. S., Tans, P. P., Teclaw, R. M., and Baumann, D. D.: 2006, Application of a differential fuel-cell analyzer for measuring atmosppheric oxygen variations, J. Atmos. Ocean Tech., 24, 82–94, https://doi.org/10.1175/JTECH1959.1, 2007.

Thoning, K. W., Tans, P. P., and Komhyr, W. D.: Atmospheric carbon dioxide at Mauna Loa Observatory 2. Analysis of the NOAA GMCC data, 1974–1985, J. Geophys. Res., 94, 8549–8565, 1989.

Tohjima, Y.: Method for measuring changes in the atmospheric O2∕N2 ratio by a gas chromatograph equipped with a thermal conductivity detector, J. Geophys. Res., 105, 14,575–14,584, 2000.

Tohjima, Y., Mukai, H., Machida, T., and Nojiri, Y.: Gas-chromatographic measurements of the atmospheric oxygen/nitrogen ratio at Hateruma Island and Cape Ochi-ishi, Japan, Geophys. Res. Lett., 30, 1653, https://doi.org/10.1029/2003GL017282, 2003.

Tohjima, Y., Machida, T., Watai, T., Akama, I., Amari, T., and Moriwaki, Y.: Preparation of gravimetric standards for measurements of atmospheric oxygen and re-evaluation of atmospheric oxygen concentration, J. Geophys. Res., 110, D11302, https://doi.org/10.1029/2004JD005595, 2005a.

Tohjima, Y., Mukai, H., Machida, T., Nojiri, Y., and Gloor, M.: First measurements of the latitudinal atmospheric O2 and CO2 distributions across the western Pacific, Geophys. Res. Lett., 32, L17805, https://doi.org/10.1029/2005GL023311, 2005b.

Tohjima, Y., Mukai, H., Nojiri, Y., Yamagishi, H., and Machida, T.: Atmospheric O2∕N2 measurements at two Japanese sites: Estimation of global oceanic and land biotic carbon sinks and analysis of the variations in atmospheric potential oxygen (APO), Tellus B, 60, 213–225, https://doi.org/10.1111/j.1600-0889.2007.00334.x, 2008.

Tohjima, Y., Mukai, H., Hashimoto, S., and Patra, P. K.: Increasing synoptic scale variability in atmospheric CO2 at Hateruma Island associated with increasing East-Asian emissions, Atmos. Chem. Phys., 10, 453–462, https://doi.org/10.5194/acp-10-453-2010, 2010.

Tohjima, Y., Minejima, C., Mukai, H., Machida, T., Yamagishi, H., and Nojiri, Y.: Analysis of seasonality and annual mean distribution of atmospheric potential oxygen (APO) in the Pacific region, Glob. Biogeochem. Cy., 26, GB4008, https://doi.org/10.1029/2011GB004110, 2012.

Tohjima, Y., Kubo, M., Minejima, C., Mukai, H., Tanimoto, H., Ganshin, A., Maksyutov, S., Katsumata, K., Machida, T., and Kita, K.: Temporal changes in the emissions of CH4 and CO from China estimated from CH4∕CO2 and CO∕CO2 correlations observed at Hateruma Island, Atmos. Chem. Phys., 14, 1663–1677, https://doi.org/10.5194/acp-14-1663-2014, 2014.

Tohjima, Y., Terao, Y., Mukai, H., Machida, T., Nojiri, Y., and Maksyutov, S., ENSO-related variability in latitudinal distribution of annual mean atmospheric potential oxygen (APO) in the equatorial Western Pacific, Tellus B, 67, 25869, https://doi.org/10.3402/tellusb.v67.25869, 2015.

Worrall, F., Clay, G. D., Masiello, C. A., and Mynheer, G.: Estimating the oxidative ratio of the global terrestrial biosphere carbon, Biogeochemistry, 115, 23–32, https://doi.org/10.1007/s10533-013-9877-6, 2013.