Mapping Global Atmospheric Co 2 Concentration at High Spatiotemporal Resolution

Satellite measurements of the spatiotemporal distributions of atmospheric CO 2 concentrations are a key component for better understanding global carbon cycle characteristics. Currently, several satellite instruments such as the Greenhouse gases Observing SATellite (GOSAT), SCanning Imaging Absorption Spectrometer for Atmospheric CHartographY (SCIAMACHY), and Orbiting Carbon Observatory-2 can be used to measure CO 2 column-averaged dry air mole fractions. However, because of cloud effects, a single satellite can only provide limited CO 2 data, resulting in significant uncertainty in the characterization of the spatiotemporal distribution of atmospheric CO 2 concentrations. In this study, a new physical data fusion technique is proposed to combine the GOSAT and SCIAMACHY measurements. On the basis of the fused dataset, a gap-filling method developed by modeling the spatial correlation structures of CO 2 concentrations is presented with the goal of generating global land CO 2 distribution maps with high spatiotemporal resolution. The results show that, compared with the single satellite dataset (i.e., GOSAT or SCIAMACHY), the global spatial coverage of the fused dataset is significantly increased (reaching up to approximately 20%), and the temporal 871 resolution is improved by two or three times. The spatial coverage and monthly variations of the generated global CO 2 distributions are also investigated. Comparisons with ground-based Total Carbon Column Observing Network (TCCON) measurements reveal that CO 2 distributions based on the gap-filling method show good agreement with TCCON records despite some biases. These results demonstrate that the fused dataset as well as the gap-filling method are rather effective to generate global CO 2 distribution with high accuracies and high spatiotemporal resolution.


