Effect of ice crystal shape and effective size on snow bidirectional reflectance

We tested the applicability of three rigorous radiative transfer computational approaches, namely, the discrete ordinates radiative transfer (DISORT) method, the adding–doubling approach, and an efficient computational technique based on Ambartsumian's nonlinear integral equation for computing the bidirectional reflectance of a semi-infinite layer. It was found that each of these three models, in a combination with the truncation of the forward peak of the bulk scattering phase functions of ice particles, can be used to simulate the bidirectional reflectance of a semi-infinite snow layer with appropriate accuracy. Furthermore, we investigate the sensitivity of the bidirectional reflectance of a homogeneous and optically infinite snow layer to ice crystal habit and effective particle size. It is shown that the bidirectional reflectance is not sensitive to the particle effective size in the visible spectrum. The sensitivity of the bidirectional reflectance in the near-infrared spectrum to the particle effective size increases with the increase of the incident wavelength. The sensitivity of the bidirectional reflectance to the effective particle size and shape is attributed fundamentally to the sensitivity of the single-scattering properties to particle size and shape. For a specific ice crystal habit, the truncated phase function used in the radiative transfer computations is not sensitive to particle effective size. Thus, the single-scattering albedo is primarily responsible for the sensitivity of the bidirectional reflectance to particle size, particularly, at a near-infrared wavelength.     

全极化微波辐射计姿态对观测亮温的影响及消除

针对全极化微波辐射计精确探测的需求, 研究了辐射计姿态对于观测亮温的影响以及亮温误差校正. 建立了姿态偏移与观测入射角以及极化旋转角的关系, 模拟了观测亮温随观测入射角以及极化旋转角波动的变化; 仿真了姿态偏移情况下的辐射计原始观测亮温;运用一种基于辐射传输模型的姿态补偿方法, 以垂直极化亮温和第三Stokes参数亮温为例, 对原始观测亮温展开误差校正. 研究显示, 该方法能够有效去除辐射计姿态偏移对观测亮温造成的影响, 校正结果满足辐射计数据预处理的误差精度要求. 关键词： 全极化微波辐射计 Stokes参数 观测入射角 极化旋转角     

Polarization difference due to nonrandomly oriented ice particles at millimeter/submillimeter waveband

This paper presents polarized signature due to oriented circular columnar and planar ice crystals at millimeter/submillimeter (mm/sub-mm) waveband. DDSCAT 6.1 and RT4 code package are employed for scattering properties and radiative transfer simulations, respectively, at the three estimated window frequencies (150, 220 and 340 GHz) of FengYun-4 (FY-4). We use empirical formulas to describe realistic sizes of planar and columnar particles and assume that ice particles are in Gamma-size distribution in this study. A “resonance” feature of polarized signals as a function of median mass diameter is notably found for horizontally oriented columns and blunt plates at the frequency of 340 GHz; however, there is no promising resonance characteristic for horizontally aligned plates with empirical sizes at the three window channels of FY-4. The position of the resonance peak is related to particle aspect ratio, frequency and ice water path (IWP), and it moves to a shorter median mass diameter when the particle aspect ratio decreases or IWP in clouds increases. Considering that particle canting angle distribution (Gaussian distribution in this study), polarization difference, as well as the brightness temperature difference between clear and cloudy sky, decreases rapidly when particles gradually change from horizontally oriented to randomly oriented. The upwelling brightness temperature is insensitive to particle size and shape but sensitive to particle orientation, the difference of brightness temperature between horizontal and random orientation up to 6 K, whereas polarized signature is quite sensitive to particle microphysics as well as orientation; polarized measurements thereby could benefit retrieval of cloud microphysical parameters.     

Optical modeling of vesicular volcanic ash particles

We present a realistic shape model for nonspherical, vesicular particles and use the model to derive single-scattering properties of volcanic fine-ash particles. Light-scattering computations with discrete-dipole approximation reveal that, qualitatively, scattering by the model particles resembles that of the measured, real volcanic ash particles. Comparison of compact and vesicular ash shows that porosity promotes positive degree of linear polarization and decreases the depolarization ratio for both large and small vesicles. Yet, the single-scattering properties of ash particles with large vesicles are found to be surprisingly similar to those of compact ash particles. A comparison with Mie computations of equal-volume spheres indicates that for small size parameters, the spherical shape underestimates the asymmetry parameter of volcanic ash particles; whereas, for larger size parameters, it is overestimated.     

