二流—四流球谐函数谱展开累加辐射传输方案在全球气候模式中的应用

本文首先构建了二流—四流球谐函数谱展开累加辐射传输的新方案,然后将其应用于国家气候中心第二代大气环流模式BCC_AGCM2.0.1的新版本中,并与模式中原有的Eddington累加方案进行了比较。由于新方案本质上是单层Eddington近似方案在四流上的推广。因此新方案在计算精度上要优于原方案。通过在全球气候模式中的应用与比较,本文发现新方案对气候模拟会产生比较大的影响。在晴空条件下,新方案计算的在南纬30°到60°区间、北大西洋东北部以及非洲北部的撒哈拉沙漠区域的地表向下年平均短波辐射通量要小于原方案结果,最大差别可以达到3.5 W/m2;同时,新方案计算的在南纬30°到60°区间和北大西洋东北部的大气顶向上年平均短波辐射通量要大于原方案结果,最大差别达到3 W/m2。在有云大气情况下,新方案计算的地表向下年平均短波辐射通量要小于原方案结果,并随着纬度的增加,新旧两种方案的差别逐渐变大,在南北极时达到最大5.5 W/m2;同时,新方案计算的在赤道区域的大气顶的年平均短波向上辐射通量要小于原方案结果,最大差别为2.5 W/m2,而在南北纬30°到60°区间,新方案计算的在大气顶的年平均短波向上辐射通量则要大于原方案结果,最大差别为1.5 W/m2。新方案计算的年平均短波加热率普遍高于原方案结果,特别是在800 hPa到地表之间的低层大气以及50 hPa到100 hPa的高层大气,最大差别可达0.03 K/d。因此,新方案有助于改善全球气候模式中普遍存在的赤道平流层中下层的温度冷偏差现象。     

陆面模拟中植被辐射传输参数化方案研究

在冠层二流辐射传输模式基础上新发展了一个描述太阳短波辐射在植被中传输的冠层四流辐射传输模式.冠层四流辐射传输模式是在大气辐射传输理论的基础上得到一组描述短波辐射在植被中传输过程的冠层辐射传输基本方程,引进大气中求解辐射传输方程的四流近似解法,并求得冠层四流辐射传输方程的解析解.方程中各项参量能够反映叶子或冠层特殊的几何和光学特征.冠层向上、向下辐射通量取决于冠层散射相函数、叶子在入射光方向投影面积、单个叶子反射率和透射率、叶面积指数以及直射光入射太阳高度角等.四流模式计算叶子水平倾角时对太阳短波辐射的反照率,与二流模式结果比较可以验证模式的理论推导和建模都是正确的:计算结果的比较,表明四流模式在水平叶角分布时计算的冠层反照率与二流模式结果一致,同时直射光从任何太阳高度角入射的冠层反照率结果也一致,从而证明发展的冠层四流辐射传输模式是成功的.模拟试验中将两种模型同时耦合到同一个陆面过程模式中进行比较试验,结果表明,冠层四流辐射传输模式能够得到更精确的植被反照率,从而使得陆面模式计算的地表吸收的净太阳辐射通量更接近于观测值.     

云三维结构对加热率廓线及通量影响的模拟

大气加热率强度与结构是驱动全球大气环流的关键因素,实际加热率分布与云三维结构密切相关。作者使用三维蒙特卡洛辐射传输模式,模拟计算了云分辨模式所得3个典型三维云场的加热率廓线及通量;定义了两个参数来同时描述加热率廓线的垂直分布和强度,通过与独立像素近似算法对比,定量统计分析了高分辨率下云三维结构对辐射的影响。结果表明,在高分辨率条件下云三维结构对加热率廓线和通量影响十分显著,且不同结构云场所体现的影响各具特点,提出需要考察现有大气模式中云三维结构对当前所用加热率计算方案的订正方向。     

中国东部大气温湿廓线特征及其对辐射收支计算影响的分析

