Articles | Volume 11, issue 3
Data description paper
04 Jul 2019
Data description paper |  | 04 Jul 2019

Satellite and ground-based measurements of XCO2 in a remote semiarid region of Australia

Voltaire A. Velazco, Nicholas M. Deutscher, Isamu Morino, Osamu Uchino, Beata Bukosa, Masataka Ajiro, Akihide Kamei, Nicholas B. Jones, Clare Paton-Walsh, and David W. T. Griffith

In this study, we present ground-based measurements of column-averaged dry-air mole fractions (DMFs) of CO2 (or XCO2) taken in a semiarid region of Australia with an EM27/SUN portable spectrometer equipped with an automated clamshell cover. We compared these measurements to space-based XCO2 retrievals from the Greenhouse Gases Observing Satellite (GOSAT). Side-by-side measurements of EM27/SUN with the Total Carbon Column Observing Network (TCCON) instrument at the University of Wollongong were conducted in 2015–2016 to derive an XCO2 scaling factor of 0.9954 relative to TCCON. Although we found a slight drift of 0.13 % over 3 months in the calibration curve of the EM27/SUN vs. TCCON XCO2, the alignment of the EM27/SUN proved stable enough for a 2-week campaign, keeping the retrieved Xair values, another measure of stability, to within 0.5 % and the modulation efficiency to within 2 %. From the measurements in Alice Springs, we confirm a small bias of around 2 ppm in the GOSAT M-gain to H-gain XCO2 retrievals, as reported by the NIES GOSAT validation team. Based on the reported random errors from GOSAT, we estimate the required duration of a future campaign in order to better understand the estimated bias between the EM27/SUN and GOSAT. The dataset from the Alice Springs measurements is accessible at (Velazco et al.2018).

1 Introduction

The Greenhouse Gases Observing Satellite (GOSAT), launched on 23 January 2009, is the first satellite in orbit dedicated to the measurements of the two major anthropogenic greenhouse gases that cause global warming: CO2 and CH4 (Yokota et al.2009; Kuze et al.2009). Since its launch, GOSAT data have been used to study and determine the global distributions of CO2 and CH4 (Reuter et al.2014; Houweling et al.2015; Parker et al.2015; Turner et al.2015). The spectrometer on board GOSAT, called Thermal and Near infrared Sensor for Carbon Observation Fourier-Transform Spectrometer (TANSO-FTS), is able to simultaneously detect shortwave infrared (SWIR: bands 1 at 0.76 µm, 2 at 1.6 µm and 3 at 2.0 µm) and thermal infrared (TIR: band 4 from 5.5 to 14.3 µm). The TANSO-FTS has a footprint of about 10.5 km in diameter on the ground viewed at nadir, which is observed for 4 s. Bands 1, 2 and 3 are measured in two linear polarizations simultaneously with three selectable gains; high (H), medium (M) and low (L) (Suto et al.2013). GOSAT employs the different gain settings to compensate for the different signal levels due to the reflective properties (albedo) of the Earth's surface, which also depend on wavelength. These gain settings are prespecified at certain locations because GOSAT does not observe both M and H gains simultaneously. For the majority of the soundings over land, H gain is used. GOSAT M-gain retrievals over land are used over surfaces that are bright in the SWIR such as deserts and semiarid regions. However, a bias between GOSAT M- and H-gain retrievals of XCO2 has been reported, and there has been a lack of M-gain validation with TCCON (Yoshida et al.2013). In Australia where surface reflectivity values are generally high, GOSAT was configured to observe using M gain for much of the land surface, at first to avoid detector saturation. However, it was found that some observations using H gain did not result in saturated signals and were still useful. Therefore from 11 February 2012, GOSAT started to perform alternate observations using both H and M gains in order to investigate the differences from retrievals between these two gain settings.

Figure 1 shows the locations of GOSAT M-gain soundings over a period of 1 year obtained from the NIES (National Institute for Environmental Studies, Japan) version 2.72 XCO2 retrieval algorithm (white squares) and data from the previous version 2.60 (black dots) for reference. Satellite-based retrievals of XCO2 are validated by the Total Carbon Column Observing Network (TCCON), a ground-based network dedicated to the precise and accurate measurements of greenhouse gases (Wunch et al.2011). TCCON data are used to validate satellite instruments such as GOSAT (Morino et al.2011), GOSAT-2 (Matsunaga et al.2018), OCO-2 (Crisp et al.2017) and Sentinel-5P (Borsdorff et al.2018). However, apart from locations in the western US (Dryden and Los Alamos TCCON sites), there are no operational TCCON stations (magenta stars in Fig. 1) that are ideally located in a clean desert or semiarid region not influenced by nearby anthropogenic emissions and have plenty of nearby GOSAT M-gain soundings. The Dryden TCCON station at the Armstrong Flight Research Center (AFRC), Edwards, CA, is located in the Mojave desert at 34.960 N, 117.881 W (700 m a.s.l.). Although not as densely populated, AFRC is only approximately 100 km north of Los Angeles (pop. ca. 17.8 million) and 100 km east of Bakersfield, CA (pop. ca. 376 380).

