Search

We present a 10-year (January 2007–December 2016) time series of continuous in situ measurements of methane (CH4), carbon monoxide (CO) and nitrous oxide (N2O) made by an in situ Fourier transform infrared trace gas and isotope analyser (FTIR) operated at Lauder, New Zealand (45.04 S, 169.68 E, 370 m a. m. s. l.). Being the longest continuous deployed operational FTIR system of this type, we are in an ideal position to perform a practical evaluation of the multi-year performance of the analyser. The operational methodology, measurement precision, reproducibility, accuracy and instrument reliability are reported. We find the FTIR has a measurement repeatability of the order of 0.37 ppb (1σ standard deviation) for CH4, 0.31 ppb for CO and 0.12 ppb for N2O. Regular target cylinder measurements provide a reproducibility estimate of 1.19 ppb for CH4, 0.74 ppb for CO and 0.27 ppb for N2O. FTIR measurements are compared to co-located ambient air flask samples acquired at Lauder since May 2009, which allows a long-term assessment of the FTIR data set across annual and seasonal composition changes. Comparing FTIR and co-located flask measurements show that the bias (FTIR minus flask) for CH4 of −1.02 ± 2.61 ppb and CO of −0.43 ± 1.60 ppb are within the Global Atmospheric Watch (GAW)-recommended compatibility goals of 2 ppb. The N2O FTIR flask bias of −0.01 ± 0.77 ppb is within the GAW-recommended compatibility goals of 0.1 ppb and should be viewed as a serendipitous result due to the large standard deviation along with known systematic differences in the measurement sets. Uncertainty budgets for each gas are also constructed based on instrument precision, reproducibility and accuracy. In the case of CH4, systematic uncertainty dominates, whilst for CO and N2O it is comparable to the random uncertainty component. The long-term instrument stability, precision estimates and flask comparison results indicate the FTIR CH4 and CO time series meet the GAW compatibility recommendations across multiple years of operation (and instrument changes) and are sufficient to capture annual trends and seasonal cycles observed at Lauder. The differences between FTIR and flask N2O measurements need to be reconciled. Trend analysis of the 10-year time series captures seasonal cycles and the secular upward trend of CH4 and N2O. The CH4 and CO time series have the required precision and accuracy at a high enough temporal resolution to be used in inversion models in a data-sparse region of the world.

We describe a method for removing systematic biases of column-averaged dry air mole fractions of CO2 (XCO2) and CH4 (XCH4) derived from short-wavelength infrared (SWIR) spectra of the Greenhouse gases Observing SATellite (GOSAT). We conduct correlation analyses between the GOSAT biases and simultaneously retrieved auxiliary parameters. We use these correlations to bias correct the GOSAT data, removing these spurious correlations. Data from the Total Carbon Column Observing Network (TCCON) were used as reference values for this regression analysis. To evaluate the effectiveness of this correction method, the uncorrected/corrected GOSAT data were compared to independent XCO2 and XCH4 data derived from aircraft measurements taken for the Comprehensive Observation Network for TRace gases by AIrLiner (CONTRAIL) project, the National Oceanic and Atmospheric Administration (NOAA), the US Department of Energy (DOE), the National Institute for Environmental Studies (NIES), the Japan Meteorological Agency (JMA), the HIAPER Pole-to-Pole observations (HIPPO) program, and the GOSAT validation aircraft observation campaign over Japan. These comparisons demonstrate that the empirically derived bias correction improves the agreement between GOSAT XCO2/XCH4 and the aircraft data. Finally, we present spatial distributions and temporal variations of the derived GOSAT biases.

An algorithm for retrieval of vertical profiles from ground-based spectra in the near IR is described and tested. Known as GFIT2, the algorithm is primarily intended for CO2, and is used exclusively for CO2 in this paper. Retrieval of CO2 vertical profiles from ground-based spectra is theoretically possible, would be very beneficial for carbon cycle studies and the validation of satellite measurements, and has been the focus of much research in recent years. GFIT2 is tested by application both to synthetic spectra and to measurements at two Total Carbon Column Observing Network (TCCON) sites. We demonstrate that there are approximately 3° of freedom for the CO2 profile, and the algorithm performs as expected on synthetic spectra. We show that the accuracy of retrievals of CO2 from measurements in the 1.61μ (6220 cm−1) spectral band is limited by small uncertainties in calculation of the atmospheric spectrum. We investigate several techniques to minimize the effect of these uncertainties in calculation of the spectrum. These techniques are somewhat effective but to date have not been demonstrated to produce CO2 profile retrievals with sufficient precision for applications to carbon dynamics. We finish by discussing ongoing research which may allow CO2 profile retrievals with sufficient accuracy to significantly improve the scientific value of the measurements from that achieved with column retrievals.

High-resolution solar absorption spectra, taken within the Network for the Detection of Atmospheric Composition Change Infrared Working Group (NDACC-IRWG) in the mid-infrared spectral region, are used to infer partial or total column abundances of many gases. In this paper we present the retrieval of a column-averaged mole fraction of carbon dioxide from NDACC-IRWG spectra taken with a Fourier transform infrared (FTIR) spectrometer at the site in Ny-Ålesund, Spitsbergen. The retrieved time series is compared to colocated standard TCCON (Total Carbon Column Observing Network) measurements of column-averaged dry-air mole fractions of CO2 (denoted by xCO2). Comparing the NDACC and TCCON retrievals, we find that the sensitivity of the NDACC retrieval is lower in the troposphere (by a factor of 2) and higher in the stratosphere, compared to TCCON. Thus, the NDACC retrieval is less sensitive to tropospheric changes (e.g., the seasonal cycle) in the column average.

Consistent validation of satellite CO2 estimates is a prerequisite for using multiple satellite CO2 measurements for joint flux inversion, and for establishing an accurate long-term atmospheric CO2 data record. Harmonizing satellite CO2 measurements is particularly important since the differences in instruments, observing geometries, sampling strategies, etc. imbue different measurement characteristics in the various satellite CO2 data products. We focus on validating model and satellite observation attributes that impact flux estimates and CO2 assimilation, including accurate error estimates, correlated and random errors, overall biases, biases by season and latitude, the impact of coincidence criteria, validation of seasonal cycle phase and amplitude, yearly growth, and daily variability. We evaluate dry-air mole fraction (XCO2) for Greenhouse gases Observing SATellite (GOSAT) (Atmospheric CO2 Observations from Space, ACOS b3.5) and SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) (Bremen Optimal Estimation DOAS, BESD v2.00.08) as well as the CarbonTracker (CT2013b) simulated CO2 mole fraction fields and the Monitoring Atmospheric Composition and Climate (MACC) CO2 inversion system (v13.1) and compare these to Total Carbon Column Observing Network (TCCON) observations (GGG2012/2014). We find standard deviations of 0.9, 0.9, 1.7, and 2.1 ppm vs. TCCON for CT2013b, MACC, GOSAT, and SCIAMACHY, respectively, with the single observation errors 1.9 and 0.9 times the predicted errors for GOSAT and SCIAMACHY, respectively. We quantify how satellite error drops with data averaging by interpreting according to error2 = a2 + b2/n (with n being the number of observations averaged, a the systematic (correlated) errors, and b the random (uncorrelated) errors). a and b are estimated by satellites, coincidence criteria, and hemisphere. Biases at individual stations have year-to-year variability of ∼ 0.3 ppm, with biases larger than the TCCON-predicted bias uncertainty of 0.4 ppm at many stations. We find that GOSAT and CT2013b underpredict the seasonal cycle amplitude in the Northern Hemisphere (NH) between 46 and 53° N, MACC overpredicts between 26 and 37° N, and CT2013b underpredicts the seasonal cycle amplitude in the Southern Hemisphere (SH). The seasonal cycle phase indicates whether a data set or model lags another data set in time. We find that the GOSAT measurements improve the seasonal cycle phase substantially over the prior while SCIAMACHY measurements improve the phase significantly for just two of seven sites. The models reproduce the measured seasonal cycle phase well except for at Lauder_125HR (CT2013b) and Darwin (MACC). We compare the variability within 1 day between TCCON and models in JJA; there is correlation between 0.2 and 0.8 in the NH, with models showing 10–50 % the variability of TCCON at different stations and CT2013b showing more variability than MACC. This paper highlights findings that provide inputs to estimate flux errors in model assimilations, and places where models and satellites need further investigation, e.g., the SH for models and 45–67° N for GOSAT and CT2013b.

