We present a novel method to implement time-delayed propagation of radiation fields in cosmo- logical radiative transfer simulations. Time-delayed propagation of radiation fields requires construction of retarded-time fields by tracking the location and lifetime of radiation sources along the corresponding light-cones. Cosmological radiative transfer simulations have, until now, ignored this “light-cone effect” or implemented ray-tracing methods that are computationally demanding. We show that radiative trans- fer calculation of the time-delayed fields can be easily achieved in numerical simulations when periodic boundary conditions are used, by calculating the time-discretized retarded-time Green’s function using the Fast Fourier Transform (FFT) method and convolving it with the source distribution. We also present a direct application of this method to the long-range radiation field of Lyman-Werner band photons, which is important in the high-redshift astrophysics with first stars.
Monte-Carlo radiative models of the diffuse Galactic light (DGL) in our Galaxy are calcu- lated using the dust radiative transfer code MoCafe, which is three-dimensional and takes full account of multiple scattering. The code is recently updated to use a fast voxel traversal algorithm, which has dramatically increased the computing speed. The radiative transfer models are calculated with the gen- erally accepted dust scale-height of 0.1 kpc. The stellar scale-heights are assumed to be 0.1 or 0.35 kpc, appropriate for far-ultraviolet (FUV) and optical wavelengths, respectively. The face-on optical depth, measured perpendicular to the Galactic plane, is also varied from 0.2 to 0.6, suitable to the optical to FUV wavelengths, respectively. We find that the DGL at high Galactic latitudes is mostly due to backward or large-angle scattering of starlight originating from the local stars within a radial distance of r < 0.5 kpc from the Earth. On the other hand, the DGL measured in the Galactic plane is mostly due to stars at a distance range that corresponds to an optical depth of ≈ 1 measured from the Earth. Therefore, the low-latitude DGL at the FUV wavelength band would be mostly caused by the stars located at a distance of r . 0.5 kpc and the optical DGL near the Galactic plane mainly originates from stars within a distance range of 1 . r . 2 kpc. We also calculate the radiative transfer models in a clumpy two-phase medium. The clumpy two-phase models provide lower intensities at high Galactic latitudes compared to the uniform density models, because of the lower effective optical depth in clumpy media. However, no significant difference in the intensity at the Galactic plane is found.
RTTOV와 CRTM은 복사관측자료에 대한 관측연산자로 수치예보에 활용되고 있는 빠른 속도의 복사전달모델이다. 본 연구에서는 두 모델의 기본구조 및 입력자료를 비교했다. 또한, 다양한 파장대를 가진 AMSU-A 마이크로파 센서에 대해 구름에 대한 정보를 포함할 때와 포함하지 않을 때 두 모델로부터 계산된 밝기온도와 관측된 밝기온도를 해양에 대해 비교했다. AMSU-A의 탐측채널(5-14)에 대해서는 두 모델로부터 계산된 밝기온도 값에 큰 차이가 존재하지 않았으나, 대기의 창 채널 및 지표근처의 탐측채널에서는 RTTOV로부터 계산된 밝기온도 값이 관측과 더 가까워 CRTM에 비해 상대적으로 작은 초기추정오차를 보였다. 한편 UM으로부터 제공된 구름물과 얼음의 정보를 추가적으로 활용하였을 때 두 모델로부터 계산된 밝기온도와 관측된 밝기온도의 차이가 감소함을 확인할 수 있었고, 특히 CRTM의 31.4 GHz와 89 GHz 채널에서 모의된 밝기온도와 관측된 밝기온도의 차이가 크게 감소했다.
본 연구에서는 에어러솔 대기 상태에서 국내 COMS 연구자들이 사용하는 세 개의 단파 복사전달모델에서 산출된 복사속 성분을 비교분석하였으며, 대기 복사 수지에서 에어러솔 역할도 분석하였다. 국내 모델들의 평가를 위하여, 15개 모델값을 평균한 Halthore et al.(2005) 결과를 기준값으로 사용하였다. 두 종류 에어러솔 농도(AOT=0.08, 0.24)에서 조사된 열대 또는 한대 대기의 에어러솔 강제력은 지표에서의 하향 일사 및 상향 산란, 대기 상부의 상향 산란, 그리고 대기 흡수도의 복사 성분들에 있어서 국외 연구에 비하여 국내 결과들에서 체계적으로 약하게 나타났다. 에어러솔 강제력은 지표에서의 하향 일사에 대하여 -10∼-40Wm-2 이었으며, 지표 및 대기 상부의 상향 산란의 경우에 상대적으로 큰 태양천정각과 고농도 에어러솔 상태에서 컸다. 두 종류 에어러솔 조건에서 지표에서의 하향 및 상향 산란값들은 대기 종류보다는 태양천정각에 더 민감하였다. 하향 산란은 상대적으로 작은 태양천정각과 고농도 에어러솔 조건에서 컸다. 에어러솔 농도 증가는 하향 산란 증가에도 불구하고 하향 직달 일사의 감소가 이를 초과함으로써 지표 냉각을 유도하는 것으로 조사되었다. 동일한 에어러솔 농도 및 태양천정각 조건에서 직달일사 소산은 열대 대기에서는 주로 수증기, 그리고 한대 대기에서는 수증기뿐만 아니라 오존에도 기인하는 것으로 조사되었다. 특히 열대 대기에서는 수증기의 역할이 오존에 비하여 3∼4배 컸다. 저농도 및 고농도 에어러솔 대기에서 대기 흡수도는 국내외 연구간에 ±10% 내에서 일치하였다.
본 연구에서는 레일리(순수 기체) 대기 조건 하에서 국내 COMS 연구자들이 사용하는 여덟 개 단파 복사전달모델에서 산출된 네 종류 복사속(flux) 성분을 상호비교함으로써 상대 오차를 조사하였다. 이들 복사속 성분은 지표에서의 직달 일사, 하향 산란, 상향 산란, 그리고 대기 상부에서의 상향 산란이다. 또한 국내 모델의 평가를 위하여, 15개모델을 평균한 Halthore et al.(2005) 결과(예, H15)를 기준값으로 사용하였다. 동일한 태양천정각에서 모델 간의 불일치는 열대 대기에서 수증기에 기인하였고, 한대 대기에서는 오존에 기인하였다. STREAMER를 제외한 국내 7개 모델의 지상에서의 하향 직달일사값은 H15에 대하여 ±4%내에서 일치하였다. 이러한 상대 오차는 태양천정각이 커질 때 증가하였으며, Halthore et al.(2005)에서의 ±3%와 근접하였다. 네 종류 복사속 분석에서 SBDART 모델이 6S 모델에 비하여 전반적으로 우수하였으나, 근적외 파장역에서는 서로 비교할만하였다. 네개 기관의 연구자들이 같은 SBDART 모델에서 산출한 지표에서의 하향 직달일사값 간에도 12.1Wm-2(1.4%)의 불일치가 존재하였다. 불일치의 원인은 부분적으로 복사속 적분에 있어서 서로 다르게 설정된 파장 분해능에도 있었다. 본 연구는 단파 영역에서 최적 모델을 선정하는 데 도움을 줄 수 있다.
Ever since the identification of 6830 and 7088 features as the Raman scattered O VI 1032, 1038 resonance doublets in symbiotic stars by Schmid (1989), Raman scattering by atomic hydrogen has been a very unique tool to investigate the mass transfer processes in symbiotic stars. Discovery of Raman scattered He II in young planetary nebulae (NGC 7027, NGC 6302, IC 5117) allow one to expect that Raman scattering can be an extremely useful tool to look into the mass loss processes in these objects. Because hydrogen is a single electron atom, their wavefunctions are known in closed form, so that exact calculations of cross sections are feasible. In this paper, I review some basic properties of Raman scattered features and present detailed and explicit matrix elements for computation of the scattering cross section of radiation with atomic hydrogen. Some astrophysical objects for which Raman scattering may be observationally pertinent are briefly mentioned.
We have written a code called QDM_sca, which numerically solves the problem of radiative transfer in an anisotropically scattering, spherical atmosphere. First we formulate the problem as a second order differential equation of a quasi-diffusion type. We then apply a three-point finite differencing to the resulting differential equation and transform it to a tri-diagonal system of simultaneous linear equations. After boundary conditions are implemented in the tri-diagonal system, the QDM_sca radiative code fixes the field of specific intensity at every point in the atmosphere. As an application example, we used the code to calculate the brightness of atmospheric diffuse light(ADL) as a function of zenith distance, which plays a pivotal role in reducing the zodiacal light brightness from night sky observations. On the basis of this ADL calculation, frequent uses of effective extinction optical depth have been fully justified in correcting the atmospheric extinction for such extended sources as zodiacal light, integrated starlight and diffuse galactic light. The code will be available on request.
Theoretical models of radiative transfer are developed to simulate the 85 GHz brightness temperature (T85) observed by the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer as a function of rain rate. These simulations are performed separately over regions of the convective and stratiform rain. TRMM Precipitation Radar (PR) observations are utilized to construct vertical profiles of hydrometeors in the regions. For a given rain rate, the extinction in 85 GHz due to hydrometeors above the freezing level is found to be relatively weak in the convective regions compared to that in the stratiform. The hydrometeor profile above the freezing level responsible for the weak extinction in convective regions is inferred from theoretical considerations to contain two layers: 1) a mixed (or mixed-phase) layer of 2 km thickness with mixed-phase particles, liquid drops and graupel above the freezing level, and 2) a layer of graupel extending from the top of the mixed layer to the cloud top. Strong extinction in the stratiform regions is inferred to result from slowly-falling, low-density ice aggregates (snow) above the freezing level. These theoretical results are consistent with the T85 measured by TMI, and with the rain rate deduced from PR for the convective and stratiform rain regions. On the basis of this study, the accuracy of the rain rate sensed by TMI is inferred to depend critically on the specification of the convective or stratiform nature of the rain.
We have developed a Monte Carlo code, which solves the problem of radiative transfer in anisotropically scattering atmosphere. The radiative code is flexible in handlings of the system geometry, the distribution of scattering particles, and the source-particle geometry. This code treats the case of highly forward throwing scattering. As performance tests, we have compared the result of Monte Carlo calculations with that of Quasi-Diffusion method for a spherically symmetric cloud model.
We calculated dose rate using radiative transfer equations to consider radiative processes distinctly. The dose rate at Pohang(36°02'N, 129°23'E) was calculated using measured ozone and meteorological data and two-stream approximations(quadrature, Eddington, delta Eddington, PIFM(practical improved flux method), discrete ordinate, delta discrete ordinate) are used in solving equation. The purpose of this study is to determine the most compatible radiative transfer approximation for simulating the radiative and photochemical processes of atmosphere through comparision between calculated and measured values. Dose rate of the biologically effective irradiance in the region 0.28-0.32μm showed the highest value when quadrature and Eddington was used and lower value on condition that delta scaling was applied. Correlation coefficient between dose rate at surface using radiation transfer equation and measured UV-B at Pohang was 0.78, 0.79 and 0.81 when detla Eddington, PIFM and delta discrete ordinate were used. Also, in case of above approximations were used, MBE(Mean Bias Error) was within -0.3MED/30min and RMBE(Relative Mean Bias Error) was below 10% between 1200 LST and 1400 LST. Approximations which are compatible in estimating radiative process are detla Eddington, PIFM and delta discrete ordinate. Especially, in case that radiative process is considered more detail, delta discrete ordinate increased the number of stream is proper.