空气动力学阻抗算法在半干旱区的应用比较和遥感反演

结合2005年中国西北半干旱区定西试验站的观测资料,使用9种空气动力学方法,估算了阻抗和感热通量,并与涡动相关阻抗结果进行了对比,分析了5～9月涡度相关测得阻抗的平均日变化特征,结果表明:(1)由于西北干旱、半干旱区和黄土高原的稀疏植被地区在植物生长季节内存在温度高、降水少、空气干燥等特点,热量和水汽输送具有特殊性;Choudh-1、Verma-R空气动力学阻抗方法对热量和动量传输的粗糙度长度、热量和动量及水汽输送的修正函数描述得较好,空气动力学阻抗估算精度较好,由此获得的感热效果也较好,较适用于该地区。(2)风速对阻抗起主要作用;植被覆盖度增加和植被密度增大的同时也增加了热量和水汽传输的粗糙度,有助于空气动力学阻抗降低。(3)采用Choudh-1方法和遥感资料反演的阻抗效果较好。     

2008年夏季中国干旱-半干旱区陆面主要物理参数的平均特征

采用2008年7-9月观测的中国干旱-半干旱区试验观测协同与集成研究资料,选取12个测站(涵盖不同气候环境区的草地、农田、果林和森林等下垫面)比较分析了干旱-半干旱区的动力、热力和水汽粗糙度长度、总体输送系数、反照率以及土壤热传导率和热扩散率的夏季平均特征,并与陆面模式的理论参考值进行了对比。结果表明,植被下垫面的水、热特性和植被高度对主要陆面参数有重要影响,但植被下垫面的土壤热传导率和扩散率与下垫面类型之间无明显关系。动力学粗糙度长度与下垫面的植被高度呈正相关关系,无论是稳定还是不稳定层结下,动量拖曳系数都随着植被高度的增大而增大,而反照率随植被高度的增加而降低。在研究区的下垫面中,草地下垫面动力学粗糙度长度的实测值比模式理论参考值偏小,但农田和果林下垫面的实测值比理论参考值偏大;草地下垫面动量拖曳系数理论的参考值比实测值高,但其他下垫面的实测值高于理论参考值;实测反照率处于可见光反照率和近红外光反照率之间,但反照率的理论参考值比实测值偏高。     

卫星遥感结合地面观测资料对中国西北干旱区地表热力输送系数的估算

本文利用黑河野外试验 (HEIFE) 地面观测资料, 采用空气动力学方法计算了干旱区内不同下垫面的地表热力输送系数CH, 结合由美国国家海洋和大气局 (NOAA) 系列卫星遥感观测的反映地表植被特征的归一化差值植被指数 (NDVI) 资料, 经拟合得到了针对我国西北干旱区不同下垫面的CH-NDVI参数化关系式, 并对此关系式进行了合理性检验。结果表明: 对于区域尺度而言, 在缺乏用其他方法获得较准确的区域CH值的情况下, 利用卫星遥感结合地面观测资料对其估算是较为可靠的方法。     

绿洲-沙漠复合地表条件下的局地和有效粗糙度

地表空气动力学粗糙度是研究水平非均匀条件下陆面过程的基础,文中对黑河试验区各不同下垫面上地表空气动力学粗糙度进行了估算和分析,分别得到了地表空气动力学粗糙度ｚ０ｍ的有效值和局地值,发现各不同典型下垫面的局地ｚ０ｍ值明显不同,尤其在有作物存在时,绿洲的局地ｚ０ｍ值比沙漠、戈壁的值明显要大;而对同一测站,有效空气动力学粗糙度比局地空气动力学粗糙度值明显要大,运用应力分解理论对这种差异进行了解释,认为有效空气动力学粗糙度与复合表面的总应力有关,其代表了较大水平尺度内的表面和较高障碍物的综合作用,当涉及对陆表通量参数化问题时,必须首先确定水平尺度及相应的地表粗糙度参数     

农牧交错带陆面过程的数值模拟研究

利用CLM(Common Land Model)模式对我国内蒙古奈曼旗农牧交错带沙漠和农田两种不同典型下垫面的陆面过程进行了数值模拟试验,并与外场试验观测结果进行了对比分析。结果表明:无论是沙漠还是农田试验,CLM都能够较好地模拟其辐射通量和土壤中的热传导特征,CLM的模拟结果能够真实地再现试验期间土壤热传导过程对天气过程的响应。相比而言,模式对沙漠地区长波辐射通量和干燥时期短波辐射通量的模拟结果好于农田,其原因可能是因为农田下垫面植被及土壤特征较沙漠复杂,有着很大的不确定性,造成了农田地表反照率和温度模拟的偏差。而对农田热传导的模拟结果好于沙漠,反映了CLM对含水量较大、持水力较强的农田下垫面的热传导模拟能力较好,而对含水量较小、持水力较弱的沙漠下垫面的热传导模拟能力相对较差。     

