紫外波段多分散系气溶胶散射相函数随机抽样方法研究

基于电磁散射与辐射传输中的基本理论,对紫外波段霾尺度范围内满足特定分布的多种气溶胶粒子的散射相函数进行了研究.提出了一种直接随机抽样拟合散射相函数的方法.比较了H-G相函数、改进的H-G相函数及随机抽样拟合的相函数与多分散系Mie相函数的偏离程度.数值计算了不同相函数拟合方法对应气溶胶的传输特性.计算结果表明,相函数的准确模拟计算对于蒙特卡罗方法等辐射传输问题的解决具有十分重要的意义.     

Solving the atmospheric scattering optical transfer function using the multi-coupled single scattering method

The atmospheric scattering optical transfer function （OTF） is solved by applying the multi-coupled single scattering （MCSS） method to the three-dimensional radiative transfer equation （RTE） under the periodic ground condition. This approach is a direct hit to the atmospheric scattering OTF using the same original context of modulation transfer function （MTF） measurement, i.e., images of sinusoidal grating at different spatial frequencies. Both the amplitude and phase shift of the OTF at various zenith and azimuth angles can be obtained at an arbitrary spatial frequency.     

大气下行长波辐射快速实验测量方法研究

采用辐射传输模式模拟了特殊角度随气溶胶光学厚度及大气温湿度的变化情况,结果显示晴空下该角度不受气溶胶影响,受大气温度影响小,但受大气水汽含量影响显著,因此晴空下热红外光谱仪测量下行辐射时特殊天顶角可只根据大气水汽含量确定。西部荒漠地区下行辐射测量实验结果表明,据大气水汽含量确定的特殊角度值与实验结果一致,且该角度处辐射作为积分辐射符合精度要求。特殊角度值可用研究区大气水汽含量确定,该快速测量法可广泛用于野外实验,方便快捷测量准确瞬时下行辐射。     

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.     

Correction factors for a total scatter/backscatter nephelometer

Total aerosol scattering and backscattering atmospheric values are typically obtained with an integrating nephelometer. Due to design limitations, measurements do not cover the full (0°–180°) angular range, and correction factors are necessary. The effect of angle cutoff is examined for a number of particle size distributions and refractive indices. Scattering data for sub-micron particles can be corrected by the use of a modified Anderson approximation, while data for larger particle distributions can be approximated by a function of the effective size parameter. Correction factors for the hemispheric backscatter ratio are found to be small if nonsphericity is assumed. Such approximations will help more accurate corrections for angle range, particularly at large size parameter values.     

地基消光测量确定大气气溶胶模型

 分别取大陆型、海洋型、城市型和Junge谱分布气溶胶模型，用6S辐射传输算法计算出对应于太阳光度计测量时的各波段大气气溶胶光学厚度。将模式计算值与测量值进行比较，确定测量地区的大气气溶胶模型。将该方法用于2004年在北京地区测量的太阳光度计数据，结果显示该地区当日实际大气在几种气溶胶模型中较为符合城市型气溶胶模型。     

Reflector-based phase calibration of ultrasound transducers

Recently, the measurement of phase transfer functions (PTFs) of piezoelectric transducers has received more attention. These PTFs are useful for e.g. coding and interference based imaging methods, and ultrasound contrast microbubble research. Several optical and acoustic methods to measure a transducer’s PTF have been reported in literature. The optical methods require a setup to which not all ultrasound laboratories have access to. The acoustic methods require accurate distance and acoustic wave speed measurements. A small error in these leads to a large error in phase, e.g. an accuracy of 0.1% on an axial distance of 10 cm leads to an uncertainty in the PTF measurement of ±97° at 4 MHz. In this paper we present an acoustic pulse-echo method to measure the PTF of a transducer, which is based on linear wave propagation and only requires an estimate of the wave travel distance and the acoustic wave speed. In our method the transducer is excited by a monofrequency sine burst with a rectangular envelope. The transducer initially vibrates at resonance (transient regime) prior to the forcing frequency response (steady state regime). The PTF value of the system is the difference between the phases deduced from the transient and the steady state regimes. Good agreement, to within 7°, was obtained between KLM simulations and measurements on two transducers in a 1-8 MHz frequency range. The reproducibility of the method was ±10°, with a systematic error of 2° at 1 MHz increasing to 16° at 8 MHz. This work demonstrates that the PTF of a transducer can be measured in a simple laboratory setting.     