The impact of ice particle roughness on the scattering phase matrix

The goal of this study is to explore the influence of ice particle habit (or shape) and surface roughness on the scattering phase matrix. As an example, reported here are the results for two wavelengths: 0.67 and 1.61 μm. For this effort, a database of single-scattering properties has been computed for a set of habits including hexagonal plates, hollow and solid columns, hollow and solid 3D bullet rosettes, droxtals, aggregates of solid columns, and aggregates of plates. The database provides properties for each of the habits at 101 wavelengths between 0.45 and 2.24 μm for smooth, moderately roughened, and severely roughened particles. At each wavelength, the scattering properties are provided at 233 discrete particle diameters ranging from 2 to 10,000 μm. A single particle size distribution from a very cold ice cloud sampled during the CRYSTAL-FACE field campaign (T_cld=–76 °C) is used to illustrate the influence of habit and roughness on the phase matrix. In all, four different habit mixtures are evaluated. The nonzero elements of the phase matrix are shown to be quite sensitive to the assumed habit, particularly in the case of ?P₁₂/P₁₁ that is associated with the degree of linear polarization of scattered radiation. Surface roughness is shown to smooth out maxima in the scattering phase function and in the other elements of the phase matrix, consistent with other studies. To compare with the theoretical simulations of the phase matrix for smooth and roughened particles, a full year of cloud-aerosol lidar with orthogonal polarization (CALIOP) data from 2008 is analyzed to provide global statistics on the values of P₁₁ and P₂₂/P₁₁ in the backscattering direction. In a comparison of two of the habit mixtures (one used for MODIS Collection 5 and another that incorporates new habits including hollow bullet rosettes and aggregates of plates) with the CALIOP data, the values for P₁₁ are higher regardless of the degree of particle surface roughness, and the values for P₂₂/P₁₁ are lower than those for CALIOP. Further investigation is warranted to better understand this discrepancy.     

A single-scattering approximation for infrared radiative transfer in limb geometry in the Martian atmosphere

We present a single-scattering approximation for infrared radiative transfer in limb geometry in the Martian atmosphere. It is based on the assumption that the upwelling internal radiation field is dominated by a surface with a uniform brightness temperature. It allows the calculation of the scattering source function for individual aerosol types, mixtures of aerosol types, and mixtures of gas and aerosol. The approximation can be applied in a Curtis-Godson radiative transfer code and is used for operational retrievals from Mars Climate Sounder measurements. Radiance comparisons with a multiple scattering model show good agreement in the mid- and far-infrared although the approximate model tends to underestimate the radiances in realistic conditions of the Martian atmosphere. Relative radiance differences are found to be about 2% in the lowermost atmosphere, increasing to ∼10% in the middle atmosphere of Mars. The increasing differences with altitude are mostly due to the increasing contribution to limb radiance of scattering relative to emission at the colder, higher atmospheric levels. This effect becomes smaller toward longer wavelengths at typical Martian temperatures. The relative radiance differences are expected to produce systematic errors of similar magnitude in retrieved opacity profiles.     

Numerical accuracy of “equivalent” spherical approximations for computing ensemble-averaged scattering properties of fractal soot aggregates

Numerical accuracy is quantitatively assessed in conjunction with the application of four “equivalent” spherical approximations in the computation of the optical properties of small aggregate soot particles. The approximations are based on equal volume, equal surface area, the radius of gyration, and a collection of independent spheres with the same volume and the same volume-to-projected area ratio as the original nonspherical particle. A diffusion-limited cluster-cluster aggregation algorithm is used to specify the geometries of soot particles. Furthermore, the Generalized Multi-particle Mie (GMM) method is utilized to compute the single-scattering properties of individual soot aggregate particles assumed to be randomly oriented in space. The ensemble-averaged single-scattering properties of the particles are obtained by accounting for the probability distribution functions (PDF) of the number of monomers per aggregate at two wavelengths, 0.628 and 1.1 μm. It is shown that all of the aforementioned equivalent-spherical approximations lead to large errors in the computation of the phase function.     