辐射平衡是气候决定性因素之一, 温湿廓线对辐射计算有重要影响。利用中国东部地区1960—2008年74个探空站资料, 将SBDART（Santa Barbara DISORT Atmospheric Radiative Transfer）辐射传输模式自带的大气温湿廓线和探空测得的实际大气温湿廓线进行对比, 并代入SBDART模式中, 分别计算晴空条件下不同温湿廓线对辐射通量的影响。结果表明:（1）中国东部地区的大气温湿廓线有明显的区域和季节特点;（2）模式地面向下辐射通量与实际相比有较明显差异, 尤其是在东北地区的冬季、中东部地区的夏季和华南地区的冬季, 差值达20—35 W/m2, 相对误差达2.01%—3.18%;（3）在东北地区的冬季、中东部地区的冬季以及华南地区的夏、冬季, 模式计算的大气顶向上辐射通量与实际相比差值达10—22 W/m2, 相对误差可达3.67%—8.94%;（4）模式与实际辐射加热率的差值在0.03—0.29 K/d。研究表明, 模式自带的大气温湿廓线区域和季节划分并不够细致, 不足以代表中国东部地区各个区域及季节的温湿特点。建立一套中国东部地区的大气标准廓线可以为辐射模拟提供更准确的输入量。      

利用TOVS资料计算我国东南地区的太阳直接辐射和散射辐射

根据多次散射理论--离散纵标法，将极轨卫星产品资料（TOVS资料）作为初始资料，利用中分辨率大气传输模式MODTRAN3分别计算出晴空和有云情况下5个层次的观测波段的太阳直接辐射和散射辐射，与国家一级探空探测的大气温、湿廓线计算出的辐射值进行比较。发现晴空状况下基本一致，有云时存在一定误差，表明由卫星的TOVS资料估算我国东部地区的晴空辐射资料，可以弥补我国辐射资料的不足。     

大气粒子散射相函数的参数化方案比较及其改进

大气粒子散射相函数的参数化是大气辐射传输参数化的重要组成部分。文中全面比较了大气粒子的HenyeyGreenstein(HG)方案和双Henyey-Greenstein(DHG)方案,并在四流球谐函数展开累加法中,应用这两种相函数参数化方案计算气溶胶、云、霾粒子的反射率、透射率或吸收率。该研究结果表明:HG方案无法表现相函数的后向峰值,因而其计算的大气粒子反射率和透射率精度较差;DHG方案能较好地表征相函数的整体特征,但是该方案计算的相函数易出现后向异常峰值或为负值,并导致计算得到的气溶胶、云、霾粒子的反射率和透射率精度甚至会低于HG方案。对DHG方案进行进一步研究,提出了改进的DHG方案(MDHG)。MDHG方案计算结果稳定,并能很好表征相函数的前向和后向峰值的特征,其计算的大气粒子的反射率和透射率精度也较高。因此,MDHG方案是一种理想的相函数参数化方案。     

南半球中高纬度区域不同类型云的辐射特性

利用CloudSat的2B-CLDCLASS-LIDAR云分类产品和2B-FLXHR-LIDAR辐射产品4 a（2007-2010年）的数据,定量分析了单层云（高云、中云、低云）和3种双层云（如:高云与中云共存、高云与低云共存以及中云与低云共存）在南半球中高纬度（40°-65°S）的云量、云辐射强迫和云辐射加热率。其中云辐射加热率定义为有云时的大气加热率廓线与晴空大气加热率廓线的差值。结果表明:研究区域盛行单层低云和单层中云,其云量分别为44.1%和10.3%。并且,中云重叠低云在双层云中云量也是最大（8.7%）。不同类型云的云量也显著影响着其云辐射强迫。单层低云在大气层顶、地表以及大气中的净云辐射强迫分别是-64.8、-56.5和-8.4 W/m2,其绝对值大于其他类型云。虽然单层的中云在大气层顶和地表的净辐射强迫也为负值,但其在大气中的净云辐射强迫为正值（2.3 W/m2）。最后,讨论了不同类型云对大气中辐射能量垂直分布的影响。所有类型云的短波（或长波）云辐射加热率都随高度升高表现为由负值转为正值（或由正值转为负值）。对于大部分云,其净云辐射加热率主要由长波云辐射加热率决定。这些研究结果旨在为模式中云重叠参数化方案在区域的适用性评估及改进提供观测依据。      

陆面过程中冠层四流辐射传输模式的模拟性能检验

