Algorithm Flowchart

Description of the CLIMCAPS Optimal Estimation inversion algorithm, a retrieval method to invert hyperspectral infrared radiance measurements into vertical atmospheric profiles.

CLIMCAPS Optimal Estimation retrieval

CLIMCAPS retrieves multiple Earth surface, cloud and atmospheric state variables from infrared and microwave measurements using Optimal Estimation inversion (Rodgers, 2000). We describe CLIMCAPS in detail elsewhere, with specific reference to uncertainty (Smith and Barnet, 2019) and information content (Smith and Barnet, 2020). Figure 1 outlines the main retrieval steps. Each box with a variable name represents an Optimal Estimation retrieval.

Flow chart of the CLIMCAPS retrieval algorithm

Figure 1: Flow diagram of the CLIMCAPS sequential retrieval algorithm. This gives a broad overview of the main retrieval steps and their logical flow towards two final products files, CLIMCAPS retrievals (CLIMCAPS RET) and cloud cleared radiances (CLIMCAPS CCR). Note that we discuss different aspects of CLIMCAPS algorithm flow also in (Smith and Barnet, 2019, 2020). See Table 1 for a description of the acronyms and symbols used here.

Table 1: List of acronyms and symbols used in Figure 1.

Table symbol

Description

\(\alpha\)

Cloud fraction retrieval

\(\varepsilon\)

Earth surface emissivity

\(\rho\)

Earth surface reflectivity

\(\delta\)T(p)***

Delta temperature: T(p) with subset of MW channels minus T(p)** from MW+IR channels. This step tests the quality of the T(p)** retrieval.

CAMEL

CLIMCAPS implementation of the Combined ASTER and MODIS Emissivity database over Land (Hook, 2019c, 2019a, 2019b)

CC

Cloud Clearing that includes retrieval of cloud fraction and cloud top pressure

CCR

Level 2 Cloud Cleared Radiance product

CD flags

Constituent Detection flags for isoprene, ethane, propylene and ammonia

CH4

Methane retrieval on 100 pressure layers using a subset of IR channels

Climatology

Global representation of atmospheric variables

CO

Carbon monoxide retrieval on 100 pressure layers using a subset of IR channels

CO2

Carbon dioxide retrieval on 100 pressure layers using a subset of IR channels

CTP

Cloud top pressure retrieval

H2O

Water vapor retrieval on 100 pressure layers using a subset of MW+IR channels

HNO3

Nitric acid retrieval on 100 pressure layers using a subset of IR channels

IR

Infrared

LAC

Local angle correction of IR radiances within 3 x 3 fields of view

Level 1

NASA geolocated, calibrated radiance products for IR and MW measurements

Level 2

NASA geophysical products retrieved from Level 1 radiance measurements.

LIQ

Liquid water path

Masuda

CLIMCAPS implementation of the Infrared sea surface emissivity model: (Masuda et al., 1988) as modified by (Wu and Smith, 1997)

MERRA2

Modern-Era Retrospective analysis for Research and Applications Version 2 (GMAO, 2015) collocated in time and space to the CLIMCAPS instrument footprints.

MW

Microwave

N2O

Nitrous oxide retrieval on 100 pressure layers using a subset of IR channels

O3

Ozone retrieval on 100 pressure layers using a subset of IR channels

PS

Surface pressure

RET

Level 2 geophysical retrieval product

SO2

Sulphur dioxide retrieval on 100 pressure layers using a subset of IR channels

T(p)*

Temperature retrieval on 100 pressure levels using a subset of MW+IR channels – first retrieval

T(p)**

Temperature retrieval on 100 pressure levels using a subset of MW+IR channels – second and final retrieval

TS

Surface skin temperature

 About the flowchart

The flowchart in Figure 1 outlines the main CLIMCAPS retrieval steps. We highlight only those steps that result in the primary retrieval products, namely cloud cleared radiances (CCR) and the nine profile retrievals – T(p), H2O, O3, CO, CO2, CH4, SO2, HNO3, N2O.

CLIMCAPS profiles are available on 100 pressure coordinates, with T(p) represented on 100 pressure levels (air_pres) and the gas species on 100 pressure layers (air_pres_lay). For the sake of simplicity, we omit reference to those steps where (i) diagnostic metrics and quality control indices are calculated, and (ii) quantities are derived from the main retrieval variables, such as outgoing long-wave radiation (OLR) and relative humidity.

