Engineering PapersSearch

SEARCH · Engineering Papers

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 55 records · Page 3

Developing a Spectral Correlation based Method for mitigating Angular Mismatch Effects in Intercalibration Using a Benchmark Hyperspectral Sensor

The CLARREO Pathfinder (CPF) mission will implement a state-of-the-art intercalibration method for transferring CPF’s in-orbit Système Internationale (SI)-traceable reference to the shortwave channel of the Clouds and the Earth’s Radiance Energy System (CERES) and the reflective solar bands of the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the NOAA-20 satellite platform with a targeted intercalibration methodology uncertainty of 0.3% (k=1). In order to achieve such a high intercalibration accuracy, the CPF intercalibration measurements to be scheduled need to closely match those from CERES and VIIRS in time, space, angles, and wavelength. Despite’s CPF’s two-axis pointing capability to match its boresight line of sight to that of a target sensor, there will be finite residual differences in angular samplings of the two instruments. The potential angular mismatch between CPF and the target sensors can introduce systematic errors in the intercalibration results, and therefore needs to be corrected using a rigorously designed algorithm. The CPF intercalibration team has been developing a correction method for mitigating the impact of these angular anisotropic effects in the CPF-VIIRS and CPF-CERES intercalibration. The method explores the spectral correlation relationship between the reflected solar radiances from the same surface target that would be measured by CPF at two adjacent angles. Our studies have shown that the spectral information based on CPF measurements can be used to accurately predict the spectral radiance or reflectance difference due to a given mismatch in the viewing and solar geometry. The hyper-spectral information can also be used to provide scene stratification without using any auxiliary data, which is critical to reduce the angular correction uncertainty. The angular correction relationship can be well established using simulated CPF-like top-of-atmosphere spectral radiances observed at all sorts of viewing geometry and solar angles and for different scenes. Intensive radiative transfer simulations have been completed using a state-of-art radiative transfer model developed by the CPF team members. The implementation of the algorithm on high-fidelity event simulation data and the characterization for the angular adjustment uncertainty will be presented.

Wan Wu

Do Fractal Models of Clouds Produces the Right 3D Radiative Effects?

Stochastic fractal models of clouds are often used to study 3D radiative effects and their influence on the remote sensing of cloud properties. Since it is important that the cloud models produce a correct radiative response, some researchers require the model parameters to match observed cloud properties such as scale-independent optical thickness variability. Unfortunately, matching these properties does not necessarily imply that the cloud models will cause the right 3D radiative effects. First, the matched properties alone only influence the 3D effects but do not completely determine them. Second, in many cases the retrieved cloud properties have been already biased by 3D radiative effects, and so the models may not match the true real clouds. Finally, the matched cloud properties cannot be considered independent from the scales at which they have been retrieved. This paper proposes an approach that helps ensure that fractal cloud models are realistic and produce the right 3D effects. The technique compares the results of radiative transfer simulations for the model clouds to new direct observations of 3D radiative effects in satellite images.

Varnai, Tamas

An Assessment of Differences Between Cloud Effective Particle Radius Retrievals for Marine Water Clouds from Three MODIS Spectral Bands