Introduction
Atmospheric carbon dioxide (CO 2 ) is the most important anthropogenic greenhouse gas, and since the industrial revolution, the CO 2 concentration in the Earth's atmosphere has increased significantly from 280 to 379 ppm in 2005 [1].Currently, global warming caused by atmospheric CO 2 has attracted the attention of scientists around the world.Predicting and mitigating climate change due to increased CO 2 depends on the accurate quantification of distribution and variability for CO 2 sources and sinks, which have been derived from atmospheric CO 2 concentration measurements by using inverse modeling [2][3][4].For this purpose, globally distributed measurements of atmospheric CO 2 concentrations with high accuracy and precision as well as high measurement density are required.In fact, the utility of CO 2 concentrations has been demonstrated in most carbon cycle-related studies.For example, Rayner and O'Brien [5] demonstrated that global column-averaged CO 2 concentrations (precision ≤1%) can help to reduce the uncertainties in regional CO 2 source and sink estimates.Alkhaled et al. [3] indicated that using atmospheric CO 2 concentration data has the potential to improve the scientific understanding of regional carbon cycle processes and budgets.
Existing ground-based CO 2 monitoring networks provide accurate measurements of atmospheric CO 2 concentration.However, these CO 2 measurements are very sparse to understand and capture the global distribution of carbon sources and sinks [3,6].Effectively grasping the global CO 2 distribution with high spatiotemporal resolution has been a long-standing problem.Satellites provide the potential to derive CO 2 column-averaged dry air mole fractions (XCO 2 ) accurately with high spatiotemporal resolutions on global scale.To date, the Orbiting Carbon Observatory-2 (OCO-2) [7], the Greenhouse gases Observing SATellite (GOSAT) [8], and the SCanning Imaging Absorption Spectrometer for Atmospheric CHartographY (SCIAMACHY) [9] are the only three satellite instruments that can retrieve XCO 2 with significant sensitivity in the boundary layer [10].Although these satellites have been used to measure global XCO 2 concentrations, because of cloud contamination and limitations due to the instrument observation modes (e.g., GOSAT observes in lattice points), the available data points for retrieved XCO 2 from any single satellite are very limited [3,11].
For instance, Morino et al. [12] showed that only approximately 10% of the GOSAT data can be used for the retrieval of XCO 2 because of cloud contamination.The limited satellite observations restrict the generation of level 3 (L3) XCO 2 maps with high spatiotemporal resolution when only a single satellite-based XCO 2 dataset is considered.Considering the weakness of single satellites, Wang et al. [13] proposed a physical fusing algorithm to generate a continuous spatiotemporal CO 2 dataset.However, in certain regions, even after combining GOSAT with SCIAMACHY retrievals at time scales of up to a few months, some gaps remain, thereby influencing the analysis of global and regional spatiotemporal characteristics of XCO 2 concentrations.Fortunately, researchers have recently suggested gap-filling methods based on the Kriging method to generate full-coverage global maps (L3 data products) using single satellite products.For example, Hammerling et al. [6,14] adopted a statistical mapping approach to create full-coverage maps (L3 data products) from GOSAT XCO 2 observations.The National Institute for Environmental Studies (NIES) GOSAT Project generates GOSAT L3 data products by interpolating and extrapolating the GOSAT level 2 (L2) XCO 2 products using ordinary Kriging [15].Liu et al. [16] and Tomasada et al. [17] also exploited methods to produce the spatial distribution of CO 2 concentrations.Creating L3 data products is one way to obtain new carbon cycle information [6].Hammerling et al. [6] indicated that these L3 XCO 2 products are not intended for direct use in inversion studies but can be used to make direct and independent comparisons with existing carbon flux and atmospheric transport models.Such comparison studies could be used to verify and track reported CO 2 emissions [6].
Nevertheless, the fact is that these single satellite-based strategies could not ensure sufficient data for stable semivariogram estimation during Kriging interpolation because the number of valid data points is limited for a single satellite [13,18], which may cause large uncertainties.Rather than using L2 XCO 2 from a single dataset (e.g., GOSAT), as performed in existing literature, fusing available CO 2 measurements derived from various space-based data would facilitate the generation of highly reliable full coverage (L3) maps with high spatiotemporal resolution.Wang et al. [13] proved that the spatial coverage of the fused data is wider than that of a single satellite.Thus, using fused data to generate a global land full-coverage map of XCO 2 distribution based on the Kriging approach is a better choice.
The overall objective of this study is to develop a global land gap-filling method and generate a CO 2 L3 map with high spatiotemporal resolution by using the new combined CO 2 dataset within 1° × 1° grids.The remainder of this paper is as follows.Section 2 describes the data used in this study.Section 3 introduces the fused algorithm and gap-filling methods based on ordinary Kriging.The results of the spatial variability analysis are presented in Section 4. In this section, comparison between single dataset-based prediction and that of fused-dataset-based prediction as well as the interpolated predictions versus ground-based measurements are investigated.Our conclusions are presented in Section 5.
SCIAMACHY is a spectrometer that measures reflected, scattered, and transmitted solar radiation in the spectral region of 214-2380 nm at moderate spectral resolution [26,27].It was onboard the Environmental Satellite (ENVISAT) [28]; however, contact with ENVISAT was lost in April 2012 [29].The primary objective of SCIAMACHY was to monitor trace gases in the troposphere and stratosphere.
NASA's OCO-2 launched in July 2014 will soon be providing approximately 100,000 high-quality daily measurements of CO 2 concentrations from around the globe [7].Currently, no XCO 2 data is publicly available from OCO-2.
In this study, version 2.9 of the GOSAT Atmospheric CO2 Observations from Space (ACOS) XCO 2 L2 data and version v02.00.08 of the SCIAMACHY Bremen optimal estimation (BESD) L2 data [30] are employed.ACOS v2.9 rather than ACOS v3.3 is used because ACOS v3.3 is still being evaluated and some deficiencies exist in this version [31].BESD v02.00.08 is the newest algorithm and latest version from IUP Bremen.The fused data used in our study are obtained by combining the XCO 2 retrievals of ACOS with those of BESD.The detailed fusing algorithm is described in the following section.
In addition, CO 2 data of Carbon Tracker (CT), which is a data assimilation system built by the National Oceanic and Atmospheric Administration's Earth System Research Laboratory, USA [32], are taken as a reference profile in the fusing algorithm.The Total Carbon Column Observing Network (TCCON, http://tccon.ipac.caltech.edu/)XCO 2 data, which are ground-based Fourier transform spectrometer measurements [33], are used for comparison with the interpolated results.In this study, Lamont (USA, 36.60°N,97.49°W), Park Falls (USA, 45.94°N, 90.27°W), and Garmisch (Germany, 47.48°N, 11.06°E) sites located in the Northern Hemisphere, and Wollongong (Australia, 34.41°S, 150.88°E) site located in the Southern Hemisphere, are used.

