scholarly journals Diffusion approximation in a radiative transfer model for astrophysical flows

2018 ◽  
Vol 64 ◽  
pp. 78-92
Author(s):  
Laurent Di Menza ◽  
Claire Michaut ◽  
Océane Saincir

In this work, we present the diffusion approximation model for radiative transfer when we deal with optically thick astrophysical flows. Since the initial model is high CPU time demanding when dealing with its numerical approximation, solving this simpler system can provide a low cost strategy for the simulation of radiative media. We then use a finite-volume algorithm coupled with an implicit scheme for radiative contributions to solve this simplified system. Numerical experiments in the one-dimensional and two dimensional cases are presented to validate our numerical strategy and to prove the relevance of this asymptotic model.

2016 ◽  
Author(s):  
Giuliano Liuzzi ◽  
Guido Masiello ◽  
Carmine Serio ◽  
Daniela Meloni ◽  
Claudia Di Biagio ◽  
...  

Abstract. In the context of the ChArMEx campaign, we present here some results concerning the quantitative comparison between simulated and observed radiances during a dust event occurred between June and July 2013 in the southern Mediterranean basin, involving the airmass above Lampedusa island. In particular, comparisons have been performed between radiances as observed by the Infrared Atmospheric Sounder Interferometer (IASI) and those simulated using the σ-IASI-as radiative transfer model, which takes into account aerosol extinction effect through a set of fast parameterizations. Simulations have been carried on with different sets of input complex refractive indices, which take into account the parent soils of the aerosols, and using the high-quality characterization of desert dust aerosol microphysical properties, achieved through direct measurements in the ChArMEx experiment; on the one hand, this comparison has offered the possibility to test the feasibility of the radiative transfer model. On the other hand, this work goes through a direct validation of different refractive indices sets for desert dust in the thermal infrared. Results show a good consistency between calculations and observations, especially in the spectral interval 800–1000 cm−1; moreover, the comparison between calculations and observations suggests that further efforts are needed to better characterize desert dust optical properties in the short wave (above 2000 cm−1). In any case, we show that it is necessary to properly tune the refractive indices according to the geographical origin of the observed aerosol.


2013 ◽  
Vol 53 (A) ◽  
pp. 832-838
Author(s):  
Smadar Bressler ◽  
Giora Shaviv ◽  
Nir J. Shaviv

We present a radiative transfer model for Earth-Like-Planets (ELP). The model allows the assessment of the effect of a change in the concentration of an atmospheric component, especially a greenhouse gas (GHG), on the surface temperature of a planet. The model is based on the separation between the contribution of the short wavelength molecular absorption and the long wavelength one. A unique feature of the model is the condition of energy conservation at every point in the atmosphere. The radiative transfer equation is solved in the two stream approximation without assuming the existence of an LTE in any wavelength range. The model allows us to solve the Simpson paradox, whereby the greenhouse effect (GHE) has no temperature limit. On the contrary, we show that the temperature saturates, and its value depends primarily on the distance of the planet from the central star. We also show how the relative humidity affects the surface temperature of a planet and explain why the effect is smaller than the one derived when the above assumptions are neglected.


2006 ◽  
Vol 6 (1) ◽  
pp. 1199-1248 ◽  
Author(s):  
F. Spada ◽  
M. C. Krol ◽  
P. Stammes

Abstract. A new multiple-scattering Monte Carlo 3-D radiative transfer model named McSCIA (Monte Carlo for SCIAmachy) is presented. The backward technique is used to efficiently simulate narrow field of view instruments. The McSCIA algorithm has been formulated as a function of the Earth's radius, and can thus perform simulations for both plane-parallel and spherical atmospheres. The latter geometry is essential for the interpretation of limb satellite measurements, as performed by SCIAMACHY on board of ESA's Envisat. The model can simulate UV-vis-NIR radiation. First the ray-tracing algorithm is presented in detail, and then successfully validated against literature references, both in plane-parallel and in spherical geometry. A simple 1-D model is used to explain two different ways of treating absorption. One method uses the single scattering albedo while the other uses the equivalence theorem. The equivalence theorem is based on a separation of absorption and scattering. It is shown that both methods give, in a statistical way, identical results for a wide variety of scenarios. Both absorption methods are included in McSCIA, and it is shown that also for a 3-D case both formulations give identical results. McSCIA limb profiles for atmospheres with and without absorption compare well with the one of the state of the art Monte Carlo radiative transfer model MCC++. A simplification of the photon statistics may lead to very fast calculations of absorption features in the atmosphere. However, these simplifications potentially introduce biases in the results. McSCIA does not use simplifications and is therefore a relatively slow implementation of the equivalence theorem. For the first time, however, the validity of the equivalence theorem is demonstrated in a spherical 3-D radiative transfer model.


2006 ◽  
Vol 6 (12) ◽  
pp. 4823-4842 ◽  
Author(s):  
F. Spada ◽  
M. C. Krol ◽  
P. Stammes