The Moderate Resolution Imaging Spectroradiometer (MODIS) cloud product provides three separate 1 km resolution retrievals of cloud particle effective radii (r (sub e)), derived from 1.6, 2.1 and 3.7 micron band observations. In this study, differences among the three size retrievals for maritime water clouds (designated as r (sub e), 1.6 r (sub e), 2.1 and r (sub e),3.7) were systematically investigated through a series of case studies and global analyses. Substantial differences are found between r (sub e),3.7 and r (sub e),2.1 retrievals (delta r (sub e),3.7-2.l), with a strong dependence on cloud regime. The differences are typically small, within +/- 2 micron, over relatively spatially homogeneous coastal stratocumulus cloud regions. However, for trade wind cumulus regimes, r (sub e),3.7 was found to be substantially smaller than r (sub e),2.1, sometimes by more than 10 micron. The correlation of delta r(sub e),3.7-2.1 with key cloud parameters, including the cloud optical thickness (tau), r (sub e) and a cloud horizontal heterogeneity index (H-sigma) derived from 250 m resolution MODIS 0.86 micron band observations, were investigated using one month of MODIS Terra data. It was found that differences among the three r (sub e) retrievals for optically thin clouds (tau <5) are highly variable, ranging from - 15 micron to 10 micron, likely due to the large MODIS retrieval uncertainties when the cloud is thin. The delta r (sub e),3.7-2.1 exhibited a threshold-like dependence on both r (sub e),2.l and H-sigma. The re,3.7 is found to agree reasonably well with re,2.! when re,2.l is smaller than about 15J-Lm, but becomes increasingly smaller than re,2.1 once re,2.! exceeds this size. All three re retrievals showed little dependence when H-sigma < 0.3 (defined as standard deviation divided by the mean for the 250 m pixels within a 1 km pixel retrieval). However, for H-=sigma >0.3, both r (sub e),1.6 and r (sub e),2.1 were seen to increase quickly with H-sigma. On the other hand, r (sub e),3.7 statistics showed little dependence on H-sigma and remained relatively stable over the whole range of H-sigma values. Potential contributing causes to the substantial r (sub e),3.7 and r (sub e),2.1 differences are discussed. In particular, based on both 1-D and 3-D radiative transfer simulations, we have elucidated mechanisms by which cloud heterogeneity and 3-D radiative effects can cause large differences between r (sub e),3.7 and r (sub e),2.l retrievals for highly inhomogeneous clouds. Our results suggest that the contrast in observed delta r (sub e)3.7-2.1 between cloud regimes is correlated with increases in both cloud r (sub e) and H-sigma. We also speculate that in some highly inhomogeneous drizzling clouds, vertical structure induced by drizzle and 3-D radiative effects might operate together to cause dramatic differences between r (sub e),3.7 and r (sub e),2.1 retrievals.

Platnick, Steven

Simulations of the effects of water vapor, cloud liquid water, and ice on AMSU moisture channel brightness temperatures

Radiative transfer simulations are performed to determine how water vapor and nonprecipitating cloud liquid water and ice particles within typical midlatitude atmospheres affect brightness temperatures T(sub B)'s of moisture sounding channels used in the Advanced Microwave Sounding Unit (AMSU) and AMSU-like instruments. The purpose is to promote a general understanding of passive top-of-atmosphere T(sub B)'s for window frequencies at 23.8, 89.0, and 157.0 GHz, and water vapor frequencies at 176.31, 180.31, and 182.31 GHz by documenting specific examples. This is accomplished through detailed analyses of T(sub B)'s for idealized atmospheres, mostly representing temperate conditions over land. Cloud effects are considered in terms of five basic properties: droplet size distribution, phase, liquid or ice water content, altitude, and thickness. Effects on T(sub B) of changing surface emissivity also are addressed. The brightness temperature contribution functions are presented as an aid to physically interpreting AMSU T(sub B)'s. Both liquid and ice clouds impact the T(sub B)'s in a variety of ways. The T(sub B)'s at 23.8 and 89 GHz are more strongly affected by altostratus liquid clouds than by cirrus clouds for equivalent water paths. In contrast, channels near 157 and 183 GHz are more strongly affected by ice clouds. Higher clouds have a greater impact on 157- and 183-GHz T(sub B)'s than do lower clouds. Clouds depress T(sub B)'s of the higher-frequency channels by suppressing, but not necessarily obscuring, radiance contributions from below. Thus, T(sub B)'s are less closely associated with cloud-top temperatures than are IR radiometric temperatures. Water vapor alone accounts for up to 89% of the total attenuation by a midtropospheric liquid cloud for channels near 183 GHz. The Rayleigh approximation is found to be adequate for typical droplet size distributions; however, Mie scattering effects from liquid droplets become important for droplet size distribution functions with modal radii greater than 20 micrometers near 157 and 183 GHz, and greater than 30-40 micrometers at 89 GHz. This is due mainly to the relatively small concentrations of droplets much larger than the mode radius. Orographic clouds and tropical cumuli have been observed to contain droplet size distributions with mode radii in the 30-40 micrometers range. Thus, as new instruments bridge the gap between microwave and infrared to frequencies even higher than 183 GHz, radiative transfer modelers are cautioned to explicitly address scattering characteristics of such clouds.

