Retrieval of geophysical parameters from moderate resolution imaging spectroradiometer thermal infrared data: evaluation of a two-step physical algorithm

A two-step physical algorithm that simultaneously retrieves geophysical parameters from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements was developed. The retrieved geophysical parameters include atmospheric temperature-humidity profile, surface skin temperature, and two surface emissivities within the shortwave (3-5-mum) and the longwave (8-14.5-mum) regions. The physical retrieval is accomplished in two steps: (i) The Tikhonov regularization method is employed to generate a regularization solution along with an optimum regularization parameter; (ii) the nonlinear Newtonian iteration algorithm is carried out with the regularization solution as a first-guess profile to obtain a final maximum probability solution for geophysical parameters. The algorithm was tested with both simulated and real MODIS Airborne Simulator (MAS) data. Sensitivity studies on simulated MAS data demonstrate that simultaneous retrievals of land and atmospheric parameters improve the accuracy of the retrieved geophysical parameters. Finally, analysis and accuracy of retrievals from real MAS data are discussed.     

Evaluation of the aerosol models for SeaWiFS and MODIS by AERONET data over open oceans

The operational atmospheric correction algorithm for Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer (MODIS) uses the predefined aerosol models to retrieve aerosol optical properties, and their accuracy depends on how well the aerosol models can represent the real aerosol optical properties. In this paper, we developed a method to evaluate the aerosol models (combined with the model selection methodology) by simulating the aerosol retrieval using the Aerosol Robotic Network (AERONET) data. Our method can evaluate the ability of aerosol models themselves, independent of the sensor performance. Two types of aerosol models for SeaWiFS and MODIS operational atmospheric correction algorithms are evaluated over global open oceans, namely the GW1994 models and Ahmad2010 models. The results show that GW1994 models significantly overestimate the aerosol optical thicknesses and underestimate the ?ngstr?m exponent, which is caused by the underestimation of the scattering phase function. However, Ahmad2010 models can significantly reduce the overestimation of the aerosol optical thickness and the underestimation of the ?ngstr?m exponent as a whole, but this improvement depends on the backscattering angle. Ahmad2010 models have a significant improvement in the retrieval of the aerosol optical thickness at a backscattering angle less than 140°. For a backscattering angle larger than 140°, GW1994 models are better at retrieving the aerosol optical thickness than the Ahmad2010 models.     

Retrieval of atmospheric-temperature and water-vapor profiles by use of combined satellite and ground-based infrared spectral-radiance measurements

A nonlinear sounding retrieval algorithm is used to produce vertical-temperature and water-vapor profiles from coincident observations taken by the airborne High-resolution Interferometer Sounder (HIS) and the ground-based Atmospheric Emitted Radiance Interferometer (AERI) during the SUbsonic Contrails and Clouds Effects Special Study (SUCCESS). Also, clear sky Geostationary Operational Environmental Satellite (GOES) and AERI radiance measurements, achieved on a daily real-time basis at the Department of Energy's Oklahoma CART (Cloud and Radiation Testbed) site, are used to demonstrate the current profiling capability by use of simultaneous geostationary satellite and ground-based remote sensing observations under clear-sky conditions. The discrepancy principle, a method to find the proper smoothing parameters from the minimum value between the normalized spectral residual norm and the a priori upper bound, is used to demonstrate the feasibility and effectiveness of on-line simultaneous tuning of the multiple weighting and smoothing parameters from the combined satellite/airborne and ground-based measurements for the temperature and water-vapor retrieval in this nonlinear-retrieval process. An objective method to determine the degrees of freedom (d.f.) of the observation signal is derived. The d.f. of the radiance signal for the combined GOES and AERI measurements is larger than that for either instrument alone; while the d.f. of the observation signal for the combined GOES and AERI measurements is larger than that for either instrument alone and of the combined GOES and AERI measurements. The use of simultaneous clear-sky AERI and GOES data now provides improved vertical temperature and moisture soundings on an hourly basis for use in the Atmospheric Radiation Measurement program [J. Appl. Meteorol. 37, 875 (1998)].     

Multiparameter retrieval of water optical properties from above-water remote-sensing reflectance using the simulated annealing algorithm

