We infer CO₂ surface fluxes using satellite observations of mid-tropospheric CO₂ from the Tropospheric Emission Spectrometer (TES) and measurements of CO₂ from surface flasks in a time-independent inversion analysis based on the GEOS-Chem model. Using TES CO₂ observations over oceans, spanning 40° S–40° N, we find that the horizontal and vertical coverage of the TES and flask data are complementary. This complementarity is demonstrated by combining the datasets in a joint inversion, which provides better constraints than from either dataset alone, when a posteriori CO₂ distributions are evaluated against independent ship and aircraft CO₂ data. In particular, the joint inversion offers improved constraints in the tropics where surface measurements are sparse, such as the tropical forests of South America. Aggregating the annual surface-to-atmosphere fluxes from the joint inversion for the year 2006 yields −1.13±0.21 Pg C for the global ocean, −2.77±0.20 Pg C for the global land biosphere and −3.90±0.29 Pg C for the total global natural flux (defined as the sum of all biospheric, oceanic, and biomass burning contributions but excluding CO₂ emissions from fossil fuel combustion). These global ocean and global land fluxes are shown to be near the median of the broad range of values from other inversion results for 2006. To achieve these results, a bias in TES CO₂ in the Southern Hemisphere was assessed and corrected using aircraft flask data, and we demonstrate that our results have low sensitivity to variations in the bias correction approach. Overall, this analysis suggests that future carbon data assimilation systems can benefit by integrating in situ and satellite observations of CO₂ and that the vertical information provided by satellite observations of mid-tropospheric CO₂ combined with measurements of surface CO₂, provides an important additional constraint for flux inversions.