半干旱草原下垫面能量平衡特征及土壤热通量对能量闭合率的影响

土壤热通量在半干旱草原下垫面能量平衡研究中极为重要,土壤热通量估计不够准确是导致地表能量不平衡的一个重要原因。利用2008年6—9月锡林郭勒草原主生长期地表辐射、通量和土壤温度梯度观测资料,研究中纬度半干旱草原下垫面地表能量平衡特征。首先,在分析能量平衡各分量月平均日变化特征的基础上,通过对土壤热流量板观测的5 cm深度土壤热通量(G)的相位前移,研究了土壤热通量相位滞后对地表能量平衡产生的影响;其次,利用谐波分析方法,通过计算地表土壤热通量(Gs),分析了地表到热流量板之间的土壤热量储存对地表能量平衡的影响。结果表明:(1)将土壤热通量相位前移30 min,湍流通量与可利用能量(Rn-G)线性回归的斜率从0.835增加到0.842,地表能量闭合率提高了0.7%,但仍有15.8%的能量不闭合;(2)考虑了地表到热流量板之间的土壤热量储存之后,湍流通量与可利用能量之间的回归斜率达到0.979,能量不闭合程度仅为2.1%。     

Effective parameters of surface energy balance in heterogeneous landscape

This paper addresses the problem of estimating surface fluxes at large scale over heterogeneous terrain, and the corresponding determination of effective surface parameters. Two kinds of formulation are used to calculate the fluxes of sensible and latent heat: the basic diffusion equations (Ohm's law type) and the Penman-Monteith equations. The strategy explored is based upon the principle of flux conservation, which stipulates that the average flux over a large area is simply the area-weighted mean of the contributions from the different patches making up the area. We show that the application of this strategy leads to different averaging schemes for the surface parameters, depending on the type of flux (latent heat, sensible heat) and on the type of formulation used to express the flux. It appears that the effective value of a given parameter must be appraised for each individual application, because it is not unique, but differs according to the magnitude being conserved and the equation used to express this magnitude. Numerical simulations are carried out to test over contrasted areas the aggregation procedures obtained. The areal fluxes estimated from these effective parameters, together with the areal fluxes calculated by means of a simple areal averaging of the parameters, are compared to the true average fluxes, calculated as area-weighted means of the elementary fluxes. The aggregation procedures obtained prove to be much more accurate for estimating areal fluxes and for closing the energy balance equation than those based upon simple areal averaging of the parameters.     

Effects of the soil heat flux estimates on surface energy balance closure over a semi-arid grassland

Soil heat flux is important for surface energy balance (SEB), and inaccurate estimation of soil heat flux often leads to surface energy imbalance. In this paper, by using observations of surface radiation fluxes and soil temperature gradients at a semi-arid grassland in Xilingguole, Inner Mongolia, China from June to September 2008, the characters of the SEB for the semi-arid grassland were analyzed. Firstly, monthly averaged diurnal variations of SEB components were revealed. A 30-min forward phase displacement of soil heat flux (G) observed by a fluxplate at the depth of 5-cm below the soil surface was conducted and its effect on the SEB was studied. Secondly, the surface soil heat flux (G _s) was computed by using harmonic analysis and the effect of the soil heat storage between the surface and the fluxplate on the SEB was examined. The results show that with the 30-min forward phase displacement of observed G, the slope of the ordinary linear regression (OLR) of turbulent fluxes (H+LE) against available energy (R _n-G) increased from 0.835 to 0.842, i.e., the closure ratio of SEB increased by 0.7%, yet energy imclosure of 15.8% still existed in the SEB. When G _s, instead of G was used in the SEB equation, the slope of corresponding OLR of (H+LE) against (R _n-G _s) reached 0.979, thereby the imclosure ratio of SEB was reduced to only 2.1%.     

Influence of heterogeneous land surfaces on surface energy and mass fluxes

Summary Land-surface heterogeneity affects surface energy fluxes. The magnitudes of selected land-surface influences are quantified by comparing observations with model simulations of the FIFE (First ISLSCP Field Experiment) domain. Several plausible heterogeneous and homogeneous initial and boundary conditions are examined, although soilmoisture variability is emphasized. It turns out that simple spatial averages of surface variation produced biased flux values. Simulated maximum latent-heat fluxes were approximately 30 to 40 W m^–2 higher, and air temperatures 0.4 °C lower (at noon), when computations were initialized with spatially averaged soil-moisture and leaf-area-index fields. The planetary boundary layer (PBL) height and turbulent exchanges were lower as well. It additionally was observed that (largely due to the nonlinear relationship between initial soil-moisture availability and the evapotranspiration rate), real latent-heat flux can be substantially less than simulated latent-heat flux using models initialized with spatially averaged soil-moisture fields. Differences between real and simulated fluxes also vary with the resolution at which real soil-moisture heterogeneity is discretized.With 8 Figures     

Observations of fluxes over heterogeneous surfaces

This study analyzes data collected from repeated aircraft runs 30 m over alternating regions of irrigated and dry nonirrigated surfaces, each region on the order of 10 km across, during the California Ozone Deposition Experiment (CODE). After studying the scale dependence of the flow, the variables and their fluxes are decomposed into means for sublegs defined in terms of irrigated and nonirrigated regions and deviations from such subleg means. Since the repeated runs were flown over the same track, compositing the eight flight legs for each of the two days allows partial isolation of the influences of surface heterogeneity and transient mesoscale motions.A variance analysis is carried out to quantify the relative importance of surface heterogeneity and transient mesoscale motions on the variability of the turbulence fluxes. The momentum and ozone fluxes are more influenced by transient mesoscale motions while fluxes of heat, moisture and carbon dioxide are more influenced by surface heterogeneity. The momentum field is also influenced by a quasi-stationary mesoscale front and larger scale velocity gradients.For the present case, the mesoscale modulation of the turbulent flux is numerically more important than the direct mesoscale flux. This spatial modulation of the turbulent fluxes leads to extra Reynolds terms which act to reduce the area-averaged turbulent momentum flux and enhance the area-averaged turbulent heat flux.     

The bulk aerodynamic formulation over heterogeneous surfaces

This interpretative literature survey examines problems with application of the bulk aerodynamic method to spatially averaged fluxes over heterogeneous surfaces. This task is approached by tying together concepts from a diverse range of recent studies on subgrid parameterization, the roughness sublayer, the roll of large inactive boundary-layer eddies, internal boundary-layer growth, the equilibrium sublayer, footprint theory and the blending height. Although these concepts are not completely compatible, qualitative scaling arguments based on these concepts lead to a tentative unified picture of the qualitative influence of surface heterogeneity for a wide spectrum of spatial scales.Generalization of the velocity scale is considered to account for nonvanishing heat and moisture fluxes in the limit of vanishing time-averaged wind speed and to account for the influence of subgrid mesoscale motions on the grid-averaged turbulent flux. The bulk aerodynamic relationship for the heat flux usually employs the surface radiation temperature or, equivalently, the temperature from the modelled surface energy budget. The corresponding thermal roughness length is quite variable and its dependence on available parameters is predictable only in special cases.An effective transfer coefficient to relate the spatially averaged surface fluxes to spatially averaged air-ground differences of temperature and other scalars can be most clearly defined when the blending height occurs below the reference level (observational level or first model level). This condition is satisfied only for surface heterogeneity occurring over horizontal scales up to a few times the boundary-layer depth, depending on the stability and height of the reference level. For surface heterogeneity on larger scales (small mesoscale), an effective transfer coefficient for the spatially averaged flow must be defined, for which predictive schemes are unavailable. For surface variations on large mesoscales, homogeneous subareas may be maintained where traditional similarity theory is locally applicable. Surface variations on these scales may generate thermally-driven mesoscale motions.     

Acoustic tomography and conventional meteorological measurements over heterogeneous surfaces

Summary Conventional micro-meteorological measuring methods are not particularly suitable for the investigation of the energy exchange under heterogeneous surface conditions. To consider the influence of the different surface properties, area-covered and spatially averaged meteorological measurements in combination with highly resolved simulations are necessary. In this context, the method of acoustic travel time tomography is introduced to provide information about the horizontal temperature and wind field.To check the applicability of the tomographic method for this problem and to provide data for the model initialisation a field experiment was carried out.The tomographic system and conventional meteorological equipment were utilised inside the investigation area with an extension of several hundred meters (300×700m²), which was arranged over areas with different surface properties: grassland and bare soil.The results of the field experiment show that differences between the measuring systems (in-situ and remote sensing) exist and the heterogeneity of the underlying surface is visible in the near surface temperature and wind field at a scale, which can be resolved with highly resolved numerical models. Depending on the incoming solar radiation and the local advection regional distinctions in the air temperature and wind field as well as in the vertical sensible heat fluxes were observed.The investigations demonstrate that the sensitivity of the acoustic tomography is sufficient to verify gradients in the meteorological fields even when the horizontal differences in sensible heat flux are small.     

北京郊区草地夏季能量收支平衡的数值模拟

北京郊区地表能量分配可能影响北京地区的天气和气候。为了进一步检验陆面过程模式对北京郊区具有代表性的稀疏草地地表能量分配的模拟能力,利用原版和改进版简化生物圈模式(SiB2,Simple Biosphere Model 2)模拟了2010年7月22日-8月5日期间北京郊区阳坊镇坦克打靶场草地的辐射平衡、能量收支以及地表热通量。并将模拟结果与实际测量的数据进行对比,结果表明:1)原版SiB2低估净辐射11.32%,改进版SiB2则低估净辐射5.81%,主要原因是改进版SiB2更新了土壤热传导率计算方法,从而提高了土壤温度(包括地表温度)模拟结果的精度,进而改善了地表向上的长波辐射模拟结果的准确性;2)改进版SiB2同时改善了感热通量和潜热通量的模拟结果,但是原版SiB2和改进版SiB2均低估了土壤热通量。     

