Total Carbon Column Observing Network Research Articles

The Total Carbon Column Observing Network's GGG2020 data version

Abstract. The Total Carbon Column Observing Network (TCCON) measures column-average mole fractions of several greenhouse gases (GHGs), beginning in 2004, from over 30 current or past measurement sites around the world using solar absorption spectroscopy in the near-infrared (near-IR) region. TCCON GHG data have been used extensively for multiple purposes, including in studies of the carbon cycle and anthropogenic emissions, as well as to validate and improve observations from space-based sensors. Here, we describe an update to the retrieval algorithm used to process the TCCON near-IR solar spectra and to generate the associated data products. This version, called GGG2020, was initially released in April 2022. It includes updates and improvements to all steps of the retrieval, including but not limited to the conversion of the original interferograms into spectra, the spectroscopic information used in the column retrieval, post hoc air mass dependence correction, and scaling to align with the calibration scales of in situ GHG measurements. All TCCON data are available through https://tccondata.org/ (last access: 22 April 2024) and are hosted on CaltechDATA (https://data.caltech.edu/, last access: 22 April 2024). Each TCCON site has a unique DOI for its data record. An archive of all the sites' data is also available with the DOI https://doi.org/10.14291/TCCON.GGG2020 (Total Carbon Column Observing Network (TCCON) Team, 2022). The hosted files are updated approximately monthly, and TCCON sites are required to deliver data to the archive no later than 1 year after acquisition. Full details of data locations are provided in the “Code and data availability” section.

Advantages of assimilating multispectral satellite retrievals of atmospheric composition: a demonstration using MOPITT carbon monoxide products

Abstract. The Measurements Of Pollution In The Troposphere (MOPITT) is an ideal instrument to understand the impact of (1) assimilating multispectral and joint retrievals versus single spectral products, (2) assimilating satellite profile products versus column products, and (3) assimilating multispectral and joint retrievals versus assimilating individual products separately. We use the Community Atmosphere Model with chemistry with the Data Assimilation Research Testbed (CAM-chem+DART) to assimilate different MOPITT carbon monoxide (CO) products to address these three questions. Both anthropogenic and fire CO emissions are optimized in the data assimilation experiments. The results are compared with independent CO observations from TROPOspheric Monitoring Instrument (TROPOMI), the Total Carbon Column Observing Network (TCCON), NOAA Carbon Cycle Greenhouse Gases (CCGG) sites, In-service Aircraft for a Global Observing System (IAGOS), and Western wildfire Experiment for Cloud chemistry, Aerosol absorption and Nitrogen (WE-CAN). We find that (1) assimilating the MOPITT joint (multispectral; near-IR and thermal IR) column product leads to better model–observation agreement at and near the surface than assimilating the MOPITT thermal-IR-only column retrieval. (2) Assimilating column products has a larger impact and improvement for background and large-scale CO compared to assimilating profile products due to vertical localization in profile assimilation. However, profile assimilation can outperform column assimilations in fire-impacted regions and near the surface. (3) Assimilating multispectral and joint products results in similar or slightly better agreement with observations compared to assimilating the single spectral products separately.

Optimizing the Atmospheric CO2 Retrieval Based on the NDACC-Type FTIR Mid-Infrared Spectra at Xianghe, China

Optimizing the Atmospheric CO2 Retrieval Based on the NDACC-Type FTIR Mid-Infrared Spectra at Xianghe, China

Ship- and aircraft-based XCH4 over oceans as a new tool for satellite validation