星载极化相关型全极化微波辐射计天线交叉极化校正技术 (II) : 校正试验

创建了地球场景数据集, 结合全极化微波辐射传输模型, 仿真了地球场景亮温. 基于自主推导的全极化天线温度方程, 通过GRASP9软件生成天线方向图, 模拟了辐射计的天线温度. 进而利用多元线性回归方法, 求取了天线交叉极化校正M矩阵, 实现了对星载极化相关型全极化微波辐射计天线交叉极化的校正. 试验结果表明: 天线温度与地球场景亮温之间具有良好的线性关系; 天线交叉极化对全极化微波辐射计正交通道亮温影响明显, 尤其以对垂直极化亮温误差的影响最为显著; 校正后各通道的天线交叉极化得到了有效的减小, 交叉极化优于-23 dB, 极化纯度大于99.5%, 采用M矩阵校正及消除天线温度中交叉极化亮温影响的方案是切实可行的. 该校正技术可以实现星载极化相关型全极化微波辐射计在轨运行后对于天线交叉极化的最终校正. 关键词： 全极化微波辐射计 天线交叉极化 天线温度方程 M矩阵')" href="#">M矩阵     

Validation of the community radiative transfer model

To validate the Community Radiative Transfer Model (CRTM) developed by the U.S. Joint Center for Satellite Data Assimilation (JCSDA), the discrete ordinate radiative transfer (DISORT) model and the line-by-line radiative transfer model (LBLRTM) are combined in order to provide a reference benchmark. Compared with the benchmark, the CRTM appears quite accurate for both clear sky and ice cloud radiance simulations with RMS errors below 0.2 K, except for clouds with small ice particles. In a computer CPU run time comparison, the CRTM is faster than DISORT by approximately two orders of magnitude. Using the operational MODIS cloud products and the European Center for Medium-range Weather Forecasting (ECMWF) atmospheric profiles as an input, the CRTM is employed to simulate the Atmospheric Infrared Sounder (AIRS) radiances. The CRTM simulations are shown to be in reasonably close agreement with the AIRS measurements (the discrepancies are within 2 K in terms of brightness temperature difference). Furthermore, the impact of uncertainties in the input cloud properties and atmospheric profiles on the CRTM simulations has been assessed. The CRTM-based brightness temperatures (BTs) at the top of the atmosphere (TOA), for both thin (τ<5) and thick (τ>30) clouds, are highly sensitive to uncertainties in atmospheric temperature and cloud top pressure. However, for an optically thick cloud, the CRTM-based BTs are not sensitive to the uncertainties of cloud optical thickness, effective particle size, and atmospheric humidity profiles. On the contrary, the uncertainties of the CRTM-based TOA BTs resulting from effective particle size and optical thickness are not negligible in an optically thin cloud.     

Attenuation of electromagnetic wave propagation in sandstorms incorporating charged sand particles

A theoretical approach for predicting the attenuation of microwave propagation in sandstorms is presented, with electric charges generated on the sand grains taken into account. It is found that the effect of electric charges distributed partially on the sand surface is notable. The calculated attenuation is in good agreement with that measured in certain conditions. The distribution of electric charges on the surface of sand grains, which is not easy to measure, can be approximately determined by measuring the attenuation value of electromagnetic waves. Some effects of sand radius, dielectric permittivity, frequency of electromagnetic wave, and visibility of sandstorms on the attenuation are also discussed quantitatively. Finally, a new electric parameter is introduced to describe the roles of scattering, absorption and effect of charges in attenuation.     

Temperature Diagnostics of a Z-Pinch Plasma Using Calculations of the Spectral Brightness of X-Ray Radiation in a Large Interval of Radiation Energies

