复杂下垫面地表土壤热通量算法改进:以黑河流域为例

地表土壤热通量是地表能量平衡的重要组成部分,对地表蒸散发的估算至关重要.利用土壤温湿度廓线观测资料基于热扩散方程计算地表土壤热通量,并通过冻土融化前后土壤液态水含量变化估算土壤含冰量,分析了土壤含冰量对土壤热通量的影响,旨在分析黑河流域典型下垫面(高山草地、农田和森林)地表土壤热通量的时空变化特征.研究结果表明:(1)黑河流域不同下垫面的地表土壤热通量有明显的日变化差异,日最大值时刻提前净辐射通量几分钟至几小时不等,这与土壤质地、湿度、热属性和植被覆盖度有关;(2)净辐射通量有显著的季节变化,一般夏季达最大值,冬季最小,地表土壤热通量也有明显的季节变化,但并不总是与净辐射通量变化保持一致,春季达最大值,夏季由于植被覆盖的原因反而降低;(3)地表土壤热通量占净辐射通量的比例因季节及下垫面不同而有差异,1月份月平均比值分别为:阿柔25.6%、盈科22.9%和关滩4.3%,7月份月平均比值分别为:阿柔2.3%、盈科1.6%和关滩0.3%;(4)冬季考虑了冰的热容量使得土壤热容量增加,土壤热存储增加,从而由热扩散方程计算的地表土壤热通量增加,使得能量平衡闭合率提高了4.3%.     

塔里木河下游荒漠河岸林蒸散规律及其关键控制机制

干旱区内陆河下游荒漠河岸林带蒸散耗水规律及其控制机制是内陆河流域下游水资源管理的重要科学基础.本研究对塔里木河下游两种典型河岸林群落胡杨群落和柽柳群落的地下水、土壤水和地表蒸散过程进行了连续观测.结果发现河岸林带生长季的蒸散过程季节趋势受到植被物候期的控制,非生长季蒸散微弱;蒸散日过程是大气要素综合影响的结果,与参考蒸散显著线性相关;蒸散的空间格局则受到植被叶面积指数的控制,植被盖度越大,蒸散量越大;地下水是地表蒸散的水分来源,浅层土壤水参与水循环微弱;胡杨和柽柳蒸散过程特征一致,但胡杨耗水能力高于柽柳.进一步分析表明,植物蒸腾是地表蒸散的主体,土壤蒸发占地表蒸散的比例很小;地下水位是控制干旱区荒漠河岸林水循环的关键因素,它通过影响植物生长和植被盖度来影响荒漠河岸林带地表蒸散过程和蒸散量;采用地下水位,而不是土壤含水量定量表达荒漠河岸林水分胁迫过程更为合适.基于揭示的水分运动和蒸散规律,系统描绘了一个表达干旱区内陆河下游河岸林带水循环过程控制机制的概化框架,并描述了一个简单实用,能满足内陆河水资源管理需要的荒漠河岸林需水量估算模型.     

北京城市及周边热岛日变化及季节特征的卫星遥感研究与影响因子分析

以EOS-MODIS遥感信息反演的地表温度、土地覆盖类型、植被覆盖、地表蒸散,结合常规气象资料,并采用GIS空间分析技术和多元统计相关,对北京城市及周边2001年城市热岛(UHI)空间分布的季节规律和日变化及影响因子进行研究.分析北京地区的土地覆盖、地形高程、植被绿地状况、城市和郊区地表蒸散与热岛时空分布状况的关系.揭示出北京UHI主要特征为:(ⅰ)北京城市下垫面的高热容和密集建筑物的多次发射,加之北京特殊的三面环山地形特征,使得北京城区一年四季均存在明显的热岛分布,并以夏季最为明显,UHI与城市结构的轮廓相一致;北京城区与地势相对平坦的近郊区的地表温度差异在4~6℃左右,与地势较高的西北远郊区的地表温度差异在8~10℃左右;(ⅱ)北京地区日间和夜晚的UHI的季节分布和程度不同,以夜间UHI明显;夏季白天郊区地表比城区蒸散量大,潜热交换明显,反映出城市与郊区的温度差异显著;(ⅲ)地表覆盖类型对UHI的效应明显,北京地区植被绿地状况与UHI呈现明显反相关分布;夏季地表NDVI与下垫面的温度散点图的回归方程的负相关系数的平方R2达到0.6481,即植被覆盖好,则UHI不明显;揭示出植被绿地对降低UHI具有重要的作用.大范围的绿地建设能有效降低UHI.     