中国西北干旱区戈壁下垫面夏季的热力输送

以敦煌戈壁站2004年6月和2008年8月的常规观测和超声观测为例,分析了西北干旱区戈壁下垫面夏季热力输送的一般过程及特征。首先评价了湍流通量的观测质量以及仪器观测的地表能量通量闭合问题,结果表明敦煌戈壁站的观测在白天总体较好。夏季地表能量通量的平均日变化显示,潜热通量整天都很小,可以忽略,白天到达地表的短波辐射以及地表向上的长波辐射非常强,地表净辐射主要转化为感热输送（敦煌戈壁站在中午时平均分别达380W·m-2以上和250W·m-2以上）;夜间土壤释放热量以平衡地表的辐射冷却,感热通量略低于0。白天时地表大气经常触发自由对流活动,影响动量通量的观测质量,并有效输送地表热力至上层大气中,有助于形成超厚大气边界层。分析了戈壁下垫面的动量粗糙度特征和热力粗糙度特征（敦煌戈壁站动量粗糙度约为0．6mm）,热力粗糙度基本小于动量粗糙度一个量级,这符合目前对干旱区戈壁下垫面热力输送特征的初步认识。     

城市下垫面空气动力学参数的估算

为了定量描述北京城市下垫面的空气动力学特征,为模式提供准确的下垫面参数.利用气象塔大气湍流观测资料,结合Martano(2000)由单层超声风速、温度资料估算非均匀下垫面空气动力学参数的方法,计算了中国科学院大气物理研究所气象塔附近的下垫面空气动力粗糙度z0和零平面位移d,即z0为1.75 m和d为40.12 m.该结果与前人研究成果的比较结果表明最近的8年内,该处的零平面位移和空气动力粗糙度明显增大,这与该塔周围城市建设状况吻合.     

城市下垫面空气动力学参数确定方法综述

地表粗糙度和零平面位移是2个重要的空气动力学参数,对于城市下垫面,粗糙元的形态非常复杂,空间分布上存在较大的非均匀性,如何确定城市下垫面地表粗糙度和零平面位移目前尚无最佳方案。概述了这2个参数适用于城市非均匀下垫面的参数化方法,包括风温观测方法和粗糙元形态学方法。探讨了各自在应用过程中存在的问题及进一步的研究工作,并且提出,通过数学模型的改进和卫星遥感技术的引入,粗糙元形态学方法会有更好的发展前景：     

一种大气污染物干沉积速率的计算方法及其应用

介绍了一种计算大气污染物干沉积速率的方法。这种方法详细地考虑了植物生理特征和冠层内湍流结构对干沉积的影响，考虑了3层覆盖物对表面阻抗的影响，通过计算7个不同的阻抗因子获得干沉积速率。利用该方法计算了农田下垫面上7种不同大气污染物沉积的表面阻抗和江西(鹰洋)红壤典型地区农田下垫面上的二氧化硫(SO2)和硫酸盐(SO4^2—)的干沉积速率，估算了大气硫输入农田生态系统的干沉积通量，并与其他方法作了定量比较。结果表明：农田下垫面上SO2与SO4^2—的干沉积速率平均值分别为0．31cm／s，0．20cm／s；干沉积速率具有明显的日变化特征，一般白天大于夜间，中午时分出现最大值；月际变化也较明显，在冬季有大值；大气硫输入农田生态系统的全年干沉积通量为7．35g/m^2。     

城市下垫面空气动力学参数确定方法研究综述

地表粗糙度和零平面位移是2个重要的空气动力学参数,对于城市下垫面,粗糙元的形态非常复杂,空间分布上存在较大的非均匀性,如何确定城市下垫面地表粗糙度和零平面位移目前尚无最佳方案。概述了这2个参数适用于城市非均匀下垫面的参数化方法,包括风温观测方法和粗糙元形态学方法。探讨了各自在应用过程中存在的问题及进一步的研究工作,并且提出,通过数学模型的改进和卫星遥感技术的引入,粗糙元形态学方法会有更好的发展前景。     

藏东南地区复杂下垫面能量收支特征分析

