CarbonTracker CT2019

CarbonTracker is a CO2 measurement and modeling system developed by NOAA to keep track of sources (emissions to the atmosphere) and sinks (removal from the atmosphere) of carbon dioxide around the world. CarbonTracker uses atmospheric CO2 observations from a host of collaborators and simulated atmospheric transport to estimate these surface fluxes of CO2. The current release of CarbonTracker, CT2019, provides global estimates of surface-atmosphere fluxes of CO2 from January 2000 through December 2018.

Global CO 2 budget

From 2000 through 2018, CO2 emissions to the atmosphere from burning of fossil fuels rose from 6.7 PgC yr-1 to 10.2 PgC yr-1 (1 petagram of carbon is 1015 gC, or 1 billion metric tons C, or 3.67 billion metric tons CO2). Global fossil fuel emissions have increased steadily year upon year, with the exception of 2009 following the global economic recession and 2014-2016 when emissions held nearly constant (Figure 1). Fossil fuel emissions are concentrated in areas with high population density and economic activity, and emissions inventory information used in CT2019 indicates that 82% of fossil fuel emissions come from the industrialized northern extratropics.

The other major source of CO2 is wildfires, which in CT2019 add an additional 2.0-2.6 PgC yr-1 to the atmosphere. In contrast to fossil fuel emissions, wildfire CO2 comes principally from tropical and southern land. 87% of wildfire emissions in CT2019 are in those regions.

Offsetting these sources are natural sinks on land and in the ocean. Together, these sinks absorb about half the anthropogenic CO2 emitted into the atmosphere. Over the period 2001-2018, the CT2019 global sum of "natural" fluxes (fire emissions, the land biosphere sink, and the ocean sink) is 49% of fossil fuel emissions over the same time period. The atmospheric CO2 growth rate would be about twice the observed rate without these sinks. CarbonTracker is designed to identify these sinks in order to better understand the mechanisms behind them.

According to CT2019, the world oceans absorb 1.4 to 4.1 PgC yr-1 . This natural sink exists as a direct result of increasing atmospheric CO2 concentrations, because dissolved carbon concentrations in the ocean increase to reach equilibrium with the atmosphere. However, the large-scale circulation of the ocean and biological, physical, and chemical carbon cycling cause there to be a source of carbon to the atmosphere in the tropics. In CT2019, this natural source of between 0.3 and 0.9 PgC yr-1 in the tropics is offset by large extratropical sinks of 2.3 - 4.4 PgC yr-1 .

The terrestrial biosphere is also a net sink of CO2, due mainly to two processes. These are CO2 fertilization, in which plants grow faster since they can more easily acquire carbon dioxide for photosynthesis, and the effects of human land-use practices, including fertilization, irrigation, fire suppression, and recovery from past land use. CarbonTracker attempts to locate these land sinks spatially and temporally, so hypotheses about their mechanisms can be tested. CT2019 finds widely-scattered terrestrial CO2 sinks, with significant absorption of carbon dioxide by northern temperate and boreal regions (1.7 - 2.8 PgC yr-1 , about 58% of the global land total of 2.5 - 4.8 PgC yr-1 ).

In this text, fluxes reported are reported as ranges when possible, to provide some context for how well CarbonTracker constrains the long-term average flux. These ranges are computed as the minimum and maximum values from the sequence of annual CT2019 flux estimates, often excluding the year 2000 as we consider it a spin-up year in which estimats may not be representative of the actual long-term mean. A large range indicates a high degree of interannual variability in the flux estimate.

CO 2 sources and sinks over North America

CT2019 results indicate that North America ecosystems have been a net sink of 0.6 (range 0.1 - 1.0) PgC yr-1 over the period 2001-2018. This natural sink offsets about one-third of the emissions of about 1.8 PgC yr-1 from the burning of fossil fuels in the U.S.A., Canada and Mexico combined.

