陆面过程模式中植被冠层短波辐射传输参数化方案研究进展

阐述植被冠层内短波辐射传输模式研究在陆面过程及陆气相互研究中的重要性,全面地回顾有关冠层短波辐射传输模式的研究进展,结合目前在陆面过程领域具有代表性的冠层短波辐射传输模式,指出当前在这一领域研究中存在并有待于解决的问题,探讨了未来的发展趋势.     

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

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

安徽省淮南森林冠层辐射传输过程的特征

本文利用淮南森林观测站2018年7月1日至2019年6月30日冠层辐射观测,分析了淮南栎树森林下垫面冠层内外辐射变化特征。结果表明:（1）从春季到夏季,栎树冠层之上向下的太阳短波辐射增加,到冬季逐渐减少。从早春开始,由于叶片生长增多,冠层中间和冠层之下向下的太阳短波辐射下降,从秋季到冬季树叶凋落,其向下的太阳辐射增加,与冠层之上的变化趋势相反;对于向上的短波辐射,无论冠层之上、冠层中间还是冠层之下,随季节的变化都与向下的短波辐射相似,只是数值小很多。（2）冠层之上、冠层中间和冠层之下向下的长波辐射,随时间的变化从春季逐渐开始增大至夏季达到最大,随后逐渐减小并在冬季达到最小;就空间变化而言,冠层中间和冠层之下向下的长波辐射值比冠层之上的辐射值高,使得冠层对长波辐射的振幅增大,晴空条件最高可达1.3倍。（3）淮南森林区冠层之上（距地面25 m）年平均反照率为0.14,比中国北方地区（35°N）温带季风气候区（混交林为主）反照率的整体水平低0.01,表明淮南的森林茂密、灌丛更多些。（4）冠层上部分和整层的短波辐射透射率主要受叶片的影响。夏季,冠层的短波透射率平均为0.1。到了冬天,叶子凋落,透射率增加并趋于一个平稳的波动。冠层的短波辐射吸收率在夏季最高,秋季逐渐降低,随着叶子凋落在冬季迅速减小,趋于一常值。     

多层四流球谐函数算法的构建及在大气辐射传输模式中的应用

为了在不大幅度增加计算成本的情况下提高大气辐射传输计算的精度,利用单层四流球谐函数结合多层二流累加法,构造了可用于多层大气的四流球谐函数算法。为了比较与其他辐射传输算法的差异,引入48流离散纵坐标算法作为比较标准,Eddington 近似、四流离散纵坐标算法作为比较对象。在真实大气廓线条件下,计算了晴空和有云大气顶向上辐射通量、地表向下辐射通量以及加热率廓线。得出以下结论：在晴空情况下,与作为标准的48流离散纵坐标法相比,Eddington 近似、四流离散纵坐标法和新构造的四流球谐函数方法加热率绝对误差都小于0?3 K／d;向上、向下辐射通量的相对误差分别小于1％和0?6％。这表明在晴空情况下,3种算法对加热率的计算精度差别不大;对辐射通量的计算精度,两种四流近似算法比传统的 Eddington 近似更为精确。在有云情况下,与48流离散纵坐标法相比,四流球谐函数和四流离散纵坐标法计算的云顶加热率相对误差小于1％,而 Eddington 近似计算的云顶加热率相对误差大于5％。结果表明：新构造的四流球谐函数算法可用于大气辐射传输模式,在不大幅度增加计算成本的同时,提高了晴空大气的整体辐射计算精度和有云大气辐射加热率的计算精度。     

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

本文首先构建了二流—四流球谐函数谱展开累加辐射传输的新方案,然后将其应用于国家气候中心第二代大气环流模式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。因此,新方案有助于改善全球气候模式中普遍存在的赤道平流层中下层的温度冷偏差现象。     

辐射传输模式中地表参数对大气长波辐射的影响

本文利用Liou-Ou一维宽带辐射传输模式，对地表热力参数取值部分作了改进，使用模式大气和青藏高原实测资料对下垫面温度与地表空气温度两者不能合二为一的问题进行了分析，同时还讨论了下垫面温度的日变化对大气长波辐射通量日变化的影响及地表比辐射率的变化对大气长波辐射通量计算结果的修正作用。     

影响植被内部辐射状况的冠层结构特征研究