Figure 1Locations of GOSAT M-gain soundings for 2016. Black dots are from version 2.60 and the white squares are from version 2.72. Magenta stars indicate the location of operational TCCON sites as of October 2018. The magenta cross marks the location of Alice Springs in Australia and the black and yellow dashed circle has a radius of approximately 1000 km centered on Alice Springs. The Darwin and Wollongong (W'gong) sites in Australia are labeled.


The world's deserts and semiarid regions encompass large areas that are mostly undisturbed by recent anthropogenic emissions and are important for understanding the carbon cycle. Recently, the importance of semiarid regions in the carbon cycle inter-annual variability has been highlighted (Poulter et al.2014; Ahlström et al.2015; Trudinger et al.2016). If the goal is to continuously improve the accuracy of carbon cycle studies, then the accuracy of satellite retrievals needs to be improved as well because if observations over high albedo are not available or biased, flux estimates would likely be biased as well and would lead to misinterpretation. Therefore, measurements over M-gain regions are needed by the satellite community (Yoshida et al.2013) and highlight the significance of this study.

Central Australia has a large semiarid region that is relatively easy to access compared to other vast semiarid sites in the world. There are plenty of uninhabited, vast, homogeneous areas that can accommodate large footprints of satellite-based sensors such as GOSAT, GOSAT-2, OCO-2, and Sentinel-5P. Central Australia is also an ideal place to measure undisturbed atmospheric conditions that can serve as a calibration point for satellite retrievals of atmospheric composition. At the same time, the desert environment has high surface reflectivity, which is a challenge for satellite retrievals because aerosols, depending on type, can lead to an effect called optical path lengthening and this effect is dominant in regions with high albedo (Yoshida et al.2013). Recently, Iwasaki et al. (2019) showed that there is an increase in the XCO2 retrievals using their PPDF-S algorithm when the albedo at 1.6 µm was high, implying that the retrieved XCO2 is strongly related to the surface albedo. This challenge leads to the improvement of satellite retrievals. Despite the importance of desert locations like Central Australia in remote sensing, infrastructure support is not available. Benchmark measurements and pilot studies for desert sites are needed to assess the benefit and feasibility of such sites because setting up a TCCON site in remote deserts will be difficult logistically and financially. We address this need by utilizing a well-established portable solar-viewing spectrometer, an EM27/SUN by Bruker Optics GmbH (Gisi et al.2012; Frey et al.2015, 2019; Hase et al.2015; Hedelius et al.2016), which was retrofitted with a protective fairing and automated solar tracker clamshell cover for operations in a harsh environment. For brevity, we will refer to this instrument as EM27. The instrument can measure spectra covering the spectral bands in the near infrared range necessary to derive column-averaged dry-air mole fractions of CO2, CH4 and CO with sufficient stability for short-term campaigns (Hedelius et al.2016). The instrument was transported to Alice Springs in Central Australia with the primary objective of validating the GOSAT XCO2 signal and making benchmark measurements in the region. The Alice Springs measurements are unique considering that they have been collected from a clean desert environment where GOSAT M-gain soundings are also abundant and close enough (within 100 km) to compare with the EM27.

This paper is organized as follows. A description of Alice Springs and desert Australia is given in Sect. 2. In Sect. 3, we briefly discuss the instruments and methods, which are already well established. Results of measurement comparisons with the TCCON station in Wollongong and comparisons with the GOSAT M-gain and H-gain soundings in Alice Springs are shown and discussed under Sect. 4. We focus on XCO2 measurements from the EM27 and GOSAT. We compare NIES GOSAT retrieval version 2.60 and the new version 2.72, which are both not bias corrected, with the retrievals from the EM27. Statistical calculations and recommendations for a future Alice Springs campaign are also under Sect. 4. We provide our conclusions in Sect. 5.

2 Alice Springs, Australia, site description

We conducted a measurement campaign with the EM27 system at the Australian Bureau of Meteorology (BOM) facility in Alice Springs (23.79 S, 133.89 E) from 29 September to 6 October 2016, with the primary objective of validating the GOSAT XCO2 signal in this semiarid region. Alice Springs is located in Central Australia, also called the “Red Centre”. Next to Darwin and Palmerston, it is the third-largest town in the Northern Territory of Australia with a population of 23 726 (2016 census). The vegetation around Alice Springs is composed mostly of dry scrubby grassland. The Alice Springs terrain, consisting mostly of sandy plains with some areas of rocky highland, is bounded by several deserts: the Tanami desert to the north, Simpson desert to the east and southeast, the Great Victoria desert to the south and the Gibson desert to the west. In September, statistical data from the Australian Bureau of Meteorology (BOM) show an average monthly rainfall of 8.7 mm, average daily high and low temperatures of 27.3 and 10.3 C, respectively, and mean monthly sunshine of 300 h. Collected GOSAT M-gain soundings over land for one whole year shown in Fig. 1 include a large part of this Red Centre including Alice Springs. Note that for version 02.72 FTS SWIR L2 retrievals, upgraded input and reference products were used. For example, there was an improvement in the spatial resolution of the cloud flagging procedure, which employs the CAI L2 (Cloud and Aerosol Imager, level 2) data. This improvement resulted in better screening of the data and may have resulted in an increased number of soundings that were passed for the version 02.72 FTS SWIR L2 retrievals. Apart from remote regions in North Africa, the Middle East and near-densely populated areas in California, US, Central Australia is the only region that provides an abundant number of M-gain soundings.

3 Instruments and methods

3.1 EM27/SUN with automated clamshell cover

The EM27 system and its characterization are thoroughly described in the works of Gisi et al. (2012), Frey et al. (2015), Hase et al. (2015) and Hedelius et al. (2016). For total column measurements of CO2, CH4, H2O and O2 spectra in the near infrared range, the instrument is fitted with an indium gallium arsenide (InGaAs) detector dedicated to 5000–12 000 cm−1. This also enables measurements of spectra covering the O2 bands necessary to derive column-averaged dry-air mole fractions of CO2 and CH4 similar to the method used by TCCON. TCCON uses a maximum optical path difference (MOPD) of 45 cm, corresponding to a spectral resolution of 0.02 cm−1. However, in contrast to TCCON, the EM27 we used has the typical MOPD of 1.8 cm, corresponding to a spectral resolution of 0.5 cm−1.

We equipped the Bruker EM27 with a weather station-controlled automated clamshell cover for the solar tracker and protective fairing (Fig. 2). Our objective was to achieve autonomous and remote measurements of greenhouse gases in harsh environments, in particular, to perform much needed measurements of XCO2 and XCH4, in desert Australia. The design and more details of this construction can be obtained by contacting the authors.

Figure 2Drawing of the Wollongong EM27 automated clamshell cover and fairing (by Steve Selby, UOW).

Recently, Hedelius et al. (2016) published a long-term assessment of errors and biases in retrievals of several TCCON gases from the EM27. In spite of a reported drift, they found that the stability of the EM27 is sufficient for short-term campaigns. Therefore, before deployment to Alice Springs, we have operated the EM27 from December 2015 to September 2016 and tested the automated clamshell cover at the University of Wollongong, where the TCCON instrument (Griffith et al.2014) is also located. Results of the collocated measurements are in Sect. 4.

3.2 Retrieval of XCO2 and Xair

To retrieve XCO2, we use the same retrieval software used by TCCON called GGG2014 (Wunch et al.2017), which was also used by Hedelius et al. (2016) for their EM27. Column-averaged dry-air mole fractions (DMFs) of gases (Xgas) are retrieved from the EM27 measurement as in Wunch et al. (2010):

(1) X gas = VC gas VC dryair = 0.2095 VC gas VC O 2 ,

where VCgas is the vertical column of the gas. VCdryair is the dry pressure column of air and 0.2095 is the known DMF of oxygen. The DMF of air (“Xair”) can be calculated using the measured oxygen column from the EM27 spectrum:

(2) X air = VC air VC O 2 0.2095 - XH 2 O m H 2 O m air dry ,

where mH2O and mairdry are the mean molecular masses of water (18.02 g mol−1) and dry air (28.964 g  mol−1), XH2O is the retrieved DMF of water vapor and VCair is calculated from the surface pressure, Ps:

(3) VC air = P s g m air dry N a ,

where g is the column-averaged acceleration due to gravity and Na is Avogadro’s constant. Xair is a good indicator of instrument stability and changes in spectrometer alignment because (1) VCair is calculated using the surface pressure, which is independently measured by a pressure sensor to better than 0.3 hPa, keeping accuracy over long periods (Wunch et al.2011), and (2) the atmospheric oxygen column is not particularly variable in dry air; hence the retrieved VCO2 by the spectrometer should be close to constant. From Eq. (2) it follows that a perfectly accurate measurement would lead to an Xair value of unity; however due to inaccuracies in the O2 spectroscopy, the actual value is approximately 0.98 for all TCCON sites. We also find that Xair values well outside of the second and 98th percentiles may indicate an obstruction in the solar beam (e.g., birds, fast moving clouds, leaves from trees). Therefore the Xair can also serve as data filtering criteria.

In Wollongong, both the TCCON and EM27 solar tracker covers were opened and closed by the same pneumatic system. The same weather station provided the meteorological data that were used to pre-filter the data (e.g., fractional variation in solar intensity, wind speed and direction, pressure). For this study, we did not filter the data according to solar zenith angles (SZAs) anymore; instead, in addition to the pre-filter we filtered out noisy retrievals by selecting only those with Xair values within 0.5 and 1.5 because anything beyond that would be unrealistic in the atmosphere but most likely be the cause of an interference or obstruction.

3.3 GOSAT specific point observations

The TANSO-FTS on GOSAT has a two-axis pointing system. Normally, TANSO-FTS follows an M-shaped grid on a five-point cross-track scan mode (, last access: 12 June 2019). By using this pointing system to vary the observation geometry, it is able to observe specific points, i.e., it can view targets with angles up to ±20 along the satellite track and by ±30 across the track. Specific point observations over Alice Springs were requested from July 2016, in preparation for the campaign in September. Five locations within 100 km of the center of Alice Springs were targeted (see Fig. 3). We only used the specific point observation data to compare with the EM27. However, to construct the time series shown in Fig. 8, we used all available GOSAT data version V2.72 from NIES spanning the years 2010 to 2017. Daily averages within 1000 km of Alice Springs are calculated for M- and H-gain retrievals separately.

Figure 3A Google Earth map (Data SIO, NOAA, U. S. Navy, NGA, GEBCO © 2018 Google) showing the GOSAT specific point observations around Alice Springs (cyan circles labeled “Obs#”) . The ground-based site is securely fenced and located at the airport, ca. 14 km south of the town center. The magenta circle shows a 60 km radius from the site, for visual reference. The inset shows the EM27 during the measurements.

4 Results

4.1 Comparison with Wollongong TCCON station

In this section, we present a side-by-side comparison of retrievals of XCO2 between the EM27 and the Wollongong TCCON station. Apart from the works by Hedelius et al. (2016) and Frey et al. (2019), studies on long-term comparisons of the EM27 with TCCON are rare. Here, we focus on comparisons of Xair and XCO2 from measurements spanning almost 1 year (November 2015 to September 2016) under varied environmental conditions. Both instruments normally measure at the same time, apart from interruptions due to occasional mid-infrared measurements with the 125HR or rare software glitches (e.g., JAVA issues).


As mentioned in Sect. 3.2, inaccuracies in the O2 spectroscopy lead to an Xair value of approximately 0.98 for all TCCON sites. This value varies by approximately 1 % with solar zenith angle, indicating an air mass dependence in the O2 retrievals (Pollard et al.2017). Deviations from the characteristic values for Xair are generally indicative of erroneous behavior in the measurement and retrieval system such as interferometer misalignment, tracking errors or fitting to an incorrect air mass due to timing errors, as recently reported in Pollard et al. (2017).

Figure 4 shows the normalized probability distribution functions (PDFs) of Xair values retrieved by both instruments in 2016. One manifestation of instrumental differences between the Wollongong TCCON and the EM27 is the difference in the Xair retrieved from both instruments. The corresponding mean values of Xair from both instruments are 0.9824 and 0.9849 for TCCON and EM27, respectively. The TCCON Xair PDF exhibits a more Gaussian pattern, while the EM27 Xair PDF is slightly skewed.

Figure 5 shows the mean Xair values calculated for each month in 2016 for TCCON (red) and EM27 (blue). The error bars represent the spread (1 standard deviation) in the monthly mean values of Xair. This gives us a measure of how much the Xair changes and therefore another measure of instrument stability. In Fig. 5, we see that the mean Xair values from the EM27 are slightly larger compared to the TCCON instrument but a slight seasonal dependence can be seen in both instruments. The difference in the mean Xair values does not appear to be purely related to sampling differently across that seasonal dependence. The SZAs seem to have an effect as well, although small. To show this in Fig. 5, we filtered Xair values for SZA ≤45 (black dots and circles) and above 45 (gray pluses and crosses). The SZA effect is within approximately 1 % for both instruments. We note that there were no realignments performed on the EM27 and TCCON instruments during this period; however the EM27 clamshell cover and fairing were fitted on the EM27 in February 2016. It is possible that this may have affected the alignment resulting in a shift in the Xair.

Figure 4Normalized probability distribution functions (PDFs) of Xair from TCCON (b) and Xair from the EM27 (a) in Wollongong for 2016. The PDFs for each month are separated by colors and the positions of the mean Xair monthly values are indicated by the month number. The mean Xair for the whole year for the EM27 is 0.9849 and 0.9824 for TCCON. The clamshell cover and fairing were installed on the EM27 in February; this may have affected the alignment resulting in a shift in the Xair.


Figure 5A measure of the spread of the of Xair values from TCCON (red) and EM27 (blue) by taking the mean (central markers) and corresponding standard deviations (bars) for each of the monthly probability distribution functions shown in Fig. 4. Filtered Xair values for angles ≤45 (black dots and circles) and above 45 (gray pluses and crosses) are also shown for reference.


4.1.2XCO2 comparison with Wollongong TCCON and apparent drift

From the EM27 and TCCON measurements in 2016, we derived a scaling factor for the EM27. We scaled the EM27 measurements to Wollongong TCCON data because the Wollongong XCO2 has been calibrated against aircraft profile measurements that are traceable to the WMO in situ scale (Wunch et al.2010). The retrieval method is predicted to be both linear and have a zero intercept (Wunch et al.2010). Therefore, we fit hourly mean data from TCCON and EM27 employing linear least-squares and force a zero intercept. The standard errors of the weighted means are used as weights in the fit. From this exercise, we arrive at a scaling factor of EM27 =0.9954 TCCON (Fig. 6). Frey et al. (2015) reported a comparable scaling factor of 0.9951⋅ TCCON for XCO2 using collocated measurements of five EM27 instruments with the TCCON instrument at the Karlsruhe Institute of Technology in Germany. However, Hedelius et al. (2016) reported a smaller bias of +0.03 % between EM27 and TCCON. We think that the slightly larger bias in this work (−0.46 %) is probably due to an imperfect alignment of our EM27, with a modulation efficiency (ME) of 96 % at maximum optical path difference (OPD). The ME calculation was performed using water lines measured at ambient laboratory/room air (Frey et al.2015). The instrument line shape (ILS) was monitored before, during and after the campaign and the ME remained stable within 94 %–96 %. Therefore, we are confident that the scaling factor relative to TCCON was consistent during the campaign.

Hedelius et al. (2016) reported a noticeable drift in the EM27 measurements over several months. We also observed a very small drift in the EM27 vs. TCCON XCO2 values shown in Fig. 6, where the calibration line has been derived. For clarity, we color-coded the EM27 vs. TCCON XCO2 values in Fig. 6 according to days after 1 January 2016. The colors progress with time, i.e., from dark blue (oldest) to dark brown (most recent). We found that the drift within 3 months before the campaign is very small, i.e., 0.13 % over 3 months as calculated from the EM27/TCCON scaling factors from July to mid-September. We did not consider the averaging kernels in this work. The averaging kernels of the EM27 have been previously presented and compared to TCCON in a study by Hedelius et al. (2016). In their study, they found that although there are differences in the TCCON and EM27 averaging kernels, the effect of the differences in averaging kernels from the top of the atmosphere cancel out the effect of differences at the bottom. Further work and more measurements may be necessary to better understand the cause of this phenomenon but this is beyond the scope of this study.

Figure 6Weighted hourly mean XCO2 from TCCON vs. XCO2 from EM27 in Wollongong for 2015–2016. To avoid noisy data, only measurements corresponding to TCCON with retrieved Xair values within 0.9783 and 0.9853 were selected (Xair values within the second to 98th percentiles). The inset shows the drift in the hourly weighted mean scaling factors over 3 months (red dots). The error bars are the sum, in quadrature, of the relative standard deviations of the weighted hourly means from each instrument.


4.2 Measurements in Alice Springs

Results of the ground-based XCO2 measurements during the Alice Springs campaign are shown in Fig. 7. The measurements started on 29 September 2016 with interruptions due to cloudy and/or rainy weather. In total, we gathered 6 d of good ground-based measurements of XCO2 under good weather conditions, which coincided with nine GOSAT specific point soundings within 100 km of the site. A slight air mass dependence in the XCO2 retrievals results in lower XCO2 values at low sun elevation, as can be observed in Fig. 7. This air mass dependence is well-known and also discussed in Wunch et al. (2011). But for the purpose of GOSAT comparisons, this air mass dependence is negligible because GOSAT passes the site at around 13:05 local time (03:35 UTC), corresponding to a high sun elevation (see Table 1 for the dates).

Figure 7Retrievals of XCO2 from the EM27 in Alice Springs, compressed along the x axis (days). The measurement days on the x axis are not continuous due to interruptions from bad weather. See Sect. 6 for data availability.


Table 1GOSAT specific point observation opportunities over Alice Springs.

Download Print Version | Download XLSX

4.3 Comparisons with GOSAT measurements in Alice Springs

To demonstrate the variability of XCO2 in the region, we plot in Fig. 8 the daily mean time series of GOSAT soundings using these two coincidence criteria: (1) all soundings within a 1000 km radius centered at the BOM facility in Alice Springs and (2) all soundings measured within the same day (local time). Each data point considered in the calculation of the daily mean is weighted by the corresponding reported retrieval error for that particular GOSAT sounding. Vertical gray lines represent the corresponding standard deviation for each daily mean calculation. The H-gain retrievals (red triangles) and M-gain retrievals (blue squares) are separated in this plot. Weighted linear least-squares fitted straight lines were calculated for the M-gain XCO2 (cyan) and H-gain XCO2 data (magenta). We confirm that the M-gain retrievals are biased high (around 2 ppm) compared to the H-gain retrievals. Nevertheless, both least-squares fitted lines to the H- and M-gain retrievals show an increase of about 2.28 ppm yr−1 with y intercept values at 133.0 ppm (H gain) and 134.92 ppm (M gain).

Compared to the TCCON sites in Wollongong and Darwin (Deutscher et al.2010, 2014), the XCO2 signal in Alice Springs is relatively smooth and undisturbed, which is expected because Alice Springs is in the middle of the Australian continent, a desert environment with no large sources to interfere with the XCO2 signal. The terrestrial biosphere is the largest driver of variability in the Southern Hemisphere column XCO2 (Deutscher et al.2014). Moreover, Deutscher et al. (2014) have shown that the magnitude of the seasonal variability in the column-average dry-air mole fraction of XCO2 is comparable in magnitude to the annual increase. In Alice Springs, only a slight seasonal cycle in the measured XCO2 can be seen. This slight cycle is probably less driven by the surrounding vegetation because the center of the Australian arid zone is not affected by seasonality and has low aggregate vegetation (Lawley et al.2011) but is more affected by meridional effect.

Figure 8Time series XCO2 from GOSAT M-gain and H-gain retrievals in Alice Springs from the NIES version 2.72 product. The timing of the Alice Springs campaign from 27 September to 6 October is indicated by the yellow arrow. Gray lines represent 1 standard error associated with each daily mean. The inset shows a map projection with the approximate locations of the soundings used for this plot. Only data within a 1000 km radius from Alice Springs (black X on the map) were considered.


Rainfall records indicate that April–August 2011 had the least amount of rainfall on record in Alice Springs after 2002, which was only broken by 2015, and then by 2017 (, last access: 12 June 2019, station number 015590). This significant absence of rain in the region could result in stunted vegetation growth or mortality, which could have led to brighter surfaces. Bright surfaces mean more M-gain measurements and this may explain the 74 % increase in M-gain soundings in April–August 2011 compared to the April–August 2010–2017 average of ca. 220 soundings per month. In contrast, the number of H-gain retrievals seems to have diminished around 2011, coinciding with the dry months starting from April 2011. The number of H-gain soundings from April–August 2011 was 23 % fewer compared to the April–August average from the years 2010–2017 (ca. 450 soundings per month).

GOSAT specific point observations that were scheduled during the campaign are shown on Table 1. The times correspond to the times when the satellite is directly above the site. The satellite normally performs five observation points across track, with an interferometric scan time lasting 4 s (Shiomi et al.2006). But specific point observations deviate from this pattern by pointing and maximizing observations near the target. M- and H-gain retrievals (6 d apart) were obtained during the short campaign. We averaged and compared the GOSAT retrievals directly with the coincident measurements from the EM27 (±0.5 h) and as in Sect. 4.1.2, we fit the data using weighted linear least squares and force a zero intercept. From here, we derived a scaling factor of 0.9927⋅ EM27 for GOSAT H gain and 0.9983⋅ EM27 for the GOSAT M gain, as shown in Fig. 9. This is only slightly different from the previous version 2.60 of the data (H gain: 0.9935⋅ EM27 and M gain: 0.9997⋅ EM27). Indeed, there is a slight bias in the M- and H-gain retrievals from GOSAT; i.e., M-gain retrievals are biased slightly low and H-gain retrievals are very close compared to the EM27. However, the M- and H-gain biases relative to each other have been improved in the version 2.72 data release.

Figure 9Mean XCO2 from the EM27 vs. GOSAT XCO2 version 2.60 (a) and version 2.72 (b) in Alice Springs. The calculated means are weighted by the 1-σ uncertainty of the individual measurements. The scaling factor derived from the TCCON comparisons, EM27 =0.9954 TCCON, has been taken into account.


4.4 Recommendations for a future Alice Springs campaign: simulating error reduction in bias estimates

In this subsection, we attempt to estimate how long a measurement campaign in Alice Springs should last in order to improve the statistical errors in the estimated bias between GOSAT and the EM27. In order to estimate the required measurements from the EM27, we make the following assumptions and conditions.

  1. There is at least one specific point overpass per week (M gain or H gain), similar to this campaign.

  2. The standard errors of the weighted hourly mean GOSAT specific point observations are 0.25 and 0.26 ppm for M gains (versions 2.6 and 2.72, respectively) and 0.29 and 0.28 ppm for H gains (versions 2.6 and 2.72, respectively). These values are taken from the averaged standard errors of the weighted hourly means from all specific point observation data falling within 100 km of the site for 2 September 2015–February 2017.

  3. 1-σ standard deviation of all EM27 measurements is 0.05 ppm at an averaging time of 1 h. Chen et al. (2016) reported a 1-σ standard deviation of 0.04 ppm at 10 min averaging.

  4. The error in the estimated bias is calculated as σT=σGOSAT2+σEM272, which then improves by the number of weeks (N), i.e., σT2/N.

  5. Any long-term drifts in the satellite and ground-based measurements are negligible over the whole duration of the campaign.

It is important to keep the statistical errors in the estimated bias as low as possible. Miller (2007) showed that a comparison of surface CO2 concentration data and XCO2 data flux inversions clearly reveals a land–ocean bias in the XCO2 retrievals, even when the bias is only 0.1 ppm. Figure 10 shows a plot of how the statistical errors in the estimated bias improve with the number of weeks in the campaign. The inset shows the normalized histogram of the reported GOSAT single sounding errors within 100 km of the site. Note that although H-gain retrievals have smaller errors on average, the number of M-gain retrievals improves the standard error of the weighted mean. A 2-week campaign would already reduce the statistical error between GOSAT M-gain and EM27 measurements to less than the GOSAT M-gain weighted mean standard error (0.26 ppm, v2.72). To achieve a statistical error of 0.1 ppm between GOSAT M-gain and EM27 measurements, it would take about 6 to 7 weeks of measurements or more. Conversely, around 8.5 weeks for H gain is needed because the standard error of the hourly weighted mean H-gain retrievals are slightly higher due to fewer soundings near the site. The ideal time window to perform this would be from March to November (autumn–winter–spring). It is best to avoid the summer months (December–February) because of the desert heat, which could be physically demanding for the operator(s) and may affect the operation of an EM27 without a temperature-regulated enclosure. A temperature-regulated enclosure would improve stability and make the EM27 suitable for summer operations even in the high heat (>40C) of Alice Springs. Note however that our error estimates are conservative and this implies that there are no drifts in the M- and H-gain measurements relative to one another that would require long-term validation.

Figure 10The statistical errors in the estimated bias of XCO2 improve with the number of weeks in the campaign. Inset: normalized histogram of the reported GOSAT single sounding errors within 100 km of the site.


5 Data availability

Column averaged dry-air mole fractions of CO2 in Alice Springs, Australia, are hosted by the University of Wollongong and can be downloaded via, last access: 12 June 2019, with the DOI: (Velazco et al.2018). GOSAT data can be obtained from the GOSAT Data Archive Service (GDAS) after registration for access via, last access: 12 June 2019.

6 Conclusions

Alice Springs would be a welcome addition to ground sites dedicated to environmental satellite sensor calibration and validation. With the exception of Lauder, New Zealand, there are no permanently dedicated ground targets for satellite calibration and validation in the Southern Hemisphere that are far enough from water (to avoid a possible land–ocean bias in retrievals) and surrounded by homogeneous topography (to avoid a possible altitude bias). Currently, vicarious radiometric calibration of passive optical instruments from satellites is carried out on the large homogeneous Railroad Valley playa (RRV) in Nevada, US. However plans for future lithium mining threaten to disrupt the playa floor, which could render this critical asset useless for large footprint sensors on satellites like GOSAT, GOSAT-2, OCO-2 and OCO-3. Alice Springs would offer unique opportunities in the context of satellite surface products and associated emissivity assessments because the environment of Alice Springs has ideal conditions to measure undisturbed atmospheric conditions that can serve as a calibration point for satellite retrievals of atmospheric composition. At the same time, the desert environment, like RRV, provides high surface reflectivity, which is a challenge for satellite retrievals and thus observations are especially needed here for validation purposes. Subject to funding and interest from the community, a repeat of the campaign is possible. In a future study, we will investigate GOSAT methane retrievals using the same approach employed here.

Author contributions

VAV conceptualized and wrote the paper, performed the calculations and statistics, made the figures, performed EM27 measurements in Alice Springs, led the EM27 automated clamshell cover project, characterized the EM27 and TCCON ILS and oversaw TCCON measurements in Wollongong during this study.

NMD helped in writing the paper, organized and co-acquired funding for the Alice Springs campaign, liaised with BOM, performed the GGG2014 retrievals of the EM27 data, performed measurements in Alice Springs and oversaw the EM27 measurements in Wollongong.

IM helped in writing the paper, provided ideas and conceptualized the need for statistical calculations (Sect. 4.4).

OU presented and highlighted the M- and H-gain bias and the need for validating GOSAT in the M-gain region as representative of the GOSAT project in the research announcement (RA) meeting, helped in writing the paper and provided advice.

DWTG supervised the EM27 automated clamshell cover project, supervised the EM27 and TCCON measurements, created the preprocessing software, helped in writing the paper and provided advice.

NBJ carried out EM27 measurements in Wollongong and helped in determining the EM27 ILS.

CPW performed TCCON measurements in Wollongong and helped improve the paper.

BB carried out TCCON measurements in Wollongong and helped improve the paper.

AK and MA enabled the specific observation mode in Alice Springs from 2016.

Competing interests

The authors declare that they have no conflict of interest.


Voltaire A. Velazco would like to thank University of Wollongong SMAH for the advancement grant scheme that funded the fabrication of the EM27 clamshell cover and fairing. We would like to thank the Australian Bureau of Meteorology (BOM), especially site manager Victoria McLean for the support in Alice Springs. We thank Graham Kettlewell for technical advice, programming the weather station and liaising with the workshop. We thank Martin Riggenbach for the design and creation of the weather station interface and for providing technical advice. We thank Steve Selby for the design and fabrication of the EM27 clamshell and cover and Phd student Neil Page for assistance with the EM27 measurements. The specific point observation modes have been made possible by the GOSAT project office. Part of this activity has been conducted under the framework of the GOSAT Research Announcement 8 (RA8) project. Nicholas M. Deutscher acknowledges funding from the Australian Research Council (ARC) via DECRA grant DE140100178.

Financial support

This research has been supported by the University of Wollongong (grant no. 2016/SPGA-S/07) and the Australian Research Council (grant nos. DP0879468, DP110103118, DP140101552, DP160101598, and LE0668470).

Review statement

This paper was edited by David Carlson and reviewed by two anonymous referees.


Ahlström, A., Raupach, M. R., Schurgers, G., Smith, B., Arneth, A., Jung, M., Reichstein, M., Canadell, J. G., Friedlingstein, P., Jain, A. K., Kato, E., Poulter, B., Sitch, S., Stocker, B. D., Viovy, N., Wang, Y. P., Wiltshire, A., Zaehle, S., and Zeng, N.: The dominant role of semi-arid ecosystems in the trend and variability of the land CO2 sink, Science, 348, 895–899,, 2015. a

Borsdorff, T., aan de Brugh, J., Hu, H., Hasekamp, O., Sussmann, R., Rettinger, M., Hase, F., Gross, J., Schneider, M., Garcia, O., Stremme, W., Grutter, M., Feist, D. G., Arnold, S. G., De Mazière, M., Kumar Sha, M., Pollard, D. F., Kiel, M., Roehl, C., Wennberg, P. O., Toon, G. C., and Landgraf, J.: Mapping carbon monoxide pollution from space down to city scales with daily global coverage, Atmos. Meas. Tech., 11, 5507–5518,, 2018. a

Chen, J., Viatte, C., Hedelius, J. K., Jones, T., Franklin, J. E., Parker, H., Gottlieb, E. W., Wennberg, P. O., Dubey, M. K., and Wofsy, S. C.: Differential column measurements using compact solar-tracking spectrometers, Atmos. Chem. Phys., 16, 8479–8498,, 2016. a

Crisp, D., Pollock, H. R., Rosenberg, R., Chapsky, L., Lee, R. A. M., Oyafuso, F. A., Frankenberg, C., O'Dell, C. W., Bruegge, C. J., Doran, G. B., Eldering, A., Fisher, B. M., Fu, D., Gunson, M. R., Mandrake, L., Osterman, G. B., Schwandner, F. M., Sun, K., Taylor, T. E., Wennberg, P. O., and Wunch, D.: The on-orbit performance of the Orbiting Carbon Observatory-2 (OCO-2) instrument and its radiometrically calibrated products, Atmos. Meas. Tech., 10, 59–81,, 2017. a

Deutscher, N. M., Griffith, D. W. T., Bryant, G. W., Wennberg, P. O., Toon, G. C., Washenfelder, R. A., Keppel-Aleks, G., Wunch, D., Yavin, Y., Allen, N. T., Blavier, J.-F., Jiménez, R., Daube, B. C., Bright, A. V., Matross, D. M., Wofsy, S. C., and Park, S.: Total column CO2 measurements at Darwin, Australia – site description and calibration against in situ aircraft profiles, Atmos. Meas. Tech., 3, 947–958,, 2010. a

Deutscher, N. M., Sherlock, V., Mikaloff Fletcher, S. E., Griffith, D. W. T., Notholt, J., Macatangay, R., Connor, B. J., Robinson, J., Shiona, H., Velazco, V. A., Wang, Y., Wennberg, P. O., and Wunch, D.: Drivers of column-average CO2 variability at Southern Hemispheric Total Carbon Column Observing Network sites, Atmos. Chem. Phys., 14, 9883–9901,, 2014. a, b, c

Frey, M., Hase, F., Blumenstock, T., Groß, J., Kiel, M., Mengistu Tsidu, G., Schäfer, K., Sha, M. K., and Orphal, J.: Calibration and instrumental line shape characterization of a set of portable FTIR spectrometers for detecting greenhouse gas emissions, Atmos. Meas. Tech., 8, 3047–3057,, 2015. a, b, c, d

Frey, M., Sha, M. K., Hase, F., Kiel, M., Blumenstock, T., Harig, R., Surawicz, G., Deutscher, N. M., Shiomi, K., Franklin, J. E., Bösch, H., Chen, J., Grutter, M., Ohyama, H., Sun, Y., Butz, A., Mengistu Tsidu, G., Ene, D., Wunch, D., Cao, Z., Garcia, O., Ramonet, M., Vogel, F., and Orphal, J.: Building the COllaborative Carbon Column Observing Network (COCCON): long-term stability and ensemble performance of the EM27  SUN Fourier transform spectrometer, Atmos. Meas. Tech., 12, 1513–1530,, 2019. a, b

Gisi, M., Hase, F., Dohe, S., Blumenstock, T., Simon, A., and Keens, A.: XCO2-measurements with a tabletop FTS using solar absorption spectroscopy, Atmos. Meas. Tech., 5, 2969–2980,, 2012. a, b

Griffith, D. W. T., Velazco, V. A., Deutscher, N. M., Murphy, C., Jones, N., Wilson, S., Macatangay, R., Kettlewell, G., Buchholz, R., and Riggenbach, M.: Total Column Carbon Observing Network (TCCON) Site: Wollongong, CaltechDATA,, 2014. a

Hase, F., Frey, M., Blumenstock, T., Groß, J., Kiel, M., Kohlhepp, R., Mengistu Tsidu, G., Schäfer, K., Sha, M. K., and Orphal, J.: Application of portable FTIR spectrometers for detecting greenhouse gas emissions of the major city Berlin, Atmos. Meas. Tech., 8, 3059–3068,, 2015. a, b

Hedelius, J. K., Viatte, C., Wunch, D., Roehl, C. M., Toon, G. C., Chen, J., Jones, T., Wofsy, S. C., Franklin, J. E., Parker, H., Dubey, M. K., and Wennberg, P. O.: Assessment of errors and biases in retrievals of XCO2, XCH4, XCO, and XN2O from a 0.5 cm–1 resolution solar-viewing spectrometer, Atmos. Meas. Tech., 9, 3527–3546,, 2016. a, b, c, d, e, f, g, h, i

Houweling, S., Baker, D., Basu, S., Boesch, H., Butz, A., Chevallier, F., Deng, F., Dlugokencky, E. J., Feng, L., Ganshin, A., Hasekamp, O., Jones, D., Maksyutov, S., Marshall, J., Oda, T., O’Dell, C. W., Oshchepkov, S., Palmer, P. I., Peylin, P., Poussi, Z., Reum, F., Takagi, H., Yoshida, Y., and Zhuravlev, R.: An intercomparison of inversemodels for estimating sources and sinks of CO2 using GOSAT measurements, J. Geophys. Res.-Atmos., 120, 5253–5266,, 2015. a

Iwasaki, C., Imasu, R., Bril, A., Oshchepkov, S., Yoshida, Y., Yokota, T., Zakharov, V., Gribanov, K., and Rokotyan, N.: Optimization of the Photon Path Length Probability Density Function-Simultaneous (PPDF-S) Method and Evaluation of CO2 Retrieval Performance Under Dense Aerosol Conditions, Sensors, 19, 1262,, 2019. a

Kuze, A., Suto, H., Nakajima, M., and Hamazaki, T.: Thermal and near infrared sensor for carbon observation Fourier-transform spectrometer on the Greenhouse Gases Observing Satellite for greenhouse gases monitoring, Appl. Optics, 48, 6716,, 2009. a

Lawley, E., Lewis, M., and Ostendorf, B.: Environmental zonation across the Australian arid region based on long-term vegetation dynamics, J. Arid Environ., 75, 576–585,, 2011. a

Matsunaga, T., Morino, I., Yoshida, Y., Saito, M., Noda, H., Ohyama, H., Kamei, A., Kawazoe, F., Maksyutov, S., Yokota, T., Kuze, A., Suto, H., Nakajima, M., and Imasu, R.: FTS and FTS-2 Onboard Japanese GOSAT Earth Observation Satellite Series, in: Light, Energy and the Environment 2018 (E2, FTS, HISE, SOLAR, SSL), OSA Technical Digest (Optical Society of America, 2018), paper FW2B.2,, 2018. a

Miller, C. E., Crisp, D., DeCola, P. L., Olsen, S. C., Randerson, J. T., Michalak, A. M., Alkhaled, A., Rayner, P., Jacob, D. J., Suntharalingam, P., Jones, D. B. A., Denning, A. S., Nicholls, M. E., Doney, S. C., Pawson, S., Boesch, H., Connor, B. J., Fung, I. Y., O'Brien, D., Salawitch, R. J., Sander, S. P., Sen, B., Tans, P., Toon, G. C., Wennberg, P. O., Wofsy, S. C., Yung, Y. L., and Law, R. M.: Precision requirements for space-based XCO2 data, J. Geophys. Res., 112, D10314,, 2007. a

Morino, I., Uchino, O., Inoue, M., Yoshida, Y., Yokota, T., Wennberg, P. O., Toon, G. C., Wunch, D., Roehl, C. M., Notholt, J., Warneke, T., Messerschmidt, J., Griffith, D. W. T., Deutscher, N. M., Sherlock, V., Connor, B., Robinson, J., Sussmann, R., and Rettinger, M.: Preliminary validation of column-averaged volume mixing ratios of carbon dioxide and methane retrieved from GOSAT short-wavelength infrared spectra, Atmos. Meas. Tech., 4, 1061–1076,, 2011. a

Parker, R. J., Boesch, H., Byckling, K., Webb, A. J., Palmer, P. I., Feng, L., Bergamaschi, P., Chevallier, F., Notholt, J., Deutscher, N., Warneke, T., Hase, F., Sussmann, R., Kawakami, S., Kivi, R., Griffith, D. W. T., and Velazco, V.: Assessing 5 years of GOSAT Proxy XCH4 data and associated uncertainties, Atmos. Meas. Tech., 8, 4785–4801,, 2015. a

Pollard, D. F., Sherlock, V., Robinson, J., Deutscher, N. M., Connor, B., and Shiona, H.: The Total Carbon Column Observing Network site description for Lauder, New Zealand, Earth Syst. Sci. Data, 9, 977–992,, 2017. a, b

Poulter, B., Frank, D., Ciais, P., Myneni, R. B., Andela, N., Bi, J., Broquet, G., Canadell, J. G., Chevallier, F., Liu, Y. Y., Running, S. W., Sitch, S., and van der Werf, G. R.: Contribution of semi-arid ecosystems to interannual variability of the global carbon cycle, Nature, 509, 600–603,, 2014. a

Reuter, M., Buchwitz, M., Hilker, M., Heymann, J., Schneising, O., Pillai, D., Bovensmann, H., Burrows, J. P., Bösch, H., Parker, R., Butz, A., Hasekamp, O., O'Dell, C. W., Yoshida, Y., Gerbig, C., Nehrkorn, T., Deutscher, N. M., Warneke, T., Notholt, J., Hase, F., Kivi, R., Sussmann, R., Machida, T., Matsueda, H., and Sawa, Y.: Satellite-inferred European carbon sink larger than expected, Atmos. Chem. Phys., 14, 13739–13753,, 2014. a

Shiomi, K., Yoshida, M., and Sekio, N.: Study for sun glint observation of GOSAT FTS using MODIS and AMSR-E data, in: International Geoscience and Remote Sensing Symposium (IGARSS), Denver, CO, USA, 31 July–4 August 2006, IEEE,, 2006. a

Suto, H., Yoshida, J., Desbiens, R., Kawashima, T., and Kuze, A.: Characterization and correction of spectral distortions induced by microvibrations onboard the GOSAT Fourier transform spectrometer, Appl. Optics, 52, 4969–4980,, 2013. a

Trudinger, C. M., Haverd, V., Briggs, P. R., and Canadell, J. G.: Interannual variability in Australia's terrestrial carbon cycle constrained by multiple observation types, Biogeosciences, 13, 6363–6383,, 2016.  a

Turner, A. J., Jacob, D. J., Wecht, K. J., Maasakkers, J. D., Lundgren, E., Andrews, A. E., Biraud, S. C., Boesch, H., Bowman, K. W., Deutscher, N. M., Dubey, M. K., Griffith, D. W. T., Hase, F., Kuze, A., Notholt, J., Ohyama, H., Parker, R., Payne, V. H., Sussmann, R., Sweeney, C., Velazco, V. A., Warneke, T., Wennberg, P. O., and Wunch, D.: Estimating global and North American methane emissions with high spatial resolution using GOSAT satellite data, Atmos. Chem. Phys., 15, 7049–7069,, 2015. a

Velazco, V. A., Deutscher, N. M., and Griffith, D. W. T.: Column averaged dry mole fractions of carbon dioxide in Alice Springs, Australia, University of Wollongong,, 2018. a, b

Wunch, D., Toon, G. C., Wennberg, P. O., Wofsy, S. C., Stephens, B. B., Fischer, M. L., Uchino, O., Abshire, J. B., Bernath, P., Biraud, S. C., Blavier, J.-F. L., Boone, C., Bowman, K. P., Browell, E. V., Campos, T., Connor, B. J., Daube, B. C., Deutscher, N. M., Diao, M., Elkins, J. W., Gerbig, C., Gottlieb, E., Griffith, D. W. T., Hurst, D. F., Jiménez, R., Keppel-Aleks, G., Kort, E. A., Macatangay, R., Machida, T., Matsueda, H., Moore, F., Morino, I., Park, S., Robinson, J., Roehl, C. M., Sawa, Y., Sherlock, V., Sweeney, C., Tanaka, T., and Zondlo, M. A.: Calibration of the Total Carbon Column Observing Network using aircraft profile data, Atmos. Meas. Tech., 3, 1351–1362,, 2010. a, b, c

Wunch, D., Toon, G. C., Blavier, J.-F. L., Washenfelder, R. A., Notholt, J., Connor, B. J., Griffith, D. W. T., Sherlock, V., and Wennberg, P. O.: The Total Carbon Column Observing Network, Philos. T. Roy. Soc. A, 369, 2087–2112,, 2011. a, b, c

Wunch, D., Toon, G. C., Sherlock, V., Deutscher, N. M., Liu, C., Feist, D. G., and Wennberg, P. O.: Documentation for the 2014 TCCON Data Release, CaltechDATA,, 2017. a

Yokota, T., Yoshida, Y., Eguchi, N., Ota, Y., Tanaka, T., Watanabe, H., and Maksyutov, S.: Global Concentrations of CO2 and CH4 Retrieved from GOSAT: First Preliminary Results, SOLA, 5, 160–163,, 2009. a

Yoshida, Y., Kikuchi, N., Morino, I., Uchino, O., Oshchepkov, S., Bril, A., Saeki, T., Schutgens, N., Toon, G. C., Wunch, D., Roehl, C. M., Wennberg, P. O., Griffith, D. W. T., Deutscher, N. M., Warneke, T., Notholt, J., Robinson, J., Sherlock, V., Connor, B., Rettinger, M., Sussmann, R., Ahonen, P., Heikkinen, P., Kyrö, E., Mendonca, J., Strong, K., Hase, F., Dohe, S., and Yokota, T.: Improvement of the retrieval algorithm for GOSAT SWIR XCO2 and XCH4 and their validation using TCCON data, Atmos. Meas. Tech., 6, 1533–1547,, 2013. a, b, c

Short summary
We present ground-based measurements of atmospheric carbon dioxide columns from a portable spectrometer taken in a semiarid region of Australia. We compared these measurements to space-based retrievals from the Greenhouse Gases Observing Satellite (GOSAT) and calibrated them against a Total Carbon Column Observing Network (TCCON) instrument to ascertain a retrieval bias. We also present the unique opportunities that Central Australia could offer in the context of satellite product validation.