利用公益性行业（气象）科研专项项目＂藏东南地区复杂下垫面地气交换观测研究＂在藏东南地区进行的地气交换观测实验数据,分析典型晴天和阴天条件下不同下垫面能量过程的特征及其差异。结果表明：在典型晴空状态下,不同下垫面的地表净辐射均具有明显的日变化特征,在典型阴天的情况下,不同下垫面地表净辐射日均值显著减小;在典型晴空天气下,4种类型下垫面上潜热均随着净辐射的增加而表现为增加的趋势,在典型阴天的情况下,潜热通量明显比晴空天气小;不同下垫面感热通量的日变化存在显著的差异,不同下垫面感热的变化特征在典型晴空和阴天条件下的差异不明显;不同下垫面土壤热通量与净辐射的变化趋势基本一致,阴天夜晚土壤热流的交换与晴天夜晚的差异不大,四种下垫面土壤均存在能量损失,土壤处于降温状态。     

A modified parameterization of flux-profile relationships in the surface layer using different roughness length values for heat and momentum

The surface-layer flux-profile formulae of Louis (1979), used in many atmospheric models, are modified in a simple way to allow for different values of the roughness lengths for heat and momentum. The modified set of formulae simplifies the calculation of surface-layer fluxes over most natural land surfaces, where the roughness length for momentum can be almost an order of magnitude greater than that for heat.     

Bulk Formulation of the Surface Heat Flux

An interpretive literature survey examines different approachesfor applying the bulk aerodynamic formulato predict the surface heat flux. The surface heat flux is often predicted in terms of the surface radiation temperature, which is also used to predict the upward longwave radiation and the heat flux into the soil. In models, the thermal roughness length based on the surface radiation temperature (radiometric roughness length) is often specified to be smaller than the roughness length for momentum for a number of distinct reasons. The definition of the radiometric roughness length depends on the way that the surface temperature is measured, the choice of stability functions and displacement height and inclusion of any additional resistances.Using airborne eddy correlation data collected over eight different sites including bare soil, crops and grassland and several types of forests, the radiometric roughness length is found to vary by orders of magnitude in a manner that is difficult to formulate. Alternatively, we evaluate the approach where the thermal roughness length is equated with the better behaved roughness length for momentum and the corresponding aerodynamic surface temperature is modelled in terms of the surface radiation temperature, solar radiation, and vegetation index. The influence of wind speed and soil moisture on the difference between the aerodynamic and surface radiation temperatures is also examined.     

Noah-MP陆面模式对干旱区不同下垫面水热通量模拟评估

陆面水热通量的准确模拟可为气候模式提供高质量的下边界条件，对气候模拟和预测具有重要意义。本研究基于干旱区张掖国家气候观象台2021年1月-2022年6月和大满灌区绿洲农田站2020年1-12月的观测数据，评估Noah-MP模式对干旱区荒漠和农田两种下垫面的水热通量的模拟性能。结果表明：Noah-MP模式模拟的干旱区荒漠下垫面辐射及感热和潜热通量的变化特征、峰谷值与观测值总体一致，模拟效果较好；模拟的5 cm土壤湿度对降水过程有明显的响应，表现出明显的冷暖季差异，但其模拟性能仍有待改进。Noah-MP模式对干旱区农田下垫面辐射通量及各层土壤温度的模拟效果较理想，但模拟的潜热通量及各层土壤湿度较观测值偏低，尤其在生长季模拟性能并不理想。总体而言，Noah-MP模式对干旱区荒漠下垫面水热通量的模拟性能优于农田下垫面，优化和发展模式水文过程的参数化方案，在模式中考虑人为作用，是提高干旱区陆面过程模式模拟能力的重要方向。     

海气界面动量通量算法的改进

COARE（Coupled Ocean-Atmosphere Response Experiment）算法是国际上较先进的计算海气界面通量的算法。最新的COARE 3.0算法包含TY01方案和O02方案两种考虑真实海浪状态的海面空气动力学粗糙度方案;当有效波高和谱峰周期缺省时,利用Taylor01风浪特征量参数化方案可对它们进行参数化。在此基础上引入自主设计的WHP（Wind-wave significant wave Height and dominant wave Period）风浪特征量参数化方案和三种海面空气动力学粗糙度方案,即SCOR方案、GW03方案、PYP07方案以及更简单的Andreas12、Vickers15摩擦速度算法,利用NDBC（National Data Buoy Center）浮标数据、解放军理工大学风浪流水槽实验数据,针对中高风速条件（10 m/s≤U₁₀≤25 m/s）比较上述方案对摩擦速度的预测效果。结果表明:WHP方案在COARE 3.0算法中的应用效果优于Taylor01方案,且新引入的SCOR、GW03、PYP07海面粗糙度方案对风浪特征量的观测误差敏感性更小、稳定度更高;水槽实测数据的对比结果表明,WHP方案结合SCOR海面粗糙度方案计算的摩擦速度与实测值最接近;引入其他三种实测数据的检验结果表明,原始COARE 3.0算法会低估摩擦速度,而WHP方案结合SCOR海面粗糙度方案能更准确地预测摩擦速度随10 m风速的增长趋势。      