Whereas North American fossil emissions are generally steady over the CarbonTracker period, ranging between 1.7 and 1.9 PgC yr-1 , the amount of CO2 taken up by the North American biosphere varies significantly from year to year (Figure 2, bottom panel). In terrestrial biosphere models, inter-annual variability in land uptake can be related to anomalies in large-scale temperature and precipitation patterns. While the CarbonTracker period of analysis is relatively short compared to the dynamics of slowly-changing pools of biospheric carbon, episodes of extremes in net ecosystem exchange (NEE) have been associated with climatic anomalies (see e.g. Peters et al., 2007). CT2019 annual estimates of the land sink over temperate N. America are clearly related to continental-scale drought intensity (Figure 2). It is interesting to note that the inferred year-to-year variabilty (the "range") of land uptake is actually as big as the average sink itself.

Widespread droughts in the U.S. west and Canada during 2002, 2006, 2008, 2011, and 2012 resulted in relatively small annual uptake by terrestrial ecosystems in temperate North America (Figure 2). In these years, land ecosystems accounted for a sink of only about 0.2 PgC yr-1 . This is about half the sink of an "average" year, in which these same land ecosystems remove about 0.4 PgC yr-1 .

Spatial distribution of North American surface fluxes

CarbonTracker flux estimates include sub-continental patterns of sources and sinks coupled to the distribution of dominant ecosystem types across the continent (Figure 3). We have greater confidence in countrywide totals than in estimates of regional sources and sinks, but we expect that such finer-scale estimates will become more robust with future expansion of the CO2 observing nework. Our results indicate that the sinks are mainly located in the agricultural regions of the U.S. and Canadian midwest, and boreal forests in Canada.

Word of caution about high-resolution biological flux maps Figure 3 shows estimated fluxes by ecoregion. While we also provide flux maps and data with a finer 1° x 1° spatial resolution, for example on our flux maps pages, these ecoregions define the actual scales at which CarbonTracker operates. With the present observing network, the detailed one-degree fluxes should not be interpreted as quantitatively meaningful for each block. Any within-ecoregion patterns come directly from results of the terrestrial biosphere model. Part of this high-resolution patterning comes from variations of temperature, precipitation, light, plant species, and soil type over each ecoregion. To spread the influence of measurements from the sparse observation network, CarbonTracker makes adjustments uniformly over an entire ecoregion. These adjustments scale the net ecosystem flux of CO2 predicted by the terrestrial biosphere model by the same factor across each ecoregion. Thus we caution that the 1° x 1° spatial detail in CarbonTracker land fluxes is based on the simulations of the terrestrial biosphere model and the assumption of large-scale ecosystem coherence. This has not been verified by observations.

The CarbonTracker observing system

CarbonTracker surface flux estimates are optimally consistent with atmospheric CO2 observations from the GLOBALVIEWplus-5.0 ObsPack, comprising 460 time series datasets from around the world using a variety of measurement techniques and platforms (Table 1, Figure 4). These observation are contributed by collaborators from 54 different laboratories. Most of the CO2 observational data can be accessed by downloading the GLOBALVIEWplusv5.0 ObsPack, or if modeled observations are also required, the CT2019 ObsPack. More information on CO2 measurements used in CT2019 can be found in the observations documentation.

Observation Number of Number of Not For Assimilation Observations Type Datasets Observations Assimilation Total Accepted Rejected Withheld surface-insitu 155 7 317 370 5 717 099 1 600 271 1 523 688 23 042 53 541 surface-flask 110 66 296 11 121 55 175 51 567 852 2 756 surface-pfp 27 34 503 17 488 17 015 16 826 189 0 tower-insitu 102 12 890 823 11 657 010 1 233 813 1 149 294 25 342 59 177 aircraft-pfp 44 70 452 10 650 59 802 57 862 1 678 262 aircraft-flask 4 6 387 3 306 3 081 2 862 67 152 aircraft-insitu 32 9 788 210 9 687951 100 259 99 994 265 0 shipboard-flask 3 3 135 501 2 634 2 374 130 130 shipboard-insitu 8 1 365 297 186 340 1 178 957 1 138 383 40 574 0 Total 485 31 542 473 27 291 466 4 251 007 4 042 850 92 139 116 018

Calculated time-dependent CO 2 fields throughout the global atmosphere