We elaborate a modern approach to the temperature diagnostics of a Z-pinch plasma. The approach is based on quantum-mechanical calculations of spectral brightness for X-ray radiation performed in a large interval of the photon energy for several temperatures and densities. In a large interval of the photon energy, a range can be found where the spectral brightness is highly sensitive to the temperature variation. This fact enables temperature diagnostics without complicated analysis of the spectral-line shape used in traditional diagnostic methods. In our calculations of the spectral brightness of X-ray radiation, we use a theoretical model known as the ion model of a plasma. We discuss important features of this model along with the other theoretical models used for calculating the radiative properties of the plasma. We calculate the spectral brightness of X-ray radiation for molybdenum plasma at temperatures of 1 and 1.2 keV and plasma densities of 1 and 2 g/cm³ and find the range of X-ray radiation energies that can be used for the temperature diagnostics.     

Direct simulation Monte Carlo ray tracing model of light scattering by a class of real particles and comparison with PROGRA experimental results

The Direct Simulation Monte Carlo (DSMC) model is presented for three-dimensional single scattering of natural light by suspended, randomly oriented, optically homogeneous and isotropic, rounded and stochastically rough cubic particles. The modelled particles have large size parameter that allows geometric optics approximation to be used. The proposed computational model is simple and flexible. It is tested by comparison with known geometric optics solution for a perfect cube and Lorenz–Mie solution for a sphere, as extreme cases of the class of rounded cubes. Scattering and polarization properties of particles with various geometrical and optical characteristics are examined. The experimental study of real NaCl crystals with new Progra² instrument in microgravity conditions is conducted. The experimental and computed polarization and brightness phase curves are compared.     

A strong ice cloud event as seen by a microwave satellite sensor: Simulations and observations

In this article, brightness temperatures observed by channels of the Advanced Microwave Sounding Unit-B (AMSU-B) instrument are compared to those simulated by a radiative transfer model, which can take into account the multiple scattering due to ice particles by using a discrete ordinate iterative solution method. The input fields, namely, the pressure, temperature, humidity, and cloud water content are taken from the short range forecast from the Met Office mesoscale model (UKMES). The comparison was made for a case study on the 25 January 2002 when a frontal system associated with significant cloud was present over the UK. It is demonstrated that liquid clouds have maximum impact on channel 16 of AMSU whereas ice clouds have maximum impact on channel 20. The main uncertainty for simulating microwave radiances is the assumptions about microphysical properties, such as size distribution, shape and orientation of the cloud particles, which are not known in the mesoscale model. The article examines the impact of these parameters on the cloud signal. The polarisation signal due to oriented ice particles at these frequencies is also discussed.     

Scattering and absorption of light by ice particles: Solution by a new physical-geometric optics hybrid method

A new physical-geometric optics hybrid (PGOH) method is developed to compute the scattering and absorption properties of ice particles. This method is suitable for studying the optical properties of ice particles with arbitrary orientations, complex refractive indices (i.e., particles with significant absorption), and size parameters (proportional to the ratio of particle size to incident wavelength) larger than ∼20, and includes consideration of the edge effects necessary for accurate determination of the extinction and absorption efficiencies. Light beams with polygon-shaped cross sections propagate within a particle and are traced by using a beam-splitting technique. The electric field associated with a beam is calculated using a beam-tracing process in which the amplitude and phase variations over the wavefront of the localized wave associated with the beam are considered analytically. The geometric-optics near field for each ray is obtained, and the single-scattering properties of particles are calculated from electromagnetic integral equations. The present method does not assume additional physical simplifications and approximations, except for geometric optics principles, and may be regarded as a “benchmark” within the framework of the geometric optics approach. The computational time is on the order of seconds for a single-orientation simulation and is essentially independent of the size parameter. The single-scattering properties of oriented hexagonal ice particles (ice plates and hexagons) are presented. The numerical results are compared with those computed from the discrete-dipole-approximation (DDA) method.     

A numerical study on the contribution of melting layer composed by coated ice spheres to the upwelling brightness temperatures in TMI channels using VDISORT algorithm