Study of Aerodynamic Parameters on Different Underling Surfaces

Aerodynamic parameters including the zero-plane displacement (d), roughness length (20), and friction velocity (u*) on the different underlying surfaces of heavy-grazing site, medium-grazing site, light-grazing site, no-grazing site, dune, inter-dune, grassland, rice paddy site, wheat site, soybean site, and maize site have been computed based on the Monin-Obukhov similarity theory by utilizing the micrometeorologically observed data of dune and vegetation in the semi-arid area at Naiman, Inner Mongolia of China, conducted jointly by the Institute of Desert Research, Chinese Academy of Sciences and the National Institute of Agro-Environmental Sciences of Japan in 1990-1994. And their relationships between wind speed and Richardson number are analyzed. The aerodynamic characteristics of different man-made disturbed grassland ecosystems are also compared. Result shows that the vegetation coverage and the above-ground biomass decrease with the increase in man-made stress of the grassland. The roughness length for different underlying surfaces is closely related to vegetation height, above-ground biomass, and ground surface undulation, and Richardson number Ri is also its influencing factor. The friction velocity varies largely on different underlying surfaces, and it is positively proportional to wind speed and roughness length. The aerodynamic parameters of various times on the same underlying surface are different, too. Above results indicate that grassland and vegetation are of significance in preventing desertification, especially in the arid and semi-arid land ecosystems. And the results of this paper are also important for constructing the land surface physical process as well as regional climate model.     

Interpretation of Crown Radiation Temperatures of a Dense Douglas fir Forest with Similarity Theory

Infrared crown radiation temperatures as observed over a dense Douglas fir forest are analyzed in the context of similarity theory and the concept of transport resistances. As such we obtain a rather high value of the roughness length for heat, which is about equal to the roughness length for momentum. This value can be explained by the more efficient transport of heat relative to momentum in the roughness sublayer of the forest. Correcting for this effect we arrive at the classic value for homogeneous terrain of about 0.1 times the roughness length for momentum. For unstable cases the presence of enhanced mixing of heat in the roughness sublayer leads to a modified integral stability function for the dimensionless potential temperature difference between the surface and the top of the roughness sublayer. The observations give some evidence for this different stability behaviour. The analysis suggests that during daytime the radiative surface temperature and the aerodynamic surface temperature are not significantly different when used to estimate fluxes. Daytime trunk space air temperature is satisfactory parameterized with the concept of gusts and with surface renewal analysis. As such it is related to the sensible heat flux and the storage heat flux. Night time radiation temperatures at times strongly deviate from the expected behaviour based on similarity theory and the roughness length for heat, suggesting that the concept of a single surface temperature is too simple for such cases.     

Comprehensive Parametrization of Surface-Layer Transfer Coefficients for Use in Atmospheric Numerical Models

A new non-iterative bulk parametrization for surface-layer transfer coefficients for momentum and heat is presented. It is applicable for a wide range of aerodynamic and thermal roughness lengths, and includes the effect of the roughness sublayer. As a consequence, the non-iterative method is suitable for every surface type, especially for urban surfaces for which existing non-iterative parametrizations fail. The analytical approximation compares very well with an iterative approach. Our method can be easily implemented in atmospheric numerical models that already employ a non-iterative approach.     

Flux–Profile Relationships Over a Fetch Limited Beech Forest

The influence of an internal boundary layer and a roughness sublayer on flux–profile relationships for momentum and sensible heat have been investigated for a closed beech forest canopy with limited fetch conditions. The influence was quantified by derivation of local scaling functions for sensible heat flux and momentum (_h and _m) and analysed as a function of atmospheric stability and fetch. For heat, the influences of the roughness sublayer and the internal boundary layer were in agreement with previous studies. For momentum, the strong vertical gradient of the flow just above the canopy top for some wind sectors led to an increase in _m, a feature that has not previously been observed. For a fetch of 500 m over the beech forest during neutral atmospheric conditions, there is no height range at the site where profiles can be expected to be logarithmic with respect to the local surface. The different influence of the roughness sublayer on _h and _m is reflected in the aerodynamic resistance for the site. The aerodynamic resistance for sensible heat is considerably smaller than the corresponding value for momentum.