Available energy and energy balance closure at four coniferous forest sites across Europe

The available energy (AE), driving the turbulent fluxes of sensible heat and latent heat at the earth surface, was estimated at four partly complex coniferous forest sites across Europe (Tharandt, Germany; Ritten/Renon, Italy; Wetzstein, Germany; Norunda, Sweden). Existing data of net radiation were used as well as storage change rates calculated from temperature and humidity measurements to finally calculate the AE of all forest sites with uncertainty bounds. Data of the advection experiments MORE II (Tharandt) and ADVEX (Renon, Wetzstein, Norunda) served as the main basis. On-site data for referencing and cross-checking of the available energy were limited. Applied cross checks for net radiation (modelling, referencing to nearby stations and ratio of net radiation to global radiation) did not reveal relevant uncertainties. Heat storage of sensible heat J _H, latent heat J _E, heat storage of biomass J _veg and heat storage due to photosynthesis J _C were of minor importance during day but of some importance during night, where J _veg turned out to be the most important one. Comparisons of calculated storage terms (J _E, J _H) at different towers of one site showed good agreement indicating that storage change calculated at a single point is representative for the whole canopy at sites with moderate heterogeneity. The uncertainty in AE was assessed on the basis of literature values and the results of the applied cross checks for net radiation. The absolute mean uncertainty of AE was estimated to be between 41 and 52 W m^?2 (10–11 W m^?2 for the sum of the storage terms J and soil heat flux G) during mid-day (approximately 12% of AE). At night, the absolute mean uncertainty of AE varied from 20 to about 30 W m^?2 (approximately 6 W m^?2 for J plus G) resulting in large relative uncertainties as AE itself is small. An inspection of the energy balance showed an improvement of closure when storage terms were included and that the imbalance cannot be attributed to the uncertainties in AE alone.     

