Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies

In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (F_T) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating F_T to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (S_T) and F_T have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if S_T and F_T are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.     

A bulk similarity approach in the atmospheric boundary layer using radiometric skin temperature to determine regional surface fluxes

Profiles of wind velocity and temperature in the outer region of the atmospheric boundary layer (ABL) were used together with surface temperature measurements, to determine regional shear stress and sensible heat flux by means of transfer parameterizations on the basis of bulk similarity. The profiles were measured by means of radiosondes and the surface temperatures by infrared radiation thermometry over hilly prairie terrain in northeastern Kansas during the First ISLSCP Field Experiment (FIFE). In the analysis, the needed similarity functions were determined and tested; the main scaling variables used for the ABL were h _i , the height of the convectively mixed layer, and V _a and _a, the wind speed and potential temperature averaged over the mixed layer. Good agreement (r = 0.80) was obtained between values of friction velocity u _* determined by this ABL bulk similarity approach and those obtained by Monin-Obukhov similarity in the surface sublayer. Similarly, values of surface flux of sensible heat H determined by this method compared well (r = 0.90) with the regional means measured at six ground stations. The corresponding regional evaporation values, determined with the energy budget equation, also compared favorably (r = 0.94).     

Radiative surface temperature and energy balance of a wheat canopy

The sensible heat loss from a stand of winter wheat was calculated from radiometric measurements of crop surface temperature, measurements of air temperature, and an atmospheric resistance to momentum transfer; corresponding latent heat flux was obtained through the energy balance equation. These estimates of sensible and latent heat were compared with fluxes from the Bowen Ratio method. When radiative temperature was derived using a measured canopy emissivity of 0.98, calculations of sensible heat flux were systematically 50–100 W m^-2 less than Bowen Ratio values. The two techniques agreed more closely when an apparent emissivity of 0.96 was used with an apparent reflectivity of 0.03. The mean difference between the estimates of latent heat flux was then -16 ± 32 W m^-2.The surface temperature method showed less systematic error in comparison with the Bowen Ratio values than did estimates using the aerodynamic method.On leave from: University of Nottingham, School of Agriculture, Loughborough LE12 5RD.     

Regional surface fluxes from satellite-derived surface temperatures (AVHRR) and radiosonde profiles

Radiometric surface temperatures, derived from measurements by the AVHRR instrument aboard the NOAA-9 and the NOAA-11 polar orbiting satellites, were used in combination with wind velocity and temperature profiles measured by radiosondes, to calculate surface fluxes of sensible heat. The measurements were made during FIFE, the First ISLSCP (International Satellite Land Surface Climatology Project) Field Experiment, in a hilly tall grass prairie area of northeastern Kansas. The method of calculation was based on turbulent similarity formulations for the atmospheric boundary layer. Good agreement (r = 0.7) was obtained with reference values of sensible heat flux, taken as arithmetic means of measurements with the Bowen ratio method at six ground stations. The values of evaporation (latent heat fluxes), derived from these sensible heat fluxes by means of the energy budget, were also in good agreement (r = 0.94) with the corresponding reference values from the ground stations.     

Heat flux parameterization for sparse and dense grasslands with the analytical land-atmosphere radiometer model (alarm)

The dependence of radiometric surface temperature (_s) on view angle and the unclear definition of the aerodynamic temperature, which is the temperature that gives the correct sensible heat flux estimate at a specified roughness length, bring about a challenge in estimating sensible heat flux from _s. An analytical-land-atmosphere-radiometer model (ALARM) has been developed to convert _s taken at any zenith view angle to a clearly defined equivalent isothermal surface temperature, _i, at a defined scalar roughness length. ALARM is an analytical model based on K-theory that links the foliage temperature profile to the radiometric surface temperature and the temperature felt by the turbulent lower atmosphere. ALARM has previously been applied with slightly different values of its parameters to several grassland sites of varying canopy density. Our objective in this study was to apply ALARM to these and to one additional dataset with a single parameterization. When compared to the reference (measured) values of sensible heat flux H, ALARM estimates of H had root mean square errors of about 35 W m^-2. These results were comparable to those from two other simple canopy models also tested with these datasets.     

A lagrangian random-walk model for simulating water vapor,CO2 and sensible heat flux densities and scalar profiles over and within a soybean canopy

