Vertical Profiles of CO 2 and the Latitudinal Partitioning of Carbon Fluxes Britton B. Stephens*, K. R. Gurney, P. P. Tans, C. Sweeney, W. Peters, L. Bruhwiler, P. Ciais, M. Ramonet, P. Bousquet, T. Nakazawa, S. Aoki, T. Machida, G. Inoue, N. Vinnichenko, J. Lloyd, A. Jordan, M. Heimann, O. Shibistova, R. L. Langenfelds, L. P. Steele, R. J. Francey, A. S. Denning *National Center for Atmospheric Research, Boulder, Colorado, USA, [email protected], see reference for co-author affiliations. Overview: Measurements of midday vertical atmospheric CO 2 distributions reveal annual-mean vertical CO 2 gradients which are inconsistent with atmospheric models that estimate a large transfer of terrestrial carbon from tropical to northern latitudes. The 3 models that most closely reproduce the observed annual-mean vertical CO 2 gradients estimate weaker northern uptake of -1.5 PgCyr -1 and weaker tropical emission of +0.1 PgCyr -1 compared to previous consensus estimates of -2.4 and +1.8 PgCyr -1 respectively. This suggests a smaller role for northern terrestrial uptake of industrial CO 2 emissions than previously thought and, after subtracting land use emissions, that tropical ecosystems may currently be strong sinks for CO 2 . Adding airborne flask data to existing inverse models without first fixing their transport biases will not improve their flux estimates. Future inverse models must be validated against available metrics, including vertical profiles of CO 2 , SF 6 , and other tracers. 1. Airborne flask sampling for CO 2 was first conducted five decades ago by Dave Keeling. Although funding limitations ended the program, these early measurements revealed the attenuation of the seasonal CO 2 cycle with altitude and latitude. Although considerable uncertainty remains, this result implies that northern countries may not be able to rely on their forests to offset their industrial emissions as much as hoped, but also that tropical countries can help to mitigate CO 2 increases by preserving their forests. 4. Observed northern-hemisphere average profiles reveal stronger vertical mixing in summer than winter. We extracted and averaged model-predicted profiles to match the observations. After optimization to marine-boundary-layer surface data, these models show large variations with generally overestimated summer mixing and both over and underestimated winter mixing. 6. Reference: Stephens, B.B., K.R. Gurney, P.P. Tans, C. Sweeney, W. Peters, L. Bruhwiler, P. Ciais, M. Ramonet, P. Bousquet, T. Nakazawa, S. Aoki, T. Machida, G. Inoue, N. Vinnichenko, J. Lloyd, A. Jordan, M. Heimann, O. Shibistova, R.L. Langenfelds, L.P. Steele, R.J. Francey, A.S. Denning, Weak northern and strong tropical land carbon uptake from vertical profiles of atmospheric CO 2 , Science, 316, 1732-1735, 2007. 20 -15 10 -10 10 -10 0 -5 20 -15 10 -10 10 -10 0 -5 Northern Hemisphere sites include Briggsdale, Colorado, USA (CAR); Estevan Point, British Columbia, Canada (ESP); Molokai Island, Hawaii, USA (HAA); Harvard Forest, Massachusetts, USA (HFM); Park Falls, Wisconsin, USA (LEF); Poker Flat, Alaska, USA (PFA); Orleans, France (ORL); Sendai/Fukuoka, Japan (SEN); Surgut, Russia (SUR); and Zotino, Russia (ZOT). Southern Hemisphere sites include Rarotonga, Cook Islands (RTA) and Bass Strait/Cape Grim, Australia (AIA). Map of airborne flask sampling locations 3. Existing airborne sampling includes these 12 sites with records over 5 years long and several with records over 20 years long, with reasonable global and continental coverage. However, the use of these data in global inverse models has been limited. Altitude-time CO 2 contour plots for all sampling locations Observations Model Predictions J ul-Aug-S ep CO 2 (ppm) -12 -8 -4 0 Pressure (mbar) 250 400 550 700 850 1000 Annual means CO 2 (ppm) -2 0 2 4 Jan-Feb-Mar CO 2 (ppm) -4 0 4 8 12 Altitude (km) 0.0 1.0 2.5 4.0 6.0 9.0 A B CAR ESP HAA HFM LEF PFA RTA AIA ORL SEN SUR ZOT NH mean C J ul-Aug-S ep CO 2 (ppm) -8 -6 -4 -2 0 Pressure (mbar) 250 400 550 700 850 1000 Annual means CO 2 (ppm) -1 0 1 2 3 4 Jan-Feb-Mar CO 2 (ppm) 2 4 6 8 Altitude (km) 0.0 1.0 2.5 4.0 6.0 9.0 A B Modeled NH mean Observed NH mean C J ul-Aug-S ep 1 km - 4 km CO 2 (ppm) -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 P os t-invers ion flux (P gC yr -1 ) -14 -12 -10 -8 -6 -4 -2 1 2 3 4 5 6 7 8 9 A B C 1 2 3 45 6 7 8 9 A B C Annual means 1 km - 4 km CO 2 (ppm) 0.0 0.5 1.0 1.5 2.0 2.5 -6 -4 -2 0 2 4 6 1 2 3 4 6 7 8 9 A B C 1 2 3 4 5 6 7 8 9 A B C Jan-Feb-Mar 1 km - 4 km CO 2 (ppm) 1.0 1.5 2.0 2.5 3.0 3.5 4.0 0 2 4 6 8 10 12 2 3 4 5 6 8 9 A B C 1 2 3 4 5 6 8 9 A B C A B r 2 = 0.11 r 2 = 0.54 r 2 = 0.54 r 2 = 0.27 r 2 = 0.76 r 2 = 0.67 1 - C Northern land 1 - C Tropical land Observed C 5 1 7 7 5. Seasonal and annual-mean comparisons show that no model accurately reproduces the observations at all times of year. We use the annual-mean 1 - 4 km gradient as a metric to select three models, which we then average to produce the new flux estimates. Observed profiels at all sites (colors) and averaged across northern hemisphere sites (black). Northern-hemisphere average modeled (red) and observed profiles (gray). Model estimated fluxes plotted versus predicted vertical profiles. The observed vertical profiles and uncertainties are indicated by gray bars. 2. TransCom 3 Level 2 estimated fluxes for northern and tropical aggregated regions vary across models by as much as 5 billion tons of carbon per year. These variations are systematically related to how the models transport CO 2 in and out of the northern-hemisphere lower troposphere. Without validation, the across-model mean is no more robust than any individual estimate. Post-inversion Flux (PgCyr -1 ) -4 -2 0 2 4 1 2 3 4 5 6 7 8 A B C Tropical Total, r 2 = 0.61 Northern Total, r 2 = 0.78 9 -4 -2 0 2 4 Pre-inversion 53°N-Eq. Terrestrial Response (ppm) -0.5 0.0 0.5 1.0 1.5 2.0 Tropical Land, r 2 = 0.63 Northern Land, r 2 = 0.60 A B Post-inversion Flux (PgCyr -1 ) TM3 C TM2 B NIRE A NIES 9 MATCH.MACCM2 8 MATCH.NCEP 7 MATCH.CCM3 6 JMA 5 UCI 4 UCB 3 GCTM 2 CSU 1 Model Name Model TM3 C TM2 B NIRE A NIES 9 MATCH.MACCM2 8 MATCH.NCEP 7 MATCH.CCM3 6 JMA 5 UCI 4 UCB 3 GCTM 2 CSU 1 Model Name Model 349.5 350.0 350.5 351.0 351.5 352.0 352.5 353.0 353.5 354.0 -90 -70 -50 -30 -10 10 30 50 70 90 Transcom 3 Neutral Biosphere Response