植被变化对西北地区陆气耦合强度的影响

西北地区地处欧亚大陆腹地,生态系统对于气候变化和人为影响十分敏感,同时该区也是湿润的东亚季风区与干燥的中亚干旱区的过渡区域,陆气相互作用比较强烈.本文对西北地区植被变化对当地的陆气耦合强度及其与之相关的地表水文过程的影响进行了分析研究,并且找出适于增加植被以缓解西北地区荒漠化趋势的最具成效的地区.本文利用美国国家大气科学研究中心(NCAR,National Center for Atmospheric Research)研制的通用大气模式CAM3(Community Atmosphere Model Version 3)对西北地区植被变化的影响进了数值模拟.本文共设计了三个试验,使用正常地表植被覆盖的参考试验,地表下垫面变为裸土的去植被试验和植被增加的生态环境好转试验.首先,本文对西北地区植被变化对于当地降水量、地表水分盈余量、径流量、地表土壤含水量等地表水文变量的影响进行了分析研究.然后对西北地区植被变化对当地的陆气耦合强度的影响进了分析研究,陆气耦合强度是衡量局地陆气相互作用强弱程度的一个新标准,基于计算年降水量与蒸散量的协方差与降水量方差之比而得到.它利用观测数据或模式输出数据,计算起来简便容易,物理意义明确清晰,陆气相互作用越强烈的地区,其陆气耦合强度也越高.最后,本文计算了一个蒸散-水汽通量散度指数来衡量植被变化对局地蒸散与大气水汽通量散度的影响,其在一定程度上反应了植被变化对局地陆气相互作用和大尺度大气环流输送作用的影响,也可以视为一个评估人为生态环境工程效果的指标.西北地区陆气耦合强度由东南向西北递增.去植被之后,西北地区降水与蒸发普遍减少,其中在东南部区域,地表径流增加约10~40mm,渗流量与地表土壤含水量分别减少约40~80mm和5~20mm3·mm-3,陆气耦合强度上升,这有可能导致水土流失,不利于当地植被的恢复.生态环境好转之后,内陆地区降水与蒸发明显增加,但地表盈余水分有所减少,主要原因是蒸散量相较于降水量增加的更多.其中在沙漠戈壁区边缘的新疆南部与内蒙西部,渗流量与地表土壤含水量分别上升约5~20mm和5~20mm3·mm-3,陆气耦合强度降低,蒸散-水汽通量散度指数较高,这可能主要是由于植被变化对局地陆气相互作用的改变而造成的.植被对于西北地区地表水文过程有着明显的影响,植被的存在能加速西北地区地表水文循环过程,减小陆面蒸散的变化,降低陆气耦合强度.在有限的人力与财力条件下,集中力量在在沙漠戈壁区边缘的新疆南部与内蒙西部适当种植灌木与青草并防止过度放牧,能有效降低当地陆气耦合强度,缓解西北地区荒漠化加剧的趋势.本文下一步还需考虑如模式地表植被数据与真实情况的差异性,海洋因素变化对于植被变化的反馈,以及进行集合实验来增加研究结果的可靠性.     

我国黄土高原地区陆面能量的空间分布规律及其与气候环境的关系