An integrated canopy micrometeorological model is described for calculating CO₂, water vapor and sensible heat exchange rates and scalar concentration profiles over and within a crop canopy. The integrated model employs a Lagrangian random walk algorithm to calculate turbulent diffusion. The integrated model extends previous Lagrangian modelling efforts by employing biochemical, physiological and micrometeorological principles to evaluate vegetative sources and sinks. Model simulations of water vapor, CO₂ and sensible heat flux densities are tested against measurements made over a soybean canopy, while calculations of scalar profiles are tested against measurements made above and within the canopy. The model simulates energy and mass fluxes and scalar profiles above the canopy successfully. On the other hand, model calculations of scalar profiles inside the canopy do not match measurements.The tested Lagrangian model is also used to evaluate simpler modelling schemes, as needed for regional and global applications. Simple, half-order closure modelling schemes (which assume a constant scalar profile in the canopy) do not yield large errors in the computation of latent heat (LE) and CO₂ (F _c ) flux densities. Small errors occur because the source-sink formulation of LE andF _c are relatively insensitive to changes in scalar concentrations and the scalar gradients are small. On the other hand, complicated modelling frames may be needed to calculate sensible heat flux densities; the source-sink formulation of sensible heat is closely coupled to the within-canopy air temperature profile.     

Intra-field variability of scalar flux densities across a transition between a desert and an irrigated potato field

This paper reports on measurements of sensible and latent heat and CO₂ fluxes made over an irrigated potato field, growing next to a patch of desert. The study was conducted using two eddy correlation systems. One measurement system was located within the equilibrium boundary layer 800 m downwind from the edge of the potato field. The other measurement system was mobile and was placed at various downwind positions to probe the horizontal transition of vertical scalar fluxes. Latent (LE) and sensible (H) heat fluxes, measured at 4 m above the surface, exhibited marked variations with downwind distance over the field. Only after the fetch to height ratio exceeded 75 to 1 didLE andH become invariant with downwind distance. When latent and sensible heat fluxes were measured upwind of this threshold, significant advection of humidity-deficit occurred, causing a vertical flux divergence ofH andLE.The measured fluxes of momentum, heat, and moisture were compared with predictions from a second-order closure two-dimensional atmospheric boundary layer model. There is good agreement between measurements and model predictions. A soil-plant-atmosphere model was used to examine nonlinear feedbacks between humidity-deficits, stomatal conductance and evaporation. Data interpretation with this model revealed that the advection of hot dry air did not enhance surface evaporation rates near the upwind edge of the potato field, because of negative feedbacks among stomatal conductance, humidity-deficits, andLE. This finding is consistent with results from several recent studies.     

Spatial Variability of Both Turbulent Fluxes and Temperature Profiles in an Urban Roughness Layer

The spatial variability of both turbulent flow statistics in the roughness sublayer (RSL) and temperature profiles within and above the canopy layer (CL) were investigated experimentally in a densely built-up residential area in Tokyo, Japan. Using five towers with measuring devices, each tower isolated from the others by at least 200 m, we collected high-frequency measurements of velocity and temperature at a height z=1.8 z _H, where z _H, the mean building height in the area, is 7.3 m. Also, temperature profiles were measured from z=0.4 to 1.8 z _H. The ‘areal mean’ geometric parameters that were obtained for the areas within 200 m of each tower were fairly homogeneous among the tower sites. The main results are as follows: (1) The spatial variability of all RSL turbulent statistics, except the sensible heat flux, was comparable to that reported in a pine forest. Also, the variability decreased with increasing friction velocity. (2) The spatial variability of the RSL sensible heat flux was larger than that reported in a pine forest. Also, the variability depended on the time of the day and became larger in the morning. The difference among the sites was well related to the areal fraction of vegetation. (3) The spatial variability of the CL temperature profile depended on the time of the day and became larger in the morning. Nevertheless, the spatial standard deviation of CL temperature was always below 0.7 K. (4) It is suggested that the “warming-up” process in the morning when heat storage is dominant increases the spatial variation of RSL sensible heat flux and CL temperature according to the local properties around each tower and the variation decreases once there is further convective mixing in the midday     

珠海凤凰山森林下垫面干季和湿季气象要素的对比分析与动量和感热交换系数的参数化研究