Muller, Bradley M.

Iterative Atmospheric Correction Scheme and the Polarization Color of Alpine Snow

Characterization of the Earth's surface is crucial to remote sensing, both to map geomorphological features and because subtracting this signal is essential during retrievals of the atmospheric constituents located between the surface and the sensor. Current operational algorithms model the surface total reflectance through a weighted linear combination of a few geometry-dependent kernels, each devised to describe a particular scattering mechanism. The information content of these measurements is overwhelmed by that of instruments with polarization capabilities: proposed models in this case are based on the Fresnel reflectance of an isotropic distribution of facets. Because of its remarkable lack of spectral contrast, the polarized reflectance of land surfaces in the shortwave infrared spectral region, where atmospheric scattering is minimal, can be used to model the surface also at shorter wavelengths, where aerosol retrievals are attempted based on well-established scattering theories. In radiative transfer simulations, straightforward separation of the surface and atmospheric contributions is not possible without approximations because of the coupling introduced by multiple reflections. Within a general inversion framework, the problem can be eliminated by linearizing the radiative transfer calculation, and making the Jacobian (i.e., the derivative expressing the sensitivity of the reflectance with respect to model parameters) available at output. We present a general methodology based on a Gauss-Newton iterative search, which automates this procedure and eliminates de facto the need of an ad hoc atmospheric correction. In this case study we analyze the color variations in the polarized reflectance measured by the NASA Goddard Institute of Space Studies Research Scanning Polarimeter during a survey of late-season snowfields in the High Sierra. This insofar unique dataset presents challenges linked to the rugged topography associated with the alpine environment and a likely high water content due to melting. The analysis benefits from ancillary information provided by the NASA Langley High Spectral Resolution Lidar deployed on the same aircraft. The results obtained from the iterative scheme are contrasted against the surface polarized reflectance obtained ignoring multiple reflections, via the simplistic subtraction of the atmospheric scattering contribution. Finally, the retrieved reflectance is modeled after the scattering properties of a dense collection of ice crystals at the surface. Confirming that the polarized reflectance of snow is spectrally flat would allow to extend the techniques already in use for polarimetric retrievals of aerosol properties over land to the large portion of snow-covered pixels plaguing orbital and suborbital observations.

Polarized BRDF

Ocean Color and Evidence of Chlorophyll Signature in the TOMS Minimum Reflectivity Data

Analysis of the TOMS minimum reflectivity data for 380 nm channel (R380) show regions of high reflectivity values (approx. 7 to 8%) over Sargasso Sea in the Northern Atlantic, anti-cyclonic region in the Southern Atlantic, and a large part of the ocean in the Southern Pacific, and low values (5 approx. 6 %) over the rest of the open ocean. Through radiative transfer simulations we show that these features are highly correlated with the distribution of chlorophyll in the ocean. Theoretical minimum reflectivity values derived with the help of CZCS chlorophyll concentration data as input into a vector ocean-atmosphere radiative transfer code developed by Ahmad and Fraser show very good agreement with TOMS minimum reflectivity data for the winter season of year 1980. For the summer season of year 1980, good qualitative agreement is observed in the equatorial and northern hemisphere but not as good in the southern hemisphere. Also, for cloud-free conditions, we find a very strong correlation between R340 minus R380 values and the chlorophyll concentration in the ocean. Results on the possible effects of absorbing and non-absorbing aerosols on the TOMS minimum reflectivity will also be presented. The results also imply that ocean color will affect the aerosol retrieval over oceans unless corrected.

Ahmad, Z.

Determination of Ice Cloud Models Using MODIS and MISR Data

