Vector radiative transfer numerical model of coupled ocean-atmosphere system using matrix-operator method

A vector radiative transfer numerical model of the coupled ocean-atmosphere system is developed based on the matrix-operator method, which is named PCOART. Using the Fourier analysis, the vector radiative transfer equation (VRTE) is separated into a set of equations depending only on the observa-tion zenith angle. Using the Gaussian-Quadrature method, VRTE is finally transferred into the matrix equation solved by the adding-doubling method. According to the reflective and refractive properties of the ocean-atmosphere interface, the vector radiative transfer numerical model of the ocean and at-mosphere is coupled in PCOART. Compared with the exact Rayleigh scattering look-up tables of MODIS (Moderate-resolution Imaging Spectroradiometer), it is shown that PCOART is an exactly numerical model, and the processing methods of the multi-scattering and polarization are correct. Also, validated with the standard problems of the radiative transfer in water, it is shown that PCOART can be used to calculate the underwater radiative transfer problems. Therefore, PCOART is a useful tool for exactly calculating the vector radiative transfer of the coupled ocean-atmosphere system, which can be used to study the polarization properties of the radiance in the whole ocean-atmosphere system and the remote sensing of the atmosphere and ocean.     

Passive microwave remote sensing contribution to hydrological variables

Physical principles governing passive microwave remote sensing of hydrological variables are outlined and illustrated by actual observations by ground-based, air-borne and space-borne microwave radiometers operating at different frequencies. Specific hydrological variables addressed in this paper are soil moisture, seasonal inundation of rivers and swamps, vegetation, snow, and rainfall. Passive remote sensing provides measurements of electromagnetic radiation emitted by the land-atmosphere system, which can be related more directly to the radiative characteristics of the system than to physical or physiological characteristics. Estimation of hydrological variables from microwave observations necessarily involves models relating the radiative to the physical characteristics, and in general more than one physical characteristics determine the microwave observations. This non-uniqueness in the relationship between microwave observations to a particular hydrological variable leads to uncertainties in the estimation of the variable. Notwithstanding this limitation, the principles and the examples given in this paper illustrate the value of passive microwave observations to regional and global hydrology at a temporal resolution of days aggregated to a decade.     

On the determination of the atmospheric ozone profile for ground based microwave measurements

It is shown how to retrieve the atmospheric ozone profile by means of ground based microwave measurements of the radiative intensity. Chahine's iteration method is used. The method is tested by a numerical experiment. The retrieval rms about the mean error is approximately 0.4 ppm. The required measurement accuracy for the brightness temperatures is ±0.01 K.     

星载微波SSM/I多时相辐射观测的特征指数监测与评估2003年7月中国淮河流域汛情

由辐射传输方程推导星载被动微波遥感SSM/I辐射观测的极化指数PI,凸现地表土壤湿度的变化与距平异常ΔnPI,用特征性指数PI、ΔnPI和散射指数SI来监测大尺度陆地地面的降雨和地区土壤湿度、洪涝汛情的变化.用中国淮河流域7年同月SSM/I观测数据的平均<PI>,洪涝期间的PI,ΔnPI和SI监测该流域2003年7月的洪涝汛情.提出距平异常值指数K分布,定量评估了洪涝汛情的范围和程度.     

Modeling Full Seismogram Envelopes Using Radiative Transfer Theory with Born Scattering Coefficients