利用珠海凤凰山陆气相互作用观测塔站2014年11月至2016年5月的观测数据,对比分析了干湿季森林下垫面能量通量和气象要素的变化特征,分析了在不同稳定度下3个风向范围（315°～45°、45°～135°和135°～225°）的动量和感热交换系数随冠层表面风速的变化特征,并对动量和感热交换系数进行了参数化研究。结果表明:干季感热和潜热通量值相当,湿季潜热远大于感热。干季和湿季的夜晚都出现负感热现象,感热从大气向森林输送。相对湿度的变化幅度大,与该地气象状况密切相关,相对湿度的垂直梯度夜晚较大,白天较小。干季的气温垂直梯度比湿季的明显。风速在冬季变化平缓,夏季变化剧烈,低层风速随高度变化梯度明显,高层较紊乱。各高度风向差异不大。中性和近中性状态下,在风向为315°～45°、45°～135°和135°～225°时,动量交换系数C_dn分别为0.05、0.0055和0.022,感热交换系数C_hn分别为0.0055、0.003和0.004。在稳定和不稳定状态下,动量交换系数C_d、感热交换系数C_h随冠层表面风速v明显发生变化,稳定条件下,C_d、C_h随v的增大而增大;不稳定条件下,C_d、C_h随v的增大而减小。分不同风向对森林冠层C_d、C_h在稳定和不稳定条件下与v的关系进行了拟合,得到了参数化公式。     

简单生物圈模式 (SiB2) 中湍流能量通量对近地层两类阻抗系数的敏感性研究

陆—气之间的能量交换是通过近地层湍流热量和水汽通量来实现的.以往的研究表明近地层阻抗对湍流能量通量有着不同程度的影响,但是关于陆面模式中阻抗系数的取值范围却始终没有统一的标准.为了深入了解简单生物圈模式 (SiB₂) 中近地层阻抗系数取值变化对湍流能量通量的影响,我们以那曲站为例,分别采用传统的逐个因子分析法和考虑参数间相互作用的部分因子分析法定量地研究了夏季该观测站近地层湍流能量通量分别对冠层阻抗系数C₁和地表阻抗系数C₂的敏感性响应.结果表明,感热通量对地表阻抗系数C₂更为敏感,而潜热通量则对冠层阻抗系数C₁较为敏感;感热通量随C₁增加而增大,随C₂增加而减小,而潜热通量则随C₁或C₂的增加而减小;不管是感热通量还是潜热通量,它们对阻抗系数的敏感度随阻抗系数的增大而减小,而对阻抗系数的相对敏感度则随阻抗系数的增大而增大.最后,结合那曲站夏季下垫面稀疏短草的分布特点分析了造成感热通量和潜热通量敏感变化各异的原因.     

Turbulent heat transfer from a sparsely vegetated surface: Two-component representation

The conventional calculation of heat fluxes from a vegetated surface, involving the coefficient of turbulent heat transfer, which increases logarithmically with surface roughness (commonly taken as about 0.12 of the plant height), appears inappropriate for highly structured surfaces such as desertscrub or open forest. An approach is developed here for computing sensible heat flux from sparsely vegetated surfaces, where the absorption of insolation and the transfer of absorbed heat to the atmosphere are calculated separately for the plants and for the soil. This approach is applied to a desert-scrub surface for which the turbulent transfer coefficient of sensible heat flux from the plants is much larger than that from the soil below, as shown by an analysis of plant, soil and air temperatures measured in an animal exclosure in the northern Sinai. The plant density is expressed as the sum of products (plant-height) x (plant-diameter) of plants per unit horizontal surface area (the dimensionless silhouette parameter of Lettau). The solar heat absorbed by the plants is assumed to be transferred immediately to the airflow. The effective turbulent transfer coefficientk _g-eff for sensible heat from the desert-scrub/soil surface computed under this assumption increases sharply with increasing solar zenith angle, as the plants absorb a greater fraction of the incoming irradiation. The surface absorptivity (the co-albedo) also increases sharply with increasing solar zenith angle, and thus the sensible heat flux from such complex surfaces (which include open forests) is a much broader function of time of day than when computed under constantk _g-eff and constant albedo assumptions. The major role that desert-fringe plants play in the genesis of convection and advection cannot be evaluated properly in the conventional calculations....you know not the way of the wind... Ecclesiastes, Ch. XI, verse 5     

