Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiative transfer equation”

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.

At least 325 records · Page 18

Using Satellite Data in Weather Forecasting: I

The GOES Product Generation System (GPGS) is a set of computer codes and scripts that enable the assimilation of real-time Geostationary Operational Environmental Satellite (GOES) data into regional-weather-forecasting mathematical models. The GPGS can be used to derive such geophysical parameters as land surface temperature, the amount of precipitable water, the degree of cloud cover, the surface albedo, and the amount of insolation from satellite measurements of radiant energy emitted by the Earth and its atmosphere. GPGS incorporates a priori information (initial guesses of thermodynamic parameters of the atmosphere) and radiometric measurements from the geostationary operational environmental satellites along with mathematical models of physical principles that govern the transfer of energy in the atmosphere. GPGS solves the radiative-transfer equation and provides the resulting data products in formats suitable for use by weather-forecasting computer programs. The data-assimilation capability afforded by GPGS offers the potential to improve local weather forecasts ranging from 3 hours to 2 days - especially with respect to temperature, humidity, cloud cover, and the probability of precipitation. The improvements afforded by GPGS could be of interest to news media, utility companies, and other organizations that utilize regional weather forecasts.

Jedlovec, Gary J.↗

SMOS/SMAP Synergy for SMAP Level 2 Soil Moisture Algorithm Evaluation

Soil Moisture Active Passive (SMAP) satellite has been proposed to provide global measurements of soil moisture and land freeze/thaw state at 10 km and 3 km resolutions, respectively. SMAP would also provide a radiometer-only soil moisture product at 40-km spatial resolution. This product and the supporting brightness temperature observations are common to both SMAP and European Space Agency's Soil Moisture and Ocean Salinity (SMOS) mission. As a result, there are opportunities for synergies between the two missions. These include exploiting the data for calibration and validation and establishing longer term L-band brightness temperature and derived soil moisture products. In this investigation we will be using SMOS brightness temperature, ancillary data, and soil moisture products to develop and evaluate a candidate SMAP L2 passive soil moisture retrieval algorithm. This work will begin with evaluations based on the SMOS product grids and ancillary data sets and transition to those that will be used by SMAP. An important step in this analysis is reprocessing the multiple incidence angle observations provided by SMOS to a global brightness temperature product that simulates the constant 40 degree incidence angle observations that SMAP will provide. The reprocessed brightness temperature data provide a basis for evaluating different SMAP algorithm alternatives. Several algorithms are being considered for the SMAP radiometer-only soil moisture retrieval. In this first phase, we utilized only the Single Channel Algorithm (SCA), which is based on the radiative transfer equation and uses the channel that is most sensitive to soil moisture (H-pol). Brightness temperature is corrected sequentially for the effects of temperature, vegetation, roughness (dynamic ancillary data sets) and soil texture (static ancillary data set). European Centre for Medium-Range Weather Forecasts (ECMWF) estimates of soil temperature for the top layer (as provided as part of the SMOS ancillary data) were used to correct for surface temperature effects and to derive microwave emissivity. ECMWF data were also used for precipitation forecasts, presence of snow, and frozen ground. Vegetation options are described below. One year of soil moisture observations from a set of four watersheds in the U.S. were used to evaluate four different retrieval methodologies: (1) SMOS soil moisture estimates (version 400), (2) SeA soil moisture estimates using the SMOS/SMAP data with SMOS estimated vegetation optical depth, which is part of the SMOS level 2 product, (3) SeA soil moisture estimates using the SMOS/SMAP data and the MODIS-based vegetation climatology data, and (4) SeA soil moisture estimates using the SMOS/SMAP data and actual MODIS observations. The use of SMOS real-world global microwave observations and the analyses described here will help in the development and selection of different land surface parameters and ancillary observations needed for the SMAP soil moisture algorithms. These investigations will greatly improve the quality and reliability of this SMAP product at launch.

Bindlish, Rajat↗

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

A Spatially Coherent Global Soil Moisture Product with Improved Temporal Resolution