The global energy balance from a surface perspective

In the framework of the global energy balance, the radiative energy exchanges between Sun, Earth and space are now accurately quantified from new satellite missions. Much less is known about the magnitude of the energy flows within the climate system and at the Earth surface, which cannot be directly measured by satellites. In addition to satellite observations, here we make extensive use of the growing number of surface observations to constrain the global energy balance not only from space, but also from the surface. We combine these observations with the latest modeling efforts performed for the 5th IPCC assessment report to infer best estimates for the global mean surface radiative components. Our analyses favor global mean downward surface solar and thermal radiation values near 185 and 342 Wm^?2, respectively, which are most compatible with surface observations. Combined with an estimated surface absorbed solar radiation and thermal emission of 161 and 397 Wm^?2, respectively, this leaves 106 Wm^?2 of surface net radiation available globally for distribution amongst the non-radiative surface energy balance components. The climate models overestimate the downward solar and underestimate the downward thermal radiation, thereby simulating nevertheless an adequate global mean surface net radiation by error compensation. This also suggests that, globally, the simulated surface sensible and latent heat fluxes, around 20 and 85 Wm^?2 on average, state realistic values. The findings of this study are compiled into a new global energy balance diagram, which may be able to reconcile currently disputed inconsistencies between energy and water cycle estimates.     