Abstract. A new multiple-scattering Monte Carlo 3-D radiative transfer model named McSCIA (Monte Carlo for SCIAmachy) is presented. The backward technique is used to efficiently simulate narrow field of view instruments. The McSCIA algorithm has been formulated as a function of the Earth's radius, and can thus perform simulations for both plane-parallel and spherical atmospheres. The latter geometry is essential for the interpretation of limb satellite measurements, as performed by SCIAMACHY on board of ESA's Envisat. The model can simulate UV-vis-NIR radiation. First the ray-tracing algorithm is presented in detail, and then successfully validated against literature references, both in plane-parallel and in spherical geometry. A simple 1-D model is used to explain two different ways of treating absorption. One method uses the single scattering albedo while the other uses the equivalence theorem. The equivalence theorem is based on a separation of absorption and scattering. It is shown that both methods give, in a statistical way, identical results for a wide variety of scenarios. Both absorption methods are included in McSCIA, and it is shown that also for a 3-D case both formulations give identical results. McSCIA limb profiles for atmospheres with and without absorption compare well with the one of the state of the art Monte Carlo radiative transfer model MCC++. A simplification of the photon statistics may lead to very fast calculations of absorption features in the atmosphere. However, these simplifications potentially introduce biases in the results. McSCIA does not use simplifications and is therefore a relatively slow implementation of the equivalence theorem.


2012 ◽  
Vol 33 (6) ◽  
pp. 1611-1624 ◽  
Author(s):  
Iñigo Mendikoa ◽  
Santiago Pérez-Hoyos ◽  
Agustín Sánchez-Lavega

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Rehman S. Eon ◽  
Charles M. Bachmann

AbstractThe advent of remote sensing from unmanned aerial systems (UAS) has opened the door to more affordable and effective methods of imaging and mapping of surface geophysical properties with many important applications in areas such as coastal zone management, ecology, agriculture, and defense. We describe a study to validate and improve soil moisture content retrieval and mapping from hyperspectral imagery collected by a UAS system. Our approach uses a recently developed model known as the multilayer radiative transfer model of soil reflectance (MARMIT). MARMIT partitions contributions due to water and the sediment surface into equivalent but separate layers and describes these layers using an equivalent slab model formalism. The model water layer thickness along with the fraction of wet surface become parameters that must be optimized in a calibration step, with extinction due to water absorption being applied in the model based on equivalent water layer thickness, while transmission and reflection coefficients follow the Fresnel formalism. In this work, we evaluate the model in both field settings, using UAS hyperspectral imagery, and laboratory settings, using hyperspectral spectra obtained with a goniometer. Sediment samples obtained from four different field sites representing disparate environmental settings comprised the laboratory analysis while field validation used hyperspectral UAS imagery and coordinated ground truth obtained on a barrier island shore during field campaigns in 2018 and 2019. Analysis of the most significant wavelengths for retrieval indicate a number of different wavelengths in the short-wave infra-red (SWIR) that provide accurate fits to measured soil moisture content in the laboratory with normalized root mean square error (NRMSE)< 0.145, while independent evaluation from sequestered test data from the hyperspectral UAS imagery obtained during the field campaign obtained an average NRMSE = 0.169 and median NRMSE = 0.152 in a bootstrap analysis.


2020 ◽  
Vol 13 (1) ◽  
pp. 116
Author(s):  
Lucie Leonarski ◽  
Laurent C.-Labonnote ◽  
Mathieu Compiègne ◽  
Jérôme Vidot ◽  
Anthony J. Baran ◽  
...  

The present study aims to quantify the potential of hyperspectral thermal infrared sounders such as the Infrared Atmospheric Sounding Interferometer (IASI) and the future IASI next generation (IASI-NG) for retrieving the ice cloud layer altitude and thickness together with the ice water path. We employed the radiative transfer model Radiative Transfer for TOVS (RTTOV) to simulate cloudy radiances using parameterized ice cloud optical properties. The radiances have been computed from an ice cloud profile database coming from global operational short-range forecasts at the European Center for Medium-range Weather Forecasts (ECMWF) which encloses the normal conditions, typical variability, and extremes of the atmospheric properties over one year (Eresmaa and McNally (2014)). We performed an information content analysis based on Shannon’s formalism to determine the amount and spectral distribution of the information about ice cloud properties. Based on this analysis, a retrieval algorithm has been developed and tested on the profile database. We considered the signal-to-noise ratio of each specific instrument and the non-retrieved atmospheric and surface parameter errors. This study brings evidence that the observing system provides information on the ice water path (IWP) as well as on the layer altitude and thickness with a convergence rate up to 95% and expected errors that decrease with cloud opacity until the signal saturation is reached (satisfying retrievals are achieved for clouds whose IWP is between about 1 and 300 g/m2).


Sign in / Sign up

Export Citation Format

Share Document