我国黄土高原地区地处夏季风边缘,分布在气候和生态过渡带,气候环境的空间差异很大,对陆面能量的空间分布格局影响非常显著.然而,由于受该地区陆面过程观测站点较少的局限,对整个黄土高原区域陆面能量的空间分布规律及其影响机制的认识十分有限.在对CLM模式模拟的陆面能量平衡分量资料进行试验验证的基础上,利用CLM模式模拟的近30年黄土高原地区陆面能量平衡分量资料,分析了该地区近30年平均陆面能量平衡分量的空间变化特征以及与最干燥年和最湿润年的差异,研究了陆面能量平衡分量空间分布与经、纬度和海拔高度等地理因素及降水和气温等气候因子之间的关系.发现,黄土高原地区陆面能量平衡分量空间差异非常显著,地表净辐射和感热通量由南至北增加,潜热通量和土壤热通量从东南向西北减少;空间最干格点和最湿格点之间的地表感热、潜热和土壤热通量几乎相差1倍左右,地表能量分配由最干格点的感热通量主导转变为最湿格点的感热和潜热平分秋色;年际干湿波动对地表能量平衡分量的影响也相当显著,对感热和潜热通量的改变幅度最大接近30%.而且,经纬度和海拔高度等地理因素及温度和降水等气候要素均与陆面能量平衡分量空间分布有一定的相关性,但地表净辐射与海拔高度和纬度的关系更密切,感热通量与降水和纬度的关系更密切,而潜热和土壤热通量只与降水的关系比较密切.     

利用NOAA/AVHRR数据获取地表特征参数的方法研究

地表特征参数的正确与否直接影响到区域陆面蒸散量估算精度,因此在区域非均匀陆面蒸散研究中,地表特征参数的获取方法是一个值得探讨的问题．与传统的方法相比,卫星遥感技术在求取地表特征参数时有其独特的优势．NOAA气象卫星AVHRR资料以其时间分辨率高、覆盖面广、价格低廉等优点广泛应用于非均匀陆面蒸散研究和应用中．本文建立了NOAA／AVHRR计算地表特征参数的参数化模型,选取中国东北松嫩平原西部地区2000年7月8日的AVHRR资料,试算了研究区地表温度、地表发射率、地表反照率、NDVI等主要地表特征参数,并且参照2000年研究区土地利用数据对各参数的空间分布特征及合理性进行了分析。     

改进的区域缺水遥感监测方法

植物水分与土壤水分状况密切相关, 提取地表缺水指数是研究区域缺水的一种有效途径. 选择不同退化程度的草地为研究对象, 通过建立基于亚像元尺度的双层蒸散模型, 成功计算了区域的地表缺水指数(SWDI); 考虑到半干旱地区植被冠层水主要由土壤水提供, 但植被供水状况与植被水分实际状况之间存在着一定的“滞后”, 鉴于此, 试图建立地表缺水指数与植被冠层水分含量、表层土壤水分含量(0~20 cm)之间的定量关系, 从而更准确、及时地反演区域的表层土壤水分含量(0~20 cm). 通过本次研究探讨一种直接应用遥感技术研究半干旱地区区域缺水的实用方法.     

热带季节雨林碳通量年变化特征及影响因子初探

为深入分析西双版纳热带季节雨林碳通量年变化特征及其与各种因子的关系,通过西双版纳主要自然生态系统(热带季节雨林)2003～2004年林冠-大气间和近地层的碳通量以及不同覆盖状况下的地表碳通量(土壤呼吸)的长期观测,并结合植物光合作用、叶面积指数、凋落物和凋落物分解速率以及温度、辐射等常规气象的测定,对热带季节雨林碳通量的年变化特征及其影响因子进行了综合的分析与研究.结果表明热带季节雨林的碳通量表现出和其它热带雨林不同的特征,在干季(11～4月)的林冠-大气间碳通量为负值,森林生态系统呈现碳汇效应;而在雨季(5～10月)表现出较弱的碳源效应;森林生态系统碳通量具有明显的日变化特征,在白天呈现碳汇效应,而夜间为碳源效应,并且干季昼间碳通量较大;雨季较小;夜间则呈现相反趋势.林冠植物在昼间具有较强的光合作用,对昼间林冠-大气间碳通量有较大的贡献;林冠和植物林内低矮植物的光合速率均与林冠-大气间碳通量有显著的相关关系,而林内低矮植物的光合速率与林内近地层碳通量仅在于热季存在显著的相关关系.林内不同覆盖状态的地表碳通量具有明显的季节变化;林冠-大气间碳通量与地表碳通量同样具有较好的相关性,地表碳通量是导致热带季节雨林生态系统碳通量呈现特殊分布的主要因子;另外,林冠-大气间碳通量与凋落物量、凋落物分解速率、降水量、土壤含水率和土壤温度均表现出较好的相关性.初步的统计表明,西双版纳热带季节雨林林冠-大气间碳通量在不同季节呈现不同的汇/源效应,在总体上表现为一个较弱的碳汇.     

