Engineering topics
Publications and source records attributed to Liu, Xu.
Different methods for retrieving atmospheric profiles in the presence of clouds from hyperspectral satellite remote sensing data will be described. We will present results from the JPSS cloud-clearing algorithm and NASA Langley cloud retrieval algorithm.
Ultraspectral resolution infrared (IR) radiances obtained from nadir observations provide information about the atmosphere, surface, aerosols, and clouds. Surface spectral emissivity (SSE) and surface skin temperature from current and future operational satellites can and will reveal critical information about the Earth s ecosystem and land-surface-type properties, which might be utilized as a means of long-term monitoring of the Earth s environment and global climate change. In this study, fast radiative transfer models applied to the atmosphere under all weather conditions are used for atmospheric profile and surface or cloud parameter retrieval from ultraspectral and/or hyperspectral spaceborne IR soundings. An inversion scheme, dealing with cloudy as well as cloud-free radiances observed with ultraspectral IR sounders, has been developed to simultaneously retrieve atmospheric thermodynamic and surface or cloud microphysical parameters. This inversion scheme has been applied to the Infrared Atmospheric Sounding Interferometer (IASI). Rapidly produced SSE is initially evaluated through quality control checks on the retrievals of other impacted surface and atmospheric parameters. Initial validation of retrieved emissivity spectra is conducted with Namib and Kalahari desert laboratory measurements. Seasonal products of global land SSE and surface skin temperature retrieved with IASI are presented to demonstrate seasonal variation of SSE.
Advanced satellite sensors are tasked with improving global observations of the Earth's atmosphere, clouds, and surface to enable enhancements in weather prediction, climate monitoring capability, and environmental change detection. Achieving such an improvement in geophysical information inferred from these observations requires optimal usage of data from current systems as well as instrument system enhancements for future sensors. This presentation addresses results of tradeoff studies evaluating the impact of spectral resolution, spectral coverage, instrument noise, and a priori knowledge on remote sensing system information content, with a specific emphasis on thermodynamic state and trace species information obtainable from advanced atmospheric sounders. Particular attention will be devoted toward information achievable from the Atmospheric InfraRed Sounder (AIRS) on the NASA EOS Aqua satellite in orbit since 2002, the Infrared Atmospheric Sounding Interferometer (IASI) aboard MetOp-A since 2006, and the Cross-track Infrared Sounder (CrIS) instrument to fly aboard the NPP and JPSS series of satellites expected to begin in late 2011. While all of these systems cover nearly the same infrared spectral extent, they have very different number of channels, instrument line shapes, coverage continuity, and instrument noise. AIRS is a grating spectrometer having 2378 discrete spectral channels ranging from about 0.4 to 2.2/cm resolution; IASI is a Michelson interferometer with 8461 uniformly-spaced spectral channels of 0.5/cm (apodized) resolution; and CrIS is a Michelson interferometer having 1305 spectral channels of 0.625, 1.250, and 2.50/cm (unapodized) spectral resolution, respectively, over its three continuous but non-overlapping bands. Results of tradeoff studies showing information content sensitivity to assumed measurement system characteristics will be presented.
Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.
General-purpose computing on graphics processing units (GPGPU) is a recent technique that allows the parallel graphics processing unit (GPU) to accelerate calculations performed sequentially by the central processing unit (CPU). To introduce GPGPU to radiative transfer, the Gauss-Seidel solution of the well-known expressions for 1-D and 3-D homogeneous, isotropic media is selected as a test case. Different algorithms are introduced to balance memory and GPU-CPU communication, critical aspects of GPGPU. Results show that speed-ups of one to two orders of magnitude are obtained when compared to sequential solutions. The underlying value of GPGPU is its potential extension in radiative solvers (e.g., Monte Carlo, discrete ordinates) at a minimal learning curve.
Infrared sounders, such as the Atmospheric Infrared Sounder (AIRS), the Infrared Atmospheric Sounding Interferometer (IASI), and the Cross-track Infrared sounder (CrIS), have a cloud-impenetrable disadvantage in observing the atmosphere and surface under opaque cloudy conditions. However, recent studies indicate that hyperspectral, infrared sounders have the ability to detect cloud effective-optical and microphysical properties and to penetrate optically thin clouds in observing the atmosphere and surface to a certain degree. We have developed a retrieval scheme dealing with atmospheric conditions with cloud presence. This scheme can be used to analyze the retrieval accuracy of atmospheric and surface parameters under clear and cloudy conditions. In this paper, we present the surface emissivity results derived from IASI global measurements under both clear and cloudy conditions. The accuracy of surface emissivity derived under cloudy conditions is statistically estimated in comparison with those derived under clear sky conditions. The retrieval error caused by the clouds is shown as a function of cloud optical depth, which helps us to understand how well infrared sounders can observe the atmosphere and surface through clouds.
Advanced satellite sensors are tasked with improving global measurements of the Earth s atmosphere, clouds, and surface to enable enhancements in weather prediction, climate monitoring capability, and environmental change detection. Achieving such measurement improvements requires instrument system advancements. This presentation focuses on the impact of spectral resolution and coverage changes on remote sensing system information content, with a specific emphasis on thermodynamic state and trace species variables obtainable from advanced atmospheric sounders such as the Infrared Atmospheric Sounding Interferometer (IASI) and Cross-track Infrared Sounder (CrIS) systems on the MetOp and NPP/NPOESS series of satellites. Key words: remote sensing, advanced sounders, information content, IASI, CrIS
A principal-component based radiative transfer model (PCRTM) is developed for simulating the infrared spectral radiance at the top of the atmosphere (TOA). The PCRTM approach removes the redundancy in radiative transfer calculation in high resolution infrared spectra, and saves significant amount of computational time with great accuracy. In PCRTM, both ice and water clouds are treated as effective transmissivity and reflectivity stored in a pre-calculated lookup tables. These quantities are calculated using cloud single scattering properties such as cloud optical depth, cloud particle size, cloud phase, etc. The cloud can be inserted into any pressure layer in the PCRTM model (up to 100 layers). The effective temperature of each cloud layer is treated as a function of its optical depth. To test the accuracy of this approximation, the results are compared with the more rigorous DISORT model, which treats cloud as a plane parallel layer. The root-mean-square error of PCRTM, with respect to DISORT results, is generally less than 0.5 K in brightness temperature. However, the CPU time required by PCRTM was approximately two orders of magnitude less than that required by DISORT.
As a part of the National Polar-orbiting Operational Environmental Satellite System (NPOESS) and the NPOESS Preparatory Project (NPP), the instruments Cross-track Infrared Sounder (CrIS) and Advanced Technology Microwave Sounder (ATMS) make up the Cross-track Infrared and Microwave Sounder Suite (CrIMSS). CrIMSS will primarily provide global temperature, moisture, and pressure profiles and calibrated radiances [1]. In preparation for the NPOESS/NPP launch, porting and testing of the CrIMSS Environmental Data Record (EDR) algorithms need to be performed.
Surface spectral emissivity derived from current and future satellites can and will reveal critical information about the Earth s ecosystem and land surface type properties, which can be utilized as a means of long-term monitoring of global environment and climate change. Hyperspectrally-resolved surface emissivities are derived with an algorithm utilizes a combined fast radiative transfer model (RTM) with a molecular RTM and a cloud RTM accounting for both atmospheric absorption and cloud absorption/scattering. Clouds are automatically detected and cloud microphysical parameters are retrieved; and emissivity is retrieved under clear and optically thin cloud conditions. This technique separates surface emissivity from skin temperature by representing the emissivity spectrum with eigenvectors derived from a laboratory measured emissivity database; in other words, using the constraint as a means for the emissivity to vary smoothly across atmospheric absorption lines. Here we present the emissivity derived under optically thin clouds in comparison with that under clear conditions.
The current and next generation satellite infrared sensors such as AIRS, IASI, CrIS, and CLARREO all have one thing in common: large number of spectral channels. In order to fully utilize the information content of these sensors, a large number of Radiative Transfer (RT) calculations through the inhomogeneous atmosphere are needed. It is also computationally intensive to invert atmospheric and surface properties using all the channels. Usually, only subsets of channels are used to perform physical inversions for atmospheric profiles. We will explore ways to speed up RT calculations and inversions and test the algorithms using data from the Infrared Atmospheric Sounding Interferometer (IASI) satellite instrument. We will describe a retrieval algorithm using a Principal Component-based Radiative Transfer Model (PCRTM) for generating atmospheric temperature/moisture/ozone profiles and surface properties. The retrieval algorithm performs both the radiative transfer calculations and inversions in the Principal Component (PC) domain. The inversion algorithm is based on a non-linear Levenberg- Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the ultraspectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels. We will also use data collected during the Joint Airborne IASI Validation Experiment (JAIVEx) field campaign to validate the algorithm and IASI retrievals.
Surface and atmospheric thermodynamic parameters retrieved with advanced ultraspectral remote sensors aboard Earth observing satellites are critical to general atmospheric and Earth science research, climate monitoring, and weather prediction. Ultraspectral resolution infrared radiance obtained from nadir observations provide atmospheric, surface, and cloud information. Presented here is the global surface IR emissivity retrieved from Infrared Atmospheric Sounding Interferometer (IASI) measurements under "clear-sky" conditions. Fast radiative transfer models, applied to the cloud-free (or clouded) atmosphere, are used for atmospheric profile and surface parameter (or cloud parameter) retrieval. The inversion scheme, dealing with cloudy as well as cloud-free radiances observed with ultraspectral infrared sounders, has been developed to simultaneously retrieve atmospheric thermodynamic and surface (or cloud microphysical) parameters. Rapidly produced surface emissivity is initially evaluated through quality control checks on the retrievals of other impacted atmospheric and surface parameters. Surface emissivity and surface skin temperature from the current and future operational satellites can and will reveal critical information on the Earth s ecosystem and land surface type properties, which can be utilized as part of long-term monitoring for the Earth s environment and global climate change.
A novel radiative transfer model and a physical inversion algorithm based on principal component analysis will be presented. Instead of dealing with channel radiances, the new approach fits principal component scores of these quantities. Compared to channel-based radiative transfer models, the new approach compresses radiances into a much smaller dimension making both forward modeling and inversion algorithm more efficient.
Ultraspectral resolution infrared spectral radiance obtained from near nadir observations provide atmospheric, surface, and cloud property information. The intent of the measurement of tropospheric thermodynamic state and trace abundances is the initialization of climate models and the monitoring of air quality. The NPOESS Airborne Sounder Testbed-Interferometer (NAST-I), designed to support the development of future satellite temperature and moisture sounders, aboard high altitude aircraft has been collecting data throughout many field campaigns. An advanced retrieval algorithm developed with NAST-I is now applied to satellite data collected with the Atmospheric InfraRed Sounder (AIRS) on the Aqua satellite launched on 4 May 2002 and the Infrared Atmospheric Sounding Interferometer (IASI) on the MetOp satellite launched on October 19, 2006. These instruments possess an ultra-spectral resolution, for example, both IASI and NAST-I have ~0.25 cm-1 and a spectral coverage from 645 to 2760 cm-1. The retrieval algorithm with a fast radiative transfer model, including cloud effects, is used for atmospheric profile and cloud parameter retrieval. The physical inversion scheme has been developed, dealing with cloudy as well as cloud-free radiance observed with ultraspectral infrared sounders, to simultaneously retrieve surface, atmospheric thermodynamic, and cloud microphysical parameters. A fast radiative transfer model, which applies to the clouded atmosphere, is used for atmospheric profile and cloud parameter retrieval. A one-dimensional (1-d) variational multi-variable inversion solution is used to improve an iterative background state defined by an eigenvector-regression-retrieval. The solution is iterated in order to account for non-linearity in the 1-d variational solution. It is shown that relatively accurate temperature and moisture retrievals can be achieved below optically thin clouds. For optically thick clouds, accurate temperature and moisture profiles down to cloud top level are obtained. For both optically thin and thick cloud situations, the cloud top height can be retrieved with relatively high accuracy (i.e., error less than 1 km). Retrievals of atmospheric soundings, surface properties, and cloud microphysical properties with the AIRS and IASI observations are obtained and presented. These retrievals are further inter-compared with those obtained from airborne FTS system, such as the NPOESS Airborne Sounder Testbed? Interferometer (NAST I), dedicated dropsondes, radiosondes, and ground based Raman Lidar. The capabilities of satellite ultra-spectral sounder such as the AIRS and IASI are investigated. These advanced satellite ultraspectral infrared instruments are now playing an important role in satellite meteorological observation for numerical weather prediction.
Atmospheric-thermodynamic parameters and surface properties are basic meteorological parameters for weather forecasting. A physical geophysical parameter retrieval scheme dealing with cloudy and cloud-free radiance observed with satellite ultraspectral infrared sounders has been developed and applied to the Infrared Atmospheric Sounding Interferometer (IASI) and the Atmospheric InfraRed Sounder (AIRS). The retrieved parameters presented herein are from radiance data gathered during the Joint Airborne IASI Validation Experiment (JAIVEx). JAIVEx provided intensive aircraft observations obtained from airborne Fourier Transform Spectrometer (FTS) systems, in-situ measurements, and dedicated dropsonde and radiosonde measurements for the validation of the IASI products. Here, IASI atmospheric profile retrievals are compared with those obtained from dedicated dropsondes, radiosondes, and the airborne FTS system. The IASI examples presented here demonstrate the ability to retrieve fine-scale horizontal features with high vertical resolution from satellite ultraspectral sounder radiance spectra.
The Infrared Atmospheric Sounding Interferometer (IASI) is a hyperspectral sensor with 8461 spectral channels and a nominal spectral resolution of 0.25 cm(sup -1). It is computationally intensive to perform radiative transfer calculations and inversions using all these channels. We will present a Principal Component-based Radiative Transfer Model (PCRTM) and a retrieval algorithm which perform all the necessary calculations in EOF domain. Since the EOFs are orthogonal to each other, only about 100 principal components are needed to represent the information content of the 8461 channels. The PCRTM provides the EOF coefficients and associated derivatives with respect to atmospheric and surface parameters needed by the inversion algorithm. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as constraints. The retrieved parameters include atmospheric temperature, moisture and ozone profiles, cloud parameters, surface skin temperature, and surface emissivities. To make the retrieval system even more compact and stable. The atmospheric vertical profiles are compressed into the EOF space as well. The surface emissivities are also compressed into EOF space.
Since modern ultra-spectral remote sensors have thousands of channels, it is difficult to include all of them in a 1D-var retrieval system. We will describe a physical inversion algorithm, which includes all available channels for the atmospheric temperature, moisture, cloud, and surface parameter retrievals. Both the forward model and the inversion algorithm compress the channel radiances into super channels. These super channels are obtained by projecting the radiance spectra onto a set of pre-calculated eigenvectors. The forward model provides both super channel properties and jacobian in EOF space directly. For ultra-spectral sensors such as Infrared Atmospheric Sounding Interferometer (IASI) and the NPOESS Airborne Sounder Testbed Interferometer (NAST), a compression ratio of more than 80 can be achieved, leading to a significant reduction in computations involved in an inversion process. Results will be shown applying the algorithm to real IASI and NAST data.
A state-of-the-art IR-only retrieval algorithm has been developed with an all-season-global EOF Physical Regression and followed by 1-D Var. Physical Iterative Retrieval for IASI, AIRS, and NAST-I. The benefits of this retrieval are to produce atmospheric structure with a single FOV horizontal resolution (approx. 15 km for IASI and AIRS), accurate profiles above the cloud (at least) or down to the surface, surface parameters, and/or cloud microphysical parameters. Initial case study and validation indicates that surface, cloud, and atmospheric structure (include TBL) are well captured by IASI and AIRS measurements. Coincident dropsondes during the IASI and AIRS overpasses are used to validate atmospheric conditions, and accurate retrievals are obtained with an expected vertical resolution. JAIVEx has provided the data needed to validate the retrieval algorithm and its products which allows us to assess the instrument ability and/or performance. Retrievals with global coverage are under investigation for detailed retrieval assessment. It is greatly desired that these products be used for testing the impact on Atmospheric Data Assimilation and/or Numerical Weather Prediction.