Consistent and accurate long-term data sets of global atmospheric concentrations of carbon dioxide (CO2) are required for carbon cycle and climate-related research. However, global data sets based on satellite observations may suffer from inconsistencies originating from the use of products derived from different satellites as needed to cover a long enough time period. One reason for inconsistencies can be the use of different retrieval algorithms. We address this potential issue by applying the same algorithm, the Bremen Optimal Estimation DOAS (BESD) algorithm, to different satellite instruments, SCIAMACHY on-board ENVISAT (March 2002–April 2012) and TANSO-FTS on-board GOSAT (launched in January 2009), to retrieve XCO2, the column-averaged dry-air mole fraction of CO2. BESD has been initially developed for SCIAMACHY XCO2 retrievals. Here, we present the first detailed assessment of the new GOSAT BESD XCO2 product. GOSAT BESD XCO2 is a product generated and delivered to the MACC project for assimilation into ECMWF's Integrated Forecasting System. We describe the modifications of the BESD algorithm needed in order to retrieve XCO2 from GOSAT and present detailed comparisons with ground-based observations of XCO2 from the Total Carbon Column Observing Network (TCCON). We discuss detailed comparison results between all three XCO2 data sets (SCIAMACHY, GOSAT and TCCON). The comparison results demonstrate the good consistency between SCIAMACHY and GOSAT XCO2. For example, we found a mean difference for daily averages of −0.60 ± 1.56 ppm (mean difference ± standard deviation) for GOSAT–SCIAMACHY (linear correlation coefficient r=0.82), −0.34 ± 1.37 ppm (r = 0.86) for GOSAT–TCCON and 0.10 ± 1.79 ppm (r = 0.75) for SCIAMACHY–TCCON. The remaining differences between GOSAT and SCIAMACHY are likely due to non-perfect collocation (± 2 h, 10° x 10° around TCCON sites), i.e. the observed air masses are not exactly identical but likely also due to a still non-perfect BESD retrieval algorithm, which will be continuously improved in the future. Our overarching goal is to generate a satellite-derived XCO2 data set appropriate for climate and carbon cycle research covering the longest possible time period. We therefore also plan to extend the existing SCIAMACHY and GOSAT data set discussed here by also using data from other missions (e.g. OCO-2, GOSAT-2, CarbonSat) in the future.

The Total Carbon Column Observing Network (TCCON) is a global ground-based network of Fourier transform spectrometers that produce precise measurements of column-averaged dry-air mole fractions of atmospheric methane (CH4). Temporal variability in the total column of CH4 due to stratospheric dynamics obscures fluctuations and trends driven by tropospheric transport and local surface fluxes that are critical for understanding CH4 sources and sinks. We reduce the contribution of stratospheric variability from the total column average by subtracting an estimate of the stratospheric CH4 derived from simultaneous measurements of hydrogen fluoride (HF). HF provides a proxy for stratospheric CH4 because it is strongly correlated to CH4 in the stratosphere, has an accurately known tropospheric abundance (of zero), and is measured at most TCCON stations. The stratospheric partial column of CH4 is calculated as a function of the zonal and annual trends in the relationship between CH4 and HF in the stratosphere, which we determine from ACE-FTS satellite data. We also explicitly take into account the CH4 column averaging kernel to estimate the contribution of stratospheric CH4 to the total column. The resulting tropospheric CH4 columns are consistent with in situ aircraft measurements and augment existing observations in the troposphere.