To examine the influence of bright band on the retrieval of precipitation rate, the performance of melting layer composed by coated ice precipitable particles on the satellite-based measurement of polarized microwave brightness temperatures is discussed in this article by a vector discrete ordinate radiative transfer model. After comparing the simulated brightness temperatures in different TMI channels with and without the melting layer, we conclude that: (1) The melting layer composed by liquid-coated ice spheres weakens the upwelling microwave brightness temperatures because of the absorption/emission effect caused by the liquid coat. This effect is more conspicuous in middle and high frequency channels (19, 37 and 85 GHz) but, in 85 GHz channel, with the increase of rain rate, the multi-scattering can weaken its effect. (2) In a specific frequency, the horizontally polarized brightness temperature is more severely weakened by the melting layer than the vertically polarized. With the “cold” background (ocean surface, for example), this character is more conspicuous than that with a warm background. That is to say, the inner structure of a cloud system is easier to be detected under a cold background. Only in the 85 GHz frequency and when the rain rate is larger than 4 mm/h can we find that the vertically polarized brightness temperature is more severely weakened than the horizontally polarized one. (3) The melting layer with the assumption of coated ice spheres can change the difference of brightness temperatures between the vertically and horizontally polarized channels in the same frequency. In general, the value of such difference with the assumption of melting layer is larger than that without it. With a warm background, this value is negative and only in middle frequency (37 GHz), it is both stable and conspicuous.     

Prediction of particle sticking efficiency for fly ash deposition at high temperatures

The tendency of ash particles to stick under high temperatures is dictated by the ash chemistry, particle physical properties, deposit surface properties and furnace operation conditions. A model has been developed in order to predict the particle sticking efficiency for fly ash deposition at high temperatures. The model incorporates the particle properties relevant to the ash chemistry, particle kinetic energy and furnace operation conditions and takes into consideration the partial sticking behaviour and the deposit layer. To test the model, the sticking behaviours of synthetic ash in a drop tube furnace are evaluated and the slagging formation from coal combustion in a down-fired furnace is modelled. Compared with the measurements, the proposed model presents reasonable prediction performance on the particle sticking behaviour and the ash deposition formation. Through a sensitivity analysis, furnace operation conditions (velocity and temperature), contact angle and particle size have been found to be the significant factors in controlling the sticking behaviours for the synthetic ash particles. The ash chemistry and furnace temperature dictate the wetting potential of the ash particles and the melting ability of the deposit surface; particle size and density not only control the particle kinetic energy, but also affect the particle temperature. The furnace velocity condition has been identified as being able to influence the selective deposition behaviour, where the maximum deposition efficiency moves to smaller particles when increasing the gas velocity. In addition, the thermophoresis effect on the arrival rate of the particles reduces with increasing the gas velocity. Further, increasing the melting degree of the deposit layer could greatly enhance the predicted deposition formation, in particular for the high furnace velocity condition.     

Optical modeling of mineral dust particles: A review

Dust particles are uniquely and irregularly shaped, they can be inhomogeneous, form agglomerates, be composed of anisotropic materials, and have a preferred orientation. As such, modeling their light scattering is very challenging. This review takes a look at the advances in dust optical modeling over the last decade. It is obvious that our ability to model the single-scattering properties of dust particles accurately depends on the size parameter. Unfortunately, our ability to account realistically for all the relevant physical properties in light-scattering modeling is the best for small particles; whereas, the realistic treatment of the particles would be most important for large size parameters. When particles are not much larger than the wavelength, even simple model shapes such as homogeneous spheroids appear to perform well; practically any reasonable shape distribution of non-spherical model particles seems superior compared to the Mie theory. Our ability to model scattering by dust particles much larger than the wavelength is very limited: no method presently exists to predict reliably and accurately the single-scattering properties of such particles, although there are models that can be tuned to agree well with the laboratory-measured reference scattering matrices. The intermediate size parameters between the resonance domain and the geometric-optics domain appear to be almost uncharted territory and, consequently, very little can be said about the impact of different physical properties on scattering in this region. Despite the challenges, the use of Mie theory should be avoided: contrary to the popular belief, the use of Mie spheres is a major source of error even in radiation-budget considerations.     

The impact of Martian aerosols on the retrieval of temperature profiles from PFS measurements