Abstract. Satellite-based estimations of dry-air column-averaged mixing ratios of methane (XCH4) contribute to a better understanding of changes in CH4 emission sources and variations in its atmospheric growth rates. High accuracy of the satellite measurements is required, and therefore, extensive validation is performed, mainly against the Total Carbon Column Observing Network (TCCON). However, validation opportunities at open-ocean areas outside the coastal regions are sparse. We propose a new approach to assess the accuracy of satellite-derived XCH4 trends and variations. We combine various ship and aircraft observations with the help of atmospheric chemistry models, mainly used for the stratospheric column, to derive observation-based XCH4 (obs. XCH4). Based on our previously developed approach for the application to XCO2, we investigated three different advancements, from a simple approach to more elaborate approaches (approaches 1, 2, and 3), to account for the higher tropospheric and stratospheric variability in CH4 as compared to CO2. Between 2014 and 2018, at 20–40° N of the western Pacific, we discuss the uncertainties in the approaches and the derived obs. XCH4 within 10° by 20° latitude–longitude boxes. Uncertainties were 22 ppb (parts per billion) for approach 1, 20 ppb for approach 2, and 16 ppb for approach 3. We analyzed the consistency with the nearest TCCON stations and found agreement of approach 3 with Saga of 1±12 ppb and -1±11 ppb with Tsukuba for the northern and southern latitude box, respectively. Furthermore, we discuss the impact of the modeled stratospheric column on the derived obs. XCH4 by applying three different models in our approaches. Depending on the models, the difference can be more than 12 ppb (0.6 %), showing the importance for the appropriate choice. We show that our obs. XCH4 dataset accurately captures seasonal variations in CH4 over the ocean. Using different retrievals of the Greenhouse Gases Observing Satellite (GOSAT) from the National Institute for Environmental Studies (NIES), the RemoTeC full-physics retrieval operated at the Netherlands Institute for Space Research (SRON), and the full-physics retrieval of the University of Leicester (UoL-OCFP), we demonstrate the applicability of the dataset for satellite evaluation. The comparison with results of approach 3 revealed that NIES showed a difference of −0.04 ± 13 ppb and strong scatter at 20–30° N, while RemoTeC and OCFP have a rather systematic negative bias of −12.1 ± 8.1 and −10.3 ± 9.6 ppb. Our new approach to derive XCH4 reference datasets over the ocean can contribute to the validation of existing and upcoming satellite missions in future.

Comparison of Atmospheric Carbon Dioxide Concentrations Based on GOSAT, OCO-2 Observations and Ground-Based TCCON Data

Carbon dioxide (CO2) is one of the most significant greenhouse gases, and its concentration and distribution in the atmosphere have always been a research hotspot. To study the temporal and spatial characteristics of atmospheric CO2 globally, it is crucial to evaluate the consistency of observation data from different carbon observation satellites. This study utilizes data from the Total Carbon Column Observing Network (TCCON) to verify the column-averaged dry air mole fractions of atmospheric CO2 (XCO2) retrieved by satellites from October 2014 to May 2016, specifically comparing the XCO2 distributions from the Greenhouse Gases Observing Satellite (GOSAT) and Orbiting Carbon Observatory 2 (OCO-2). Our analysis indicates a strong correlation between the TCCON and both the GOSAT (correlation coefficient of 0.85) and OCO-2 (correlation coefficient of 0.91). Cross-validation further reveals that the measurements of the GOSAT and OCO-2 are highly consistent, with an average deviation and standard deviation of 0.92 ± 1.16 ppm and a correlation coefficient of 0.92. These differences remain stable over time, indicating that the calibration in the data set is reliable. Moreover, monthly averaged time-series and seasonal climatology comparisons were also performed separately over the six continents, i.e., Asia, North America, Europe, Africa, South America, and Oceania. The investigation of monthly XCO2 values across continents highlights greater consistency in Asia, North America, and Oceania (standard deviation from 0.15 to 0.27 ppm) as compared to Europe, South America, and Africa (standard deviation from 0.45 to 0.84 ppm). A seasonal analysis exhibited a high level of consistency in spring (correlation coefficient of 0.97), but lower agreement in summer (correlation coefficient of 0.78), potentially due to cloud cover and aerosol interference. Although some differences exist among the datasets, the overall findings demonstrate a strong correlation between the satellite measurements of XCO2. These results emphasize the importance of continued monitoring and calibration efforts to ensure the accurate assessment and understanding of atmospheric CO2 levels.

Global Evaluation and Intercomparison of XCO2 Retrievals from GOSAT, OCO-2, and TANSAT with TCCON