A "byproduct" of the data assimilation system, once sources and sinks have been estimated, is that the mole fraction of CO2 is calculated everywhere in the model domain and over the entire 2000-2016 time period, based on the optimized source and sink estimates (Figure 1). As a check on model transport properties and CarbonTracker inversion performance, calculated CO2 mole fractions are regularly compared with measurements from aircraft campaign datasets taken by NOAA/ESRL and collaborators. These independent samples are not used to estimate fluxes in CarbonTracker, but rather set aside for cross-validation.

Since CarbonTracker simulates CO2 throughout the entire atmospheric column, the model atmosphere can be sampled like satellite (GOSAT and OCO-2) and ground-based remote sensing instrument (TCCON) retrievals of CO2. Examples of our agreement with the latter can be found on our TCCON page.

Flux uncertainties

It is important to note that at this time the uncertainty estimates for CarbonTracker sources and sinks are themselves quite uncertain. They have been derived from the mathematics of the ensemble data assimilation system, which requires several educated guesses for initial uncertainty estimates. The paper describing CarbonTracker (Peters et al. (2007), Proc. Nat. Acad. Sci. vol. 104, p. 18925-18930) presents different uncertainty estimates based on the sensitivity of the results to 14 alternative yet plausible ways to construct the CarbonTracker system. For example, the 14 realizations produce a range of the net annual average terrestrial emissions in North America of -0.40 to -1.01 PgC yr-1 (negative emissions indicate a sink). The procedure is described in the Supporting Information Appendix to that paper, which is freely downloadable from the PNAS web site.

Furthermore, the estimates do not take into account several additional factors noted below. The calculation is set up for sources and sinks to slowly revert, in the absence of observational data, to first guesses of net ecosystem exchange, which are close to zero on an annual basis. This set-up may result in a bias. Also due to the sparseness of measurements, we have had to assume coherence of ecosystem processes over large distances, giving existing observations perhaps an undue amount of weight. The process model for terrestrial photosynthesis and respiration was very basic, and will likely be greatly improved in future releases of CarbonTracker. Easily the largest single annual average source of CO2 is emissions from fossil fuel burning, which are currently not estimated by CarbonTracker. We use estimates from emissions inventories (economic accounting) and subtract the CO2 mole fraction signatures of those fluxes from observations. As a result, the biosphere and ocean fluxes estimated by CarbonTracker inherit error from the assumed fossil fuel emissions. While these emissions inventories may have a small relative error on global scales (perhaps 5 or 10%), any such bias translates into a larger relative error in the annual average ecosystem sources and sinks, since those fluxes have smaller magnitudes. We expect to add a process model of fossil fuel combustion in future releases of CarbonTracker. Finally, additional measurement sites are expected to lead to the greatest improvements, especially to more robust and specific source/sink results at smaller spatial scales.

Consistency of modeled and observed atmospheric CO 2 growth rates

Global atmospheric CO2 growth rates inferred directly from observed carbon dioxide at marine surface sites are consistent with those modeled by CarbonTracker, both in their average values and in their year-to-year variations (Figure 6). These global growth rates hovered at around 4 PgC yr-1 , or around 1.9 ppm yr-1 (Figure 6) in the decade of the 2000s, but have increased in recent years. The observed global growth rates in 2015 and 2016 are particularly high, at 6.3 PgC yr-1 and 6.0 PgC yr-1 (just under 3.0 ppm yr-1 ). This anomaly is believed to be related to the large 2015-16 El Niño (see below).

The 2015-2016 El Niño

The large El Niño of 2015 and 2016 was responsible for two years of sustained record global growth rates in atmospheric CO2 (Figure 6). With its extended assimilation window and expanded network of measurement data, CT2019 simulates these growth rate anomalies quite well. In terms of surface fluxes, CT2019 finds that land ecosystems in the tropics were responsible for a CO2 emissions anomaly of about 1.2 PgC during this event (see Figure 7).

2007 CarbonTracker PNAS publication

CarbonTracker is a NOAA contribution to the North American Carbon Program

Other CarbonTrackers