非球形气溶胶粒子散射相函数经验公式

散射相函数是研究电磁波传输特性的重要参数，直接影响电磁波传输方程的简化程度和解的精度。基于电磁散射与辐射传输中的基本理论，对非球形粒子散射相函数的经验公式进行了研究。为了很好的模拟非球形粒子的后向散射峰值，提高辐射传输方程的简化程度和解的精度，提出了一种新的相函数经验公式。分析新的相函数对非球形粒子的适用性，以单个沙尘性气溶胶为例，计算了不同形状粒子的Henyey-Greenstein*相函数和新的相函数随角度的变化，并与T矩阵法的计算结果进行了对比，发现椭球形粒子的长短轴比和有限长圆柱形粒子的径长比大于0.5时，新的相函数在大角度后向散射部分与T矩阵法的吻合程度较高。考虑波长变化，对比了尺寸谱满足对数正态分布的四种气溶胶粒子的Henyey-Greenstein*相函数和新的相函数与T矩阵法的计算结果。研究表明，对于椭球形粒子和有限长圆柱形粒子，在大角度(大于90°)后向散射部分，除了0.694时的椭球形海洋性气溶胶，新的相函数均方根差较小的占100%，证明了新的相函数可以较好的模拟非球形粒子的后向散射特征。新的相函数对准确模拟辐射传输过程具有重要意义。     

Parameterization of land surface albedo

Remote measurements of Earth’s surface from ground, airborne, and spaceborne instruments show that its albedo is highly variable and is sensitive to solar zenith angle（SZA） and atmospheric opacity. Using a validated radiative transfer calculating toolbox, DISORT and a bidirectional reflectance distribution function library, AMBRALS, a land surface albedo（LSA） lookup table（LUT） is produced with respect to SZA and aerosol optical depth. With the LUT, spectral and broadband LSA can be obtained at any given illumination geometries and atmospheric conditions. It provides a fast and accurate way to simulate surface reflectance over large temporal and spatial scales for climate study.     

Light scattering properties of sea-salt aerosol particles inferred from modeling studies and ground-based measurements

Direct climate radiative forcing depends on the aerosol optical depth τ, the single scattering albedo ?, and the up-scatter fraction β; these quantities are functions of the refractive index of the particles, their size relative to the incident wavelength, and their shape. Sea-salt aerosols crystallize into cubic shapes or in agglomerates of cubic particles under low relative humidity conditions. The present study investigates the effects of the shape of dried sea-salt particles on the detection of light scattering from the particles. Ground-based measurements of scattering and backscattering coefficients have been performed with an integrating nephelometer instrument for a wavelength . The measurements are compared to two models: the Mie theory assuming a spherical shape for the particles and the Discrete Dipole Approximation (DDA) model for the hypothesis of cubic shape of the sea-salt aerosols. The comparison is made accurately by taking into account the actual range of the scattering angles measured by the nephelometer in both models that is from 7° to 170° for the scattering coefficient and from 90° to 170° for the backscattering coefficient. Modeled scattering and backscattering coefficients increase for nonspherical particles compared to spherical shape of particles with diameter larger than about 1 μm. However, the comparison of the modeling results with the measurements gives best agreement for particles diameter less than about 1 μm. The size distribution of the particles is measured with two instruments with different size bins: an electrical low-pressure impactor (ELPI) and an aerodynamic particle sizer (APS). It is found that the size of the bins of the instruments to determine the number concentration of the particles in accordance with their diameter is critical in the comparison of measurements with modeling.     

Application of SH surface acoustic waves for measuring the viscosity of liquids in function of pressure and temperature

Viscosity measurements were carried out on triolein at pressures from atmospheric up to 650 MPa and in the temperature range from 10 °C to 40 °C using ultrasonic measuring setup. Bleustein–Gulyaev SH surface acoustic waves waveguides were used as viscosity sensors. Additionally, pressure changes occurring during phase transition have been measured over the same temperature range. Application of ultrasonic SH surface acoustic waves in the liquid viscosity measurements at high pressure has many advantages. It enables viscosity measurement during phase transitions and in the high-pressure range where the classical viscosity measurement methods cannot operate. Measurements of phase transition kinetics and viscosity of liquids at high pressures and various temperatures (isotherms) is a novelty. The knowledge of changes in viscosity in function of pressure and temperature can help to obtain a deeper insight into thermodynamic properties of liquids.     