Radiative surface temperature and convective flux calculation over crop canopies

The analysis presented in this paper aims at a better understanding of the potential role of radiative temperature, as measured by a radiometer over crops, in sensible heat flux calculation. Defining radiative temperature as the mean temperature of the surfaces viewed by the radiometer (leaves and soil surface) and assuming that an Ohm's law type formula can be used to express sensible heat flux as a function of the difference between air temperature and radiative temperature, the aerodynamic resistance which divides this temperature difference has been analytically defined. The parameters which appear in the resistance expression depend essentially on wind velocity and canopy structure but also on the inclination angle of the radiometer. Finally an experimental verification is presented with data obtained over a potato crop.     

Surface exchange of water vapour between an open sphagnum fen and the atmosphere

Water loss by evapotranspiration (ET) is a principal component of the hydrologic cycle in wetlands. Using micrometeorological techniques, we measured ET from a Sphagnum-dominated open fen in northcentral Minnesota (U.S.A.) from May to October in 1991 and 1992. The daily ET rate ranged from 0.2–4.8 mm d^-1 with a growing season average of 3.0 mm d^-1. The evapotranspiration rate of the fen was near the potential rate of open water evaporation when the vascular plants were actively growing and the water table level was within or above the rooting zone. Using a dual-source modification of the Penman-Monteith equation (Massman, 1992), we partitioned the measured ET into evaporation from the non-vascular Sphagnum surfaces and transpiration from vascular plants. The analysis indicated that about two thirds of the water vapour flux to the atmosphere was from evaporation when the Sphagnum surface was wet. Such an evaporative flux was expected because of vertical distribution of vascular plant leaves which had a small leaf area index (0.4–0.7) and intercepted only about 30% of net radiation (R _n ) during the day. The remainder of R _n was thus available for evaporation from Sphagnum. Evaporation significantly decreased as the Sphagnum surface dried out. When the water table was within the rooting zone (0–0.4 m), the vascular plants absorbed Sphagnum-generated sensible heat, which amounted up to one third of their transpiration energy flux. Under these conditions, the total water vapour flux remained near its potential rate owing to the enhanced transpiration from vascular plants. A drop in water table of 0.15–0.2 m below the hollow bottom during vascular plant senescence resulted in ET rates lower than the potential rates by 5–65%.     

A Vegetated Urban Canopy Model for Meteorological and Environmental Modelling

An urban canopy model is developed for use in mesoscale meteorological and environmental modelling. The urban geometry is composed of simple homogeneous buildings characterized by the canyon aspect ratio (h/w) as well as the canyon vegetation characterized by the leaf aspect ratio (σ _l ) and leaf area density profile. Five energy exchanging surfaces (roof, wall, road, leaf, soil) are considered in the model, and energy conservation relations are applied to each component. In addition, the temperature and specific humidity of canopy air are predicted without the assumption of thermal equilibrium. For radiative transfer within the canyon, multiple reflections for shortwave radiation and one reflection for longwave radiation are considered, while the shadowing and absorption of radiation due to the canyon vegetation are computed by using the transmissivity and the leaf area density profile function. The model is evaluated using field measurements in Vancouver, British Columbia and Marseille, France. Results show that the model quite well simulates the observations of surface temperatures, canopy air temperature and specific humidity, momentum flux, net radiation, and energy partitioning into turbulent fluxes and storage heat flux. Sensitivity tests show that the canyon vegetation has a large influence not only on surface temperatures but also on the partitioning of sensible and latent heat fluxes. In addition, the surface energy balance can be affected by soil moisture content and leaf area index as well as the fraction of vegetation. These results suggest that a proper parameterization of the canyon vegetation is prerequisite for urban modelling.     

Further Development of the Vegetated Urban Canopy Model Including a Grass-Covered Surface Parametrization and Photosynthesis Effects