Accurate global monitoring of carbon dioxide (CO2) is essential for understanding climate change and informing policy decisions. This study compares column-averaged dry-air mole fractions of CO2 (XCO2) between ACOS_L2_Lite_FP V9r for Japan’s Greenhouse Gases Observing Satellite (GOSAT), OCO-2_L2_Lite_FP V10r for the USA’s Orbiting Carbon Observatory-2 (OCO-2), and IAPCAS V2.0 for China’s Carbon Dioxide Observation Satellite (TANSAT) collectively referred to as GOT, with data from the Total Carbon Column Observing Network (TCCON). Our findings are as follows: (1) Significant data quantity differences exist between OCO-2 and the other satellites, with OCO-2 boasting a data volume 100 times greater. GOT shows the highest data volume between 30–45°N and 20–30°S, but data availability is notably lower near the equator. (2) XCO2 from GOT exhibits similar seasonal variations, with lower concentrations during June, July, and August (JJA) (402.72–403.74 ppm) and higher concentrations during December, January, and February (DJF) (405.74–407.14 ppm). XCO2 levels are higher in the Northern Hemisphere during March, April, and May (MAM) and DJF, while slightly lower during JJA and September, October, and November (SON). (3) The differences in XCO2 (ΔXCO2) reveal that ΔXCO2 between OCO-2 and TANSAT are minor (−0.47 ± 0.28 ppm), whereas the most significant difference is observed between GOSAT and TANSAT (−1.13 ± 0.15 ppm). Minimal differences are seen in SON (with the biggest difference between GOSAT and TANSAT: −0.84 ± 0.12 ppm), while notable differences occur in DJF (with the biggest difference between GOSAT and TANSAT: −1.43 ± 0.17 ppm). Regarding latitudinal variations, distinctions between OCO-2 and TANSAT are most pronounced in JJA and SON. (4) Compared to TCCON, XCO2 from GOT exhibits relatively high determination coefficients (R2 > 0.8), with GOSAT having the highest root mean square error (RMSE = 1.226 ppm, <1.5 ppm), indicating a strong relationship between ground-based observed and retrieved values. This research contributes significantly to our understanding of the spatial characteristics of global XCO2. Furthermore, it offers insights that can inform the analysis of differences in the inversion of carbon sources and sinks within assimilation systems when incorporating XCO2 data from satellite observations.

Inter-comparison and evaluation of global satellite XCO2 products

ABSTRACT Carbon dioxide (CO2) is one of the main greenhouse gases and has become a major concern as its concentration has been growing in recent years. Satellite remote sensing is an efficient way to monitor CO2 in the atmosphere, and several satellites are already used for CO2 monitoring. It is imperative to investigate the spatial coverage and spatio-temporal trends of satellite products, as well as identify the satellites with higher levels of accuracy. Additionally, examining the disparities between the older and new generations of satellites would be meaningful. Therefore, this paper provides a comprehensive evaluation and inter-comparison for the commonly used satellite column-averaged dry-air mole fraction of CO2 (XCO2) products. Specifically, the temporal trends and monthly coverage of the Greenhouse Gases Observing SATellite (GOSAT), Greenhouse Gases Observing SATellite-2 (GOSAT-2), Orbiting Carbon Observatory-2 (OCO-2), Orbiting Carbon Observatory-3 (OCO-3), and SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) are investigated. The accuracy of these satellite products is evaluated and analyzed based on Total Carbon Column Observing Network (TCCON) data. The results indicate that the XCO2 of all the satellite products show a year-by-year increase, with seasonal periodicity. In terms of overall accuracy, the OCO series satellites exhibit a slightly higher level of accuracy compared to the GOSAT series. The products of the new generation of satellites are less stable than those of the older generation, probably due to the impacts of the inversion algorithm and platforms.

A retrieval of xCO2 from ground-based mid-infrared NDACC solar absorption spectra and comparison to TCCON

Abstract. Two global networks of ground-based Fourier transform spectrometers are measuring abundances of atmospheric trace gases that absorb in the near infrared and mid-infrared: the Network for the Detection of Atmospheric Composition Change (NDACC) and the Total Carbon Column Observing Network (TCCON). The first lacks a CO2 product; therefore, this study focuses on developing an xCO2 retrieval method for NDACC from a spectral window in the 4800 cm−1 region. This retrieval will allow extending ground-based measurements back in time, which we will demonstrate with historical data available from Ny-Ålesund, Svalbard. At this site, both TCCON and NDACC measurements are routinely performed, which is an advantage for collocated comparisons. The results are compared with collocated TCCON measurements of column-averaged dry-air mole fractions of CO2 (denoted by xCO2) in Ny-Ålesund, Svalbard, and only TCCON in Burgos, Philippines. We found that it is possible to retrieve xCO2 from NDACC spectra with a precision of 0.2 %. The comparison between the new retrieval and TCCON showed that the sensitivity of the new retrieval is high in the troposphere and lower in the upper stratosphere, similar to TCCON, as seen in the averaging kernels, and that the seasonality is well captured as seen in the retrieved time series. Additionally, we have included a retrieval strategy suggestion to improve the quality of the xCO2 product.

