Jan 24, 2017

New inverse algorithm for CO2 retrieval from satellite observations

Atmospheric Carbon Dioxide (CO2) is a primary greenhouse gas that has been persistently increasing over the past few decades. It is a major driver of regional and global climate change (IPCC, 2007).

Estimated Article Reading Time: 1 min.

Most CO2 sources and sinks are located near the surface. Observations from shortwave infrared sounders loaded on satellite, such as Greenhouse Gases Observing SATellite (GOSAT) and Orbiting Carbon Observatory (OCO-2) can provide accurate measurements of the column-averaged atmospheric CO2 concentration.

The atmospheric radiative transfer equation (RTE) is a Fredholm integral equation of the first kind, which is recognized as likely to be ill-conditioned. Thus the inverse problem based on the RTE is not well-posed. The inverse method for atmospheric gas profiles retrieval is usually based on optimization theory.

In an paper coauthored with Liangfu Chen, Shenshen Li, Meng Fan, Jinhua Tao, Ying Zhang, scholars at The State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, an improved constraining method for satellite CO2 remote sensing in the (short-wave infrared) SWIR band was stated, which combines a pre-processing step applied to the a priori state vector prior to retrievals, with the modified damped Newton method (MDNM).

MDNM contains two constraining factors that stabilize the retrieval iterations. The Levenberg-Marquardt parameter (?) is used to ensure a positive Hessian matrix, and a scale factor (a) is used to adjust the step size. The algorithm iteratively searches for an optimized solution using observed spectral radiances, and parameters (y and a) are appropriately adjusted.

A pre-processor for initializing the first guess (X0) prior to the retrievals, when the algorithm detects that X0 is far from the true state vector. A new data-screening method for detecting cloud scenes is also presented based on the different spectrum shapes in the oxygen-A and 2 micron bands.

This research aims to stabilize the retrieval iterations. And preliminary validations indicate that the quality of the MDNM-based retrieval results is relatively stable.