一种基于土壤温湿资料计算地表土壤热通量的温度预报校正法

分析大气边界层观测站地表能量平衡需要估计其地表土壤热通量。发展了一种由多层土壤温度和湿度观测资料估算土壤热通量的新方法。该方法首先求解一维热扩散方程得到土壤温度的基本廓线,然后校正所求温度廓线与观测值的偏差,最后积分温度廓线得到土壤各层的热通量。与众多的方法不同,该方法不需要事先给定不易准确测量和推求的热传导（或热扩散）系数值。通过与实测资料对比、模型合成数据试验、以及敏感性分析等,表明该方法的计算结果稳定可靠,对土壤表层数厘米深度内有无观测资料也不敏感。此外,指出热流板可准确测量热通量的方向和相位,但所测通量值的误差常较大。     

甘肃岷县漳县MS6.6地震前后岷县植被覆盖度动态变化

利用Landsat8遥感影像数据,应用像元二分模型估算了岷县植被覆盖度,分析了岷县漳县地震前后岷县植被覆盖度动态变化,并结合高程和坡度数据,研究了地形因素对植被受损和恢复的影响,结果表明:3个时段岷县植被覆盖度以中高和高植被覆盖度为主,反映出研究区植被覆盖状况良好。①空间分布上来看,中高和高植被覆盖度主要分布在岷县西南和东南部的高海拔地区,低植被覆盖区由岷县县城所在地向西北方向两侧沿居民地延展;②2013-2014年平均植被覆盖度下降了0.09,受损面积达到96.53%,主要为高植被覆盖度植被受损,植被受损在高程上主要分布在2 344~3 283 m,坡度5°~35°范围是植被受损的主要区域;③2014-2017年平均植被覆盖度回升至0.43,恢复面积达2 901.93 km^2,以高植被覆盖度植被恢复为主,其次为低植被覆盖度;植被恢复在高程2 031~2 970 m,呈增长趋势,随后开始下降,在坡度0°~25°范围内植被恢复面积呈增加趋势,随后降低。     

Modeling Water and Heat Balance Components for Large Agricultural Region Utilizing Information from Meteorological Satellites

The method has been developed to evaluate water and heat balance components for vegetation covered area of regional scale based on the refined physical-mathematical model of vertical water and heat exchange between land surface and atmosphere (Land Surface Model, LSM) for vegetation season adapted to satellite information on land surface and meteorological conditions. The LSM is accommodated for utilizing satellite-derived estimates of vegetation and meteorological characteristics as model parameters and input variables. Estimates of these characteristics presented as distributions of their values over the study area have been obtained from AVHRR/NOAA, MODIS/EOS Terra and Aqua, SEVIRI/Meteosat-9, -10 data. To build such estimates methods and technologies have been developed and refined using results of thematic processing measurement data from these sensors. Among them the original Multi Threshold Method (MTM) has been developed and tested to calculate daily precipitation sums using rainfall intensity estimates retrieved from AVHRR and SEVIRI data with subsequent replacement of ground-measured rainfall amounts by these daily rainfalls. All technologies have been adapted to the study area with square of 227300 km² being the part of the Central Black Earth Region of European Russia. Developed earlier procedures of utilizing satellitederived estimates of vegetation and meteorological characteristics (including precipitation) in the model have been refined and verified. Final result of modeling is the fields of soil water content, evapotranspiration and other water and heat balance components of the region under study for years 2012–2014 vegetation seasons.     