植被冠层结构特征是影响植被内部辐射状况的重要因子之一,对于植被群体的光合特性来说,冠层结构特性影响了光合有效辐射的截获、吸收和透射。本文着重分析了冠层结构的数量特征（如叶面积指数）和冠层几何结构特征,综合考虑叶倾角、叶方位角和Ｇ函数等群体几何特征,以便直接用于遥感光合植被生物量机理估测模型的建立。     

CLM4.0土壤水分传输方案改进在青藏高原陆面过程模拟中的效应

利用2010年5月25日-12月31日玛曲高寒草原的气象观测资料和陆面过程模式(CLM4.0)对玛曲高寒草原陆面过程进行了数值模拟。通过评估模式的模拟性能、模式对含砂量的敏感程度以及模式土壤水分传输方案改进对青藏高原地区陆面过程模拟的影响,发现CLM4.0模式能较好地再现观测站土壤温、湿度、地表辐射、湍流通量等的变化趋势,但土壤温度模拟偏低,感热通量模拟偏大;含砂量增多会减弱土壤的持水能力,使得夏季感热通量增大而潜热通量减小;CLM4.0模式中新引入的有机质对土壤温、湿度模拟均有重要影响,Richards方程和径流计算的修改则对土壤含水量模拟影响较大,这对其他陆面模式的改进具有一定的指导意义。     

大气NLTE辐射传输中分子光谱参数间的关系

讨论Ｅｉｎｓｔｅｉｎ参数、谱振子强度、权重跃迁矩的平方、线强及吸收系数等分子光谱参数的定义，导出了它们之间的关系，为计算上层大气辐射提供了物理基础。     

用云和辐射传输模式研究对流性降水云微物理及辐射特性

利用中国气象科学研究院的三维对流云模式模拟了 1 998年发生在我国武汉地区的特大暴雨过程中的对流云。在模拟云的成熟阶段 ,模拟结果和地面最大降水量实况及热带测雨卫星TRMM (TropicalRainfallMeasureMission)的降水雷达PR(Pre cipitationRadar)的零度层高度和降水粒子达到的高度较接近。将该云物理模式的输出结果输入到Liu的三维辐射模式的输入场 ,模拟了以TRMM的微波成像仪TMI(TRMMMicrowaveImager)的扫描角为入射角 ,频率为 85 5GHz和 37 0GHz的星载微波辐射仪接受到的亮温 ,模拟的结果表明云中各微物理量的分布对上行微波辐射的影响较大。     

A Comparison of Two Canopy Radiative Models in Land Surface Processes

This paper compares the predictions by two radiative transfer models-the two-stream approximation model and the generalized layered model （developed by the authors） in land surface processes -for different canopies under direct or diffuse radiation conditions. The comparison indicates that there are significant differences between the two models, especially in the near infrared （NIR） band. Results of canopy reflectance from the two-stream model are larger than those from the generalized model. However, results of canopy absorptance from the two-stream model are larger in some cases and smaller in others compared to those from the generalized model, depending on the cases involved. In the visible （VIS） band, canopy reflectance is smaller and canopy absorptance larger from the two-stream model compared to the generalized model when the Leaf Area Index （LAI） is low and soil reflectance is high. In cases of canopies with vertical leaf angles, the differences of reflectance and absorptance in the VIS and NIR bands between the two models are especially large. Two commonly occurring cases, with which the two-stream model cannot deal accurately, are also investigated. One is for a canopy with different adaxial and abaxial leaf optical properties; and the other is for incident sky diffuse radiation with a non-uniform distribution. Comparison of the generalized model within the same canopy for both uniform and non-uniform incident diffuse radiation inputs shows smaller differences in general. However, there is a measurable difference between these radiation inputs for a canopy with high leaf angle. This indicates that the application of the two-stream model to a canopy with different adaxial and abaxial leaf optical properties will introduce non-negligible errors.     

A simplified scheme of the generalized layered radiative transfer model