利用实际观测资料检验新发展的植被冠层四流辐射传输模式的模拟性能.将浙江大学田间观测数据用于模式试验,其结果表明:在可见光波段,因叶片散射率值较小,四流模式和二流模式模拟结果差别不大,模拟的冠层反照率都接近观测值.这表明在可见光波段,二流模式已经能够较好模拟芯层反照率,四流模式能够提高的精度范围有限;在近红外波段,因叶片散射率值较大,两个模式模拟结果差别较明显,四流模式模拟的冠层反照率相对二流模式的模拟结果更接近观测值.利用第2次雪模式比较计划SNOWMIP资料和改进的10层陆面模式BATS进行耦合试验,将四流模式以及二流模式均耦合到该陆面模式中,其目的是为考察四流模式对陆面模式模拟地表反射太刚短波辐射通量的影响.耦合后得到的结果与SNOWMIP中加拿大BERMS草地站和森林站的观测资料进行了对比.对比结果表明,采用网流模式、二流模式、BATS原辐射传输模式后,耦合四流模式的陆面模式模拟地表反射太阳短波辐射通量最接近观测值.说明采用四流模式能够改善陆面模式对地表反射太阳短波辐射通啦的模拟.     

基于MCARaTS模式的三维云辐射通量模拟分析

云通过辐射反馈机制影响着云系的演变过程,对天气和气候系统起重要的调节作用。目前在研究辐射传输问题时普遍采用基于一维平面平行辐射传输的独立柱近似方法,忽略水平云大气柱之间的辐射传输过程,无法考虑三维云辐射效应,由此造成偏差可能影响天气和气候预报的准确性。为此,该文基于Linux系统,采用MCARaTS模式,基于大涡模拟的大气和地表光学性质参数,通过模拟光子随机路径,得到0.67 μm、2.13 μm、11 μm 3个波段下三维蒙特卡罗辐射传输和独立柱近似的向上反射和向下透射的三维辐射通量以及大气加热率空间分布场,通过Matlab进行数据预处理和图表绘制。三维辐射传输模式考虑了在水平方向上的光子传输,得出的辐射场具有真实的三维辐射现象,比独立柱近似模拟结果具有更高的准确性。     

两种计算长波辐射传输的近似方法

在二流近似的基础上,提出了两种计算晴空大气中长波辐射传输的方法,整个长波辐射区域被分成两个宽的光谱带,它们分别对相应的宽带光学厚度和发射率作参数化处理。这两种方法,所需的计算时间很少,但有相当好的精度,因而适于在大气环流模式中进行辐射传输的计算。     

计算非均一大气条件下太阳辐射通量的一个简单而又精确的模式

发展了一个计算非均一大气条件下太阳辐射通量的一个简单而又精确的模式，其中包括关于大气球反射率与透过率的一个参数化表达式，并引入加权一次散射反照率和加权不对称因子，用于拟合非均一大气条件下计算辐射通量的四个经验订正因子。对清洁和浑浊的两类大气，都具有120060组的辐射通量模拟试验，以检验本模式的精度。这些模拟试验覆盖0－50的云光学厚度、0－0．8地表反射率、Junge和对数正态的气溶胶谱分布、－0．05气溶胶折射率虚部。在均一大气条件下，由本模式计算的120060组向上通量的标准差对清洁和浑浊两类大气分别为1．08％和1．04％；而向下通量的标准差分别为4.12%和3．31％。在非均一大气条件下，由本模式计算的向上通量的标准差对清洁和浑浊两类大气分别为3．01％和3．48％；而向下通量的标准差分别为4．54％和4．89％，其精度远优于均一假设下的计算结果。     

A simple yet more accurate model to calculate solar radiative flux in the inhomogeneous atmosphere

A simple yet more accurate semiempirical model is developed to calculate solar radiative flux in the optically inhomogeneous atmosphere. In the model a parameterized expression of spherical reflectance and transmitance of the atmosphere is confirmed, and the weighted single scatter albedo and weighted asymmetric factor are introduced to fit four empirical correction factors responsible for radiative fluxes in the inhomogeneous atmosphere. For both clean and turbid models, there are 120060 sets of radiative flux simulations for accuracy checks of the model, which cover 0-50 cloud optical depths, 0-0.8 surface reflectance, Junge and Log-normal aerosol size distributions, and 0-0.05 imaginary parts of aerosol refractive indexes. In case of the homogeneous atmosphere, standard errors of the 120060 upward fluxes from the present model are 1.08% and 1.04% for clean and turbid aerosol models, respectively; and those of the downward fluxes are 4.12% and 3.31%. In case of the inhomogeneous atmosphere, standard errors of the upw ard fluxes from the present model are 3.01% and 3.48% for clean and turbid aerosol models.respectively; and those of the downward fluxes are 4.54% and 4.89%, showing a much better accuracy than the results calculated by using an assumption of the homogeneous atmosphere.     