Our purpose was a qualitative assessment of the impact of dust and water ice aerosols on the retrieved temperature profiles and the retrieval process itself in the Martian atmosphere. It aims to quantify the related uncertainties in the atmospheric temperature profiles derived from radiance measurements of the Planetary Fourier Spectrometer (PFS), currently operating on the Mars Express orbiter. In this study the effects of aerosol opacities on simulated data and retrieved temperature profiles were also investigated.From the analysis of the model atmosphere including dust and water ice with different size distributions it results that the dust component affects weighting functions and brightness temperatures less than water ice. A similar situation is also observed when different vertical distributions are considered. Unlike dust, water ice with different sizes of crystals evidently influences weighting functions and brightness temperatures. The impact of the considered water ice vertical distributions on brightness temperatures is noticeable only near 840 cm⁻¹.Considering different dust opacities, the largest differences—5 K maximum—between retrieved temperature profiles were observed close to the surface, regardless assumptions on a size distribution or the refractive indices. Contrary to dust, the different sizes of water ice particles assumed during retrieval stronger affected the retrieved temperature profiles than water ice opacities. Moreover, the differences in the retrieved temperature profiles were amplified while wrong optical properties for dust as well as for water ice aerosol were assumed instead of the nominal case. This means that the wrong assumption can induce an additional source of the retrieval error and lead to unreasonable temperature profiles. In the cases of expected heavily loads water ice crystals, their size distribution in the Martian atmosphere should be known from other observations before the retrieval of the temperature profile is attempted.For the analyzed examples of real PFS measurements the impact of different dust vertical distributions on the retrieval of temperature profile is prominent only in layers close to the surface. However, these differences remain comparable with retrieval errors. All influences of dust on weighting functions, brightness temperatures and during retrieval can be neglected if the noise equivalent radiance (NER) of PFS is taken into account.     

Electron spin resonance studies on quantum tunneling in spinel ferrite nanoparticles

The electron spin resonance (ESR) spectrometer, a very sensitive instrument with fast detecting window to explore quantum phase transitions for magnetic nanoparticles, was exploited to study the fascinating interplay between thermal and quantum fluctuations in the vicinity of a quantum critical point. We have measured ESR in ferrofluid samples containing nanosize particles of Fe₂O₃. The evolution of the ESR spectrum with temperature suggests that quantum tunneling of spins occurs in single domain magnetic particles in the low temperature regime. The effects of various microwave fields, particle sizes, and temperatures on the magnetic states of single domain spinel ferrite nanoparticles are investigated. We can consistently explain experimental data assuming that, as the temperature decreases, the spectrum changes from superparamagnetic (SPR) to blocked SPR and finally evolves quantum superparamagnetic behaviour as the temperature lowers down further. A nanoparticle system of a highly anisotropic magnetic material can be qualitatively specified by a simple quantum spin model, or by the Heisenberg model with strong easy-plane anisotropy.Received: 29 August 2003, Published online: 15 October 2003PACS: 76.30.-v Electron paramagnetic resonance and relaxation - 75.40.Cx Static properties (order parameter, static susceptibility, heat capacities, critical exponents, etc.) - 05.30.-d Quantum statistical mechanics - 75.50.Dd Nonmetallic ferromagnetic materials     

切比雪夫雾霾粒子的紫外光散射偏振特性研究

雾霾天气已严重影响人们的日常生活,通过检测雾霾粒子的紫外光散射偏振特性可以有效分析雾霾成因。矿物粒子、灰尘粒子等雾霾颗粒均有小规模表面粗糙度的形态学特征,因此可用切比雪夫粒子作为模型分析。“日盲”紫外光与切比雪夫雾霾粒子相互作用发生散射,散射光偏振特性可反演切比雪夫雾霾颗粒物理性质(如粒子尺寸参数、复杂折射率、粒子形变、波纹参数)。采用紫外光单次散射模型和T矩阵方法,仿真分析切比雪夫雾霾粒子物理参数与散射光偏振特性(Stokes矢量和偏振度)之间的关系,结果表明：粒径对散射光Stokes矢量Is和Qs随散射角的变化趋势影响很大,粒子的粒径和复杂折射率虚部的变化会造成散射光偏振度随散射角的变化趋势的改变;具体分析散射角度为10°时,得到粒径对Is和Qs的数值影响最大,当粒径r<1 μm时,Is随粒径呈现抛物线趋势;切比雪夫粒子形变的增大,Is呈现先增大后减小的趋势。