A blended TROPOMI+GOSAT satellite data product for atmospheric methane using machine learning to correct retrieval biases

Abstract. Satellite observations of dry-column methane mixing ratios (XCH4) from shortwave infrared (SWIR) solar backscatter radiation provide a powerful resource to quantify methane emissions in service of climate action. The TROPOspheric Monitoring Instrument (TROPOMI), launched in October 2017, provides global daily coverage at a 5.5 × 7 km2 (nadir) pixel resolution, but its methane retrievals can suffer from biases associated with SWIR surface albedo, scattering from aerosols and cirrus clouds, and across-track variability (striping). The Greenhouse gases Observing SATellite (GOSAT) instrument, launched in 2009, has better spectral characteristics and its methane retrieval is much less subject to biases, but its data density is 250 times sparser than TROPOMI. Here, we present a blended TROPOMI+GOSAT methane product obtained by training a machine learning (ML) model to predict the difference between TROPOMI and GOSAT co-located measurements, using only predictor variables included in the TROPOMI retrieval, and then applying the correction to the complete TROPOMI record from April 2018 to present. We find that the largest corrections are associated with coarse aerosol particles, high SWIR surface albedo, and across-track pixel index. Our blended product corrects a systematic difference between TROPOMI and GOSAT over water, and it features corrections exceeding 10 ppb over arid land, persistently cloudy regions, and high northern latitudes. It reduces the TROPOMI spatially variable bias over land (referenced to GOSAT data) from 14.3 to 10.4 ppb at a 0.25∘ × 0.3125∘ resolution. Validation with Total Carbon Column Observing Network (TCCON) ground-based column measurements shows reductions in variable bias compared with the original TROPOMI data from 4.7 to 4.4 ppb and in single-retrieval precision from 14.5 to 11.9 ppb. TCCON data are all in locations with a SWIR surface albedo below 0.4 (where TROPOMI biases tend to be relatively low), but they confirm the dependence of TROPOMI biases on SWIR surface albedo and coarse aerosol particles, as well as the reduction of these biases in the blended product. Fine-scale inspection of the Arabian Peninsula shows that a number of hotspots in the original TROPOMI data are removed as artifacts in the blended product. The blended product also corrects striping and aerosol/cloud biases in single-orbit TROPOMI data, enabling better detection and quantification of ultra-emitters. Residual coastal biases can be removed by applying additional filters. The ML method presented here can be applied more generally to validate and correct data from any new satellite instrument by reference to a more established instrument.

Approximation of multi-year time series of XCO2 concentrations using satellite observations and statistical interpolation methods

Long-term time series of CO2 concentrations are required for studying e.g., the effects of rising CO2 concentrations on plants and ecosystems. Different methods for approximating column-averaged dry-air mole fractions of CO2 (XCO2) have been developed to substitute missing local CO2 measurements.This study presents a novel method (XCO2SAT+) for determining monthly XCO2 concentrations for selected sites, based entirely on observations from multiple satellites as well as statistical methods in the case of missing satellite data. The product currently covers a period of 16 years from 2003 to 2018. XCO2SAT+ results were validated against data of five European sites of the Total Carbon Column Observing Network (TCCON) and compared to two established XCO2 approximation methods. All three approximation methods were compared to CO2 measurements of German towers from the Integrated Carbon Observation System (ICOS).The mean absolute error (MAE) between the XCO2SAT+ method and TCCON site Karlsruhe, Germany ranges between 0.5 and 1.1 ppm (0.1–0.3%) for years 2014–2018. Mean differences between XCO2SAT+ and XCO2 approximation method CAMS-GCM shows a slightly better MAE of 0.1–0.3 ppm for CAMS-GCM at two German TCCON sites. Multi-year comparisons with ground-level ICOS measurements revealed a comparable seasonality of CO2 concentrations during the vegetation period, but underestimated high winter concentrations, a pattern also observed when comparing with CAMS-GCM. XCO2SAT+ is a suitable alternative and simple method for providing long-term time series of monthly averaged XCO2 concentrations for selected sites, expanding the methodological toolbox for in-situ CO2 research. The new method, which so far provides XCO2 concentrations, will thus serve as a first step towards providing further CO2 data sources, which could be used to further improve numerical models such as atmospheric models or crop growth models and to study CO2 impacts in agriculture and natural ecosystems.

A XCO2 Retrieval Algorithm Coupled Spatial Correlation for the Aerosol and Carbon Detection Lidar