Global soil moisture products that are completely independent of any type of ancillary data and solely rely on satellite observations are presented. Additionally, we further develop an existing downscaling technique that enhances the spatial resolution of such products to approximately 11 km. These products are based on internal modules of the Land Parameter Retrieval Model (LPRM), an algorithm that uses the radiative transfer equation to link soil moisture, vegetation optical depth and land surface temperature to observed brightness temperatures. The soil moisture product that is independent of any type of ancillary data uses the internally calculated dielectric constant as a soil moisture proxy. This data product is not influenced by errors associated with coarse-scale global soil property maps or by any other type of forcing (e.g. re-analysis) data and is therefore solely based on satellite microwave observations. The second step builds upon recent developments to increase the spatial resolution of the LPRM retrievals using a smoothing filter downscaling method. With this method we can attain a spatial resolution that can be more useful at the scale of local and regional hydrological studies as well. The steps presented in this paper were applied to observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). The newly derived data sets were validated using ground-based observations from the International Soil Moisture Network (ISMN). The internally calculated dielectric constant product results in significantly more days with valid retrievals than the original soil moisture data products, in particular over arid regions. The dielectric constant product resulted in similar correlations with in situ data as the original soil moisture data product. Together, these findings demonstrate the usefulness of this new dielectric constant product for the hydrological modeling community and climate studies. A case study on the Australian Fitzroy catchment demonstrated that the downscaled data product has a more detailed spatial description of soil moisture, especially during wet and dry conditions with more pronounced dry and wet regions within the catchment. The increased resolution data products could therefore improve runoff predictions and this study demonstrated the potential added value of a transitioning from a spatial resolution of 56 km toward a higher resolution of 11 km. The hydrological implications of these newly developed data records are not only linked to AMSR-E satellite data, but also to the next generation Soil Moisture Active and Passive (SMAP) mission where a 9 km spatial resolution is the target resolution for satellite soil moisture products. The new data products will not replace the current LPRM products, but will be added to the existing array of data products and will become publicly available through our data portals.

spatial↗

Far-Field Lorenz-Mie Scattering in an Absorbing Host Medium: Theoretical Formalism and FORTRAN Program

In this paper we make practical use of the recently developed first-principles approach to electromagnetic scattering by particles immersed in an unbounded absorbing host medium. Specifically, we introduce an actual computational tool for the calculation of pertinent far-field optical observables in the context of the classical Lorenzâ€"Mie theory. The paper summarizes the relevant theoretical formalism, explains various aspects of the corresponding numerical algorithm, specifies the input and output parameters of a FORTRAN program available at https://www.giss.nasa.gov/staff/mmishchenko/Lorenz-Mie.html, and tabulates benchmark results useful for testing purposes. This public-domain FORTRAN program enables one to solve the following two important problems: (i) simulate theoretically the reading of a remote well-collimated radiometer measuring electromagnetic scattering by an individual spherical particle or a small random group of spherical particles; and (ii) compute the single-scattering parameters that enter the vector radiative transfer equation derived directly from the Maxwell equations.

Far-field electromagnetic scattering; Absorbing ho↗

A Neural Network Correction to the Scalar Approximation in Radiative Transfer

The next generation of advanced high-resolution sensors in geostationary orbit will gather detailed information for studying the Earth system. There is an increasing desire to perform observing system simulation experiments (OSSEs) for new sensors during the development phase of the mission in order to better leverage information content from the new and existing sensors. Forward radiative transfer calculations that simulate the observing characteristics of a new instrument are the first step to an OSSE, and they are computationally intensive. The scalar approximation to the radiative transfer equation, a simplification of the vector representation, can save considerable computational cost, but produces errors in top of the atmosphere (TOA) radiance as large as 10% due to neglecting polarization effects. This article presents an artificial neural network technique to correct scalar TOA radiance over both land and ocean surfaces to within 1% of vector-calculated radiance. A neural network was trained on a database of scalar-vector TOA radiance differences at a large range of solar and viewing angles for several thousand realistic atmospheric vertical profiles that were sampled from a high resolution (7 km) global atmospheric transport model. The profiles include Rayleigh scattering and aerosol scattering and absorption. Training and validation of the neural network was demonstrated for two wavelengths in the ultraviolet-visible (US-Vis) spectral range (354 nm and 670 nm). The significant computational savings accrued from using a scalar approximation plus neural network correction approach to simulating TOA radiance will make feasible hyperspectral forward simulations of high-resolution sensors on geostationary satellites, such as TEMPO, GOES-R, GEMS, and SENTINEL-4.

TOA↗

Matrix Exponential in C/C++ Version of Vector Radiative Transfer Code IPOL

We use only left eigenvectors to evaluate the matrix exponential in the method of discrete ordinates for the vector radiative transfer equation, which neglects circular polarization, in a plane-parallel atmosphere. This is contrary to a common practice of using the right eigenvectors to evaluate the matrix exponential combined with the left eigenvectors to avoid the inversion of the matrix of the right ones. Two numerical tests for Rayleigh and Aerosol scattering confirm our idea. For better explanation of our approach and for independent crosscheck of our results, we distribute an example in C/C++.

Sergey Korkin↗

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

Fitting Strategies of Accretion Column Models and Application to the Broadband Spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT data-set of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies

Philipp Thalhammer↗

Applicaton of radiative transfer theory to microwave transmission medium calibrations

Precise determinations of the transmission medium loss and noise temperature contribution which are important to the performance characterization of low noise microwave receiving systems and thermal noise standards are discussed. Tropospheric loss is frequently inferred from microwave radiometer noise temperature measurements. Interpretation of these measurements requires an inversion of the radiative transfer integral equation. This is inconvenient even with computer techniques. Solutions of a rapidly convergent power series of the radiative transfer equations are presented. This solution is applicable to a low loss medium with either uniform or nonuniform loss distributions. A four layer atmosphere model is investigated to demonstrate the accuracy of the solution relative to the model. Applications include thermal noise standards and single- and dual-frequency water radiometers.

Stelzried, C. T.↗

A new formulation for anisotropic radiative transfer problems. I - Solution with a variational technique

The equations of radiative transfer in anisotropically scattering media are reformulated as linear operator equations in a single independent variable. The resulting equations are suitable for solution by a variety of standard mathematical techniques. The operators appearing in the resulting equations are in general nonsymmetric; however, it is shown that every bounded linear operator equation can be embedded in a symmetric linear operator equation and a variational solution can be obtained in a straightforward way. For purposes of demonstration, a Rayleigh-Ritz variational method is applied to three problems involving simple phase functions. It is to be noted that the variational technique demonstrated is of general applicability and permits simple solutions for a wide range of otherwise difficult mathematical problems in physics.

Cheyney, H., III↗

Thermal protection for hypervelocity flight in earth's atmosphere by use of radiation backscattering ablating materials

A heat-shield-material response code predicting the transient performance of a material subject to the combined convective and radiative heating associated with the hypervelocity flight is developed. The code is dynamically interactive to the heating from a transient flow field, including the effects of material ablation on flow field behavior. It accomodates finite time variable material thickness, internal material phase change, wavelength-dependent radiative properties, and temperature-dependent thermal, physical, and radiative properties. The equations of radiative transfer are solved with the material and are coupled to the transfer energy equation containing the radiative flux divergence in addition to the usual energy terms.

Howe, John T.↗

Heat Transfer In High-Temperature Multilayer Insulation

The combined radiation/conduction heat transfer in high-temperature multilayer insulations for typical reentry of reusable launch vehicles from low Earth orbit was investigated experimentally and numerically. The high-temperature multilayer insulation investigated consisted of gold-coated reflective foils separated by alumina fibrous insulation spacers. The steady-state heat transfer through four multilayer insulation configurations was investigated experimentally over the temperature range of 300-1300 K and environmental pressure range of 1.33 10(exp -5)-101.32 kPa. It was shown that including the reflective foils reduced the effective thermal conductivity compared to fibrous insulation sample at 1.5 times the density of the multilayer sample. A finite volume numerical model was developed to solve the governing combined radiation/conduction heat transfer equations. The radiation heat transfer in the fibrous insulation spacers was modeled using the modified two-flux approximation assuming anisotropic scattering and gray medium. The numerical model was validated by comparison with steady-state experimental data. The root mean square deviation between the predicted and measured effective thermal conductivity of the samples was 9.5%.

Daryabeigi, Kamran↗

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗

Radiative transfer in a sphere illuminated by a parallel beam - An integral equation approach

The problem of multiple scattering of nonpolarized light in a planetary body of arbitrary shape illuminated by a parallel beam is formulated using the integral equation approach. There exists a simple functional whose stationarity condition is equivalent to solving the equation of radiative transfer and whose value at the stationary point is proportional to the differential cross section. The analysis reveals a direct relation between the microscopic symmetry of the phase function for each scattering event and the macroscopic symmetry of the differential cross section for the entire planetary body, and the interconnection of these symmetry relations and the variational principle. The case of a homogeneous sphere containing isotropic scatterers is investigated in detail. It is shown that the solution can be expanded in a multipole series such that the general spherical problem is reduced to solving a set of decoupled integral equations in one dimension. Computations have been performed for a range of parameters of interest, and illustrative examples of applications to planetary problems as provided.

Shia, R.-L.↗