The equation of radiative transfer is used to model the transport of seismic energy in 2-D and 3-D acoustic random media. Monte-Carlo solutions of this equation using non-isotropic Born scattering coefficients are compared to three analytical solutions: Markov approximation, radiative transfer theory with isotropic scattering coefficients, and diffusion approximation. Additionally, we compare to finite differences solutions of the full wave equation in 2-D. We find a good correspondence of radiative transfer theory to Markov approximation for the case of multiple forward scattering. The comparison to radiative transfer theory with isotropic scattering coefficients, a model frequently used in data analysis, demonstrates that in the case of forward scattering the isotropic scattering model is not better than a diffusion approach. To compare radiative transfer theory with non-isotropic scattering coefficients to finite differences solutions of the full wave equation, the finite source duration and the bandpass filter process as well as the normalization of absolute amplitudes are explicitely taken into account. We find a good coincidence of both theories for scattering parameters, which are realistic for usual Earth crust. The theory correctly describes the unscattered direct wavefront, the envelope broadening caused by multiple forward scattering, as well as the late coda caused by multiple wide angle scattering. For strong scattering, which can be expected for very heterogeneous media such as strato volcanoes, the solutions of radiative transfer differ from the more complete solutions of the full wave equation.     

利用OI 135.6 nm夜气辉辐射探测电离层峰值电子密度及电子总含量的研究

测量夜间135.6 nm大气气辉辐射强度是目前有效的电离层探测方式之一,我国即将在风云三号卫星上搭载仪器,利用该波段夜气辉辐射测量来反演电子总含量(本文所指电子总含量表示卫星高度以下大气柱的电子含量)及峰值电子密度,因此非常有必要开展相关的气辉发光模型及反演研究.在介绍氧原子135.6 nm波段夜气辉激发机制基础上,考虑辐射在传输过程中受到大气氧原子的散射及氧气分子的吸收,采用迭代法求解包含多次散射及大气吸收衰减的辐射传输方程,得到该波段的体发射率,最终通过考虑包含辐射传输的路径积分计算得到135.6 nm气辉辐射强度值.对结果的分析表明:该气辉模型能较好地描述体发射率随高度的分布特征,计算得到的135.6 nm夜气辉辐射强度在不同时空及太阳活动的分布与相应条件下峰值电子密度(N_mF₂)及电子总含量(TEC)的分布基本一致.相同的时空及太阳活动输入条件下,模式计算的135.6 nm夜气辉辐射强度与国外同类模式结果的值平均偏差约为3%.文中最后介绍了通过135.6 nm夜气辉的辐射强度探测来反演电离层峰值电子密度N_mF₂及电子总含量TEC的反演方法.     

GPU加速的水体辐射传输Monte Carlo模拟模型研究

水体辐射传输方程是复杂的微积分方程,只能利用数值方法求解,如Monte Carlo光线追踪法、不变嵌入法、离散坐标法等,其中,Monte Carlo方法是目前解决水体水下光场三维问题的唯一有效方法.根据辐射传输理论,开发了水下光场的Monte Carlo模拟模型,主要包含大气、水-气界面、层化水体和水底边界4个模块.实现了模拟任意太阳角度、不同水体固有光学属性和任意深度条件下,考虑大气、粗糙水面和水底边界的水下光场,能够获取辐亮度、辐照度等辐射量的空间分布.该模型暂不考虑Raman散射、偏振、内部光源的影响.实现了GPU加速水下光场Monte Carlo模拟,并用Mobley等提出的海洋光学标准问题中的问题1~6进行验证.在两种计算环境下,通过对不同边界条件下的CPU、GPU运行时间及加速比的对比,发现GPU计算可以达到几百至上千倍的加速比.     

The Role of Radiative Transfer in Maintaining the Indian Summer Monsoon Circulation

—The radiative-convective feedback and land-sea thermal forcing play significant roles in maintenance of the summer monsoon circulation over the Indian sub-continent. In this study, the role of radiative transfer in maintaining the monsoon circulation is examined with numerical sensitivity experiments. For this purpose, a sixteen layer primitive equation limited area model is used to perform numerical simulations with and without atmospheric radiative transfer processes parameterized in the model. The initial values and boundary conditions for the numerical integrations of the model are derived from operational analyses of the ECMWF, UK. The results show that the radiative transfer is essential in maintaining the intensity of the low level Somali Jet as well as the upper level Tropical Easterly Jet (TEJ) over the Indian sub-continent and adjoining seas. The meridional circulation over the region is also well simulated. As a result, enough moisture transports from the warm equatorial region to simulate more realistic orographic precipitation in the windward side of the mountains along the West coast of India. Without radiative transfer processes in the model atmosphere the simulated monsoon circulation weakens, moisture transport decreases and the precipitation lessens.     

A model for accurately calculating hyper-spectral,middle-shortwave infrared radiative transfer for remote sensing

Although the calculation of radiative transfer in the middle-shortwave infrared band is important in the field of optical remote sensing, studies in this area of research are rare in China. Both solar reflection and atmospheric emission should be considered when calculating radiative transfer in the middle-shortwave infrared band. This paper presents a new radiative transfer model based on the doubling and adding method. The new model uses approximate calculations of direct solar reflection,multiple scattering, and thermal emissions for a finitely thin atmospheric layer and considers both the solar and thermal sources of radiation. To verify its accuracy, the calculation results produced by the model for four typical scenarios(single layer at night,multi-layer aerosols, double-layer with ice and water clouds, and multi-layer with clouds and aerosols) were compared with those of the DISORT model. With the exception of a few channels, the absolute deviation between the two models was less than2×10~(-6) K. For the same calculation, the computation speed of the new model was approximately two to three times faster than that of the DISORT model. Sensitivity studies were performed to evaluate the error resulting from using simplified calculation methods in the new model. The results obtained in this study indicated that atmospheric thermal emission made a significant contribution to the measured radiance in the strong-absorption band(2230–2400 cm~(-1)), whereas solar radiation could be neglected in this region. However, neglecting solar radiation in the window region(2400–2580 cm~(-1)) introduced error on the order of dozens of K. Employing the average-layer temperature method simplified the calculation of thermal radiation but caused a larger error in the strong-absorption band than in the window region. In the doubling and adding method, the calculation error decreased as the value used for minimum optical thickness decreased. Under the condition of satisfying the requirement of calculation precision, we can consider using the layer-average temperature radiation method and selecting a relative larger minimum optical thickness value to improve the calculation efficiency. The new radiative calculation model proposed herein can be used in the simulation, inversion, and assimilation of middle-shortwave infrared measurements by hyper-spectral satellite instruments.     

The distribution of infrared radiative flux densities and their divergence along the symmetry line of an idealized valley

Summary This investigation concerns itself with the determination of infrared radiative flux divergence along the symmetry line of a cone. Such a cone, imbedded in the plane earth, is thought to idealize a large valley. The radiative transfer equations were set up in such a manner as to deal with the special geometry involved. It was found that radiative temperature changes within the valley are a strong function of the steepness of the valley walls.     

Radiative properties of cirrus clouds in NOAA 4 VTPR channels: Some explorations of cloud scenes from satellites

A one-dimensional spectral infrared radiative transfer model has been developed for atmospheres containing cirrus clouds and absorbing gases above, below and within the cloud. The transfer model takes into consideration the inhomogeneity of the cloudy atmosphere, the gaseous absorption in scattering cloud layers and the wavenumber dependence of radiative transfer. In addition, the cirrus cloud is further divided into a number of sub-layers to account for the non-isothermal and inhomogeneous cloud characteristics. Single-scattering properties for ice crystals are calculated assuming ice cylinders 200 and 60 m in lenght and width, respectively, randomly oriented in a horizontal plane. The spectral infrared transfer program is applied to VTPR channels of the NOAA 4 satellite to simulate upward radiances in cirrus cloud conditions.Comparisons between satellite observed and theoretically simulated upward radiances are carried out for selected cirrus cloud cases. Incorporating atmospheric profiles obtained from radiosonde and the observed cloud information into the spectral transfer program, we show a systematic agreement between observed and computed upward radiances. Systematic reduction patterns of the upward radiance caused by the increase of the cloud ice content are clearly demonstrated for VTPR channels employing tropical and midlatitude atmospheric profiles. Having the quantitative relationships between upward radiances and ice contents, procedures are described for the inference of the cloud ice content and cloud amount. The proposed method has been successfully applied to the three cirrus cloud cases.     