In this paper, firstly, a simplified version （SGRTM） of the generalized layered radiative transfer model （GRTM） within the canopy, developed by us, is presented. It reduces the information requirement of inputted sky diffuse radiation, as well as of canopy morphology, and in turn saves computer resources. Results from the SGRTM agree perfectly with those of the GRTM. Secondly, by applying the linear superposition principle of the optics and by using the basic solutions of the GRTM for radiative transfer within the canopy under the condition of assumed zero soil reflectance, two sets of explicit analytical solutions of radiative transfer within the canopy with any soil reflectance magnitude are derived： one for incident diffuse, and the other for direct beam radiation. The explicit analytical solutions need two sets of basic solutions of canopy reflectance and transmittance under zero soil reflectance, run by the model for both diffuse and direct beam radiation. One set of basic solutions is the canopy reflectance αf （written as α1 for direct beam radiation） and transmittance βf （written as β1 for direction beam radiation） with zero soil reflectance for the downward radiation from above the canopy （i.e. sky）, and the other set is the canopy reflectance （αb） and transmittance βb for the upward radiation from below the canopy （i.e., ground）. Under the condition of the same plant architecture in the vertical layers, and the same leaf adaxial and abaxial optical properties in the canopies for the uniform diffuse radiation, the explicit solutions need only one set of basic solutions, because under this condition the two basic solutions are equal, i.e., αf = αb and βf = βb. Using the explicit analytical solutions, the fractions of any kind of incident solar radiation reflected from （defined as surface albedo, or canopy reflectance）, transmitted through （defined as canopy transmittance）, and absorbed by （defined as canopy absorptance） the canopy and other properties per     

A generalized layered radiative transfer model in the vegetation canopy

In this paper, a generalized layered model for radiation transfer in canopy with high vertical resolution is developed. Differing from the two-stream approximate radiation transfer model commonly used in the land surface models, the generalized model takes into account the effect of complicated canopy morphology and inhomogeneous optical properties of leaves on radiation transfer within the canopy. In the model, the total leaf area index （LAI） of the canopy is divided into many layers. At a given layer, the influences of diffuse radiation angle distributions and leaf angle distributions on radiation transfer within the canopy are considered. The derivation of equations serving the model are described in detail, and these can deal with various diffuse radiation transfers in quite broad categories of canopy with quite inhomogeneons vertical structures and uneven leaves with substantially different optical properties of adaxial and abaxial faces of the leaves. The model is used to simulate the radiation transfer for canopies with horizontal leaves to validate the generalized model. Results from the model are compared with those from the two-stream scheme, and differences between these two models are discussed.     

Improving the vegetation dynamic simulation in a land surface model by using a statistical-dynamic canopy interception scheme

Canopy interception of incident precipitation, as a critical component of a forest＇s water budget, can affect the amount of water available to the soil, and ultimately vegetation distribution and function. In this paper, a statistical-dynamic approach based on leaf area index and statistical canopy interception is used to parameterize the canopy interception process. The statistical-dynamic canopy interception scheme is implemented into the Community Land Model with dynamic global vegetation model （CLM-DGVM） to improve its dynamic vegetation simulation. The simulation for continental China by the land surface model with the new canopy interception scheme shows that the new one reasonably represents the precipitation intercepted by the canopy. Moreover, the new scheme enhances the water availability in the root zone for vegetation growth, especially in the densely vegetated and semi-arid areas, and improves the model＇s performance of potential vegetation simulation.     

A Comparison of Two Canopy Radiative Models in Land Surface Processes

This paper compares the predictions by two radiative transfer models-the two-stream approximation model and the generalized layered model (developed by the authors) in land surface processes-for different canopies under direct or diffuse radiation conditions. The comparison indicates that there are significant differences between the two models, especially in the near infrared (NIR) band. Results of canopy reflectance from the two-stream model are larger than those from the generalized model. However, results of canopy absorptance from the two-stream model are larger in some cases and smaller in others compared to those from the generalized model, depending on the cases involved. In the visible (VIS) band, canopy reflectance is smaller and canopy absorptance larger from the two-stream model compared to the generalized model when the Leaf Area Index (LAI) is low and soil reflectance is high. In cases of canopies with vertical leaf angles, the differences of reflectance and absorptance in the VIS and NIR bands between the two models are especially large.Two commonly occurring cases, with which the two-stream model cannot deal accurately, are also investigated. One is for a canopy with different adaxial and abaxial leaf optical properties; and the other is for incident sky diffuse radiation with a non-uniform distribution. Comparison of the generalized model within the same canopy for both uniform and non-uniform incident diffuse radiation inputs shows smaller differences in general. However, there is a measurable difference between these radiation inputs for a canopy with high leaf angle. This indicates that the application of the two-stream model to a canopy with different adaxial and abaxial leaf optical properties will introduce non-negligible errors.     