Fusing CO 2 Measurements from GOSAT and SCIAMACHY
Considering the current availability of multiple XCO 2 datasets from space, combining these measurements is a promising way to obtain more data points within a certain time compared with any single satellite-based dataset.Both GOSAT and SCIAMACHY can provide XCO 2 information, and we can consider combining these two XCO 2 datasets to generate a new dataset at high spatiotemporal resolution.However, these two datasets cannot be combined directly because of differences in sampling methods, overpass times, averaging kernels, and the a priori CO 2 profile involved during retrievals [34].Therefore, a necessary adjustment is performed to adjust the two observations for a common a priori profile (i.e., the newly generated fused profile).Wang et al. [13] employed the CT CO 2 profiles (2° × 3°) to assist data fusion.They indicated that a priori CO 2 profile of both the ACOS and BESD are first interpolated or extrapolated to the level of the CT CO 2 profile according to their pressure layers.After interpolation, the a priori CO 2 profiles for both ACOS and BESD have the same dimension as the CT profile.Then equation ( 1) can be used to adjust these two XCO 2 data to a unified level.In this study, a new CO 2 profile database with 1° × 1° resolution is used by integrating the CO 2 profiles of CT, ACOS, BESD, and NIES v02.xx products.Thus, the adjustment equation can be expressed as follows [13]: In the above equation, XCO 2_adj refers to the adjusted XCO 2 for ACOS or BESD, XCO 2_ret is the retrieved XCO 2 of ACOS or BESD, h denotes a pressure weighting function, a is the column averaging kernel of ACOS or BESD, I is an identity matrix, and X CT and X a are the newly generated fused CO 2 profile and the corresponding a priori CO 2 profile for ACOS or BESD, respectively; Here, the unit of XCO 2_adj , XCO 2_ret , and CO 2 profile is ppm.
In view of different spatial samplings, to account for retrieval uncertainties, the XCO 2 data with fine spatial scale are aggregated to a new value according to the uncertainty weight of individual XCO 2 retrievals to suit the relatively coarse scale during fusion [11].The detailed "aggregated" method can be found in [13].In addition, the time difference between these two datasets is also considered by interpolating the newly generated fused profile (1° × 1°) data at a reference time (taking BESD overpass time as a reference).
On the basis of the steps described above, the globe is divided into numerous 1° × 1° latitude/longitude grid boxes (180 × 360 in total).For each grid, all CO 2 observations located within that grid are averaged as the fused value weighted by the individual XCO 2 retrieval accuracy.

Gap-Filling Method for the Fused Data
In this study, a gap-filling method based on the ordinary Kriging [6,15] technique is applied to map global land full-coverage XCO 2 distribution by using fused ACOS and BESD data.This gap-filling method employs the spatial correlation of CO 2 observations between different locations.The spatial correlation structure of the CO 2 observations is derived using a semivariogram.Then, the derived spatial correlation structure and CO 2 observations are used to estimate XCO 2 values.The XCO 2 values are predicted within global 1° × 1° grids.The detailed process is described in the following.
First, abnormal XCO 2 data points should be screened as they may have significant impact on the interpolation results when no other nearby data point exists [15].In this study, skewness and kurtosis is calculated to filter these extreme values that make the L2 data distribution differ from normal distribution.Our threshold number for skewness is 0.01 and that for kurtosis is 5.By using these two measures, the screening process is performed, and the detailed screening process is similar to that adopted by the NIES GOSAT Project [15].
Second, calculating and analyzing experimental semivariograms is indispensable to quantify spatial variability between CO 2 observations; this can be calculated by the following equation [16]: whereγ(h) denotes the experimental semivariogram, Z(X i ) is the XCO 2 value at X i , X i and X i + h are the spatial locations on the earth's surface, N(h) is the number of sample pairs separated by h, and h is the spatial distance between sample pairs of X i and X i + h, which is calculated by [15,17]: 1 cos (sin sin cos cos cos( )) where r is the Earth's radius and and are the latitude and longitude of location X i respectively.Next, to quantify the spatial correlation structure between the pairs of XCO 2 measurements, it is crucial to choose a theoretical semivariogram model.In this study, depending on the feature of the spatial variability of CO 2 observations, an exponential semivariogram model with a nugget effect component was selected to model the XCO 2 spatial variability [18].
where N is the nugget effect component, C is the sill value, and R is the range value [18].These three model parameters, which fit the experimental semivariogram, are estimated by a nonlinear least-square method.
The subsequent prediction step is an ordinary Kriging approach.A distinct feature of Kriging is that an observation is not only weighted as a function of its distance to the prediction location but also as a function of its location relative to those of other observations [6].Following the algorithm of GOSAT L3 products [15], the spatial correlation structure for each location is estimated by using local Kriging within a certain search range according to the property of a semivariogram (i.e., the farther the distance, the lesser the observation points contribute).Here, the search range of observation points adopts a moving ellipse similar to the method used by GOSAT L3 products [15]; however, the actual ellipse radii used were half that of the GOSAT L3 products.The detailed ellipse radii are listed in Table 1.
The XCO 2 at each grid point x 0 is subsequently predicted by the following equation [17]: Where Z * (X 0 ) refers to predicted XCO 2 , n is the number of observation points, and Z(X i ) is the value of the variable Z at X i . in ( 5) is calculated by the following equation [17]: where = ( − ) ⋯ ( − ) ] and = { ( − )} .
Another significant characteristic of Kriging is that it can quantify the uncertainties in the predicted value [15].The mean square prediction error can be calculated by the following equation [15]: In the interpolation processes, because of the differences between experimental semivariograms, the distribution tendency of global CO 2 concentrations varies significantly between land and sea [15].Therefore, in our study, the experiments are restricted to global land areas.