Lee生物光学模型在不同水体组分特性下的适用性

辐射传输模型和生物光学模型均可用于模拟水体遥感反射率.前者模拟精度高,但计算复杂.不利于水质参数的反演;后者简便易反演,但在浑浊水体中的模拟精度还有待进一步检验.本文通过设计大量不同组分浓度组成的水体,以辐射传输模型(即Hydrolisht模型)模拟结果为真值,对生物光学模型(即Lee模型)模拟二类水体遥感反射率的精度...     

COSMIC数据验证AMSU平流层低层观测的初步分析结果

基于Global Positioning System (GPS)掩星数据在平流层具有较高准确性、稳定性的优势,本文尝试用新一代GPS掩星观测——the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC)资料验证不同卫星平台上先进的微波探测仪（AMSU）的平流层观测结果.通过COSMIC大气温度廓线与AMSU辐射传输模式结合,得到模拟亮温,然后与AMSU平流层观测进行匹配比较.分析表明GPS掩星数据能够作为一个相对独立的参量检验NOAA15、16、18卫星平台内部的偏差.通过一年数据的比较验证,初步显示不同卫星平台的AMSU观测亮温在平流层低层都偏低,并且NOAA18平台的亮温偏低程度明显大于NOAA15、16.AMSU亮温偏差在极地冬季较为显著,尤其南极地区NOAA18的偏差幅度达到1.8 K.结合24小时内AMSU观测亮温偏差变化及其样本分布特征,可以看到明显的太阳辐射差异可能是导致AMSU观测亮温在极地偏差显著的主要原因.     

A soil moisture assimilation scheme based on the ensemble Kalman filter using microwave brightness temperature

This study presents a soil moisture assimilation scheme, which could assimilate microwave brightness temperature directly, based on the ensemble Kalman filter and the shuffled complex evolution method (SCE-UA). It uses the soil water model of the land surface model CLM3.0 as the forecast operator, and a radiative transfer model (RTM) as the observation operator in the assimilation system. The assimilation scheme is implemented in two phases: the parameter calibration phase and the pure soil moisture assimilation phase. The vegetation optical thickness and surface roughness parameters in the RTM are calibrated by SCE-UA method and the optimal parameters are used as the final model parameters of the observation operator in the assimilation phase. The ideal experiments with synthetic data indicate that this scheme could significantly improve the simulation of soil moisture at the surface layer. Furthermore, the estimation of soil moisture in the deeper layers could also be improved to a certain extent. The real assimilation experiments with AMSR-E brightness temperature at 10.65 GHz (vertical polarization) show that the root mean square error (RMSE) of soil moisture in the top layer (0–10 cm) by assimilation is 0.03355 m³ · m⁻³, which is reduced by 33.6% compared with that by simulation (0.05052 m³ · m⁻³). The mean RMSE by assimilation for the deeper layers (10–50 cm) is also reduced by 20.9%. All these experiments demonstrate the reasonability of the assimilation scheme developed in this study.     

基于辐射传递模拟及人工神经网络技术的二类水体光学组分的反演

介绍了一种基于辐射传递模拟和人工神经网络技术(ANN)的二类水体水色要素(CHL,SPM,CDOM)的反演算法.在辐射传递模拟计算中,纯海水吸收和散射、浮游植物吸收的数据或模型是已发表的被广泛采用的结果.黄色物质和非浮游植物颗粒吸收以及海洋颗粒物的散射模型从COASTLOOC数据中导出.另外,还利用了一个新的海洋颗粒物后向散射概率模型,在该模型中颗粒物后向散射概率是颗粒有机物与SPM比值和波长的甬数.把上述定义的固有光学性质作为输入,经过辐射传递模拟得到海表面以下辐照度反射比数据集,然后将该模拟数据集用于训练不同的人工神经网络,获取水色和水色要素浓度之间甬数关系的最佳近似.利用以上建立的基于人工神经网络的算法,把COASTLOOC数据集和PMNS数据集的辐照度反射比作为输入进行水色要素反演,通过比较反演值和真实测量值来评价算法性能.结果显示,建立的基于ANN的二类水体水色要素反演算法具有很好的性能.     