Response of water and energy exchange to the environmental variable in a desert‐oasis wetland of Northwest China

A case study on a desert‐oasis wetland ecosystem in the arid region of Northwest China measured the seasonal and interannual variation in energy partitioning and evapotranspiration to analyse the response of water and energy exchange on soil moisture, groundwater, and environmental variables. Energy partitioning showed a clear seasonal and interannual variability, and the process of water and energy exchange differed significantly in the monthly and interannual scales. The net radiation was 7.31 MJ m^?2· day^?1, and sensible heat flux accounted for 50.42% of net radiation in energy fluxes, 40.56% for latent heat flux, and 9.02% for ground heat flux. The parameters in energy fluxes were best described by a unimodal curve, whereas sensible heat flux followed a bimodal curve. Variations in the daily evapotranspiration and crop evapotranspiration also exhibited a single peak curve with annual values of 569.84 and 644.47 mm, respectively. Canopy conductance averaged 20.77 ± 13.75 mm s^?1 and varied from 0.16 to 83.96 mm s^?1 during the two hydrological years. The variation in water and energy exchange reflected environmental conditions and depended primarily on vapour pressure deficit, net radiation, soil moisture, and water depth. Although the effects of precipitation on evapotranspiration showed that the response of this ecosystem to climate changes was not obvious, the variation of air temperatures had a strong influence on evapotranspiration, resulting in a significant increase in evapotranspiration (R = 0.730; P < 0.01). Copyright © 2013 John Wiley & Sons, Ltd.     

Estimation of seasonal crop water consumption in a vineyard using Bowen ratio‐energy balance method

Information about seasonal crop water consumption is useful to develop the appropriate irrigation scheme. Measurements of energy balance components using the Bowen ratio method were made for a complete growing season at a vineyard in the arid region of northwest China. Vine in the experiment was furrow‐irrigated using a trellis system. The measured evapotranspiration was compared with estimates using the soil water balance method. It is shown that the Bowen ratio method provided accurate estimates of evapotranspiration from the vineyard and this requires that the Bowen ratio system is appropriately installed. The energy balance components showed typical diurnal pattern with peaks that occurred around the midday, except for the ground heat flux which delayed its peak by 2–3 h. The sensible heat flux was greater than the latent heat flux and followed the net radiation closely. The ratio of the latent heat flux to net radiation was low in the early growing season and increased over time. Under the limited irrigation experienced in the vineyard, the latent heat flux was controlled by available soil moisture and the total evapotranspiration in the growing season was 253 mm. The seasonal progression of the crop coefficient is similar to that reported in the literature, with the maximum occurring during the month of September. The crop coefficient can be estimated as a non‐linear function of day of year (DOY) and used to estimate evapotranspiration from vineyards in the region. Copyright © 2007 John Wiley & Sons, Ltd.     

Verification of land surface evapotranspiration estimation from remote sensing spatial contextual information

Estimation of evapotranspiration (ET) is of great significance in modeling the water and energy interactions between land and atmosphere. Negative correlation of surface temperature (T_s) versus vegetation index (VI) from remote sensing data provides diagnosis on the spatial pattern of surface soil moisture and ET. This study further examined the applicability of T_s–VI triangle method with a newly developed edges determination technique in estimating regional evaporative fraction (EF) and ET at MODIS pixel scale through comparison with large aperture scintillometer (LAS) and high‐level eddy covariance measurements collected at Changwu agro‐ecological experiment station from late June to late October, 2009. An algorithm with merely land and atmosphere products from MODIS onboard Terra satellite was used to estimate the surface net radiation (R_n) and soil heat flux. In most cases, the estimated instantaneous R_n was in good agreement with surface measurement with slight overestimation by 12 W/m². Validation results from LAS measurement showed that the root mean square error is 0.097 for instantaneous EF, 48 W/m² for instantaneous sensible heat flux, and 30 W/m² for daily latent heat flux. This paper successfully presents a miniature of the overall capability of T_s–VI triangle in estimating regional EF and ET from limited number of data. For a thorough interpretation, further comprehensive investigation needs to be done with more integration of remote sensing data and in‐situ surface measurements. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulation of crop growth and energy and carbon dioxide fluxes at different time steps from hourly to daily