Representation of ice clouds in radiative transfer simulations is subject to uncertainties associated with the shapes and sizes of ice crystals within cirrus clouds. In this study, we examined several ice cloud models consisting of smooth, roughened, homogeneous and inhomogeneous hexagonal ice crystals with various aspect ratios. The sensitivity of the bulk scattering properties and solar reflectances of cirrus clouds to specific ice cloud models is investigated using the improved geometric optics method (IGOM) and the discrete ordinates radiative transfer (DISORT) model. The ice crystal habit fractions in the ice cloud model may significantly affect the simulations of cloud reflectances. A new algorithm was developed to help determine an appropriate ice cloud model for application to the satellite-based retrieval of ice cloud properties. The ice cloud particle size retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) data, collocated with Multi-angle Imaging Spectroradiometer (MISR) observations, is used to infer the optical thicknesses of ice clouds for nine MISR viewing angles. The relative differences between view-dependent cloud optical thickness and the averaged value over the nine MISR viewing angles can vary from −0.5 to 0.5 and are used to evaluate the ice cloud models. In the case for 2 July 2009, the ice cloud model with mixed ice crystal habits is the best fit to the observations (the root mean square (RMS) error of cloud optical thickness reaches 0.365). This ice cloud model also produces consistent cloud property retrievals for the nine MISR viewing configurations within the measurement uncertainties.

MODIS

The influence of charge exchange on the velocity distribution of hydrogen in the Venus exosphere

The simulation of the exosphere of Venus by a Monte Carlo technique has been extended to provide velocity distribution profiles that can be used with radiative transfer simulation to model the Lyman-alpha emission seen by spacecraft. The line profiles show a narrow core due to the exobase thermal source of hydrogen superimposed on the broad profile of hot hydrogen from charge exchange with hot ions in the nighttime ionosphere. Nightside radial profiles show long tails of upward flowing, escaping atoms. There are some downward and lateral superescape signatures, because the charge exchange source extends well above the exobase. At higher altitudes, flattening and even shallow central valleys appear in the transverse profiles. The planetary average escape rate due to charge exchange was found to be 2.8 x 10 to the 7th/sq cm per s, which is several times larger than escape rates for other candidate mechanisms.

Hodges, Richard R., Jr.

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

stars: atmospheres

The Influence of Thermodynamic Phase on the Retrieval of Mixed-Phase Cloud Microphysical and Optical Properties in the Visible and Near Infrared Region

Cloud microphysical and optical properties are inferred from the bidirectional reflectances simulated for a single-layered cloud consisting of an external mixture of ice particles and liquid droplets. The reflectances are calculated with a rigorous discrete ordinates radiative transfer model and are functions of the cloud effective particle size, the cloud optical thickness, and the values of the ice fraction in the cloud (i.e., the ratio of ice water content to total water content). In the present light scattering and radiative transfer simulations, the ice fraction is assumed to be vertically homogeneous; the habit (shape) percentage as a function of ice particle size is consistent with that used for the Moderate Resolution Imaging Spectroradiometer (MODIS) operational (Collection 4 and earlier) cloud products; and the surface is assumed to be Lambertian with an albedo of 0.03. Furthermore, error analyses pertaining to the inference of the effective particle sizes and optical thicknesses of mixed-phase clouds are performed. Errors are calculated with respect to the assumption of a cloud containing solely liquid or ice phase particles. The analyses suggest that the effective particle size inferred for a mixed-phase cloud can be underestimated (or overestimated) if pure liquid phase (or pure ice phase) is assumed for the cloud, whereas the corresponding cloud optical thickness can be overestimated (or underestimated).

Lee, Joonsuk

Estimation of the radiation budget during MOSAiC based on ground-based and satellite remote sensing observations