Global Spatial Distribution of XCO 2 for Fused ACOS and BESD Data
In this study, as an example, the monthly averaged XCO 2 data for October 2010 are used to investigate the global spatial coverage of the fused ACOS and BESD data.The fused data performed at other months are similar to that of October 2010.For better visualization, all data are aggregated in 2° × 2° bins by averaging data points for every grid cell in October 2010.
As shown in Figure 1, the XCO 2 measurements from ACOS are distributed over both land and ocean areas.The XCO 2 data points from BESD are restricted to land regions because of the low signal-to-noise ratio over the ocean.Although the data points from ACOS are slightly wider than those from BESD, both ACOS and BESD show poor global coverage.These results agree with those presented in the study conducted by Wang et al. [11].
From Figure 1, we can see that, on a global scale, the spatial distribution of the fused data is more extensive compared to that of other single-satellite data.We also investigated the global land coverage percentage for fused, ACOS, and BESD data (Figure 2) within 1° × 1° grids.Figure 2 shows that the global land spatial coverage of the fused data can reach up to 20.04% for time periods up to 30 days.For 30 days, the average global coverage of ACOS and BESD within 1° × 1° grids is approximately 8.86% and 14.60%, respectively.The global land spatial coverage of the fused data reached up to 15.57% within 15 days and 13.03% within 10 days.These results illustrate the overwhelming advantages of using fused data for both space and time scales.Although the fused data generate a wider global coverage, there are still some gaps in some regions.In this case, a proper gap-filling method is required to generate a continuous full-coverage map.To address this, an ordinary Kriging interpolation approach based on the three datasets (i.e., fused, ACOS, and BESD) was attempted and is described in the next section.

Estimating and Modeling Experimental Semivariograms of ACOS, BESD, and Fused XCO 2 Data
In practice, initially a semivariogram model should be designated to fit a semivariogram obtained from actual observed data [15].Semivariogram models include Gaussian, spherical, exponential, linear, and power models.After many tests, it turned out that the exponential semivariograms can fit the characteristics of the actually observed data more accurately than other semivariogram models for all three CO 2 datasets (i.e., fused, ACOS, and BESD).Therefore, here, the exponential semivariogram model with a nugget effect component is selected.The experimental semivariograms and fitted results of fused, ACOS, and BESD data for October 2010 are shown in Figure 3.
In this study, on the basis of the relative distance, all semivariograms derived from CO 2 observations are classified on a 100 km scale.The experimental semivariograms, which are denoted as blue dots in Figure 3, represent the average value of the semivariograms for each classification.
As shown by Figure 3, there are significant spatial correlations within these CO 2 datasets, although the spatial correlation becomes weaker as the distance increases.Evidently, a potential advantage of the fused data is that it can ensure more data points for stable semivariogram estimation.