A new approach to the retrieval of surface properties from earthshine measurements

Instruments such as the MODIS and MISR radiometers on EOS AM-1, and POLDER on ADEOS have been deployed for the remote sensing retrieval of surface properties. Typically, retrieval algorithms use linear combinations of semi-empirical bidirectional reflectance distribution function (BRDF) kernels to model surface reflectance. The retrieval proceeds in two steps; first, an atmospheric correction relates surface BRDF to top-of-atmosphere (TOA) reflectances, then regression is used to establish the linear coefficients used in the kernel combination. BRDF kernels may also depend on a number of physical or empirical non-linear parameters (e.g. ocean wind speed for a specular BRDF); such parameters are usually assumed known. A major source of error in this retrieval comes from lack of knowledge of planetary boundary layer (PBL) aerosol properties.In this paper, we present a different approach to surface property retrieval. For the radiative transfer simulations, we use the discrete ordinate LIDORT model, which has the capability to generate simultaneous fields of radiances and weighting functions in a multiply scattering multi-layer atmosphere. Surface-atmosphere coupling due to multiple scattering and reflection effects is treated in full; the use of an atmospheric correction is not required. Further, it is shown that sensitivities of TOA reflectances to both linear and non-linear surface BRDF parameters may be established directly by explicit analytic differentiation of the discrete ordinate radiative transfer equations. Surface properties may thus be retrieved directly and conveniently from satellite measurements using standard non-linear fitting methods. In the fitting for BRDF parameters, lower-boundary aerosol properties can either be retrieved as auxiliary parameters, or they can be regarded as forward model parameter errors. We present examples of simulated radiances and surface/aerosol weighting functions for combinations of multi-angle measurements at several different wavelengths, and we perform some examples of self-consistent non-linear fitting to demonstrate feasibility for this kind of surface property retrieval.     

Integrating nephelometer performance in fog and thick Haze

Summary In fogs and thick hazes integrating-nephelometer measurements are affected by the angular truncation error due to the integration of the atmospheric scattering function on a truncated angular range. A method of experimental investigation on the actual consequences of this error is exposed. Equations that take into account both truncation error and photometric systematic errors are proposed and validated through a field comparison by nephelometer and transparency measurements. It is shown that the integrating nephelometer can substitute the transmissometer in all visibility conditions if photometric errors are avoided. Paper presented at the 1° Congresso del Gruppo Nazionale per la Fisica dell'Atmosfera e dell'Oceano, June 19–22, 1984, Rome.     

Sensitivity of multi-angle photo-polarimetry to vertical layering and mixing of absorbing aerosols: Quantifying measurement uncertainties

The impact of tropospheric aerosols on climate can vary greatly based upon relatively small variations in aerosol properties, such as composition, shape and size distributions, as well as vertical layering. Polarimetric measurements have been advocated in recent years as an additional tool to better understand and retrieve the aerosol properties needed for improved predictions of aerosol radiative forcing on climate. The goal of this study is to introduce a formal approach to assessing the sensitivity of both intensity and polarization signals to absorbing aerosol layering, explicitly accounting for instrument measurement uncertainties. If ignored, sensitivity to aerosol height can introduce biases in aerosol property retrievals at short (ultraviolet or blue) wavelengths; if properly exploited, it may enable the extraction of some basic information on aerosol profiles. Employing a vector successive-orders-of-scattering (SOS) radiative transfer code, we conducted modeling experiments to determine how the measured Stokes vector elements are affected at 446 nm (blue band) by the vertical distribution, mixing and layering of smoke and dust aerosols under the assumption of a simple Lambertian surface and predefined aerosol microphysical properties. We find that smoke and dust vertical layering, if ignored, can introduce biases in radiometric and polarimetric aerosol property retrievals for aerosol optical depth (AOD) above 0.3 (polarimetric) and AOD above 0.5 (radiometric), and should, therefore, be accounted for in retrievals at high aerosol loadings.     

Aerosol light absorption and its measurement: A review

