Heritage of the MARC and REFIR

retrieval codes

The KLIMA algorithm, which has been adapted for processing IASI data, is the Millimetre-Wave Atmospheric Retrieval Code (MARC) (Carli et al., 2007) originally written to analyze the Millimetre-wave Airborne Receiver for Spectroscopic CHaracterization of Atmospheric Limb-Sounding (MARSCHALS) measurements (Oldfield et al., 2001). The MARC software was designed for the analysis of limb-sounding observations in the millimetre-waves and substantial modifications were required to extend its capabilities to the processing of data acquired by passive sensors that operate at higher frequencies and with a different viewing geometry. An upgraded version of the code was implemented and tested for application to the inversion of nadir-sounding observations performed in the middle/far infrared spectral region by the balloon-borne FT spectrometer REFIR-PAD (Palchetti et al., 2007 and Bianchini et al., 2007). The code was used to retrieve temperature and water vapour information by simultaneously fitting the water vapour profile, the temperature profile and the Earth skin brightness temperature using the spectrum from 100 to 1000 cm-1 (Del Bianco et al., 2007). The new features introduced in the MARC code for the analysis of the REFIR-PAD measurements were mostly driven by the need to extend the forward model capabilities to the frequency range (100-1400 cm-1) and to viewing mode (nadir sounding) of the instrument. Accordingly, the MARC code was modified to read the HITRAN database (Rothman et al., 2005 and Gordon et al., 2007), and the related atmospheric continuum was implemented (Clough et al., 2005), while to take into account the CO2 line-mixing a dedicated database and line-by-line cross section routines was used (Niro et al., 2005a,b).

Moreover, the new observation geometry required the implementation of the retrieval of surface temperature. As for the MARC code, the surface was modelled as a black body multiplied by a frequency independent emissivity. Different instrumental features had to be modelled: the side band effects included in the model of the MARSCHALS heterodyne spectrometer were neglected and new pertinent ILS and FOV functions were taken into account. Spectral radiance and frequency units commonly adopted in the millimetre-wave range (Kelvin and GHz, respectively) and used therefore in the original version of the code were modified in the REFIR version, where frequency is expressed in cm-1 and spectral radiance in nW/(cm2 sr cm−1) . The evolution of the MARC code into the version used for the analysis of REFIR spectra constituted the basis for the development of the retrieval code labelled as KLIMA and aimed at processing Level 1 IASI data. The KLIMA-IASI code is the result of further upgrading of both the forward model, to perform accurate and efficient RT calculation taking into account the IASI instrument characteristics, and of the inverse model to retrieve CO2 concentration.

In Table 1 the main differences between the REFIR and IASI instrument are reported.

 

REFIR

IASI

Type of Instrument

Balloonborne FT spectrometer

Spaceborne FT spectrometer

Spectral  Range

100 - 1400 cm-1

645 - 2700 cm-1

Spectral Resolution

0.475 cm-1

0.250 cm-1

Viewing Geometry

Nadir

Nadir ±48.3°

Viewing Altitude

35 km

(stratospheric balloon)

817 km

(MetOp-A satellite)

IFOV

133 mrad

(5 km)

14.65 mrad

(12 km at nadir)

NESR (at 800 cm-1)

200 nW/cm2/cm-1/sr

40 nW/cm2/cm-1/sr

Interferogram

acquisition time

32 s

151 ms

Table 1 - Main differences between REFIR and IASI instruments

In Table 2 the target products retrieved from the two instruments using the KLIMA inversion algorithm are listed.

 

Retrieval Product

REFIR

IASI

Temperature Profile

YES

YES

Water Vapour Profile

YES

YES

Ozone Column

NO

YES

Earth skin temperature

YES

YES

Carbone Dioxide Profile

NO

YES

Table 2 - Main differences between REFIR and IASI target products