Comparison of the Interpolated Map of the Fused CO 2 with that of Single Satellite CO 2
In this study, the ordinary Kriging method is used to fill in the gap of global land CO 2 distributions for the fused and two single satellite datasets.A feature of this gap-filling method is that each predicted value has an associated uncertainty, which reflects the number of observations surrounding an estimation location and the spatial variability in the XCO 2 field [14].The predicted values with their associated uncertainties are used to implement comparisons between the fused data and the two single satellite observations.Figure 4 shows an example of the comparison of predicted results for these three datasets in October 2010.
Global standard deviation for interpolated XCO 2 products (ACOS, BESD, and fused) in April, July, and October 2010 was also investigated and compared (Table 2).It can be seen from Table 2 that the interpolated fused data show the smallest mean standard deviation and the minimum value.
As shown in Figure 4, although the spatial coverage of the interpolated ACOS shows a similar trend with that of the interpolated fused data, it shows a small gap in the south of Central America relative to the fused data.From the corresponding uncertainties of ACOS, we also see that the gaps exist in the north of Southern America.Therefore, the spatial coverage of the interpolated fused data is somewhat wider over the global land region compared to the interpolated CO 2 value of single satellite products.
In particular, interpolated values of the fused data show lower uncertainties than those of the two satellites.These results indicate that our fused CO 2 products are more effective for generating global full-coverage CO 2 concentrations based on the ordinary Kriging method.In addition, with regard to the predictions, the spatial coverage over land in the Southern Hemisphere, e.g., southern Africa and southern South America, is relatively good for each of the three datasets.The predicted uncertainties are also low for these regions.However, for the interpolated results of the fused, ACOS, and BESD data, regions with relatively weak constraint exist in high northern latitudes, e.g., northern Canada, Greenland, and northern Russia.This is thought to be due to limited data points.At the same time, it can be seen from Figure 4 that, in eastern Africa and the northeast part of South America, high uncertainties also exist in the fused, ACOS, and BESD product, which may also be due to limited data points.Furthermore, we can see from Figure 4 that the interpolated fused results are similar to those of ACOS in the northern part of North America and the Eurasian area.This is because the fused data points are primarily from ACOS observation points, and limited BESD data points are available in these regions.
Overall, the fused CO 2 data based on the ordinary Kriging is more helpful to map the global land CO 2 distribution at high spatiotemporal resolution, and in the following section, we further investigate monthly variability for the interpolated fused data.

Monthly Variability for Predicting the Fused Data
In this study, monthly global land distributions of the predicted CO 2 with the associated uncertainties are presented for April, July, September, and December 2010.As shown in Figure 5, the number of predicted CO 2 values decreases slightly in April and December for the Northern Hemisphere.In addition, poor coverage is found in the high northern latitudes for April and December, which is probably due to solar zenith angle restrictions and limited valid observations in these regions [14].
Furthermore, the monthly variations are captured well in the prediction maps of the fused data, especially in the Northern Hemisphere with relatively higher CO 2 values in April and lower CO 2 values in September, thus reflecting the effect of the atmospheric CO 2 seasonal cycle.The cause of these significant seasonal variations in CO 2 concentration in the Northern Hemisphere has been analyzed by many researchers.They have indicated that, in the Northern Hemisphere, XCO 2 values are higher in spring because of coal combustion as well as strong respiration of plants and soil, and XCO 2 values are lower in autumn because of strong photosynthesis activity.For example, Liu et al. [16] proposed that smoke and dust from coal combustion is a main reason for higher XCO 2 concentration in winter and spring.Bai et al. [35] stated that significant seasonal variations of CO 2 in the Northern Hemisphere were closely related to human activities and green vegetation.Wang et al. [11] also found that seasonal cycles exist in the retrieved XCO 2 concentrations from satellites in the Northern Hemisphere and that there are no significant seasonal variations in the Southern Hemisphere.From Figure 4, we can also see that there is no noticeable seasonal variation in the Southern Hemisphere, particularly for Australia and the southern part of South America.

Comparison of Kriging Interpolation Results with Total Carbon Column Observing Network (TCCON) Measurements
To investigate the effectiveness of the interpolation result based on the fused dataset, ground-based measurements are compared with the predictions for monthly mean values in 2010.The ground measurements taken between 10 a.m. and 2 p.m. are used to calculate monthly averaged values to match the overpass time of the satellites involving GOSAT and SCIAMACHY.Monthly means typically comprise n ≈ 100-1000 individual measurements for TCCON and typically n ≈ 25 data for fused Kriging.
From Figure 6 it can be seen that the monthly average predicted CO 2 data in 2010 are in good agreement with those of the TCCON sites on the whole, especially in the Lamont site located in the Northern Hemisphere.However, the monthly average predicted CO 2 data are lower than those of the TCCON sites by 0.5% in the Lamont site, which may stem from underestimated observation values; these biases were also found by Wunch et al. [36].At the same time, there are also some differences between the interpolated fused results and site data in some months for the Wollongong site.This is probably due to the lower number of data points from this site [36] for these months, which could result in inaccurate interpolated values.Note that the seasonal variations of the predicted and ground-based CO 2 data in Lamont, with higher values in spring and winter and lower values in autumn, are also sufficiently evident.However, this tendency is not found at the Wollongong site located in the Southern Hemisphere.