Light absorption by aerosols contributes to solar radiative forcing through absorption of solar radiation and heating of the absorbing aerosol layer. Besides the direct radiative effect, the heating can evaporate clouds and change the atmospheric dynamics. Aerosol light absorption in the atmosphere is dominated by black carbon (BC) with additional, significant contributions from the still poorly understood brown carbon and from mineral dust. Sources of these absorbing aerosols include biomass burning and other combustion processes and dust entrainment.For particles much smaller than the wavelength of incident light, absorption is proportional to the particle volume and mass. Absorption can be calculated with Mie theory for spherical particles and with more complicated numerical methods for other particle shapes.The quantitative measurement of aerosol light absorption is still a challenge. Simple, commonly used filter measurements are prone to measurement artifacts due to particle concentration and modification of particle and filter morphology upon particle deposition, optical interaction of deposited particles and filter medium, and poor angular integration of light scattered by deposited particles. In situ methods measure particle absorption with the particles in their natural suspended state and therefore are not prone to effects related to particle deposition and concentration on filters. Photoacoustic and refractive index-based measurements rely on the heating of particles during light absorption, which, for power-modulated light sources, causes an acoustic signal and modulation of the refractive index in the air surrounding the particles that can be quantified with a microphone and an interferometer, respectively. These methods may suffer from some interference due to light-induced particle evaporation. Laser-induced incandescence also monitors particle heating upon absorption, but heats absorbing particles to much higher temperatures to quantify BC mass from the thermal radiation emitted by the heated particles. Extinction-minus-scattering techniques have limited sensitivity for measuring aerosol light absorption unless the very long absorption paths of cavity ring-down techniques are used. Systematic errors can be dominated by truncation errors in the scattering measurement for large particles or by subtraction errors for high single scattering albedo particles. Remote sensing techniques are essential for global monitoring of aerosol light absorption. While local column-integrated measurements of aerosol light absorption with sun and sky radiometers are routinely done, global satellite measurements are so far largely limited to determining a semi-quantitative UV absorption index.     

大气模式与气溶胶模型对辐射传输计算的影响

大气模式与气溶胶模型选择是影响定量遥感应用的辐射传输计算的重要因素。人们一般凭感性认识去进行选择,有一定随意性,对其带来的影响程度关注甚少。以太阳辐射计测量为依据,对大气模式与气溶胶模型选择方法作了研究,并利用经辐射定标过的光谱辐射计地面测量对辐射传输计算精度进行了验证。在昆明进行的试验表明,在0.50~0.68μm范围内,选择中纬度冬季大气模式和大陆型气溶胶,经辐射传输计算后得到的光谱辐射亮度与光谱辐射计测量结果一致性很好,差别在3.3%以内;变换大气模式对辐射传输计算产生明显影响,差别达10%左右;选择不同的气溶胶模型对辐射传输计算影响也很大,差别达11%左右。基于辐射计测量的大气模式与气溶胶模型选择避免了主观选择的不可靠性,有益于减少辐射传输计算或卫星遥感大气订正的误差。     

Radiation field in a semi-infinite homogeneous atmosphere with internal sources

The equation of radiative transfer in a semi-infinite homogeneous atmosphere with different internal sources is solved by the method of kernel approximation—the kernel in the equation for the Sobolev resolvent function is approximated by a Gauss-Legendre sum. Then the obtained approximate equation can be solved exactly and the solution is a weighted sum of exponentials. All the necessary coefficients of the solutions may be easily found. Since the resolvent function is closely connected with the Green function of the integral radiative transfer equation, the radiation field for different internal sources can be found by simple integration. For the considered cases the formulas for the radiation field are obtained and the respective accuracy estimated. The package of codes in Fortran-77 is given at http://www.aai.ee/∼viik/homogen.for.     

大气折射对可见光波段辐射传输特性的影响