Water-leaving radiance, measured just above the ocean surface, contains important information about near-surface or subsurface processes that occur on or below the deep ocean and coastal water. As such, retrieving seawater inherent optical properties (IOPs) is an important step to determining water type, subsurface light field, turbidity, pigment concentration, and sediment loading. However, the retrieval (or inversion) of seawater IOPs from just above water radiance measurements is a multiparameter nonlinear problem that is difficult to solve by conventional optimization methods. The applicability of the simulated annealing algorithm (SA) is explored as a nonlinear global optimizer to solve this multiparameter retrieval problem. The SA algorithm is combined with widely known semianalytical relations for seawater's IOPs to parameter invert these properties from simulated and measured water-leaving reflectance spectra. Furthermore, given the versatility of the SA algorithm, the scheme is extended to retrieve water depth from input reflectance data. Extensive tests and comparisons with in situ and simulated data sets compiled by the International Ocean-Color Coordinating Group are presented. Field data include reflectance spectra acquired with a handheld GER 1500 spectroradiometer and absorption measurements, performed with the AC-9 instrument on waters around Singapore's nearby islands.     

Radiative transfer. I. Atmospheric transmission monitoring with modeling and ground-based multispectral measurements

Existing solar radiative codes such as lowtran allow us to model the radiative properties of the atmosphere and its constituents for standard atmospheric conditions. The increase in urbanization and air pollution has led to changes in the distribution, type, and concentration of the atmospheric constituents, affecting spectral atmospheric transmission and modifying weather and climate. This requires knowledge of the real optical properties of atmospheric transmission. We have developed a model for the radiative properties of atmospheric transmission, with ground-based multispectral measurements of direct solar radiation in the 310-830-nm range. An application of this model to Athens' urban atmosphere is described. The radiative properties of a U.S. Standard Atmosphere are also simulated by use of the lowtran 7 code; simulations and calculations are compared. The total ozone retrieval scheme and the algorithm for retrieving the spectral transmission function and optical thickness, considering multiple scattering, are given. Results for the spectral atmospheric transmission and aerosol and gas transmission functions as well as optical-thickness measurements for the Athens area are also presented as an application of the proposed methodology.     

Simultaneous retrieval of the complex refractive indices of the core and shell of coated aerosol particles from extinction measurements using simulated annealing

Simultaneously retrieving the complex refractive indices of the core and shell of coated aerosol particles given the measured extinction efficiency as a function of particle dimensions (core diameter and coated diameter) is much more difficult than retrieving the complex refractive index of homogeneous aerosol particles. Not only must the minimization be performed over a four-parameter space, making it less efficient, but in addition the absolute value of the difference between the measured extinction and the calculated extinction does not have an easily distinguished global minimum. Rather, there are a number of local minima to which almost all conventional retrieval algorithms converge. In this work, we develop a new (to our knowledge) retrieval algorithm that employs the numerical method known as simulated annealing with an innovative "temperature" schedule. This study is limited only to spherical particles with a concentric shell and to cases in which the diameter of both the core and the coated particle are known. We find that when the top ranking particle sizes according to their information content are combined from separate experiments to make up the particle size distribution, the simulated annealing retrieval algorithm is quite robust and by far superior to a greedy random perturbation approach often used.     

Retrieval of atmospheric profiles from satellite sounder measurements by use of the discrepancy principle

It is known that an infrared or a microwave remote-sensing equation is an integral equation of the first kind. As a result, it is ill-posed, the solution is unstable, and difficulties arise in its retrieval. To make the solution stable, either an a priori error covariance matrix or a smoothing factor gamma is necessary as a constraint. However, if the error covariance matrix is not known or if it is estimated incorrectly, the solution will be suboptimal. The smoothing factor gamma depends greatly on the observations, the observation error, the spectral coverage of channels, and the initial state or the first guess of the atmospheric profile. It is difficult to determine this factor properly during the retrieval procedure, so the factor is usually chosen empirically. We have developed a discrepancy principle (DP) to determine the gamma in an objective way. An approach is formulated for achieving an optimal solution for the atmospheric profile together with the gamma from satellite sounder observations. The DP method was applied to actual Geostationary Operational Environment Satellite (GOES-8) sounder data at the Southern Great Plains Cloud and Radiation Testbed site. Results show that the DP method yields a 21.7% improvement for low-level temperature and a 23.9% improvement for total precipitable water (TPW) retrievals compared with the traditional minimum-information method. The DP method is also compared with the Marquardt-Levenberg algorithm used in current operational GOES data processing. Results of the comparison show significant improvement, 6.5% for TPW and 11% for low-level water-vapor retrievals, in results obtained with the DP method compared with the Marquardt-Levenberg approach.     