We describe the approach to estimating the atmospheric carbon dioxide (CO2) for the Aerosol and Carbon Detection Lidar (ACDL) onboard the Atmospheric Environment Monitoring Satellite (AEMS). The method estimates the optimal state vector by maximizing the measurement posterior probability under a given prior state vector probability distribution. A priori constraint considering the spatial correlation is used as regularization to solve the ill-posed problem. We ran a series of observing system simulation experiments to demonstrate the critical outcome and character percentage uncertainty reduction. The results show that the state vector uncertainty can be reduced by ∼ 10% near the surface for the single sounding. The CO2 column-averaged dry air mole fraction (XCO2) derived by this algorithm is more stable than that obtained by the conventional algorithm and enables the monitoring of concentration changes for the multiple soundings. Similar to the Total Carbon Column Observing Network (TCCON), the averaging kernel is also provided for the subsequent flux inversion. Our simulation experiments demonstrate that the structure of the prior error covariance plays an important role in revealing vertical information from observations. In addition, we applied this algorithm to an airborne ACDL experiment for the retrieval of atmospheric CO2 over Bohai Bay on March 14, 2019. AEMS’s observations with a small footprint will yield important information on the carbon cycle, especially for small but strong emission sources.

Evaluating the consistency between OCO-2 and OCO-3 XCO2 estimates derived from the NASA ACOS version 10 retrieval algorithm

Abstract. The version 10 (v10) Atmospheric Carbon Observations from Space (ACOS) Level 2 full-physics (L2FP) retrieval algorithm has been applied to multiyear records of observations from NASA's Orbiting Carbon Observatory 2 and 3 sensors (OCO-2 and OCO-3, respectively) to provide estimates of the carbon dioxide (CO2) column-averaged dry-air mole fraction (XCO2). In this study, a number of improvements to the ACOS v10 L2FP algorithm are described. The post-processing quality filtering and bias correction of the XCO2 estimates against multiple truth proxies are also discussed. The OCO v10 data volumes and XCO2 estimates from the two sensors for the time period of August 2019 through February 2022 are compared, highlighting differences in spatiotemporal sampling but demonstrating broad agreement between the two sensors where they overlap in time and space. A number of evaluation sources applied to both sensors suggest they are broadly similar in data and error characteristics. Mean OCO-3 differences relative to collocated OCO-2 data are approximately 0.2 and −0.3 ppm for land and ocean observations, respectively. Comparison of XCO2 estimates to collocated Total Carbon Column Observing Network (TCCON) measurements shows root mean squared errors (RMSEs) of approximately 0.8 and 0.9 ppm for OCO-2 and OCO-3, respectively. An evaluation against XCO2 fields derived from atmospheric inversion systems that assimilated only near-surface CO2 observations, i.e., did not assimilate satellite CO2 measurements, yielded RMSEs of 1.0 and 1.1 ppm for OCO-2 and OCO-3, respectively. Evaluation of uncertainties in XCO2 over small areas, as well as XCO2 biases across land–ocean crossings, also indicates similar behavior in the error characteristics of both sensors. Taken together, these results demonstrate a broad consistency of OCO-2 and OCO-3 XCO2 measurements, suggesting they may be used together for scientific analyses.

Spatiotemporal variation analysis of global XCO2 concentration during 2010–2020 based on DINEOF-BME framework and wavelet function

Combining with Carbon dioxide column concentration (XCO2) remote sensing data, it is of great scientific significance to obtain XCO2 long time series data with high precision and high spatio-temporal coverage. In this study, the combination framework of DINEOF and BME were employed to integrate the XCO2 data of GOSAT, OCO-2 and OCO-3 satellites for generating global XCO2 data from January 2010 to December 2020, with the average monthly space coverage rate of more than 96 %. Through cross-validation and comparison of The Total Carbon Column Observing Network (TCCON) XCO2 data and DINEOF-BME interpolation XCO2 products, it is verified that DINEOF-BME method has better interpolation accuracy, and the coefficient of determination of interpolated XCO2 products and TCCON data is 0.920. The long time series of global XCO2 products showed a wave rising trend, with a total increase of ~23 ppm; obviously seasonal characteristics were also detected with the highest XCO2 value in spring and the lowest in autumn. According to the zonal integration analysis, the values of XCO2 in the northern hemisphere is higher than the southern hemisphere during January–May and October–December, while the values of XCO2 in the southern hemisphere is higher than the northern hemisphere during June–September, which accords with the seasonal law. Through EOF mapping, the first mode accounted for 88.93 % of the total variability, and its variation trend is consistent with that of XCO2 concentration, which verifies the variation rule of XCO2 from the time and space pattern. Through wavelet analysis, the time scale corresponding to the first main cycle of XCO2 change is 59-month, which has obvious regularity on the time scale. DINEOF-BME technology framework has good generality, while XCO2 long time series data products and the spatio-temporal variation of XCO2 revealed by the research provide a solid theoretical basis and data support for related research.