The vegetated urban canopy model (VUCM), which includes parametrizations of urban physical processes for artificial surfaces and vegetated areas in an integrated system, has been further developed by including physical processes associated with grass-covered surfaces in urban pervious surfaces and the photosynthesis effects of urban vegetation. Using measurements made from three urban/suburban sites during the BUBBLE field campaign in 2002, the model’s performance in modelling surface fluxes (momentum flux, net radiation, sensible and latent heat fluxes and storage heat flux) and canopy air conditions (canopy air temperature and specific humidity) was critically evaluated for the non-precipitation and the precipitation days. The observed surface fluxes at the urban/suburban sites were significantly altered by precipitation as well as urban vegetation. Especially, the storage heat at urban surfaces and underlying substrates varied drastically depending on weather conditions while having an important role in the formation of a nocturnal urban surface layer. Unlike the nighttime canopy air temperature that was largely affected by the storage-heat release, the daytime canopy air conditions were highly influenced by the vertical turbulent exchange with the overlying atmosphere. The VUCM well reproduced these observed features in surface fluxes and canopy air conditions at all sites while performing well for both the non-precipitation and the precipitation days. The newly implemented parametrizations clearly improved the model’s performance in the simulation of sensible and latent heat fluxes at the sites, more noticeably at the suburban site where the vegetated area fraction is the largest among the sites. Sensitivity analyses for model input parameters in VUCM elucidated the relative importance of the morphological, aerodynamic, hydrological and radiative/thermal properties in modelling urban surface fluxes and canopy air conditions for daytime and nighttime periods. These results suggest that the VUCM has great potential for urban atmospheric numerical modelling for a range of cities and weather conditions in addition to having a better physical basis in the representation of urban vegetated areas and associated physical processes.     

Estimating sensible heat flux from radiometric temperature over crop canopy

The model devised by Lhommeet al. (1988) allows one to calculate the sensible heat flux over a homogeneous crop canopy from radiometric surface temperature by adding a so-called canopy aerodynamic resistance to the classical aerodynamic resistance calculated above the canopy. This model is reformulated in order to simplify the mathematical procedure needed to calculate this additional resistance. Analytical expressions of micrometeorological profiles within the canopy are introduced. Assuming a constant leaf area density, an analytical expression of canopy aerodynamic resistance is inferred, which is a function of wind velocity, inclination angle of the radiometer and crop characteristics such as crop height, leaf area index, inclination index of the foliage and leaf width. Sensitivity of this resistance to the different parameters is investigated. The most significant are wind velocity and LAI. Finally, the predictions of the model are tested against two sets of measurements obtained for two different crops, potato and maize.     

On the subarctic North Atlantic cooling due to global warming

Various ocean reanalysis data reveal that the subarctic Atlantic sea surface temperature (SST) has been cooling during the twentieth century. A similar cooling pattern is found in the doubling CO₂ experiment obtained from the CMIP3 (coupled model intercomparison project third phase) compared to the pre-industrial experiment. Here, in order to investigate the main driver of this cooling, we perform the heat budget analysis on the subarctic Atlantic upper ocean temperature. The net surface heat flux associated with the increased concentration of greenhouse gases heats the subarctic ocean surface. In the most of models, the longwave radiation, latent heat flux, and sensible heat flux exert a warming effect, and the shortwave radiation exerts a cooling effect. On the other hand, the thermal advection by the meridional current reduces the subarctic upper ocean temperature in all models. This cold advection is attributed to the weakening of the meridional overturning circulation, which is related to the reduction in the ocean surface density. In particular, greater warming of the surface air than of the sea surface results in the reduction of surface evaporation and thereby enhanced freshening of the ocean surface water, while precipitation change was smaller than evaporation change. The thermal advections by both the wind-driven Ekman current and the density-driven geostrophic current contribute to cooling in most of the models, where the heat transport by the geostrophic current tends to be larger than that by the Ekman current.     

Regional-Scale Heat and Water Vapour Fluxes in an Agricultural Landscape: An Evaluation of CBL Budget Methods at OASIS