Bandpass-resampling effects for the retrieval of surface emissivity

The retrieval of surface emissivity in the 8-14-microm region from remotely sensed thermal imagery requires channel-averaged values of atmospheric transmittance, path radiance, and downwelling sky flux. Band-pass resampling introduces inherent retrieval errors that depend on atmospheric conditions, spectral region, bandwidth, flight altitude, and surface temperature. This simulation study is performed for clear sky conditions and moderate atmospheric water vapor contents. It shows that relative emissivity retrieval errors can reach as much as 3% for broadband sensors (1-2-microm bandwidth) and 0.8% for narrowband instruments (0.15 microm), even for constant surface emissivity. For spectrally varying surface emissivities the relative retrieval error increases for the broadband instrument by approximately 2% in channels with strong emissivity changes of 0.05-0.1. The corresponding retrieval errors for narrowband sensors increase by approximately 3-4%. The channels in the atmospheric window regions with lower transmittance, i.e., 8-8.5 and 12.5-14 microm, are most sensitive to retrieval errors.     

Volume emission rate tomography from a satellite platform

The possibility of retrieving horizontal atmospheric structure from a series of limb images taken aboard a satellite is discussed and a maximum likelihood expectation maximization algorithm is developed. Examples of the retrieval of horizontal structure with this algorithm, for different S/N (signal-to-noise) ratios and different structures, are presented. It is shown that with this algorithm and even in the presence of substantial observational noise, a S/N equal to 10 for a single observation, it is possible to retrieve both horizontal and vertical atmospheric structure.     

利用MODIS影像反演沙漠化地区地表温度的方法研究

针对当前利用MODIS数据反演地表温度过于复杂和缺乏简便易行的实地验证方法的不足,尝试了改进反演方法的研究,提出了反演结果的地面准同步测量验证方法和在沙漠化地区分类计算地表比辐射率的方法,简化了原劈窗算法所需参数的计算过程。利用改进后的方法对新疆沙漠化地区的地表温度进行了反演,并同步进行了该区域在各种天气条件和地表条件下的地表温度反演结果的验证。经分析表明,在天空晴朗无云,无风,植被稀疏,地表类型均一,面积较大,沙质土壤且观测时刻与卫星过境时刻接近的条件下,该方法精度较高。研究表明改进的利用MODIS影像数据反演沙漠化地区地表温度方法和地面验证方法具有一定的可行性,适合于对沙化地区地表温度大范围且快速的遥感监测。     

Atmospheric correction over case 2 waters with an iterative fitting algorithm: relative humidity effects

In algorithms for the atmospheric correction of visible and near-IR satellite observations of the Earth's surface, it is generally assumed that the spectral variation of aerosol optical depth is characterized by an Angstr?m power law or similar dependence. In an iterative fitting algorithm for atmospheric correction of ocean color imagery over case 2 waters, this assumption leads to an inability to retrieve the aerosol type and to the attribution to aerosol spectral variations of spectral effects actually caused by the water contents. An improvement to this algorithm is described in which the spectral variation of optical depth is calculated as a function of aerosol type and relative humidity, and an attempt is made to retrieve the relative humidity in addition to aerosol type. The aerosol is treated as a mixture of aerosol components (e.g., soot), rather than of aerosol types (e.g., urban). We demonstrate the improvement over the previous method by using simulated case 1 and case 2 sea-viewing wide field-of-view sensor data, although the retrieval of relative humidity was not successful.     

Interferometer for ground-based observations of emitted spectral radiance from the troposphere: evaluation and retrieval performance

We evaluate the spectral quality, radiometric noise, and retrieval performance of a Fourier transform infrared spectrometer, which has been developed for recording spectrally resolved observations in a region of the spectrum which is important both for the science of Earth's climate and applications, such as the remote sensing of temperature and atmospheric gas species. This spectral region extends from 100 to 1600 cm(-1) and encompasses the two fundamental, rotation and vibration, absorption bands of water vapor. The instrument is a customized version of a Bomem AERI (Atmospheric Emitted Radiance Interferometer) spectrometer, whose spectral coverage has been extended in the far infrared with the use of uncooled pyroelectric detectors. Retrieval examples for water vapor and temperature profiles are shown, which also allow us to intercompare the retrieval performance of both H(2)O vibration and rotation bands.     

