Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiative transfer simulations”

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 307 records · Page 17

Coherent Effects in Microwave Backscattering Models for Forest Canopies

In modeling forest canopies, several scattering mechanisms are taken into account: 1) volume scattering, 2) surface-volume interaction, and 3) surface scattering from forest floor. Depending on the structural and dielectric characteristics of forest canopies, the relative contribution of each mechanism in the total backscatter signal of an imaging radar can vary. In this paper, two commonly used first order discrete scattering models, Distorted Born Approximation (DBA) and Radiative Transfer (RT) are used to simulate the backscattered power received by polarimetric radars at P-, L-, and C-bands over coniferous and deciduous forests. The difference between the two models resides on the coherent effect in the surface-volume interaction terms.

Microwave Backscattering Models↗

3-D GRMHD Simulations of Disk-Jet Coupling and Emission

We have performed a fully 3-D GRMHD simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The simulation results show that a bipolar jet is initially created. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk. Recently we have developed two new codes: 3 -D GRMHD: RelAtivIStic magnetoHydrodynamica1 sImulatioN (RAISHIN) code constructed by modern high-resolution shock-capturing (HRSC) techniques and 3-D GRPIC code. We have calculated free-free and synchrotron emission from the disks and jet/outflows obtained from our GRMHD simulations using a fully covariant radiative transfer formulation.

Nishikawa, K.-I.↗

High Performance Simulations of Accretion Disk Dynamics and Jet Formations Around Kerr Black Holes

We investigate jet formation in black-hole systems using 3-D General Relativistic Particle-In-Cell (GRPIC) and 3-D GRMHD simulations. GRPIC simulations, which allow charge separations in a collisionless plasma, do not need to invoke the frozen condition as in GRMHD simulations. 3-D GRPIC simulations show that jets are launched from Kerr black holes as in 3-D GRMHD simulations, but jet formation in the two cases may not be identical. Comparative study of black hole systems with GRPIC and GRMHD simulations with the inclusion of radiate transfer will further clarify the mechanisms that drive the evolution of disk-jet systems.

Nishikawa, Ken-Ichi↗

A Global Model Simulation of Aerosol Effects of Surface Radiation Budget- Toward Understanding of the "Dimming to Brightening" Transition

We present a global model study on the role aerosols play in the change of solar radiation at Earth's surface that transitioned from a decreasing (dimming) trend to an increasing (brightening) trend. Our primary objective is to understand the relationship between the long-term trends of aerosol emission, atmospheric burden, and surface solar radiation. More specifically, we use the recently compiled comprehensive global emission datasets of aerosols and precursors from fuel combustion, biomass burning, volcanic eruptions and other sources from 1980 to 2006 to simulate long-term variations of aerosol distributions and optical properties, and then calculate the multi-decadal changes of short-wave radiative fluxes at the surface and at the top of the atmosphere by coupling the GOCART model simulated aerosols with the Goddard radiative transfer model. The model results are compared with long-term observational records from ground-based networks and satellite data. We will address the following critical questions: To what extent can the observed surface solar radiation trends, known as the transition from dimming to brightening, be explained by the changes of anthropogenic and natural aerosol loading on global and regional scales? What are the relative contributions of local emission and long-range transport to the surface radiation budget and how do these contributions change with time?

Chin, Mian↗

Uncertainty Quantification of GEOS-5 L-band Radiative Transfer Model Parameters Using Bayesian Inference and SMOS Observations

Uncertainties in L-band (1.4 GHz) radiative transfer modeling (RTM) affect the simulation of brightness temperatures (Tb) over land and the inversion of satellite-observed Tb into soil moisture retrievals. In particular, accurate estimates of the microwave soil roughness, vegetation opacity and scattering albedo for large-scale applications are difficult to obtain from field studies and often lack an uncertainty estimate. Here, a Markov Chain Monte Carlo (MCMC) simulation method is used to determine satellite-scale estimates of RTM parameters and their posterior uncertainty by minimizing the misfit between long-term averages and standard deviations of simulated and observed Tb at a range of incidence angles, at horizontal and vertical polarization, and for morning and evening overpasses. Tb simulations are generated with the Goddard Earth Observing System (GEOS-5) and confronted with Tb observations from the Soil Moisture Ocean Salinity (SMOS) mission. The MCMC algorithm suggests that the relative uncertainty of the RTM parameter estimates is typically less than 25 of the maximum a posteriori density (MAP) parameter value. Furthermore, the actual root-mean-square-differences in long-term Tb averages and standard deviations are found consistent with the respective estimated total simulation and observation error standard deviations of m3.1K and s2.4K. It is also shown that the MAP parameter values estimated through MCMC simulation are in close agreement with those obtained with Particle Swarm Optimization (PSO).