Conclusions
In this study, a physical fusion and gap-filling method was proposed to map global land CO 2 distribution by combining GOSAT with SCIAMACHY data at high spatiotemporal resolution.Initially, spatial coverage of the monthly averaged CO 2 in October 2010 was investigated for the fused, GOSAT, and SCIAMACHY data products.At the same time, the global land coverage percentage was compared among the fused, ACOS, and BESD datasets.The comparison results show that the global land spatial coverage of the fused data could reach up to 20.04% within 30 days, while the average global coverage of ACOS and BESD was approximately 8.86% and 14.60%, respectively.However, the global land spatial coverage percentage of the fused data within 15 days and 10 days reached up to 15.57% and 13.03%, respectively.Compared to ACOS or BESD, the fused data showed two or three times higher temporal resolution.These results indicate that the fused dataset is very effective for mapping the global distribution of XCO 2 concentration in either space or time scales.Based on this, the global land spatial correlation structure was evaluated and modeled using an exponential semivariogram model with a nugget effect component for fused, ACOS, and BESD CO 2 databases for October 2010.Subsequently, gap-filling maps of these CO 2 datasets with associated uncertainties were constructed on the basis of ordinary Kriging interpolation.The interpolated results imply that the interpolated fused data have the lowest standard deviation (mean value: 1.61 ppm and minimum value: 1.29 ppm) among these three databases and the largest spatial coverage.The monthly variations of the predictions for the fused data were also investigated.The results show that the monthly variations are relatively noticeable with higher values in April and December and lower values in July and September in the Northern Hemisphere, which is coincident with the CO 2 seasonal cycle.
In addition, the interpolated monthly averaged fused data were compared with data from the TCCON Lamont, Park Falls and Garmisch sites in the Northern Hemisphere and the Wollongong site in the Southern Hemisphere.The result revealed that the interpolated fused data is in good agreement with the TCCON sites on the whole.Furthermore, the seasonal cycle of the monthly averaged predicted fused data are also consistent with the TCCON measurements.These findings prove that predicted XCO 2 based on the fused dataset and gap-filling method is very useful to map global full-coverage XCO 2 distribution.
OCO-2 was launched in June 2014, and its data is still not available publicly.Once its data is released publicly, we will incorporate it and then generate improved fused XCO 2 datasets.In the current fusion of the proposed method, we only use the CO 2 profile of CT and Satellite (including GOSAT and SCIAMACHY); a new ground-based CO 2 profile and other Satellite CO 2 profiles, such as OCO-2, will be further considered to improve our method in the future.Considering the generality of our strategies shown in this study, they can be easily adapted to process other trace gases.

Figure 4 .
Figure 4. Comparison of the spatial distribution of the Kriging predictions and associated standard deviations: (A) Interpolated fused CO 2 data; (B) Standard deviation in (A); (C) Interpolated ACOS CO 2 data; (D) Standard deviation in (B); (E) Interpolated BESD CO 2 data; and (F) Standard deviation in (E).

Figure 5 .
Figure 5. Monthly mean maps of the predicted fused data with associated standard deviation (gridded in 2° × 2° bins): (A) Predicted fused data for April 2010; (B) Associated standard deviation for the predicted value (April 2010); (C) Predicted fused data for July 2010; (D) Associated standard deviation for the predicted value (July 2010); (E) Predicted fused data for September 2010; (F) Associated standard deviation for the predicted value (September 2010); (G) Predicted fused data for December 2010; and (H) Associated standard deviation for the predicted value (December 2010).

Figure 6 .
Figure 6.Error bars are pure statistical uncertainties of the monthly means (±3 sigma/sqrt(n) for TCCON site and fused Kriging data within ±2.5° in 2010; for TCCON sites, the monthly data between 10 am and 2 pm are averaged according to the overpass time of two satellites being combined.(A) Lamont site in 2010; (B) Park Falls site in 2010; (C) Garmisch site in 2010; and (D) Wollongong site in 2010.

Table 1 .
Search range used in the interpolated processes (unit: km).

Table 2 .
Basic global standard deviation statistics for interpolated XCO 2 products in April, July, and October 2010.Unit: ppm.