Optical remote sensing of marine constituents in coastal waters: a feasibility study

Optical remote sensing of ocean color is a well-established technique for inferring ocean properties. However, most retrieval algorithms are based on the assumption that the radiance received by satellite instruments is affected only by the phytoplankton pigment concentration and correlated substances. This assumption works well for open ocean water but becomes questionable for coastal waters. To reduce uncertainties associated with this assumption, we developed a new algorithm for the retrieval of marine constituents in a coastal environment. We assumed that ocean color can be adequately described by a three-component model made up of chlorophyll a, suspended matter, and yellow substance. The simultaneous retrieval of these three marine constituents and of the atmospheric aerosol content was accomplished through an inverse-modeling scheme in which the difference between simulated radiances exiting the atmosphere and radiances measured with a satellite sensor was minimized. Simulated radiances were generated with a comprehensive radiative transfer model that is applicable to the coupled atmosphere-ocean system. The method of simulated annealing was used to minimize the difference between measured and simulated radiances. To evaluate the retrieval algorithm, we used simulated (instead of measured) satellite-received radiances that were generated for specified concentrations of aerosols and marine constituents, and we tested the ability of the algorithm to retrieve assumed concentrations. Our results require experimental validation but show that the retrieval of marine constituents in coastal waters is possible.     

Atmospheric correction algorithm for hyperspectral remote sensing of ocean color from space

Existing atmospheric correction algorithms for multichannel remote sensing of ocean color from space were designed for retrieving water-leaving radiances in the visible over clear deep ocean areas and cannot easily be modified for retrievals over turbid coastal waters. We have developed an atmospheric correction algorithm for hyperspectral remote sensing of ocean color with the near-future Coastal Ocean Imaging Spectrometer. The algorithm uses lookup tables generated with a vector radiative transfer code. Aerosol parameters are determined by a spectrum-matching technique that uses channels located at wavelengths longer than 0.86 mum. The aerosol information is extracted back to the visible based on aerosol models during the retrieval of water-leaving radiances. Quite reasonable water-leaving radiances have been obtained when our algorithm was applied to process hyperspectral imaging data acquired with an airborne imaging spectrometer.     

利用AMSR-E被动微波数据反演地表温度的神经网络算法

结合对地观测卫星AQUA多传感器/多分辨率的特点,研究了利用AMSR-E被动微波数据反演地表温度的神经网络算法.MODIS地表温度(LST)产品被作为地表温度实测数据,对应的平均温度被用作对应AMSR-E像元的实际地表温度,从而克服由于AMSR-E像元尺度太大和云的影响而难以获得地表实测数据的难点.反演结果分析表明,利用神经网络能够精确地由AMSR-E数据反演地表温度.当使用5个频率10个通道反演时,反演精度最高,说明使用更多的通道能更好地消除土壤水分、粗糙度、大气和其它因素的影响.相对于MODIS温度产品,用此算法反演的平均误差约低于2K.     

Columnar aerosol properties over oceans by combining surface and aircraft measurements: simulations

We report an algorithm that can be used to invert the radiance exiting the top and bottom of the atmosphere to yield the columnar optical properties of atmospheric aerosol under clear sky conditions over the oceans. The method is an augmentation of a similar algorithm presented by Wang and Gordon [Appl. Opt. 32, 4598 (1993)] that used only sky radiance, and therefore was incapable of retrieving the aerosol phase function at the large scattering angles that are of critical importance in remote sensing of oceanic and atmospheric properties with satellites. Well-known aerosol models were combined with radiative transfer theory to simulate pseudodata for testing of the algorithm. The tests suggest that it should be possible to retrieve the aerosol phase function and the aerosol single-scattering albedo accurately over the visible spectrum at aerosol optical thicknesses as large as 2.0. The algorithm is capable of retrievals with such large optical thicknesses because all significant orders of multiple scattering are included. We believe that combining an algorithm of this type with surface-based and high-altitude aircraft-based radiance measurements could be useful for studying aerosol columnar optical properties over oceans and large lakes. The use of the retrieval method is possible over the ocean because, unlike the land surface, the albedo of the ocean is low and spatially uniform.     