Understanding the exchange processes of energy and carbon dioxide (CO₂) in the soil–vegetation–atmosphere system is important for assessing the role of the terrestrial ecosystem in the global water and carbon cycle and in climate change. We present a soil–vegetation–atmosphere integrated model (ChinaAgrosys) for simulating energy, water and CO₂ fluxes, crop growth and development, with ample supply of nutrients and in the absence of pests, diseases and weed damage. Furthermore, we test the hypotheses of whether there is any significant difference between simulations over different time steps. CO₂, water and heat fluxes were estimated by the improving parameterization method of the coupled photosynthesis–stomatal conductance–transpiration model. Soil water evaporation and plant transpiration were calculated using a multilayer water and heat‐transfer model. Field experiments were conducted in the Yucheng Integrated Agricultural Experimental Station on the North China Plain. Daily weather and crop growth variables were observed during 1998–2001, and hourly weather variables and water and heat fluxes were measured using the eddy covariance method during 2002–2003. The results showed that the model could effectively simulate diurnal and seasonal changes of net radiation, sensible and latent heat flux, soil heat flux and CO₂ fluxes. The processes of evapotranspiration, soil temperature and leaf area index agree well with the measured values. Midday depression of canopy photosynthesis could be simulated by assessing the diurnal change in canopy water potential. Moreover, the comparisons of simulated daily evapotranspiration and net ecosystem exchange (NEE) under different time steps indicated that time steps used by a model affect the simulated results. There is no significant difference between simulated evapotranspiration using the model under different time steps. However, simulated NEE produces large differences in the response to different time steps. Therefore, the accurate calculation of average absorbed photosynthetic active radiation is important for the scaling of the model from hourly steps to daily steps in simulating energy and CO₂ flux exchanges between winter wheat and the atmosphere. Copyright © 2007 John Wiley & Sons, Ltd.     

Examination of Morton's CRAE model for estimating daily evaporation from field-sized areas

Morton's complementary relationship areal evapotranspiration (CRAE) model was originally designed to provide regional estimates of monthly evapotranspiration. Often, however, hydrologists and others require estimates of evapotranspiration for field-sized land units under a specific land use, for shorter intervals of time. This paper examines CRAE with respect to the algorithms used to describe different terms and its applicability to reduced spatial and temporal scales.

Daily estimates by CRAE of atmospheric radiation fluxes during the summer months are compared with monitored values. It is shown that errors in estimation of the extra-terrestrial flux, the transmittancy of clouds to short-wave radiation, the surface albedo and the net long-wave flux result in standard deviations of the difference between ‘modelled’ and ‘measured’ net all-wave radiation for 1-, 5- and 10-day periods of 2.58, 1.8 and 1.50 MJm⁻² day⁻¹ respectively.

The assumption in CRAE that the vapour transfer coefficient is independent of wind speed may lead to appreciable error in computing evapotranspiration. A procedure for incorporating a wind correction factor is described and the improvement in estimating regional evaporation is illustrated.

Comparisons of evapotranspiration estimates by CRAE and measurements obtained from soil moisture and precipitation observations in the semi-arid, cold-climate Prairie region of western Canada demonstrate that the assumptions that the soil heat flux and storage terms are negligible, lead to large overestimation by the model during periods of soil thaw.     

Controls on evapotranspiration from jack pine forests in the Boreal Plains Ecozone