Column-averaged dry-air mole fractions of carbon dioxide and methane have been retrieved from spectra acquired by the TANSO-FTS (Thermal And Near-infrared Sensor for carbon Observations-Fourier Transform Spectrometer) and SCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric Cartography) instruments on board GOSAT (Greenhouse gases Observing SATellite) and ENVISAT (ENVIronmental SATellite), respectively, using a range of European retrieval algorithms. These retrievals have been compared with data from ground-based high-resolution Fourier transform spectrometers (FTSs) from the Total Carbon Column Observing Network (TCCON). The participating algorithms are the weighting function modified differential optical absorption spectroscopy (DOAS) algorithm (WFMD, University of Bremen), the Bremen optimal estimation DOAS algorithm (BESD, University of Bremen), the iterative maximum a posteriori DOAS (IMAP, Jet Propulsion Laboratory (JPL) and Netherlands Institute for Space Research algorithm (SRON)), the proxy and full-physics versions of SRON's RemoTeC algorithm (SRPR and SRFP, respectively) and the proxy and full-physics versions of the University of Leicester's adaptation of the OCO (Orbiting Carbon Observatory) algorithm (OCPR and OCFP, respectively). The goal of this algorithm inter-comparison was to identify strengths and weaknesses of the various so-called round- robin data sets generated with the various algorithms so as to determine which of the competing algorithms would proceed to the next round of the European Space Agency's (ESA) Greenhouse Gas Climate Change Initiative (GHG-CCI) project, which is the generation of the so-called Climate Research Data Package (CRDP), which is the first version of the Essential Climate Variable (ECV) "greenhouse gases" (GHGs). For XCO2, all algorithms reach the precision requirements for inverse modelling (< 8 ppm), with only WFMD having a lower precision (4.7 ppm) than the other algorithm products (2.4–2.5 ppm). When looking at the seasonal relative accuracy (SRA, variability of the bias in space and time), none of the algorithms have reached the demanding < 0.5 ppm threshold. For XCH4, the precision for both SCIAMACHY products (50.2 ppb for IMAP and 76.4 ppb for WFMD) fails to meet the < 34 ppb threshold for inverse modelling, but note that this work focusses on the period after the 2005 SCIAMACHY detector degradation. The GOSAT XCH4 precision ranges between 18.1 and 14.0 ppb. Looking at the SRA, all GOSAT algorithm products reach the < 10 ppm threshold (values ranging between 5.4 and 6.2 ppb). For SCIAMACHY, IMAP and WFMD have a SRA of 17.2 and 10.5 ppb, respectively.

The Total Carbon Column Observing Network (TCCON) has been established to provide ground-based remote sensing measurements of the column-averaged dry air mole fractions (DMF) of key greenhouse gases. To ensure network-wide consistency, biases between Fourier transform spectrometers at different sites have to be well controlled. Errors in interferogram sampling can introduce significant biases in retrievals. In this study we investigate a two-step scheme to correct these errors. In the first step the laser sampling error (LSE) is estimated by determining the sampling shift which minimises the magnitude of the signal intensity in selected, fully absorbed regions of the solar spectrum. The LSE is estimated for every day with measurements which meet certain selection criteria to derive the site-specific time series of the LSEs. In the second step, this sequence of LSEs is used to resample all the interferograms acquired at the site, and hence correct the sampling errors. Measurements acquired at the Izaña and Lauder TCCON sites are used to demonstrate the method. At both sites the sampling error histories show changes in LSE due to instrument interventions (e.g. realignment). Estimated LSEs are in good agreement with sampling errors inferred from the ratio of primary and ghost spectral signatures in optically bandpass-limited tungsten lamp spectra acquired at Lauder. The original time series of Xair and XCO2 (XY: column-averaged DMF of the target gas Y) at both sites show discrepancies of 0.2–0.5% due to changes in the LSE associated with instrument interventions or changes in the measurement sample rate. After resampling, discrepancies are reduced to 0.1% or less at Lauder and 0.2% at Izaña. In the latter case, coincident changes in interferometer alignment may also have contributed to the residual difference. In the future the proposed method will be used to correct historical spectra at all TCCON sites.

The column-averaged dry-air mole fractions of carbon dioxide and methane (XCO2 and XCH4) have been retrieved from Greenhouse gases Observing SATellite (GOSAT) Short-Wavelength InfraRed (SWIR) observations and released as a SWIR L2 product from the National Institute for Environmental Studies (NIES). XCO2 and XCH4 retrieved using the version 01.xx retrieval algorithm showed large negative biases and standard deviations (−8.85 and 4.75 ppm for XCO2 and −20.4 and 18.9 ppb for XCH4, respectively) compared with data of the Total Carbon Column Observing Network (TCCON). Multiple reasons for these error characteristics (e.g., solar irradiance database, handling of aerosol scattering) are identified and corrected in a revised version of the retrieval algorithm (version 02.xx). The improved retrieval algorithm shows much smaller biases and standard deviations (−1.48 and 2.09 ppm for XCO2 and −5.9 and 12.6 ppb for XCH4, respectively) than the version 01.xx. Also, the number of post-screened measurements is increased, especially at northern mid- and high-latitudinal areas.