Atmospheric correction for NO2 absorption in retrieving water-leaving reflectances from the SeaWiFS and MODIS measurements

The absorption by atmospheric nitrogen dioxide (NO2) gas in the visible has been traditionally neglected in the retrieval of oceanic parameters from satellite measurements. Recent measurements of NO2 from spaceborne sensors show that over the Eastern United States the NO2 column amount often exceeds 1 Dobson Unit (approximately 2.69x10(16) molecules/cm2). Our radiative transfer sensitivity calculations show that under high NO2 conditions (approximately 1x10(16) molecules/cm2) the error in top-of-atmosphere (TOA) reflectance in the blue channels of the sea-viewing wide field-of-view sensor (SeaWiFS) and moderate-resolution imaging spectroradiometer (MODIS) sensors is approximately 1%. This translates into approximately 10% error in water-leaving radiance for clear waters and to higher values (>20%) in the coastal areas. We have developed an atmospheric-correction algorithm that allows an accurate retrieval of normalized water-leaving radiances (nLws) in the presence of NO2 in the atmosphere. The application of the algorithm to 52 MODIS scenes over the Chesapeake Bay area show a decrease in the frequency of negative nLw estimates in the 412 nm band and an increase in the value of nLws in the same band. For the particular scene reported in this paper, the mean value of nLws in the 412 nm band increased by 17%, which is significant, because for the MODIS sensor the error in nLws attributable to the digitization error in the observed TOA reflectance over case 2 waters is approximately 2.5%.     

Retrieval of CO from nadir remote-sensing measurements in the infrared by use of four different inversion algorithms

Four inversion schemes based on various retrieval approaches (digital gas correlation, nonlinear least squares, global fit adjustment, and neural networks) developed to retrieve CO from nadir radiances measured by such downward-looking satelliteborne instruments as the Measurement of Pollution in the Troposphere (MOPITT), the Tropospheric Emission Spectrometer (TES), and the Infrared Atmospheric Sounding Interferometer (IASI) instruments were compared both for simulated cases and for atmospheric spectra recorded by the Interferometric Monitor for Greenhouse Gases (IMG). The sensitivity of the retrieved CO total column amount to properties that may affect the inversion accuracy (noise, ancillary temperature profile, and water-vapor content) was investigated. The CO column amounts for the simulated radiance spectra agreed within 4%, whereas larger discrepancies were obtained when atmospheric spectra recorded by the IMG instrument were analyzed. The assumed vertical temperature profile is shown to be a critical parameter for accurate CO retrieval. The instrument's line shape was also identified as a possible cause of disagreement among the result provided by the groups of scientist who are participating in this study.     

SZ-4飞船辐射模态沙漠地区数据模拟验证

对SZ-4飞船辐射模态获取的亮温数据,进行了模拟研究。辐射模态的工作频率范围为6．6～37GHz,对应的地面像元尺寸为6～32km。以新疆塔中地区为研究对象,通过分析被动微波遥感中裸露地表的辐射特性,基于裸露地表的辐射模型（AIEM模型）和地表的实测属性数据来模拟该地区的辐射亮温值,并将其与SZ-4辐射模态和AMSR-E实际的亮温观测值进行了比较,通过对正演结果的比较分析检验了SZ-4亮温数据的质量及其在沙漠地区反演地表参数的潜力。     

GMTR: two-dimensional geo-fit multitarget retrieval model for michelson interferometer for passive atmospheric sounding/environmental satellite observations

We present a new retrieval model designed to analyze the observations of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), which is on board the ENVironmental SATellite (ENVISAT). The new geo-fit multitarget retrieval model (GMTR) implements the geo-fit two-dimensional inversion for the simultaneous retrieval of several targets including a set of atmospheric constituents that are not considered by the ground processor of the MIPAS experiment. We describe the innovative solutions adopted in the inversion algorithm and the main functionalities of the corresponding computer code. The performance of GMTR is compared with that of the MIPAS ground processor in terms of accuracy of the retrieval products. Furthermore, we show the capability of GMTR to resolve the horizontal structures of the atmosphere. The new retrieval model is implemented in an optimized computer code that is distributed by the European Space Agency as "open source" in a package that includes a full set of auxiliary data for the retrieval of 28 atmospheric targets.