SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Radiative transfer between two infinite parallel metallic surfaces separated by nonconducting thick film ideal dielectric based on electromagnetic wave theory
TOA (top-of-atmosphere) radiance from high-spatial-resolution satellite imagery systems is important for a wide variety of research and applications. Many research initiatives require data with absolute radiometric accuracy better than a few percent. The conversion of satellite digital numbers to radiance depends on accurate radiometric calibration. A common method for determining and validating radiometric calibrations is to rely upon vicarious calibration approaches. Historically, vicarious calibration methods use radiative transfer codes with ground-based atmosphere and surface reflectance or radiance inputs for estimating TOA radiance values. These TOA radiance values are compared against the satellite digital numbers to determine the radiometric calibration. However, the radiative transfer codes used depend on many assumptions about the aerosol properties and the atmospheric point spread function. A measurement-based atmospheric radiance estimation approach for high-spatial-resolution, multispectral, visible/near-infrared sensors is presented that eliminates the use of radiative transfer codes and many of the underlying assumptions. A comparison between the radiative transfer and non-radiative transfer approaches is made.
Scattering calculations using the discrete dipole approximation and vector radiative transfer calculations were performed to model multiparameter radar return and passive microwave emission for a simple model of a winter storm. The issue of dendrite riming was addressed by computing scattering properties of thin ice disks with varying bulk density. It was shown that C-band multiparameter radar contains information about particle density and the number concentration of the ice particles. The radiative transfer modeling indicated that polarized multifrequency passive microwave emission may be used to infer some properties of ice hydrometers. Detailed radar modeling and vector radiative transfer modeling is in progress to enhance the understanding of simultaneous radar and radiometer measurements, as in the case of the proposed TRMM field program. A one-dimensional cloud model will be used to simulate the storm structure in detail and study the microphysics, such as size and density. Multifrequency polarized radiometer measurements from the SSMI satellite instrument will be analyzed in relation to dual-frequency and dual-polarization radar measurements.
Measurements from the FIRE 1991 cirrus cloud field experiment in the central United States are presented and analyzed. The first part focuses on cirrus microphysical properties. Aircraft 2D-probe in situ data at different cloud altitudes were evaluated for cirrus cases on four different days. Also presented are simultaneous data samples from balloonborne videosondes. Only these balloonsondes could detect the smaller crystals. Their data suggest (at least for midlatitude altitudes below 10 km) that ice crystals smaller than 15 microns in size are rare and that small ice crystals not detected by 2D-probe measurements are radiatively of minor importance, as overlooked 2D-probe crystals account for about 10% of the total extinction. The second part focuses on the link between cirrus cloud properties and radiation. With cloud macrophysical properties from surface remote sensing added to the microphysical data and additional radiation measurements at the surface, testbeds for radiative transfer models were created. To focus on scattering processes, model evaluations were limited to the solar radiative transfer by comparing calculated and measured transmissions of sunlight at the surface. Comparisons under cloud-free conditions already reveal a model bias of about +45 W/sq m for the hemispheric solar downward broadband flux. This discrepancy, which is (at least in part) difficult to explain, has to be accounted for in comparisons involving clouds. Comparisons under cirrus cloud conditions identify as the major obstacle in cirrus solar radiative transfer modeling the inability of one-dimensional radiative transfer models to account for horizontal inhomogeneities. The successful incorporation of multidimensional radiative transfer effects will depend not only on better models but critically on the ability to measure and to define characteristic inhomogeneity scales of cloud fields. The relative minor error related to the microphysical treatment is in part a reflection of the improved understanding on solar scattering on ice crystals over the last decade and of the available wealth on ice-crystal size and shape data for this study. In absence of this information, uncertainties from microphysical cirrus model assumptions will remain high.
This presentation will discuss the development of three-dimensional (3D) radiative transfer in cloudy atmospheres for the last two decades. The main emphasis will be put on the remote sensing problems. Observational evidence of 3D radiative effects that influence retrievals of cloud and aerosol properties will be shown using examples from MODIS and ASTER. The international Intercomparison of 3D Radiation Codes (13RC) will be also discussed. Finally, I will briefly review the results and advanced radiative transfer tools recently developed by our team at Goddard Space Flight Center (Cahalan, Chiu, Marshak, Oreopoulos, Varnai, Wen and Wiscombe). Note that the field of 3D radiative transfer is now sufficiently mature and its first book has just been published by Springer, edited by Marshak and Davis.
The second step in development of MODEL III is summarized. It combines the four radiative transfer equations of the first step with the equations for a geostrophic and hydrostatic atmosphere. This step is intended to bring radiance into a three dimensional balance with wind, height, and temperature. The use of the geostrophic approximation in place of the full set of primitive equations allows for an easier evaluation of how the inclusion of the radiative transfer equation increases the complexity of the variational equations. Seven different variational formulations were developed for geostrophic, hydrostatic, and radiative transfer equations. The first derivation was too complex to yield solutions that were physically meaningful. For the remaining six derivations, the variational method gave the same physical interpretation (the observed brightness temperatures could provide no meaningful input to a geostrophic, hydrostatic balance) at least through the problem solving methodology used in these studies. The variational method is presented and the Euler-Lagrange equations rederived for the geostrophic, hydrostatic, and radiative transfer equations.
The direct method for the solution of the spherical harmonics approximation to the equation of transfer of radiation is applied to the cases of (1) scattering of the solar radiation in the atmosphere with the Lambertian boundary and (2) thermal radiation transfer.
Infrared radiative transfer through a regular array of cuboidal clouds is studied and the interaction of the sides of the clouds with each other and the ground is considered. The theory is developed for black clouds and is extended to scattering clouds using a variable azimuth two-stream approximation. It is shown that geometrical considerations often dominate over the microphysical aspects of radiative transfer through the clouds. For example, the difference in simulated 10 micron brightness temperature between black isothermal cubic clouds and cubic clouds of optical depth 10, is less than 2 deg for zenith angles less than 50 deg for all cloud fractions when viewed parallel to the array. The results show that serious errors are made in flux and cooling rate computations if broken clouds are modeled as planiform. Radiances computed by the usual practice of area-weighting cloudy and clear sky radiances are in error by 2 to 8 K in brightness temperature for cubic clouds over a wide range of cloud fractions and zenith angles. It is also shown that the lapse rate does not markedly affect the exiting radiances for cuboidal clouds of unit aspect ratio and optical depth 10.
Inverse problems in radiative transfer applied to remote measurements of IR radiation from planetary atmospheres, noting temperature and water vapor inversions
We perform 3D radiative transfer calculations in non-local thermodynamic equilibrium (NLTE) in the simple two-level atom approximation on the Mas-Par MP-1, which contains 8192 processors and is a single instruction multiple data (SIMD) machine, an example of the new generation of massively parallel computers. On such a machine, all processors execute the same command at a given time, but on different data. To make radiative transfer calculations efficient, we must re-consider the numerical methods and storage of data. To solve the transfer equation, we adopt the short characteristic method and examine different acceleration methods to obtain the source function. We use the ALI method and test local and non-local operators. Furthermore, we compare the Ng and the orthomin methods of acceleration. We also investigate the use of multi-grid methods to get fast solutions for the NLTE case. In order to test these numerical methods, we apply them to two problems with and without periodic boundary conditions.
The Joint Center for Satellite Data Assimilation (JCSDA) Community Radiative Transfer Model (CRTM) is a fast, 1-D radiative transfer model used in numerical weather prediction, calibration/validation, etc. across multiple federal agencies and universities. The key benefit of the CRTM is that it is a satellite simulator. It provides a highly accurate representation of satellite radiances by using the specific sensor response functions convolved with a line-by-line radiative transfer model (LBLRTM). CRTM covers the spectral ranges consistent with all present operational and most research satellites, from visible to microwave. The capability to simulate ultraviolet radiances and support space-based radar sensors is being added over the next two years in CRTM Version 3.0. In addition to simulated radiances, the CRTM also provides Jacobian outputs needed to interpret satellite observations for numerical weather prediction. The Jacobian estimates how changes in geophysical parameters affect simulated measurements from satellite sensors. Using the Jacobian in modeling and weather prediction improves the accuracy and efficiency of data analysis, leading to better weather predictions. The CRTM model's success and growth depend on community contributions and evaluation. To facilitate this, we have made the CRTM highly accessible through modular programming, clear documentation and tutorials, public domain licensing, unfettered public access via Github, and a clear path to operational implementation for innovative research. We encourage and welcome contributions from the community to help us continue to improve the CRTM.
Modern infrared satellite sensors such as Atmospheric Infrared Sounder (AIRS), Cosmic Ray Isotope Spectrometer (CrIS), Thermal Emission Spectrometer (TES), Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) and Infrared Atmospheric Sounding Interferometer (IASI) are capable of providing high spatial and spectral resolution infrared spectra. To fully exploit the vast amount of spectral information from these instruments, super fast radiative transfer models are needed. This paper presents a novel radiative transfer model based on principal component analysis. Instead of predicting channel radiance or transmittance spectra directly, the Principal Component-based Radiative Transfer Model (PCRTM) predicts the Principal Component (PC) scores of these quantities. This prediction ability leads to significant savings in computational time. The parameterization of the PCRTM model is derived from properties of PC scores and instrument line shape functions. The PCRTM is very accurate and flexible. Due to its high speed and compressed spectral information format, it has great potential for super fast one-dimensional physical retrievals and for Numerical Weather Prediction (NWP) large volume radiance data assimilation applications. The model has been successfully developed for the National Polar-orbiting Operational Environmental Satellite System Airborne Sounder Testbed - Interferometer (NAST-I) and AIRS instruments. The PCRTM model performs monochromatic radiative transfer calculations and is able to include multiple scattering calculations to account for clouds and aerosols.
Moderate Resolution Imaging Spectroradiometer (MODIS) on Aqua has been in operation providing continuous global observations for science research and applications since 2002. The long-term stability of thermal emissive bands (TEBs) of Aqua MODIS was monitored through inter-comparisons with measurements by hyperspectral infrared sensors, such as AIRS on Aqua and IASI on MetOp satellite, or through long term vicarious monitoring over cold targets, such as Dome- C and deep convective clouds (DCC). In this paper, a radiative transfer modeling-based simulation model using Community Radiative Transfer Model (CRTM) is developed to perform the long-term monitoring of the stability of TEBs of Aqua MODIS. CRTM is a fast-radiative transfer model for calculations of radiances for satellite infrared radiometers and is able to output infrared radiance and brightness temperature at spectral bands of MODIS. Long term European Centre for Medium-Range Weather Forecasts (ECMWF) global atmospheric reanalysis data, such as temperature and humidity profiles, are used as inputs to the CRTM simulation. By confining the area of interest to be over low to middle latitude ocean, the long-term stabilities of selected Aqua MODIS TEBs are monitored through Observation–Background (O-B) brightness temperature (BT) bias between MODIS measurements and BT retrieval from CRTM simulation. The consistency and relative stability between Aqua MODIS and ECMWF reanalysis data for surface channels of MODIS are evaluated. In addition, the radiative transfer modeling with CRTM enables us to evaluate the impacts of long term variation of global CO2 distribution on the O-B BT biases for CO2 channels of Aqua MODIS through comparison of simulations with constant or long term variable CO2 as inputs. The O-B analysis with RTM show that Aqua MODIS surface channels are all radiometrically stable with yearly BT bias drift less than 0.004K/year for B20, B22, B23 and ~0.01K/year for B31 and B32. The CO2 absorption channels of Aqua MODIS, e.g. B33-B36, are stable with BT bias drift < 0.005K/year.
An important contribution to coral reef research is to improve spectral distinction between various health states of coral species in areas subject to harmful anthropogenic activity and climate change. New insights into radiative transfer properties of corals under healthy and stressed conditions can advance understandings of ecological processes on reefs and allow better assessments of the impacts of large-scale bleaching and disease events, Our objective was to examine the spectral and spatial properties of hyperspectral sensors that may be used to remotely sense changes in reef community health. We compare in situ reef environment spectra (healthy coral, stressed coral, dead coral, algae, and sand) with airborne hyperspectral data to identify important spectral characteristics and indices. Additionally, spectral measurements over a range of water depths, relief, and bottom types are compared to help quantify bottom-water column influences. In situ spectra were collected in July and August 2002 at the Long Rock site in the Andros Island, Bahamas coastal zone coral reef. Our primary emphasis was on Acropora palmata (or elkhorn coral), a major reef building coral, which is prevalent in the study area, but is suffering from white band disease. A. palmata is currently being, proposed as an endangered species because its populations have severely declined in many areas of the Caribbean. In addition to the A. palmata biotope, we have collected spectra of at least seven other coral biotopes that exist within the study area, each with different coral community composition, density of corals, relief, and size of corals. Coral spectral reflectance was then input into a radiative transfer model, CORALMOD (CM1), which is based on a leaf radiative transfer model. In CM1, input coral reflectance measurements produce modeled reflectance through an inversion at each visible wavelength to provide the absorption spectrum. Initially, we imposed a scattering baseline that is the same regardless of the coral condition and that coral is optically thick and no light is transmitted through coral. Here we will focus on methodology, experimental design, and initial findings of the in situ spectral measurements and preliminary output from the radiative transfer model.
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.
The radiative transfer theory is applied to interpret polarimetric radar backscatter from pine forest with clustered vegetation structures. The scattering function of each cluster is calculated by incorporating the phase interference of scattered fields from each component. The resulting phase matrix is used in the radiative transfer equations to evaluate the polarimetric backscattering coefficients from random medium layers embedded with vegetation clusters. Upon including multiscale structures (trunks, primary and secondary branches, and needles), polarimetric radar responses from pine forest for different frequencies and looking angles are interpreted and simulated. Preliminary results are shown to be in good agreement with the measured backscattering coefficients at the Landes maritime pine forest during the MAESTRO-1 experiment.
This program implements a spectral line, radiative transfer tool for interpreting Spitzer Space Telescope observations by matching them with models of protostellar disks for improved understanding of planet and star formation. The Spitzer Space Telescope detects gas phase molecules in the infrared spectra of protostellar disks, with spectral lines carrying information on the chemical composition of the material from which planets form. Input to the software includes chemical models developed at JPL. The products are synthetic images and spectra for comparison with Spitzer measurements. Radiative transfer in a protostellar disk is primarily affected by absorption and emission processes in the dust and in molecular gases such as H2, CO, and HCO. The magnitude of the optical absorption and emission is determined by the population of the electronic, vibrational, and rotational energy levels. The population of the molecular level is in turn determined by the intensity of the radiation field. Therefore, the intensity of the radiation field and the population of the molecular levels are inter-dependent quantities. To meet the computational challenges of solving for the coupled radiation field and electronic level populations in disks having wide ranges of optical depths and spatial scales, the tool runs in parallel on the JPL Dell Cluster supercomputer with C++ and Fortran compiler with a Message Passing Interface. Because this software has been developed on a distributed computing platform, the modeling of systems previously beyond the reach of available computational resources is possible.
A computationally efficient method has been developed for highly nonlinear problems in which radiative transfer is an important aspect of the heating and cooling of the medium. An approximate probabilistic radiative transfer equation is derived for one-dimensional plane-parallel atmospheres of finite or semi-infinite extent, for both spectral lines and bound-free continua. Boundary conditions, accuracy, escape probabilities, and practical aspects of complete linearization, are discussed. The method is accurate to a few tens of percent of a wide variety of realistic problems in which frequency redistribution of scattered photons dominates the transfer and escape of radiation.