Stratified flow over non-uniform surfaces: Turbulent energy model

A turbulent energy model is developed to simulate the response of a neutrally stratified atmospheric boundary layer to sudden changes in surface roughness. A mechanism of turbulent energy transfer is proposed, based upon the results of numerical experiments, that explains the distribution of shear stress and hence the distribution of velocity profiles in the atmospheric surface layer. Two length scales associated with the turbulent energy equation are obtained from experimental data and the law of the wall. Turbulent energy is also predicted.The predicted growth of the internal boundary layer is slower than that obtained from mixing-length models. Also, the predicted surface shear stress obtained from the turbulent energy model is in better agreement with field data than that obtained from mixing-length models.     

Surface energy balance complexity in GCM land surface models

 A large number of land surface models (LSMs) have been designed for use in atmospheric general circulation models (GCMs) and GCM modellers therefore have a large number of options when selecting an LSM for their GCM. This study provides information to aid LSM design choices. A framework within which sensitivity to LSM design can be tested is presented and a series of experiments carried out to investigate how general aspects of surface energy balance parametrisation affect land-atmosphere evaporation. Firstly, it is shown that a combination of surface energy balance complexity and aerodynamic parametrisation can be used to explain the gross simulation differences obtained in the Project for Intercomparison of Land-surface Parameterization Schemes (PILPS). Secondly, a simple surface energy balance parametrisation with a constant surface resistance is found to be as appropriate as more complex method for simulating annual, monthly and seasonally averaged diurnal cycles of evaporation. However, complex aspects of surface energy balance parametrisation (canopy interception, bare ground evaporation and canopy resistance) are shown to contain substantial geographic and daily functionality that is not present in the simpler parametrisation. Received: 12 June 1998/Accepted: 24 December 1998     

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.     

The energy balance of water drops on a leaf surface

The amount of heat conducted to an isolated drop of water on a leaf is described by an equation which includes an effective thermal conductivity coefficient. Measurements of the dimensions of water drops on a horizontal wheat leaf give relationships which allow the volume and exposed surface areas to be obtained from drop diameter. These relationships are used in the experimental determination of the drop boundary-layer resistance and the effective thermal conductivity coefficient for drops on a leaf in a chamber. The boundary-layer resistance of the drop appeared to be independent of drop size and the mean value was about 60% of that for one side of the leaf. For drops with diameters less than 1 mm, conduction of heat to the drop reduced the leaf-to-drop temperature difference by over 50% of the value expected without conduction. Conduction of heat to drops will significantly affect the evaporation rate of surface water from cereal canopies after rain.     

Components of surface energy balance in a temperate grassland ecosystem