中国地区夏季平均加热率的时空分布特征

The latitude-altitude distributions of radiative fluxes and heating rates are investigated by utilizing CloudSat satellite data over China during summer. The Tibetan Plateau causes the downward shortwave fluxes of the lower atmosphere over central China to be smaller than the fluxes over southern and northern China by generating more clouds. The existence of a larger quantity of clouds over central China reflects a greater amount of solar radiation back into space. The vertical gradients of upward shortwave radiative fluxes in the atmosphere below 8 km are greater than those above 8 km. The latitudinal-altitude distributions of downward longwave radiative fluxes show a slantwise decreasing trend from low latitudes to high latitudes that gradually weaken in the downward direction. The upward longwave radiative fluxes also weaken in the upward direction but with larger gradients. The maximum heating rates by solar radiation and cooling rates by longwave infrared radiation are located over 28-40°N at 7-8 km mean sea level (MSL), and they are larger than the rates in the northern and southern regions. The heating and cooling rates match well both vertically and geographically.     

Thermal radiative fluxes through inhomogeneous cloud fields: a sensitivity study using a new stochastic cloud generator

We analyze the effects of flat and bumpy top, fractional and internally inhomogeneous cloud layers on large area-averaged thermal radiative fluxes. Inhomogeneous clouds are generated by a new stochastic model: the tree-driven mass accumulation process (tdMAP). This model is able to provide stratocumulus and cumulus cloud fields with properties close to those observed in real clouds. A sensitivity study of cloud parameters is done by analyzing differences between 3D fluxes simulated by the spherical harmonic discrete ordinate method and three “standard” models likely to be used in general circulation models: plane-parallel homogeneous cloud model (PPH), PPH with fractional cloud coverage model (FCPPH) and independent pixel approximation model (IPA). We show that thermal fluxes are strong functions of fractional cloud coverage, mean optical depth, mean geometrical thickness and cloud base altitude. Fluctuations of “in-cloud” horizontal variability in optical depth and cloud-top bumps have negligible effects in the whole. We also showed that PPH, FCPPH and IPA models are not suitable to compute thermal fluxes of flat top fractional inhomogeneous cloud layer, except for completely overcast cloud. This implies that horizontal transport of photon at thermal wavelengths is important when cloudy cells are separated by optically thin regions.     

Radiative effects of Pinatubo aerosols in the middle latitudes of the Northern Hemisphere

Summary Using scattering coefficient profiles of the Pinatubo aerosols derived from the observation of skylight polarization and lidar backscattering ratio in Beijing, the radiative effect of Pinatubo aerosols in middle latitudes is assessed by a delta-four-stream radiative transfer model. It is shown that the Pinatubo aerosols significantly change the radiation field. Due to the presence of the volcanic aerosols, the downward short wave flux at the surface decreases with a maximum of 8 W/m² while the upward short wave flux at the top of the atmosphere increases with a maximum of 6.5 W/m². The volcanic aerosols are injected into the region bounded below by the tropopause and up by the 25 km level. The upward and downward radiative fluxes are changed in opposite directions at those two boundaries. Downward short wave fluxes below the tropopause are 7–9 W/m² less than background values and downward long wave fluxes below the tropopause are 2 W/m² more than background values. Upward short wave fluxes above 25 km level is 5–7 W/m² more and upward fluxes above there are about 3 W/m² less.The effects of the Pinatubo aerosols on heating rates are also significant. The maximum increase in the short wave heating rate can be as large as 0.2 K/day at 22 km. The increase in the long wave heating rate is less with a maximum amplitude of about 0.15 K/day. The maximum increase of the total heating rate is about 0.35 K/day, which is comparable with the heating rate caused by the ozone 9.6 µm band in this region.Results of this study are compared with studies of Lacis et al. (1992) and Russell et al. (1993) as well as ERBE measurements. The results generally agree well. Causes for the differences are analyzed.Based on the numerical study, it is also found that the LOWTRAN fresh volcanic model is not representative for the Pinatubo aerosols.With 9 Figures     