Inferring the vertical distribution of CO and CO2 from TCCON total column values using the TARDISS algorithm

Abstract. We describe an approach for determining limited information about the vertical distribution of carbon monoxide (CO) and carbon dioxide (CO2) from total column ground-based Total Carbon Column Observation Network (TCCON) observations. For CO and CO2, it has been difficult to retrieve information about their vertical distribution from spectral line shapes because of the errors in the spectroscopy and the atmospheric temperature profile that mask the effects of variations in their mixing ratio with altitude. For CO2 the challenge is especially difficult given that these variations are typically 2 % or less. Nevertheless, if sufficient accuracy can be obtained, such information would be highly valuable for evaluation of retrievals from satellites and more generally for improving the estimate of surface sources and sinks of these trace gases. We present here the Temporal Atmospheric Retrieval Determining Information from Secondary Scaling (TARDISS) retrieval algorithm. TARDISS uses several simultaneously obtained total column observations of the same gas from different absorption bands with distinctly different vertical averaging kernels. The different total column retrievals are combined in TARDISS using a Bayesian approach where the weights and temporal covariance applied to the different retrievals include additional constraints on the diurnal variation in the vertical distribution for these gases. We assume that the near-surface part of the column varies rapidly over the course of a day (from surface sources and sinks, for example) and that the upper part of the column has a larger temporal covariance over the course of a day. Using measurements from the five North American TCCON sites, we find that the retrieved lower partial column (between the surface and ∼ 800 hPa) of the CO and CO2 dry mole fractions (DMFs) have slopes of 0.999 ± 0.002 and 1.001 ± 0.003 with respect to lower column DMF from integrated in situ data measured directly from aircraft and in AirCores. The average error for our lower column CO retrieval is 1.51 ppb (∼ 2 %) while the average error for our CO2 retrieval is 5.09 ppm (∼ 1.25 %). Compared with classical line-shape-derived vertical profile retrievals, our algorithm reduces the influence of forward model errors such as imprecision in spectroscopy (line shapes and intensities) and in the instrument line shape. In addition, because TARDISS uses the existing retrieved column abundances from TCCON (which themselves are computationally much less intensive than profile retrieval algorithms), it is very fast and processes years of data in minutes. We anticipate that this approach will find broad application for use in carbon cycle science.

Potential improvement of XCO2 retrieval of the OCO-2 by having aerosol information from the A-train satellites

ABSTRACT Near-real time observations of aerosol properties could have a potential to improve the accuracy of XCO2 retrieval algorithm in operational satellite missions. In this study, we developed a retrieval algorithm of XCO2 (Yonsei Retrieval Algorithm; YCAR) based on the Optimal Estimation (OE) method that used aerosol information at the location of the Orbiting Carbon Observatory-2 (OCO-2) measurement from co-located measurement of the Afternoon constellation (A-train) such as the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard the Cloud-Aerosol Lidar and Infrared Pathfinder Observation (CALIPSO) and the MODerate-resolution Imaging Spectrometer (MODIS) onboard the Aqua. Specifically, we used optical depth, vertical profile, and optical properties of aerosol from MODIS and CALIOP data. We validated retrieval results to the Total Carbon Column Observing Network (TCCON) ground-based measurements and found general consistency. The impact of observed aerosol information and its constraint was examined by retrieval tests using different settings. The effect of using additional aerosol information was analyzed in connection with the bias correction process of the operational retrieval algorithm. YCAR using a priori aerosol loading parameters from co-located satellite measurements and less constraint of aerosol optical properties made comparable results with operational data with the bias correction process in three of the four cases subject to this study. Our work provides evidence supporting the bias correction process of operational algorithms and quantitatively presents the effectiveness of synergic use of multiple satellites (e.g. A-train) and better treatment of aerosol information.

Characteristics of interannual variability in space-based XCO2 global observations