Eddy correlation measurements were made of fluxes of moisture, heat and momentum at a tallgrass prairie site near Manhattan, Kansas, U.S.A. during the First ISLSCP Field Experiment (FIFE) in 1987. The study site is dominated by three C₄ grass species: big bluestem (Andropogon gerardii), indiangrass (Sorghastrum nutans), and switchgrass (Panicum virgatum). The stomatal conductance and leaf water potential of these grass species were also measured.In this paper, daily and seasonal variations in the components of the surface energy balance are examined. The aerodynamic and canopy surface conductances for the prairie vegetation are also evaluated.Published as Paper No. 8987, Journal Series, Nebraska Agricultural Research Division.

---

     

Diagnosing the diurnal surface energy balance over the summer tundra at Princess Marie Bay from simple short-period measurements

Summary Using existing physical parameterizations, a new mathematical model is formulated to diagnose the diurnal variation of the energy fluxes and temperature on the snow-free tundra surface at Princess Marie Bay, Ellesmere Island, Canada. The input to the model consists of three meteorological variables which can be readily measured by an automatic weather station: incoming short-wave radiation, windspeed and screen level temperature. The model is based on the one-dimensional heat conduction equation for unfrozen soil, with surface heat exchange by short- and long-wave radiation and by convection and evaporation. A permafrost surface is used as a lower boundary condition. The model is formulated and tuned using a series of data from the Princess Marie Bay site. It is then tested using a separate data set from the same site and an independent data set from a nearby site.With 10 Figures     

An objective method to achieve closure of overdetermined surface energy budgets

The analysis of an overdetermined energy budget provides an opportunity to investigate the existence of advective effects at the site and in the atmospheric layer represented by the energy budget. The present contribution presents a set of hourly averaged energy budget data representing the atmospheric surface layer over a suburban surface. The measured variables result in an overdetermined budget. A statistical analysis shows that the residual energy has a near-zero mean with no apparent relation to wind-speed or wind-direction, thus eliminating the possibilty of advection. It is then shown how the residual may be distributed between the two turbulent fluxes on the presumption that it represents measurement errors which are the largest in those fluxes. This method is expected to provide more reliable flux estimates in environments characterized by large turbulent fluxes subject to significant errors.     

A new device for continuous recording of the energy balance of natural surfaces

A new device for continuous computing and recording of the energy balance terms is described with some of its principal advantages and characteristics. Hourly and daily printed mean results are compared with those of a sensible weighing lysimeter. Satisfactory performances were obtained during the two last past years, so it is actually used on any natural surfaces, instead of weighing lysimeter, with its simple package which only needs power.     

Effects of vegetative heterogeneity and patch-scale harvest on energy balance closure and flux measurements

Based on the micrometeorological measurements at a heterogeneous farmland in South China, this work detects the effects of vegetative heterogeneity and patch-scale harvest on the energy balance closure and turbulent fluxes. As a quality control, the integral turbulent characteristics are analyzed in the framework of Monin-Obukhov similarity theory. Spatial representativeness of the measurements is studied in terms of footprint and source area. Firstly, in two wind sectors, the nondimensional standard deviations of turbulent quantities generally agree with some foregoing studies. Discrepancies exist in the other sectors due to the instrument-induced flow distortion. Secondly, energy balance closure is examined with two types of linear regression, which confirms that mismatching source areas between the available energy and turbulent fluxes have no preference to either energy “deficit” or “surplus”. Thirdly, turbulent fluxes exhibit greater variability when they represent smaller source areas. The patch-scale harvest adjacent to the flux mast causes notable increase and decrease in the sensible heat and latent heat fluxes, respectively, while the CO₂ exchange almost vanishes after the harvest. Interestingly, energy balance closure is less influenced despite the notable effects on the turbulent fluxes and Bowen ratio, implying that the energy balance closure check may mask some variability in the turbulent fluxes. Thus, to adjust the heat fluxes with a single “closure factor” for a perfect closure is dangerous at a patchy site.