Climatological mean and decadal change in surface ocean pCO2, and net sea–air CO2 flux over the global oceans

A climatological mean distribution for the surface water pCO2 over the global oceans in non-El Nino conditions has been constructed with spatial resolution of 4°(latitude) x 5°(longitude) for a reference year 2000 based upon about 3 million measurements of surface water pCO2 obtained from 1970 to 2007. The database used for this study is about 3 times larger than the 0.94 million used for our earlier paper [Takahashi et al., 2002. Global sea-air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects. Deep-Sea Res. 11, 49, 1601-1622]. A time-trend analysis using deseasonalized surface water pCO2 data in portions of the North Atlantic, North and South Pacific and Southern Oceans (which cover about 27% of the global ocean areas) indicates that the surface water pCO2 over these oceanic areas has increased on average at a mean rate of 1.5 µatm y-1 with basin-specific rates varying between 1.2 +/- 0.5 and 2.1 +/- 0.4 µatm y-1. A global ocean database for a single reference year 2000 is assembled using this mean rate for correcting observations made in different years to the reference year. The observations made during El Nino periods in the equatorial Pacific and those made in coastal zones are excluded from the database. Seasonal changes in the surface water pCO2 and the sea-air pCO2 difference over four climatic zones in the Atlantic, Pacific, Indian and Southern Oceans are presented. Over the Southern Ocean seasonal ice zone, the seasonality is complex. Although it cannot be thoroughly documented due to the limited extent of observations, seasonal changes in pCO2 are approximated by using the data for under-ice waters during austral winter and those for the marginal ice and ice-free zones. The net air-sea CO2 flux is estimated using the sea-air pCO2 difference and the air-sea gas transfer rate that is parameterized as a function of (wind speed)2 with a scaling factor of 0.26. This is estimated by inverting the bomb 14C data using Ocean General Circulation models and the 1979-2005 NCEP-DOE AMIP-II Reanalysis (R-2) wind speed data. The equatorial Pacific (14°N-14°S) is the major source for atmospheric CO2, emitting about +0.48 Pg-C y-1, and the temperate oceans between 14° and 50° in the both hemispheres are the major sink zones with an uptake flux of -0.70 Pg-C y-1 for the northern and -1.05 Pg-C y-1 for the southern zone. The high-latitude North Atlantic, including the Nordic Seas and portion of the Arctic Sea, is the most intense CO2 sink area on the basis of per unit area, with a mean of -2.5 tons-C month-1 km-2. This is due to the combination of the low pCO2 in seawater and high gas exchange rates. In the ice-free zone of the Southern Ocean (50°-62°S), the mean annual flux is small (-0.06 Pg-C y-1) because of a cancellation of the summer uptake v flux with the winter release Of CO2 caused by deepwater upwelling. The annual mean for the contemporary net v uptake flux over the global oceans is estimated to be -1.6+/-0.9 Pg-C y-1, which includes an undersampling correction to the direct estimate of -1.4+/-0.7 Pg-C y-1. Taking the pre-industrial steady-state ocean source of 0.4+/-0.2 Pg-C y-1 into account, the total ocean uptake flux including the anthropogenic CO2 is estimated to be -2.0+/- 1.0 Pg-C y-1 in 2000.