Abstract. Atmospheric carbon dioxide (CO2) accounts for the largest radiative forcing among anthropogenic greenhouse gases. There is, therefore, a pressing need to understand the rate at which CO2 accumulates in the atmosphere, including the interannual variations (IAVs) in this rate. IAV in the CO2 growth rate is a small signal relative to the long-term trend and the mean annual cycle of atmospheric CO2, and IAV is tied to climatic variations that may provide insights into long-term carbon–climate feedbacks. Observations from the Orbiting Carbon Observatory-2 (OCO-2) mission offer a new opportunity to refine our understanding of atmospheric CO2 IAV since the satellite can measure over remote terrestrial regions and the open ocean, where traditional in situ CO2 monitoring is difficult, providing better spatial coverage compared to ground-based monitoring techniques. In this study, we analyze the IAV of column-averaged dry-air CO2 mole fraction (XCO2) from OCO-2 between September 2014 and June 2021. The amplitude of the IAV, which is calculated as the standard deviation of the time series, is up to 1.2 ppm over the continents and around 0.4 ppm over the open ocean. Across all latitudes, the OCO-2-detected XCO2 IAV shows a clear relationship with El Niño–Southern Oscillation (ENSO)-driven variations that originate in the tropics and are transported poleward. Similar, but smoother, zonal patterns of OCO-2 XCO2 IAV time series compared to ground-based in situ observations and with column observations from the Total Carbon Column Observing Network (TCCON) and the Greenhouse Gases Observing Satellite (GOSAT) show that OCO-2 observations can be used reliably to estimate IAV. Furthermore, the extensive spatial coverage of the OCO-2 satellite data leads to smoother IAV time series than those from other datasets, suggesting that OCO-2 provides new capabilities for revealing small IAV signals despite sources of noise and error that are inherent to remote-sensing datasets.

Spatial-temporal variation in XCH4 during 2009–2021 and its driving factors across the land of the Northern Hemisphere

Methane (CH4) is the second most significant driver of global warming, following carbon dioxide. However, the spatial-temporal variation of CH4 and its driving factors largely remain unclear. Here we selected the Northern Hemisphere as the study area. We used the data from the Total Column Carbon Observing Network (TCCON) to assess the accuracy of the Greenhouse Gases Observing Satellite (GOSAT) Proxy XCH4 (column-averaged dry air mixing ratio of CH4) data. We then analyzed the spatial-temporal distribution of XCH4 in the Northern Hemisphere, and further quantified the influencing factors using geographic detectors. The results showed that during 2009–2021, the annual mean XCH4 increased from 2009 (1775.19 ppb) to 2021 (1872.71 ppb), with an increasing rate of 7.50 ppb/year. The monthly average value was the lowest in May (1805.65 ppb) and the highest in September (1825.63 ppb). The XCH4 in the low-latitude region was higher than that in the high-latitudinal region. The geographic detector showed that anthropogenic activities were the main factors affecting the XCH4. Our results revealed the spatial-temporal patterns XCH4 and their driving factors in the Northern Hemisphere, and thus provided a scientific basis for the management of this greenhouse gas in the future.

A Coupled CH4, CO and CO2 Simulation for Improved Chemical Source Modeling

Understanding greenhouse gas–climate processes and feedbacks is a fundamental step in understanding climate variability and its links to greenhouse gas fluxes. Chemical transport models are the primary tool for linking greenhouse gas fluxes to their atmospheric abundances. Hence, accurate simulations of greenhouse gases are essential. Here, we present a new simulation in the GEOS-Chem chemical transport model that couples the two main greenhouse gases—carbon dioxide (CO2) and methane (CH4)—along with the indirect greenhouse gas carbon monoxide (CO) based on their chemistry. Our updates include the online calculation of the chemical production of CO from CH4 and the online production of CO2 from CO, both of which were handled offline in the previous versions of these simulations. In the newly developed coupled (online) simulation, we used consistent hydroxyl radical (OH) fields for all aspects of the simulation, resolving biases introduced by inconsistent OH fields in the currently available uncoupled (offline) CH4, CO and CO2 simulations. We compare our coupled simulation with the existing v12.1.1 GEOS-Chem uncoupled simulations run the way they are currently being used by the community. We discuss differences between the uncoupled and coupled calculation of the chemical terms and compare our results with surface measurements from the NOAA Global Greenhouse Gas Reference Network (NOAA GGGRN), total column measurements from the Total Carbon Column Observing Network (TCCON) and aircraft measurements from the Atmospheric Tomography Mission (ATom). Relative to the standard uncoupled simulations, our coupled results suggest a stronger CO chemical production from CH4, weaker production of CO2 from CO and biases in the OH fields. However, we found a significantly stronger chemical production of CO2 in tropical land regions, especially in the Amazon. The model–measurement differences point to underestimated biomass burning emissions and secondary production for CO. The new self-consistent coupled simulation opens new possibilities when identifying biases in CH4, CO and CO2 source and sink fields, as well as a better understanding of their interannual variability and co-variation.