折射是影响辐射传输的重要因素. 为分析大气折射对辐射传输的影响, 基于Monte Carlo方法, 给出了考虑大气折射的矢量辐射传输模型, 实现了均匀气层和耦合面处光子随机运动过程的模拟, 实现了直射光及漫射光Stokes矢量、偏振度和辐射通量等参数的计算. 在考虑和不考虑大气折射两种条件下, 验证了模型的准确性; 在纯瑞利散射条件下, 讨论了大气折射对不同方向漫射光Stokes矢量的影响; 在不同太阳天顶角、大气廓线、气溶胶及含云大气条件下, 分析了大气折射对辐射传输过程的影响. 结果表明: 大气折射对漫射光Stokes矢量的影响主要体现在天顶角70°–110°区间, 且随着太阳入射角增大, 其影响更为显著; 不同大气廓线情形下, 大气折射对Stokes矢量的影响不一致, 其原因是不同大气廓线对应的折射率廓线存在差异. 含云及含气溶胶大气条件下, 大气折射对辐射传输的影响变弱, 沙尘型及海盐型气溶胶条件下, 折射对辐射传输的影响强于可溶型气溶胶情形; 不同形状气溶胶条件下, 大气折射对辐射传输的影响也存在显著差异; 不同云高条件下, 大气折射对漫射光Stokes矢量的影响无显著差异, 但随着云光学厚度增大, 大气折射的影响减弱.     

基于多角度多光谱Stokes参数Q和U测量反演气溶胶偏振相函数

偏振相函数是气溶胶重要的光学参数之一,它对气溶胶复折射指数、粒子尺度和形状都十分敏感。多角度多光谱偏振遥感可以有效获取气溶胶偏振相函数信息。新一代CIMEL太阳-天空偏振辐射计CE318-DP作为高精度地基气溶胶偏振遥感仪器已被引入全球气溶胶自动观测网AERONET(AErosol RObotic NETwork),并作为扩展多波长偏振测量的太阳-天空辐射计观测网SONET(Sun/sky-radiometer Observation NETwork)的主要仪器,已在不同类型气溶胶观测站点积累了多年的偏振数据。但目前偏振反演仅能利用线偏振度或偏振辐亮度。与线偏振度和偏振辐亮度相比,Stokes参数Q和U不仅包含天空光线偏振强度信息还包含偏振方向信息。利用CE318-DP多光谱多角度测量的天空光Stokes参数Q和U反演气溶胶偏振相函数的方法。针对CE318-DP标准主平面偏振观测模式PPP(Polarized Principal Plane)下Stokes参数U对气溶胶特性变化不敏感、信息难以利用的不足,测试了新的平纬圈偏振扫描模式ALMP(ALMucantar Polarization)获取Stokes参数Q和U,并成功应用于偏振相函数的反演。系统分析了340～1 640 nm多光谱通道上典型生物质燃烧型气溶胶和水溶性气溶胶的-P₁₂/P₁₁反演结果并测试了反演方法在晴朗和灰霾不同大气条件下的适用性。无论在主平面还是平纬圈观测几何下,反演结果在可见光和近红外通道上均与真实值具有较好的一致性。进一步讨论了模型中基于气溶胶参数初始值和大气气溶胶参数真实值计算的“大气单次散射/大气散射”的比值近似相等的假设条件在短波通道不能很好地满足是造成紫外波段反演结果偏差较大的原因之一。后续有待进一步提高反演模型在短波通道的适用性,为利用不同光谱通道上-P₁₂/P₁₁的变化特征改进气溶胶微物理参数反演奠定了基础。     

Simultaneous retrieval of aerosols and ocean properties: A classic inverse modeling approach. I. Analytic Jacobians from the linearized CAO-DISORT model

In this paper and the sequel, we investigate the application of classic inverse methods based on iterative least-squares cost-function minimization to the simultaneous retrieval of aerosol and ocean properties from visible and near infrared spectral radiance measurements such as those from the SeaWiFS and MODIS instruments. Radiance measurements at the satellite are simulated directly using an accurate coupled atmosphere-ocean-discrete-ordinate radiative transfer (CAO-DISORT) code as the main component of the forward model. For this kind of cost-function inverse problem, we require the forward model to generate weighting functions (radiance partial derivatives) with respect to the aerosol and marine properties to be retrieved, and to other model parameters which are sources of error in the retrievals.In this paper, we report on the linearization of the CAO-DISORT model. This linearization provides a complete analytic differentiation of the coupled-media radiative transfer theory, and it allows the model to generate analytic weighting functions for any atmospheric or marine parameter. For high solar zenith angles, we give an implementation of the pseudo-spherical (P-S) approach to solar beam attenuation in the atmosphere in the linearized model. We summarize a number of performance enhancements such as the use of an exact single-scattering calculation to improve accuracy. We derive inherent optical property inputs for the linearized CAO-DISORT code for a simple 2-parameter bio-optical model for the marine environment coupled to a 2-parameter bimodal atmospheric aerosol medium.