多层城市冠层模式的建立及数值试验研究

为在城市气象数值模拟中更好地体现由城市发展引起的下垫面土地利用改变及人为活动对大气过程的影响,建立了基于建筑物三维分布的多层城市冠层模式,冠层内动力方程组考虑了建筑物冠层拖曳力的作用及雷诺应力的影响,通过引入建筑物宽度、间距以及垂直分布密度指数等建筑物形态特征参数,以更好地体现城市复杂地表对大气温度、湿度及动量方程的影响.同时,该模式分屋顶、4个侧壤及地面分别考虑辐射及能量平衡求解表面温度,计算各表面与大气的通量交换,并考虑辐射阴影效应、冠层内部各个面之间的可视因子、以及与冠层内建筑物密度指数、可视因子等相关的多重反射辐射导致的辐射截陷作用.模式的离线检验结果表明:(1)冠层模式计算风廓线与风洞实验测量数据吻合良好;(2)离线冠层模式能够模拟实际小区的风速、温度垂直廓线,并能够较好地体现小区内气温日变化.冠层模式与区域边界层模式耦合检验结果表明:(1)耦合模拟的近地面(2 m处)气温及地表温度的结果明显优于传统的水泥平板方案,尤其是在夜间,水泥平板方案与实测气温最大偏差4 K左右,耦合模拟方案为1-2 K;(2)耦合模拟方案考虑了建筑物对冠层之上的拖曳力影响以及建筑物形态结构对雷诺应力的影响,风速(10 m处)计算结果与观测值相差约在1 m/s,水泥平板方案偏差3 m/s左右.     

土壤热异常影响地表能量平衡的个例分析和数值模拟

The statistical relationship between soil thermal anomaly and short-term climate change is presented based on a typical case study. Furthermore, possible physical mechanisms behind the relationship are revealed through using an off-line land surface model with a reasonable soil thermal forcing at the bottom of the soil layer.In the first experiment, the given heat flux is 5 W m-2 at the bottom of the soil layer (in depth of 6.3 m)for 3 months, while only a positive ground temperature anomaly of 0.06℃ can be found compared to the control run. The anomaly, however, could reach 0.65℃ if the soil thermal conductivity was one order of magnitude larger. It could be even as large as 0.81℃ assuming the heat flux at bottom is 10 W m-2. Meanwhile, an increase of about 10 W m-2 was detected both for heat flux in soil and sensible heat on land surface, which is not neglectable to the short-term climate change. The results show that considerable response in land surface energy budget could be expected when the soil thermal forcing reaches a certain spatial-tem poral scale. Therefore, land surface models should not ignore the upward heat flux from the bottom of the soil layer. Moreover, integration for a longer period of time and coupled land-atmosphere model are also necessary for the better understanding of this issue.     

陆面模型Noah-MP的不同参数化方案在沙漠区域的适用性研究

为研究陆面模型Noah-MP在沙漠下垫面的最优参数化方案组合,本文利用中国气象局塔克拉玛干沙漠气象野外科学试验基地观测数据,根据沙漠环境特征进行不同参数化方案组合的三组模拟实验,利用观测数据对10 cm土壤温湿度、感热、潜热通量模拟值对比分析得出最优组合。研究表明:第三组对10 cm的土壤温度模拟效果最好,主要原因是Chen97感热交换系数和全网格二流近似(gap=0)辐射传输方案比较符合沙漠的环境特征。三组试验对土壤湿度模拟效果差,其主要原因是沙漠的土壤信息未能体现在模式中,第二组选择CLM方案对土壤类型影响蒸发方面有一定考虑,其模拟结果相对较好。对于感热通量,第一、二组模拟值在波峰存在高估,尤其是第二组模拟值在降水后出现了明显低估情况,第三组模拟效果最好,主要得益于选择了感热交换系数Chen97方案,能够较为真实的刻画Ch变化特征。潜热通量在四个特征量中模拟效果最差,主要原因是沙漠土壤水分极低,观测降水和实际进入土壤的水量有差异,另外没有植被和植物根系,模式无法准确计算土壤蒸发和植被蒸散。根据统计分析和泰勒图可知,第三组能够更好地还原沙漠区域的陆面过程。