This paper evaluates convective boundary layer (CBL) budget methods as a tool for estimating regionally averaged sensible and latent heat fluxes for the study region used in OASIS (Observations at Several Interacting Scales). This is an agricultural region of mixed cropping and grazing extending about 100 km west of the town of Wagga Wagga, NSW, Australia.The analysis proceeds in three stages: first, a simpleone-dimensional model of the well-mixed layer (the CBL slab model), forced with measurements of the surface heat and evaporation fluxes, is evaluated by comparing measured and modelled CBL temperature, humidity and depths. A comparison of several entrainment schemes shows that a simple model, where the entrainment kinetic energy is parameterised as a fraction (₃) of the surface sensible heat flux, works well if is set to 0.5. Second, the slab model is coupled to a Penman–Monteith model of surface evaporation to predict regional scale evaporation and thence heat fluxes. Finally, the integral CBL budget approach, which is an inverse method using theone-dimensional slab model, is used to infer regional heat and evaporation fluxes from measured time series of CBL temperature and humidity.We find that the simple CBL slab model works reasonably well for predicting CBL depth and very well for CBL temperature, especially if approximate estimates of subsidence velocity and warming due to advection are included. Regional sensible heat fluxes estimated from the integral CBL method match those measured, although the method is very sensitive to measurement errors. Measurement-model differences were larger for short integration times, because the well-mixed assumptions are violated at particular times of the day. The corollary is that `whole-day' (0530–1530 h) estimates are in reasonable agreement with measured values. Integral methods could not be used to infer the regional evaporation flux directly because CBL humidity profiles were complex and often not well mixed until mid-afternoon. We recommend that regional evaporation fluxes be predicted either from a coupled Penman–Monteith – CBL slab model, or inferred as a residual term from estimates of the regionally averaged available energy and sensible heat flux. Furthermore, we show that inferring fluxes via integral methods will always be difficult when the scalar concentrations have either a large surface source and free atmosphere sink (in the case of water vapour and methane), or a large surface sink and upper level source (in the case of CO₂).     

Thermal effect of the sea breeze on the structure of the boundary layer and the heat budget over land

The daytime boundary-layer heating process and the air-land heat budget were investigated over the coastal sea-breeze region by means of observations over the Sendai plain in Japan during the summer. In this area, the onset of the sea breeze begins at the coast around 0900 LST, intruding about 35 km inland by late afternoon. The cold sea breeze creates a temperature difference of over 10°C between the coastal and inland areas in the afternoon. On the other hand, warm air advection due to the combination of the counter-sea breeze and land-to-sea synoptic wind occurs in the layer above the cold sea breeze in the coastal region. Owing to this local warm air advection, there is no significant difference in the daytime heating rate over the entire atmospheric boundary layer between the coastal and inland areas. The sensible heat flux from the land surface gradually decreases as distance from the coastline increases, being mainly attributed to the cold sea breeze. The daytime mean cold air advection due to the sea breeze is estimated asQ _adv ^local =–29 W m^–2 averaged over the sea breeze region (035 km from the coastline). This value is 17% of the surface sensible heat fluxH over the same region. The results of a two-dimensional numerical model show that the value ofQ _adv ^local /H is strongly affected by the upper-level synoptic wind direction. The absolute value ofQ _adv ^local /H becomes smaller when the synoptic wind has the opposite direction of the sea breeze. This condition occurred during the observations used in the present study.     

The Surface Energy Budget and Its Impact on the Freeze-thaw Processes of Active Layer in Permafrost Regions of the Qinghai-Tibetan Plateau

The surface energy budget is closely related to freeze-thaw processes and is also a key issue for land surface process research in permafrost regions.In this study,in situ data collected from 2005 to 2015 at the Tanggula site were used to analyze surface energy regimes,the interaction between surface energy budget and freeze-thaw processes.The results confirmed that surface energy flux in the permafrost region of the Qinghai-Tibetan Plateau exhibited obvious seasonal variations.Annual average net radiation(R_n)for 2010 was 86.5 W m^-2,with the largest being in July and smallest in November.Surface soil heat flux(G₀)was positive during warm seasons but negative in cold seasons with annual average value of 2.7 W m^-2.Variations in R_n and G₀ were closely related to freeze-thaw processes.Sensible heat flux(H)was the main energy budget component during cold seasons,whereas latent heat flux(LE)dominated surface energy distribution in warm seasons.Freeze-thaw processes,snow cover,precipitation,and surface conditions were important influence factors for surface energy flux.Albedo was strongly dependent on soil moisture content and ground surface state,increasing significantly when land surface was covered with deep snow,and exhibited negative correlation with surface soil moisture content.Energy variation was significantly related to active layer thaw depth.Soil heat balance coefficient K was>1 during the investigation time period,indicating the permafrost in the Tanggula area tended to degrade.