An accurate representation of the radiation budget is essential for investigating the impact of clouds on the climate system, especially in the Arctic, an environment highly sensitive to complex and rapid environmental changes. In this study, we analyse a unique dataset of observations from the central Arctic made during the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition in conjunction with state-of-the-art satellite products from CERES (Clouds and the Earth's Radiant Energy System) to investigate the radiative effect of clouds and radiative closure at the surface and the top of the atmosphere (TOA). We perform a series of radiative transfer simulations using derived cloud macro- and microphysical properties as inputs to the simulations for the entire MOSAiC period, comparing our results to collocated satellite products and ice-floe observations. The radiative closure biases were generally within the instrumental uncertainty, indicating that the simulations are sufficiently accurate to reproduce the radiation budget during MOSAiC. Comparisons of the simulated radiation budget relative to CERES show similar values in the terrestrial flux but relatively large differences in the solar flux, which are attributed to a lower surface albedo and a possible underestimation of atmospheric opacity by CERES. While the simulation results were consistent with the observations, more detailed analyses reveal an overestimation of simulated cloud opacity for cases involving geometrically thick ice clouds. In the annual mean, we found that, during the MOSAiC expedition, the presence of clouds leads to a loss of 5.2 W m −2 of the atmosphere–surface system to space, while the surface gains 25.0 W m −2 and the atmosphere is cooled by 30.2 W m −2 .

54 ENVIRONMENTAL SCIENCES

Optically thin cirrus clouds - Remote sensing, and geophysical significance

The region of the IR spectrum that is ideally suited for detecting optically thin cirrus clouds is in the window between 10 and 13 microns. Here relatively weak absorption due to the water vapor lines and continuum is present and hence the extinction characteristic of the cloud particles is readily discernible. In order to demonstrate these properties, two IR spectra are presented, one with clear skies and one with an optically thin cirrus. As a result of the cloud particle extinction, an appreciable increase in the brightness temperature from 10 to 13 microns is observed. This decrease is found to be nearly linear in the case of the tropical thin cirrus, which is also geometrically thin. On the basis of radiative transfer simulations, it is inferred that the cloud particle size that can produce this spectral character has an effective diameter of about 12 microns, which is comparable to the wavelength of the radiation.

Prabhakara, C.

An Algorithm for the Retrieval of 30-m Snow-Free Albedo from Landsat Surface Reflectance and MODIS BRDF

We present a new methodology to generate 30-m resolution land surface albedo using Landsat surface reflectance and anisotropy information from concurrent MODIS 500-m observations. Albedo information at fine spatial resolution is particularly useful for quantifying climate impacts associated with land use change and ecosystem disturbance. The derived white-sky and black-sky spectral albedos maybe used to estimate actual spectral albedos by taking into account the proportion of direct and diffuse solar radiation arriving at the ground. A further spectral-to-broadband conversion based on extensive radiative transfer simulations is applied to produce the broadband albedos at visible, near infrared, and shortwave regimes. The accuracy of this approach has been evaluated using 270 Landsat scenes covering six field stations supported by the SURFace RADiation Budget Network (SURFRAD) and Atmospheric Radiation Measurement Southern Great Plains (ARM/SGP) network. Comparison with field measurements shows that Landsat 30-m snow-free shortwave albedos from all seasons generally achieve an absolute accuracy of +/-0.02 - 0.05 for these validation sites during available clear days in 2003-2005,with a root mean square error less than 0.03 and a bias less than 0.02. This level of accuracy has been regarded as sufficient for driving global and regional climate models. The Landsat-based retrievals have also been compared to the operational 16-day MODIS albedo produced every 8-days from MODIS on Terra and Aqua (MCD43A). The Landsat albedo provides more detailed landscape texture, and achieves better agreement (correlation and dynamic range) with in-situ data at the validation stations, particularly when the stations include a heterogeneous mix of surface covers.

Shuai, Yanmin

Quantifying the Relative Impact of Model Microphysics Parameterizations and Scattering Models in Simulating Synthetic Radar and Microwave Radiometer Observations

Output from numerical weather models is often used to simulate observations from remote sensing instruments, for purposes ranging from data assimilation, synthetic retrievals of geophysical quantities, and optimization of observing systems. However, when hydrometeors are present, the level of detail provided by the weather model is generally insufficient to fully constrain the input to the radiative transfer model (RTM), and further assumptions must be made by the RTM user in order to produce synthetic observations. Using a hierarchy of models including cloud-resolving, double-moment, bin microphysical, and ice-habit predicting models, along with scattering properties from the OpenSSP, Atmospheric Radiative Transfer Simulator (ARTS) databases, as well as relatively simple geometries (e.g., cylindrical plates and columns), we demonstrate the spread in synthetic observation output and the extent to which it is reduced when microphysics is more strongly constrained by the model. As an intermediate step, an error budget for the RTM simulations was derived and from that we developed and will describe best practices for observation simulation (e.g., optimal number of hydrometeor size bins, truncation of the particle size distribution, angular resolution of scattering phase function). Some statistical comparisons with observations will also be presented.

Munchak, S. Joseph

Case Studies of the Vertical Structure of the Direct Shortwave Aerosol Radiative Forcing During TARFOX

The vertical structure of aerosol-induced radiative flux changes in the Earth's troposphere affects local heating rates and thereby convective processes, the formation and lifetime of clouds, and hence the distribution of chemical constituents. We present observationally based estimates of the vertical structure of direct shortwave aerosol radiative forcing for two case studies from the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX) which took place on the U.S. east coast in July 1996. The aerosol radiative forcings are computed using the Fu-Liou broadband radiative transfer model. The aerosol optical properties used in the radiative transfer simulations are calculated from independent vertically resolved estimates of the complex aerosol indices of refraction in two to three distinct vertical layers, using profiles of in situ particle size distributions measured aboard the University of Washington research aircraft. Aerosol single-scattering albedos at 450 nm thus determined range from 0.9 to 0.985, while the asymmetry factor varies from 0.6 to 0.8. The instantaneous shortwave aerosol radiative forcings derived from the optical properties of the aerosols are of the order of -36 Wm(exp -2) at the top of the atmosphere and about -56 Wm(exp -2) at the surface for both case studies.

Redemann, J.

Case Studies of the Vertical Structure of the Direct Shortwave Aerosol Radiative Forcing During TARFOX

The vertical structure of aerosol-induced radiative flux changes in the Earth's troposphere affects local heating rates and thereby convective processes, the formation and lifetime of clouds, and hence the distribution of chemical constituents. We present observationally-based estimates of the vertical structure of direct shortwave aerosol radiative forcing for two case studies from the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX) which took place on the US East coast in July 1996. The aerosol radiative forcings are computed using the Fu-Liou broadband radiative transfer model. The aerosol optical properties used in the radiative transfer simulations are calculated from independent vertically-resolved estimates of the complex aerosol indices of refraction in two to three distinct vertical layers, using profiles of in situ particle size distributions measured aboard the University of Washington research aircraft. Aerosol single-scattering albedos at 450 nm thus determined range from 0.9 to 0.985, while the asymmetry factor varies from 0.6 to 0.8. The instantaneous shortwave aerosol radiative forcings derived from the optical profiles of the aerosols are of the order of -36 W/sq m at the top of the atmosphere and about -56 W/sq m at the surface for both case studies.

Redemann, J.

Dust and gas absorption in the High Mass X-ray Binary IGR J16318−4848

Context. With an absorption column density on the order of 10(exp 24)per sq. cm, IGR J16318-4848 is one of the most extreme cases of a highly obscured High Mass X-ray Binary. Besides the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. Aims. We aim at a self-consistent description of the continuum absorption, the absorption edges and the fluorescent lines to constrain properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescent emission from dust grains. Methods. We use XMM-Newton and NuSTAR spectra to first constrain empirically the incident continuum and fluorescence lines. Next, we use XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In a third step, we model the absorption and fluorescence in a dusty olivine absorber and employ both, a simple analytical model for the fluorescent line emission and Monte Carlo radiative transfer, spectral shapes and line fluxes that are very close to the data are generated. Results. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is well able to produce the observed continuum absorption and edges. Conclusions. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318􀀀4848. In particular, for low optical depth of individual grains the dust will contribute significantly to the fluorescent emission, which is the case for typical densities and grain size distribution models.

Ralf Ballhausen

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante