Recorded presentation5.2
Variational retrievals of atmospheric profiles including clouds in principal component space using the Havemann-Taylor Fast Radiative Transfer Code (HT-FRTC)
Stephan Havemann, Met Office, Exeter, United Kingdom; and J. C. Thelen, A. Baran, and J. P. Taylor
The Havemann-Taylor Fast Radiative Transfer Code allows fast and exact calculations and is ideally suited for the simulation of hyperspectral instruments with hundreds or thousands of channels. The code works in any part of the spectrum and has so far been used in the shortwave and in the infrared. During a training phase, the code is trained on a wide variety of on the order of a 1000 different atmospheric and surface conditions. The line-by-line calculations carried out during this training phase are based on the recent spectroscopy of LBLRTM 11.7. From the radiance spectra obtained during the training, up to 100 Principal Components are derived, which contain the information of all the partly correlated channel radiances in a much more compact form. With a k-means clustering algorithm, a few hundred frequencies are selected. Monochromatic radiances obtained from line-by-line calculations at these centroid frequencies are then used as the predictors for the Principal Components (PC) Scores. The HT-FRTC can simulate a spectrum for any atmosphere/surface within a fraction of second, because it does only monochromatic calculations and it saves computer memory since it works in the low-dimensional PC space.
For the simulation of scattering by clouds and aerosols and Rayleigh scattering in the short-wave, a spherical harmonics line-by-line code has been integrated into the HT-FRTC, which is very similar to the Edwards-Slingo radiation code. Scattering calculations with the HT-FRTC for a complete spectrum take on the order of a second.
The HT-FRTC has been incorporated into a one-dimensional variational retrieval system that also works solely in PC space, which keeps the dimensions of the matrices involved small. The solution of the full non-linear problem is achieved with an iterative Levenberg-Marquardt minimization procedure. The variational retrieval system creates a lot of diagnostics in addition to the solution, like the averaging kernel, the Jacobians, the gain matrix and the degrees of signal as a measure of information gained. The state vector includes the vertical profiles of atmospheric temperature, water vapour and ozone, as well as the surface temperature and surface emissivity (which is represented as principal components). For a scattering atmosphere, cloud parameters and aerosol parameters can be added to the state vector. All state vector elements can be retrieved simultaneously or alternatively some elements of the state vector can be kept fixed and not retrieved. The cloud part of the state vector for cirrus cloud includes cloud top pressure, water content, cloud fraction and cloud optical depth. For water cloud there is also an effective droplet size.
For cirrus cloud a new description of the optical properties has been applied. The single-scattering properties (i.e. volume extinction coefficient, single-scattering albedo and the moments of the phase function), rather than being linked to the ice crystal size distribution function (PSD) through the concept of an effective size are instead described in terms of a parametrized PSD. The parametrized PSD is based on in situ measurements of midlatitude and tropical PSD's and is generated from just two variables, the ice water content (IWC) and the cloud temperature.
Case studies of measurements with the Atmospheric Research Interferometer Evaluation System (ARIES) on board the Met Office BAe-146 research aircraft coincident with overpasses of the Infrared Atmospheric Sounding Interferometer (IASI) above ice/water cloud fields will be used to retrieve atmospheric and cloud parameters from ARIES and IASI. Retrieved cloud parameters will be compared to cloud microphysical aircraft measurements. The effect of cloud on the retrieval of the temperature and humidity profiles will be investigated. Also, it will be studied how much information on the vertical profile of IWC can be obtained from the measurements.
Session 5, Remote Sensing I: Passive
Tuesday, 29 June 2010, 8:30 AM-10:00 AM, Pacific Northwest Ballroom
Previous paper Next paper