The solution of the vector radiative transfer equation using the discrete ordinates technique: Selected applications

The paper is devoted to the presentation of a newly developed radiative transfer code SCIAPOL_1.0 for a plane-parallel turbid slab illuminated by the monodirectional wide beam. The SCIAPOL_1.0 is based on the discrete-ordinates solution of the vector radiative transfer equation.The code is applied to a number of problems, including studies of the applicability of the scalar approximation for the calculation of light reflectance from aerosols, clouds, and molecular atmospheres. We also study the accuracy of the single scattering approximation as applied to the calculation of light reflection from molecular and cloudy atmospheres and propose new approximations for the calculation of the reflection function and the degree of light polarization under unpolarized light illumination conditions.     

Four-stream Radiative Transfer Parameterization Scheme in a Land Surface Process Model

Accurate estimates of albedos are required in climate modeling. Accurate and simple schemes for radiative transfer within canopy are required for these estimates, but severe limitations exist. This paper developed a four-stream solar radiative transfer model and coupled it with a land surface process model. The radiative model uses a four-stream approximation method as in the atmosphere to obtain analytic solutions of the basic equation of canopy radiative transfer. As an analytical model, the four-stream radiative transfer model can be easily applied efficiently to improve the parameterization of land surface radiation in climate models. Our four-stream solar radiative transfer model is based on a two-stream short wave radiative transfer model. It can simulate short wave solar radiative transfer within canopy according to the relevant theory in the atmosphere. Each parameter of the basic radiative transfer equation of canopy has special geometry and optical characters of leaves or canopy. The upward or downward radiative fluxes are related to the diffuse phase function, the G-function, leaf reflectivity and transmission, leaf area index, and the solar angle of the incident beam. The four-stream simulation is compared with that of the two-stream model. The four-stream model is proved successful through its consistent modeling of canopy albedo at any solar incident angle. In order to compare and find differences between the results predicted by the four- and two-stream models, a number of numerical experiments are performed through examining the effects of different leaf area indices, leaf angle distributions, optical properties of leaves, and ground surface conditions on the canopy albedo. Parallel experiments show that the canopy albedos predicted by the two models differ significantly when the leaf angle distribution is spherical and vertical. The results also show that the difference is particularly great for different incident solar beams. One additional experiment is carried out to evaluate the simulations of the BATS land surface model coupled with the two- and four-stream radiative transfer models. Station observations in 1998 are used for comparison. The results indicate that the simulation of BATS coupled with the four-stream model is the best because the surface absorbed solar radiation from the four-stream model is the closest to the observation.     

The Effects of Different HITRAN Versions on Calculated Long-Wave Radiation and Uncertainty Evaluation

Four editions of the High Resolution Transmission (HITRAN) databases (HITRAN96, HITRAN2K, HITRAN04, and HITRAN08) are compared by using a line-by-line (LBL) radiative model in the long-wave calculation for six typical atmospheres. The results show that differences in downward radiative fluxes between HITRAN96 and HITRAN08 at the surface can reach a maximum of 1.70 W m 2 for tropical atmospheres. The largest difference in heating rate between HITRAN96 and HITRAN08 can reach 0.1 K day 1 for midlatitude summer atmosphere. Uncertainties caused by line intensity and air-broadened halfwidths are also evaluated in this work using the uncertainty codes given in HITRAN08. The uncertainty is found to be 1.92 W m 2 for upward fluxes at the top of the atmosphere (TOA) and 1.97 W m 2 for downward fluxes at the surface. The largest heating rate caused by the uncertainty of line intensity and air-broadened half-width can reach 0.5 K day 1 . The differences in optical depths between 1300 and 1700 cm 1 caused by different HITRAN versions are larger than those caused by the uncertainties in intensity and air-broadened half-width. This paper suggests that there is inaccurate representation of line parameters over some spectral ranges in HITRAN and more attention should be paid to these ranges in fields such as remote sensing.