Influence of (FeO+TiO2) abundance on the microwave thermal emissions of lunar regolith

One of the essential controls on the microwave thermal emissions (MTE) of the lunar regolith is the abundance of FeO and TiO₂, known as the (FeO+TiO₂) abundance (FTA). In this paper, a radiative transfer simulation is employed first to study the change in the brightness temperature (T_B) with FTA under a range of frequencies and surface temperatures. Then, we analyze the influence of FTA on the MTE of the lunar regolith using microwave sounder (CELMS) data from the Chang’E-2 lunar orbiter, Clementine UV-VIS data, and lunar samples recovered from the Apollo and Surveyor projects. We conclude that: (1) FTA strongly influences the MTE of the lunar regolith, but it is not the decisive control, and (2) FTA decreases slightly with depth. This research plays an essential role in appropriately inverting CELMS data to obtain lunar regolith parameters.     

Oxygen isotope measurements of phosphate from fish teeth and bones

In situ measurements of lunar surface brightness temperatures made as a part of the Apollo Lunar Surface Experiments Package at the Apollo 15 Hadley Rille landing site are reported. Data derived from 5 thermocouples of the Heat Flow Experiment, which are lying on or just above the surface, are used to examine the thermal properties of the upper 15 cm of the lunar regolith using eclipse and nighttime cool-down temperatures. Application of finite-difference techniques in modeling the lunar soil shows the thermocouple data are best fit by a model consisting of a low-density and low-thermal conductivity surface layer approximately 2 cm thick overlying a region increasing in conductivity and density with depth. Conductivities on the order of 1 × 10^?5 W/cm-°K are postulated for the upper layer, with conductivity increasing to the order of 1 × 10^?4 W/cm-°K at depths exceeding 20 cm. An increase in mean temperature with depth indicates that the ratio of radiative to conductive transfer at 350°K is 2.7 for at least the upper few centimeters of lunar soil; this value is nearly twice that measured for returned lunar fines. The thermal properties model deduced from Apollo 15 surface temperatures is consistent with earth-based microwave observations if electrical properties measured on returned lunar fines are assumed.     

The effect of the pollution of the ground layer of the atmosphere on the albedo

Summary A simple model is proposed suitable for studying the effect of the ground layer of the atmosphere, polluted by aerosol, on the albedo. This model is founded on solving the equation of transfer of radiative energy. The numerical results are discussed, particular attention being paid to the analysis of errors due to neglecting the multiple reflection of solar radiation on the aerosol particles. A method which would also include the multiple reflection is proposed, and the conditions under which the presence of the aerosol is responsible for an increase or decrease of the solar radiation balance on the Earth's surface, are analysed.     

Radiative forcing by contrails

A parametric study of the instantaneous radiative impact of contrails is presented using three different radiative transfer models for a series of model atmospheres and cloud parameters. Contrails are treated as geometrically and optically thin plane parallel homogeneous cirrus layers in a static atmosphere. The ice water content is varied as a function of ambient temperature. The model atmospheres include tropical, mid-latitude, and subarctic summer and winter atmospheres. Optically thin contrails cause a positive net forcing at top of the atmosphere. At the surface the radiative forcing is negative during daytime. The forcing increases with the optical depth and the amount of contrail cover. At the top of the atmosphere, a mean contrail cover of 0.1% with average optical depth of 0.2 to 0.5 causes about 0.01 to 0.03 Wm⁻² daily mean instantaneous radiative forcing. Contrails cool the surface during the day and heat the surface during the night, and hence reduce the daily temperature amplitude. The net effect depends strongly on the daily variation of contrail cloud cover. The indirect radiative forcing due to particle changes in natural cirrus clouds may be of the same magnitude as the direct one due to additional cover.