MCMC↗

Principle Component Analysis of AIRS and CrIS Data

Synthetic Eigen Vectors (EV) used for the statistical analysis of the PC reconstruction residual of large ensembles of data are a novel tool for the analysis of data from hyperspectral infrared sounders like the Atmospheric Infrared Sounder (AIRS) on the EOS Aqua and the Cross-track Infrared Sounder (CrIS) on the SUOMI polar orbiting satellites. Unlike empirical EV, which are derived from the observed spectra, the synthetic EV are derived from a large ensemble of spectra which are calculated assuming that, given a state of the atmosphere, the spectra created by the instrument can be accurately calculated. The synthetic EV are then used to reconstruct the observed spectra. The analysis of the differences between the observed spectra and the reconstructed spectra for Simultaneous Nadir Overpasses of tropical oceans reveals unexpected differences at the more than 200 mK level under relatively clear conditions, particularly in the mid-wave water vapor channels of CrIS. The repeatability of these differences using independently trained SEV and results from different years appears to rule out inconsistencies in the radiative transfer algorithm or the data simulation. The reasons for these discrepancies are under evaluation.

infrared↗

TPSAS-NF1676L-10726-DND

The principal component radiative transfer model PCRTM (Xu 2006) simulates spectral longwave radiance at a 0.5cm-1 resolution over the wavelengths (50-2760cm-1) using 280 principal components. The CERES product SSF containing CERES field-of-view (FOV) cloud properties based on MODIS pixel information are used as inputs. NASA GMAO(Global Modeling and Assimilation Office)Geos-4 temperature and humidity profiles, as well as SMOBA (Stratosphere Monitoring Ozone Blended Analysis) ozone profiles allow computations for five years (Jan 2001- Dec 2006) of FOV nadir spectral radiance. The approximately twice daily FOV PCRTM simulations along with CERES broadband measurements are time-space averaged to 10deg zonal monthly data. Analysis of monthly seasonal and deseasonalized interannual variability in terms of brightness temperature are presented. Comparison of interannual variability in the portions of the spectrum seen by the AIRS instrument are shown. Comparisons of spectral integrated PCRTM simulations to CERES broadband longwave radiance measurements are made. In order to study the sensitivity of the spatially and temporally averaged TOA spectrum to various atmospheric properties, we perturbed atmospheric and surface temperatures, humidity, cloud height, cloud fraction, and particle size. The sensitivity of the spectrum to these properties depends on region and seasons. Month to month variability of the longwave spectra to these perturbations is 10 to 20% of the global monthly mean perturbation response. In addition, it depends on time period and size of the area used for averaging spectra. In this poster, the magnitude of sensitivity variability and how the sensitivity changes depending on temporal and spatial scale of averaging will be discussed.

Fred G. Rose↗

A Novel Atmospheric Correction Algorithm to Exploit the Diurnal Variability in Hypertemporal Geostationary Observations

This study developed a new atmospheric correction algorithm, GeoNEX-AC, that is independent from the traditional use of spectral band ratios but dedicated to exploiting information from the diurnal variability in the hypertemporal geostationary observations. The algorithm starts by evaluating smooth segments of the diurnal time series of the top-of-atmosphere (TOA) reflectance to identify clear-sky and snow-free observations. It then attempts to retrieve the Ross-Thick–Li-Sparse (RTLS) surface bi-directional reflectance distribution function (BRDF) parameters and the daily mean atmospheric optical depth (AOD) with an atmospheric radiative transfer model (RTM) to optimally simulate the observed diurnal variability in the clear-sky TOA reflectance. Once the initial RTLS parameters are retrieved after the algorithm’s burn-in period, they serve as the prior information to estimate the AOD levels for the following days and update the surface BRDF information with the new clear-sky observations. This process is iterated through the full time span of the observations, skipping only totally cloudy days or when surface snow is detected. We tested the algorithm over various Aerosol Robotic Network (AERONET) sites and the retrieved results well agree with the ground-based measurements. This study demonstrates that the high-frequency diurnal geostationary observations contain unique information that can help to address the atmospheric correction problem from new directions.

atmospheric correction↗

Improved Peatland Hydrology and L-Band Microwave Radiative Transfer Modeling in Version 7 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture from April 2015 to present with a mean latency of 2.5 days from the time of observation. The L4_SM estimates are derived from the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model. The recently released Version 7 of L4_SM comprises two key advances. First, the Catchment model now includes the recently developed PEATCLSM hydrology module for peatlands and uses an updated global map of peatlands. Second, revised parameters are used in the L-band radiative transfer model that converts the simulated soil moisture and temperature estimates into Tb predictions for use in the radiance-based L4_SM analysis. Specifically, the L-band scattering albedo, soil roughness, and vegetation opacity climatology are obtained from the SMAP Level-2 dual-channel soil moisture retrieval product. In this presentation, we examine the improvements gained with the model and parameter changes. In peatlands, the dynamics of water table depth, surface soil moisture and evapotranspiration are considerably improved in L4_SM Version 7 over Version 6 when evaluated against in situ measurements. These improvements are also manifested in smaller Tb observation-minus-forecast residuals.

SMAP↗

Principal Component-Based Radiative Transfer Model (PCRTM) for Hyperspectral Remote Sensors for UV, VVIS, NIR, IR, and FIR Spectral Regions

Fast and accurate radiative transfer models are needed to efficiently process satellite hyperspectral remote sensors. We have developed a Principal Component-based radiative transfer model (PCRTM) which can simulate TOA radiance or reflectance spectra in the cloudy atmosphere for far IR, IR, NIR, VIS, and UV spectral regions quickly and accurately. Multi-scattering of multiple layers of clouds/aerosols is included in the model. Polarization components can also be calculated. The computation speed is 3 to 4 orders of magnitude faster than MODTRAN5, LBLRTM, and VLIDORT. The PCRTM calculated radiance spectra agree with reference RTM calculated radiance spectra very well (0.03 K in IR and 0.05% in solar). Comparisons of the PCRTM model calculations with observed AIRS, CrIS, IASI, NAST-I, and SCIAMACHY data will be presented. The solar-PCRTM has been developed for CLARREO PathFinder (CPF) and has been extensively used for CPF science algorithms. The PCRTM can also be used for many current and future NASA future missions such as EMIT, SBG, TEMPO, and PACE.

pcrtm↗

A Principal-Component-Based Radiative Transfer Model (PCRTM) for Hyperspectral Shortwave and Longwave Satellite Sensors and Its Applications

The radiative transfer model (RTM) or forward model is an essential component in satellite remote sensing. For modern hyperspectral remote sensors, fast and accurate RTMs are needed due to a large number of spectral dimensions and high spatial resolutions. We will describe a Principal Component-based radiative transfer model (PCRTM) which can simulate the top-of-atmosphere (TOA) radiance or reflectance spectra 250 nm to 2000 micrometers quickly and accurately. The PCRTM has been demonstrated to be extremely accurate, compared to the line-by-line RTM benchmarks, and the former is several orders more computationally efficient than the latter. We will demonstrate how the PCRTM and the associated inversion algorithms are used to infer atmospheric temperature, moisture, and trace gas profiles, as well as cloud and surface properties from hyperspectral IR sounders such as Atomspheric Infrared Souder (AIRS) and Cross-track Infrared Sounder (CrIS). High-quality climate records for a 20-year duration have been derived from these IR hyperspectral data. Finally, we will show some examples of using PCRTM to retrieve cloud properties from Earth Surface Mineral Dust Source Investigation (EMIT) and its applicability of PCRTM to future missions such as the CLARREO (Climate Absolute Radiance and Refractivity Observatory) Pathfinder (CPF) CPF and the Surface Biology and Geology (SBG).

Xu Liu↗

Enabling the Assimilation of CrIS Shortwave Infrared Observations in Global NWP at NOAA. Part I: Background and Methods

Radiance observations from Earth-observing satellites have a significant positive impact on numerical weather prediction (NWP) forecasts, but some spectral regions are not fully exploited. Observations from hyperspectral infrared (IR) sounders in the longwave region (650-1100 cm -1 ), for instance, are routinely assimilated in many NWP models, but observations in the shortwave region (2155-2550 cm -1 ) are not. Each of these regions provides information on the temperature structure of the atmosphere, but the shortwave IR (SWIR) region is considered challenging to assimilate due to noise equivalent delta temperature (NEDT) that is highly variable depending on scene brightness temperature and to phenomena that are difficult to model, like non-Local Thermodynamic Equilibrium (NLTE) and solar reflectance. With recent advances in small-satellite technology, SWIR temperature sounders may provide an agile and cost-effective complement to the current constellation of IR sounders. Therefore, a better understanding of the use and impact of SWIR observations in data assimilation for NWP is warranted. In part one of this study, as presented here, the amount of unique information (as determined by Empirical Orthogonal Decomposition (EOD)) made available to a data assimilation system by Cross-track Infrared Sounder (CrIS) SWIR observations is reviewed, recent advancements to the Community Radiative Transfer Model (CRTM) for the simulation of CrIS shortwave radiances are tested, and enhancements to NOAA’s Global Data Assimilation System (GDAS) for the assimilation of CrIS SWIR observations are implemented and evaluated. Part two of this study, which seeks to assess the value of assimilating shortwave IR observations in global NWP, is also introduced.

Erin Jones↗

Radiative interactions in chemically reacting compressible nozzle flows using Monte Carlo simulations

The two-dimensional spatially elliptic Navier-Stokes equations have been used to investigate the radiative interactions in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The radiative heat transfer term in the energy equation is simulated using the Monte Carlo method (MCM). The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The spectral correlation has been considered in the Monte Carlo formulations. Results obtained demonstrate that the effect of radiation on the flow field is minimal but its effect on the wall heat transfer is significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and the nozzle size on the radiative and conductive wall fluxes.

Liu, J.↗

Vector Radiative Transfer Code SORD: Performance Analysis and Quick Start Guide

We present a new open source polarized radiative transfer code SORD written in Fortran 9095. SORD numerically simulates propagation of monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface using the method of successive orders of scattering (hence the name). Thermal emission is ignored. We did not improve the method in any way, but report the accuracy and runtime in 52 benchmark scenarios. This paper also serves as a quick start users guide for the code available from ftp:maiac.gsfc.nasa.govpubskorkin, from the JQSRT website, or from the corresponding (first) author.

polarized radiative transfer↗

A Fast Hybrid (3-D/1-D) Model for Thermal Radiative Transfer in Cirrus via Successive Orders of Scattering

We investigate the impact of cirrus cloud heterogeneity on the direct emission by cloud or surface and on the scattering by ice particles in the thermal infrared (TIR). Realistic 3-D cirri are modeled with the 3DCLOUD code, and top-of-atmosphere radiances are simulated by the 3-D Monte Carlo radiative transfer (RT) algorithm 3DMCPOL for two (8.65 micrometers and 12.05 micrometers) channels of the Imaging Infrared Radiometer on CALIPSO. At nadir, comparisons of 1-D and 3-D RT show that 3-D radiances are larger than their 1-D counterparts for direct emission but smaller for scattered radiation. For our cirrus cases, 99% of the 3-D total radiance is computed by the third scattering order, which corresponds to 90% of the total computational effort, but larger optical thicknesses need more scattering orders. To radically accelerate the 3-D RT computations (using only few percent of 3-D RT time with a Monte Carlo code), even in the presence of large optical depths, we develop a hybrid model based on exact 3-D direct emission, the first scattering order from 1-D in each homogenized column, and an empirical adjustment linearly dependent on the optical thickness to account for higher scattering orders. Good agreement is found between the hybrid model and the exact 3-D radiances for two very different cirrus models without changing the empirical parameters. We anticipate that a future deterministic implementation of the hybrid model will be fast enough to process multiangle thermal imagery in a practical tomographic reconstruction of 3-D cirrus fields.

radiative transfer↗

Simulated directional radiances of vegetation from satellite platforms

A computer modeling and simulation investigated the effects of off-nadir viewing, canopy geometry and density, solar zenith angle, and atmospheric conditions on the radiance of several vegetated targets. Effects on the normalized difference (ND) index were also investigated, and two wavelengths of 0.68 and 0.80 micron were simulated at nine sensor zenith angles, five sensor azimuths, and nine solar zenith angles. Seven canopies of low, medium, and high densities were observed, and bidirectional reflectance functions of other canopies were discussed. A modified Turner's (1974) atmospheric radiative transfer model was used for 2300 simulation runs. Off-nadir viewing effects were more pronounced in the red than the infrared, and the variability was found to be a function of canopy geometry. The ND index eliminated much of the variability seen in the individual bands, while the overall variability decreased with increasing biomass, and increased with increasing solar zenith angles.

Kirchner, J. A.↗

Fluorescence Correction Vegetation Index (FCVI): A physically based reflectance index to separate physiological and non-physiological information in far-red sun-induced chlorophyll fluorescence

Sun-induced chlorophyll fluorescence (SIF) has been used to track vegetation photosynthetic activity for improving estimation of gross primary productivity (GPP) and detecting plant stress. There are both physical and physiological controls of SIF measured at the surface and retrieved from remote sensing including satellite observations. In order to accurately use SIF for monitoring of plant physiology, the effects of physically-based radiation processes related to leaf and canopy structure, notably photosynthetically active radiation (PAR) absorption and SIF scattering and re-absorption, must be characterized. In this study, we investigate both PAR absorption and SIF scattering processes and find that although it is difficult to quantify their effects individually by using just reflectance, the combined effects of the two processes can be well approximated by a reflectance index. This index, referred to as FCVI (Fluorescence Correction Vegetation Index), is defined as the difference between near-infrared (NIR) and broad-band visible (VIS, 400–700 nm) reflectance acquired under identical sun-canopy-observer geometry of the SIF measurements. The development of the index was based on the physical connection between reflectance and far-red SIF, which was revealed by using the spectral invariant theory. The utility of FCVI to correct far-red SIF for PAR absorption and scattering effects, thus improving the link to photosynthesis, was tested with data from: (i) a field experiment for a growing season; and (ii) a numerical experiment which included a number of scenarios generated by a radiative transfer model. For both the observations and simulations, the FCVI provided a promising estimate of the impact of the physically-based radiation processes on far-red SIF of moderately dense canopies (i.e., FCVI ≥ 0.18). Normalizing the TOC far-red SIF by both the incident PAR (iPAR) and the FCVI provided a good estimate of the far-red fluorescence emission efficiency of the canopies examined. This approach enhances our ability to generalize retrievals for vegetation processes as they change through natural growth phases and seasons. Taken together, far-red SIF and FCVI may enable the assessment of the light partitioning of vegetation canopies, an essential step to facilitate the use of far-red SIF for tracking physiological processes.

Peiqi Yang↗

Machine learning based noise reduction for satellite products: application to solar-induced fluorescence retrievals using simulated and real data

In the past two decades, global satellite measurements of terrestrial chlorophyll solar-induced fluorescence (SIF) have been used widely for a number of different applications related to physiology, phenology, and productivity of plants. However, SIF retrievals are inherently noisy due to the relatively small SIF spectral signature in comparison with observational noise. In this work, we examine how a spectral-based approach that employs principal component analysis along with a relatively shallow artificial neural network can be used to reduce noise and other artifacts in satellite level 2 (L2) products. We first apply the approach in a controlled environment in which radiance spectra are simulated with a full atmospheric and surface radiative transfer model for different scenarios including various SIF values that are known. Various levels of noise can be added to the simulated spectra. Resulting noisy and noise-reduced SIF retrievals are compared with the true values to assess performance. We then apply the noise reduction approach to real SIF derived from instruments flying on meteorological satellites. The results are evaluated by comparing SIF retrievals from different platforms with each other and with other independent data sets, showing enhanced capability to capture seasonal and interannual variability in SIF.

Chlorophyll fluorescence↗