The exchanges of water, energy and carbon between the land surface and the atmosphere are tightly coupled, so that errors in simulating evapotranspiration lead to errors in simulating both the water and carbon balances. Areas with seasonally frozen soils present a particular challenge due to the snowmelt-dominated hydrology and the impact of soil freezing on the soil hydraulic properties and plant root water uptake. Land surface schemes that have been applied in high latitudes often have reported problems with simulating the snowpack and runoff. Models applied at the Boreal Ecosystem Research and Monitoring Sites in central Saskatchewan have consistently over-predicted evapotranspiration as compared with flux tower estimates. We assessed the performance of two Canadian land surface schemes (CLASS and CLASS-CTEM) for simulating point-scale evapotranspiration at an instrumented jack pine sandy upland site in the southern edge of the boreal forest in Saskatchewan, Canada. Consistent with past reported results, these models over-predicted evapotranspiration, as compared with flux tower observations, but only in the spring period. Looking systematically at soil properties and vegetation characteristics, we found that the dominant control on evapotranspiration within these models was the canopy conductance. However, the problem of excessive spring ET could not be solved satisfactorily by changing the soil or vegetation parameters. The model overestimation of spring ET coincided with the overestimation of spring soil liquid water content. Improved algorithms for the infiltration of snowmelt into frozen soils and plant-water uptake during the snowmelt and soil thaw periods may be key to addressing the biases in spring ET.     

An improved algorithm to estimate the surface soil heat flux over a heterogeneous surface: A case study in the Heihe River Basin

Surface soil heat flux(G0) is an indispensable component of the surface energy balance and plays an important role in the estimation of surface evapotranspiration(ET). This study calculated G0 in the Heihe River Basin based on the thermal diffusion equation, using the observed soil temperature and moisture profiles, with the aim to analyze the spatial-temporal variations of G0 over the heterogeneous area(with alpine grassland, farmland, and forest). The soil ice content was estimated by the difference in liquid soil water content before and after the melting of the frozen soil and its impact on the calculation of G0 was further analyzed. The results show that:(1) the diurnal variation of G0 is obvious under different underlying surfaces in the Heihe River Basin, and the time when the daily maximum value of G0 occurs is a few minutes to several hours earlier than that of the net radiation flux, which is related to the soil texture, soil moisture, soil thermal properties, and the vegetation coverage;(2) the net radiation flux varies with season and reaches the maximum in summer and the minimum in winter, whereas G0 reaches the maximum in spring rather than in summer, because more vegetation in summer hinders energy transfer into the soil;(3) the proportions of G0 to the net radiation flux are different with seasons and surface types, and the mean values in January are 25.6% at the Arou site, 22.9% at the Yingke site and 4.3% at the Guantan site, whereas the values in July are 2.3%, 1.6% and 0.3%, respectively; and(4) G0 increases when the soil ice content is included in thermal diffusion equation, which improves the surface energy balance closure by 4.3%.     

Modelling plant canopy effects on annual variability of evapotranspiration and heat fluxes for a semi‐arid grassland on the southern periphery of the Eurasian cryosphere in Mongolia

The ability to predict vegetation cover effects on thermal/water regimes can enhance our understanding of canopy controls on evapotranspiration. The Simultaneous Heat and Water (SHAW) model is a detailed process model of heat and water movement in a snow–residue–soil system. This paper describes provisions added to the SHAW model for vegetation cover and simulation of heat and water transfer through the soil–plant–air continuum. The model was applied to four full years (May 2003–April 2007) of data collected on sparse grassland at Nalaikh in north‐eastern Mongolia. Simulated soil temperature and radiation components agreed reasonably well with measured values. The absolute differences between simulated and measured soil temperatures were larger at both the surface layer and deeper layer, but relatively smaller in the layer from 0·8 to 2·4 m. Radiation components were mimicked by the SHAW model with model efficiency (ME) reaching 0·93–0·72. Latent and sensible heat fluxes were simulated well with MEs of 0·93 and 0·87, respectively. The vegetation control on evapotranspiration was investigated by sensitivity experiments of model performance with changing leaf area index (LAI) values but constant of other variables. The results suggest that annual evapotranspiration ranged from 16 to ? 22% in response to extremes of doubled and zero LAI. Copyright © 2010 John Wiley & Sons, Ltd.