CLIMCAPS output variables are saved to two product files in netCDF format:

  1. *.L2_CLIMCAPS_CCR.*: Cloud cleared radiances.
  2. *.L2_CLIMCAPS_RET.*: Retrieved and derived geophysical variables. This file is organized into two tiers, with the main group of variables as tier 1 (T(p), surface and cloud variables, relative humidity, etc.) and the remaining retrieved variables, derived quantities and diagnostic metrics organized into four subgroups on tier 2. We refer to three of these subgroups in Figure 1:
    • AUX: collection of auxiliary variables, diagnostic metrics as well as minor gas detection flags.
    • MOL_LAY: collection of the eight primary trace gas profile retrievals, H2O, O3, CO, CO2, CH4, SO2, HNO3, N2 With the exception of CO2 that is retrieved as volume mixing ratio, all trace gases are retrieved as layer column densities. The averaging kernel matrices for all nine retrieval variables (eight gases and T(p)) are contained in the AVE_KERN subgroup.
    • MW: collection of MW-only retrieval variables using the algorithm developed by (Rosenkranz, 2001, 2006)

CLIMCAPS uses several diagnostic metrics to derive a single, final quality control flag (0 = best, 1 = good, 2 = bad). Some reasons for rejection include:

  • The footprint is covered with precipitating clouds such that TMW(p) fails
  • The retrieved cloud fraction exceeds 80%
  • Large differences exist between the clear radiance estimate (from the MW-only step) and the cloud clear estimate (from the cloud clearing step)
  • The RMS of the observed minus calculated brightness temperature exceeds 1.75 K for specific window channels.
  • The boundary layer estimate of dT(p)*** > 1.5 K

None of the above checks alter the retrieved values, but instead are used to derive a final quality flag. CLIMCAPS, therefore, does not have quality control flags tailored for each individual retrieved variable, but instead applies the logic that if T(p) and H2O fail, all subsequent trace gas retrievals also fail.

CLIMCAPS and Microwave

CLIMCAPS does not use any of the MW channels in cloud clearing or cloud parameter retrievals (cloud fraction and cloud top pressure).

The MW-only retrievals of liquid water path (LIQ) and emissivity (\(\varepsilon\)MW) are used to define the background atmospheric state in subsequent IR+MW retrievals, while the MW-only temperature and water vapor retrievals are written to the output RET file as is; they are not used in subsequent IR+MW retrievals.

References

  • GMAO: MERRA-2 inst3_3d_asm_Nv: 3d,3-Hourly,Instantaneous,Model-Level,Assimilation,Assimilated Meteorological Fields V5.12.4, https://doi.org/10.5067/WWQSXQ8IVFW8, 2015.

  • Hook, S.: Combined ASTER and MODIS Emissivity database over Land (CAMEL) Coefficient Monthly Global 0.05Deg V002, https://doi.org/10.5067/MEASURES/LSTE/CAM5K30CF.002, 2019a.

  • Hook, S.: Combined ASTER and MODIS Emissivity database over Land (CAMEL) Emissivity Monthly Global 0.05Deg V002, https://doi.org/10.5067/MEASURES/LSTE/CAM5K30EM.002, 2019b.

  • Hook, S.: Combined ASTER and MODIS Emissivity database over Land (CAMEL) Uncertainty Monthly Global 0.05Deg V002, https://doi.org/10.5067/MEASURES/LSTE/CAM5K30UC.002, 2019c.

  • Rodgers, C. D.: Inverse methods for atmospheric sounding: theory and practice, World Scientific, Singapore; Hackensack, N.J., 2000.

  • Rosenkranz, P. W.: Retrieval of temperature and moisture profiles from AMSU-A and AMSU-B measurements, 39, 2429–2435, https://doi.org/10.1109/36.964979, 2001.

  • Rosenkranz, P. W.: Cloud liquid-water profile retrieval algorithm and validation, 111, https://doi.org/10.1029/2005JD005832, 2006.

  • Smith, N. and Barnet, C. D.: Uncertainty Characterization and Propagation in the Community Long-Term Infrared Microwave Combined Atmospheric Product System (CLIMCAPS), 11, 1227, https://doi.org/10.3390/rs11101227, 2019.

  • Smith, N. and Barnet, C. D.: CLIMCAPS Observing Capability for Temperature, Moisture and Trace Gases from AIRS/AMSU and CrIS/ATMS, Gases/Remote Sensing/Data Processing and Information Retrieval, https://doi.org/10.5194/amt-2020-71, 2020.

  • Wu, X. and Smith, W. L.: Emissivity of rough sea surface for 8–13 µm: modeling and verification, 36, 2609, https://doi.org/10.1364/AO.36.002609, 1997.