Here, we report preliminary estimates of the column averaged carbon dioxide (CO2) dry air mole fraction, XCO2, retrieved from spectra recorded over land by the Greenhouse gases Observing Satellite, GOSAT (nicknamed "Ibuki"), using retrieval methods originally developed for the NASA Orbiting Carbon Observatory (OCO) mission. After screening for clouds and other known error sources, these retrievals reproduce much of the expected structure in the global XCO2 field, including its variation with latitude and season. However, low yields of retrieved XCO2 over persistently cloudy areas and ice covered surfaces at high latitudes limit the coverage of some geographic regions, even on seasonal time scales. Comparisons of early GOSAT XCO2 retrievals with XCO2 estimates from the Total Carbon Column Observing Network (TCCON) revealed a global, −2% (7–8 parts per million, ppm, with respect to dry air) XCO2 bias and 2 to 3 times more variance in the GOSAT retrievals. About half of the global XCO2 bias is associated with a systematic, 1% overestimate in the retrieved air mass, first identified as a global +10 hPa bias in the retrieved surface pressure. This error has been attributed to errors in the O2 A-band absorption cross sections. Much of the remaining bias and spurious variance in the GOSAT XCO2 retrievals has been traced to uncertainties in the instrument's calibration, oversimplified methods for generating O2 and CO2 absorption cross sections, and other subtle errors in the implementation of the retrieval algorithm. Many of these deficiencies have been addressed in the most recent version (Build 2.9) of the retrieval algorithm, which produces negligible bias in XCO2 on global scales as well as a ~30% reduction in variance. Comparisons with TCCON measurements indicate that regional scale biases remain, but these could be reduced by applying empirical corrections like those described by Wunch et al. (2011b). We recommend that such corrections be applied before these data are used in source sink inversion studies to minimize spurious fluxes associated with known biases. These and other lessons learned from the analysis of GOSAT data are expected to accelerate the delivery of high quality data products from the Orbiting Carbon Observatory-2 (OCO-2), once that satellite is successfully launched and inserted into orbit.

Column-averaged volume mixing ratios of carbon dioxide and methane retrieved from the Greenhouse gases Observing SATellite (GOSAT) Short-Wavelength InfraRed observation (GOSAT SWIR XCO2 and XCH4) were compared with the reference calibrated data obtained by ground-based high-resolution Fourier Transform Spectrometers (g-b FTSs) participating in the Total Carbon Column Observing Network (TCCON). Preliminary results are as follows: the GOSAT SWIR XCO2 and XCH4 (Version 01.xx) are biased low by 8.85 ± 4.75 ppm (2.3 ± 1.2 %) and 20.4 ± 18.9 ppb (1.2 ± 1.1 %), respectively. The standard deviation of the GOSAT SWIR XCO2 and XCH4 is about 1 % (1 σ) after correcting the negative biases of XCO2 and XCH4 by 8.85 ppm and 20.4 ppb, respectively. The latitudinal distributions of zonal means of the GOSAT SWIR XCO2 and XCH4 show similar features to those of the g-b FTS data except for the negative biases in the GOSAT data.

The Total Carbon Column Observing Network (TCCON) produces precise measurements of the column average dry-air mole fractions of CO2, CO, CH4, N2O and H2O at a variety of sites worldwide. These observations rely on spectroscopic parameters that are not known with sufficient accuracy to compute total columns that can be used in combination with in situ measurements. The TCCON must therefore be calibrated to World Meteorological Organization (WMO) in situ trace gas measurement scales. We present a calibration of TCCON data using WMO-scale instrumentation aboard aircraft that measured profiles over four TCCON stations during 2008 and 2009. These calibrations are compared with similar observations made in 2004 and 2006. The results indicate that a single, global calibration factor for each gas accurately captures the TCCON total column data within error.