Spatiotemporal Analysis of XCO2 and Its Relationship to Urban and Green Areas of China’s Major Southern Cities from Remote Sensing and WRF-Chem Modeling Data from 2010 to 2019

Monitoring CO2 concentrations is believed to be an effective measure for assisting in the control of greenhouse gas emissions. Satellite measurements compensate for the sparse and uneven spatial distribution of ground observation stations, allowing for the collection of a wide range of CO2 concentration data. However, satellite monitoring’s spatial coverage remains limited. This study fills the knowledge gaps of column-averaged dry-air mole fraction of CO2 (XCO2) products retrieved from the Greenhouse Gases Observing Satellite (GOSAT) and Orbiting Carbon Observatory Satellite (OCO-2) based on the normalized output of atmospheric chemical models, WRF-Chem, in Southern China during 2010–2019. Hefei (HF)/Total Carbon Column Observing Network (TCCON), Lulin (LLN)/World Data Centre for Greenhouse Gases (WDCGG) station observations were used to validate the results of void filling with an acceptable accuracy for spatiotemporal analysis (R = 0.96, R2 = 0.92, RMSE = 2.44 ppm). Compared to the IDW (inverse distance weighting) and Kriging (ordinary Kriging) interpolation methods, this method has a higher validation accuracy. In addition, spatiotemporal distributions of CO2, as well as the sensitivity of CO2 concentration to the urban built-up areas and urban green space areas in China’s major southern cities during 2010–2019, are discussed. The approximate annual average concentrations have gradually increased from 388.56 to 414.72 ppm, with an annual growth rate of 6.73%, and the seasonal cycle presents a maximum in spring and a minimum in summer or autumn from 2010 to 2019. CO2 concentrations have a strong positive correlation with the impervious area to city area ratio, while anomaly values of the impervious area to urban green area ratio occurred in individual cities. The experimental findings demonstrate the viability of the study hypothesis that combines remote sensing data with the WRF-Chem model to produce a local area dataset with high spatial resolution and an extracted urban unit from statistical data.

Accounting for surface reflectance spectral features in TROPOMI methane retrievals

Abstract. Satellite remote sensing of methane (CH4) using the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Copernicus Sentinel-5 Precursor (S5-P) satellite is key to monitor and quantify emissions globally. Overall, the S5-P methane data are of satisfying quality, demonstrated by the product validation with ground-based measurements from the Total Carbon Column Observing Network (TCCON). However, analysis of TROPOMI dry-air column mixing ratio (XCH4) data has pointed to a few false methane anomalies up to 20–40 ppb that can potentially be misinterpreted as enhancements due to strong emission sources. These artefacts are caused by spectral features of the underlying surfaces of specific materials (e.g. carbonate rocks), which are not well represented in the forward model of the retrieval algorithm. In this study we show that the observed anomalies are due to the surface model which describes the spectral dependence of the Lambertian albedo by a second-order polynomial in wavelength. By analysing the ECOSTRESS library that contains laboratory spectra for different types of surfaces, we find that a quadratic function might not be the most optimal representation of the surface reflectance spectral dependencies in the short-wave infrared (SWIR) range. Already the use of a third-order polynomial improves the methane data such that the anomalies disappear at several locations (e.g. Siberia, Australia and Algeria) without affecting the data quality elsewhere, and the quality of the fit significantly improves. We also found that the known bias in retrieved methane for low-albedo scenes slightly improves, but still, a posterior correction needs to be applied, leaving open the question about the root cause of the albedo bias. After applying the adjusted surface model globally, we perform the routine validation with TCCON and Greenhouse gases Observing SATellite (GOSAT) data. GOSAT comparison does not significantly improve, while TCCON validation results show a small improvement in some stations of 2–4 ppb, up to a factor of 10 smaller than the artificial XCH4 enhancements. This reflects that TCCON stations are not close to any of the corrected artefacts, hinting at a limitation of the current validation approach of the S5-P XCH4 data product.