Contrasting Characteristics of the Surface Energy Balance between the Urban and Rural Areas of Beijing

A direct comparison of urban and rural surface energy balances, as well as a variety of other variables including incoming shortwave/longwave radiation and aerosol optical depth, is conducted for the Beijing metropolitan area. The results indicate that, overall, the urban area receives a smaller amount of incoming shortwave radiation but a larger amount of incoming longwave radiation. However, comparisons in the aerosol optical depth and cloud fraction at the two locations suggest that neither aerosol optical depth nor cloud fraction alone can explain the difference in the incoming shortwave radiation. The urban–rural differences in the incoming longwave radiation are unlikely to be caused by the presence of more abundant greenhouse gases over the urban area, as suggested by some previous studies, given that water vapor is the most dominant greenhouse gas and precipitable water is found to be less in urban areas. The higher incoming longwave radiation observed over the urban area is mostly likely due to the higher temperatures of the ambient air. The urban area is also found to always produce higher sensible heat fluxes and lower latent heat fluxes in the growing season. Furthermore, the urban area is associated with a larger amount of available energy（the sum of sensible and latent heat fluxes） than the rural area, except in May and October when evapotranspiration in the rural area significantly exceeds that in the urban area. This study provides observational evidence of urban–rural contrasts in relevant energy-balance components that plausibly arise from urban–rural differences in atmospheric and land-surface conditions.

An Assessment of Storage Terms in the Surface Energy Balance of a Subalpine Meadow in Northwest China

The heat storage terms in the soil-vegetation-atmosphere system may play an important role in the surface energy budget.In this paper,we evaluate the heat storage terms of a subalpine meadow based on a ficld experiment conducted in the complex terrain of the eastern Qilian Mountains of Northwest China and their impact on the closure of the surface energy balance under such non-ideal conditions.During the night, the average sum of the storage terms was -5.5 W m,which corresponded to 10.4%of net radiation.The sum of the terms became positive at 0730 LST and negative again at about 1500 LST,with a maximum value of 19 W mobserved at approximately 0830 LST.During the day,the average of the sum of the storage terms was 6.5 W m,which corresponded to 4.0%of net radiation.According to the slopes obtained when linear regression of the net radiation and partitioned fluxes was forced through the origin,there is an imbalance of 14.0%in the subalpine meadow when the storage terms are not considered in the surface energy balance.This imbalance was improved by 3.4%by calculating the sum of the storage terms.The soil heat storage flux gave the highest contribution(1.59%),while the vegetation enthalpy change and the rest of the storage terms were responsible for improvements of 1.04%and 0.77%,respectively.

Evapotranspiration Estimation Based on MODIS Products and Surface Energy Balance Algorithms for Land(SEBAL) Model in Sanjiang Plain,Northeast China

In this study,the Surface Energy Balance Algorithms for Land(SEBAL) model and Moderate Resolution Imaging Spectroradiometer(MODIS) products from Terra satellite were combined with meteorological data to estimate evapotranspiration(ET) over the Sanjiang Plain,Northeast China.Land cover/land use was classified by using a recursive partitioning and regression tree with MODIS Normalized Difference Vegetation Index(NDVI) time series data,which were reconstructed based on the Savitzky-Golay filtering approach.The MODIS product Quality Assessment Science Data Sets(QA-SDS) was analyzed and all scenes with valid data covering more than 75% of the Sanjiang Plain were selected for the SEBAL modeling.This provided 12 overpasses during 184-day growing season from May 1st to October 31st,2006.Daily ET estimated by the SEBAL model was misestimaed at the range of-11.29% to 27.57% compared with that measured by Eddy Covariance system(10.52% on average).The validation results show that seasonal ET from the SEBAL model is comparable to that from ground observation within 8.86% of deviation.Our results reveal that the time series daily ET of different land cover/use increases from vegetation on-going until June or July and then decreases as vegetation senesced.Seasonal ET is lower in dry farmland(average(Ave):491 mm) and paddy field(Ave:522 mm) and increases in wetlands to more than 586 mm.As expected,higher seasonal ET values are observed for the Xingkai Lake in the southeastern part of the Sanjiang Plain(Ave:823 mm),broadleaf forest(Ave:666 mm) and mixed wood(Ave:622 mm) in the southern/western Sanjiang Plain.The ET estimation with SEBAL using MODIS products can provide decision support for operational water management issues.

Surface energy balance of Keqicar Glacier, Tianshan Mountains, China, during ablation period

The meteorological data of ablation season in 2005 were recorded by two automatic weather stations on Keqicar Glacier, in the southwest Tianshan Mountains of China. One is operated on the glacier near the equilibrium line with an altitude of 4,265 m （Site A） and another is operated on the glacier ablation area with an altitude of 3,700 m （Site B）. These data were used to analyze the meteorological conditions and the surface energy balance （SEB） of Keqicar Glacier. Net radiation was directly measured, and turbulent heat fluxes were calculated using the bulk aerodynamic approach, including stability correction. The ablation value of 0.68 m w.e. derived from four ablation stakes is in close correspondence to the modeled value of 0.71 m w.e. During the observation period, net radiation accounts for 81.4% of the total energy with its value of 63.3 W/m2. The rest energy source is provided by the sensible heat flux with a value of 14.4 W/m2. Energy is consumed mainly by melting and evaporation, accotmting for 69.5% and 29.7% of the total energy with their values of 54.0 and 23.0 W/m2, respectively. Radiative energy dominates energy exchanges at the glacier-atmosphere interface, governed by the variation in net shortwave radiation. Net short-wave radiation varies significantly due to the effects of cloudiness and the high albedo caused by solid precipitation. Wind speed influences the turbulent heat fluxes distinctively and sensible heat flux and latent heat flux are much larger in July with high wind speed.

Estimations of Land Surface Characteristic Parameters and Turbulent Heat Fluxes over the Tibetan Plateau Based on FY-4A/AGRI Data

Accurate estimates of land surface characteristic parameters and turbulent heat fluxes play an important role in the understanding of land-atmosphere interaction. In this study, Fengyun-4A (FY-4A) Advanced Geostationary Radiation Imager (AGRI) satellite data and the China Land Data Assimilation System (CLDAS) meteorological forcing dataset CLDAS-V2.0 were applied for the retrieval of broadband albedo, land surface temperature (LST), radiation flux components, and turbulent heat fluxes over the Tibetan Plateau (TP). The FY-4A/AGRI and CLDAS-V2.0 data from 12 March 2018 to 30 April 2018 were first used to estimate the hourly turbulent heat fluxes over the TP. The time series data of in-situ measurements from the Tibetan Observation and Research Platform were divided into two halves-one for developing retrieval algorithms for broadband albedo and LST based on FY-4A, and the other for the cross validation. Results show the root-mean-square errors (RMSEs) of the FY-4A retrieved broadband albedo and LST were 0.0309 and 3.85 K, respectively, which verifies the applicability of the retrieval method. The RMSEs of the downwelling/upwelling shortwave radiation flux and downwelling/upwelling longwave radiation flux were 138.87/32.78 W m^(−2) and 51.55/17.92 W m^(−2), respectively, and the RMSEs of net radiation flux, sensible heat flux, and latent heat flux were 58.88 W m^(−2), 82.56 W m^(−2) and 72.46 W m^(−2), respectively. The spatial distributions and diurnal variations of LST and turbulent heat fluxes were further analyzed in detail.

Characteristics of land-atmosphere energy and turbulentfluxes over the plateau steppe in central Tibetan Plateau

The land-atmosphere energy and turbulence exchange is key to understanding land surface processes on the Tibetan Plateau(TP). Using observed data for Aug. 4 to Dec. 3, 2012 from the Bujiao observation point(BJ) of the Nagqu Plateau Climate and Environment Station(NPCE-BJ), different characteristics of the energy flux during the Asian summer monsoon(ASM) season and post-monsoon period were analyzed. This study outlines the impact of the ASM on energy fluxes in the central TP. It also demonstrates that the surface energy closure rate during the ASM season is higher than that of the post-monsoon period. Footprint modeling shows the distribution of data quality assessments(QA) and quality controls(QC) surrounding the observation point. The measured turbulent flux data at the NPCE-BJ site were highly representative of the target land-use type. The target surface contributed more to the fluxes under unstable conditions than under stable conditions. The main wind directions(180° and 210°) with the highest data density showed flux contributions reaching 100%, even under stable conditions. The lowest flux contributions were found in sectors with low data density, e.g., 90.4% in the 360° sector under stable conditions during the ASM season. Lastly, a surface energy water balance(SEWAB) model was used to gap-fill any absent or corrected turbulence data. The potential simulation error was also explored in this study. The Nash-Sutcliffe model efficiency coefficients(NSEs) of the observed fluxes with the SEWAB model runs were 0.78 for sensible heat flux and 0.63 for latent heat flux during the ASM season, but unrealistic values of-0.9 for latent heat flux during the post-monsoon period.

The Representation of Soil Moisture−Atmosphere Feedbacks across the Tibetan Plateau in CMIP6

Thermal processes on the Tibetan Plateau(TP)influence atmospheric conditions on regional and global scales.Given this,previous work has shown that soil moisture−driven surface flux variations feed back onto the atmosphere.Whilst soil moisture is a source of atmospheric predictability,no study has evaluated soil moisture−atmosphere coupling on the TP in general circulation models(GCMs).In this study,we use several analysis techniques to assess soil moisture−atmosphere coupling in CMIP6 simulations including:instantaneous coupling indices;analysis of flux and atmospheric behaviour during dry spells;and a quantification of the preference for convection over drier soils.Through these metrics we partition feedbacks into their atmospheric and terrestrial components.Consistent with previous global studies,we conclude substantial inter-model differences in the representation of soil moisture−atmosphere coupling,and that most models underestimate such feedbacks.Focusing on dry spell analysis,most models underestimate increased sensible heat during periods of rainfall deficiency.For example,the model-mean bias in anomalous sensible heat flux is 10 W m−2(≈25%)smaller compared to observations.Deficient dry-spell sensible heat fluxes lead to a weaker atmospheric response.We also find that most GCMs fail to capture the negative feedback between soil moisture and deep convection.The poor simulation of feedbacks in CMIP6 experiments suggests that forecast models also struggle to exploit soil moisture−driven predictability.To improve the representation of land−atmosphere feedbacks requires developments in not only atmospheric modelling,but also surface processes,as we find weak relationships between rainfall biases and coupling indexes.

Analysis of observations on the urban surface energy balance in Beijing

The 1-year（2009-2010） measurements are analyzed of the urban surface energy balance（SEB） obtained from the sensors located at three vertical layers of a 325-m tower in downtown Beijing.Results show that：（1） The measurements from the 325-m tower represent the SEB characteristics of the cities located in semi-humid warm-temperate continental monsoon climate zone.In a typical hot and rainy summer,cold and dry winter,the measured Bowen ratio is minimum in summer and maximum in winter.The Bowen ratio measured at 140 m for spring,summer,autumn,and winter are 2.86,0.82,1.17,and 4.16 respectively.（2） At the height of 140-m（in the constant flux layer）,the noontime albedo is ~0.10 for summer,~0.12 for spring and autumn,and ~0.14 for winter.The ratios of daytime sensible heat flux,latent heat flux,and storage heat flux to net radiation are 0.25,0.16,and 0.59 for clear-sky days,and 0.33,0.19,and 0.48 for cloudy days respectively.（3） Under clear-sky days,the nighttime sensible heat flux is almost zero,but the latent heat flux is greater than zero.For cloudy days,the nighttime sensible heat flux is slightly greater than the latent heat flux in winter.The nighttime upward heat flux is presumably due to the anthropogenic release（mainly latent heat for summer,while latent and sensible heat for winter）.

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-em below the soil surface was conducted and its effect on the SEB was studied. Secondly, the surface soil heat flux （Gs） 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 （Rn 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 Gs, instead of G was used in the SEB equation, the slope of corresponding OLR of （H＋LE） against （Rn-Gs） reached 0.979, thereby the imelosure ratio of SEB was reduced to only 2.1%.

Spatial distribution of surface energy fluxes over the Loess Plateau in China and its relationship with climate and the environment

China's Loess Plateau is located at the edge of the Asian summer monsoon in a transition zone of climate and ecology. In the Loess Plateau, climate and environments change along with space, which has an obvious impact on the spatial distribution of surface energy fluxes. Because of scarce land-surface observation sites and short observation time in this area, previous studies have failed to fully understand the land-surface energy balance characteristics over the entire the Loess Plateau and their effect mechanisms. In this paper, we first test the simulation ability of the Community Land Model(CLM) model by comparing its simulated data with observed data. Based on the simulation data for the Loess Plateau over the past thirty years, we then analyze the spatial distribution of surface energy fluxes and compare the pattern differences between the area averages for the driest year and wettest year. Furthermore, we analyze the relationship between the spatial distribution of the components of the surface energy balance with longitude, latitude, altitude, precipitation and temperature. The main results are as follows: the spatial distribution of surface energy fluxes are significantly different, with the surface net radiation and sensible heat flux increasing from south to north and latent heat flux and soil heat flux decreasing from southeast to northwest. The sensible heat flux at the driest point is nearly twice as high as that at the wettest point, whereas the latent heat flux and soil heat flux at the driest point are half as much as that at the wettest point. The impact of variations of annual precipitation on the components of the surface energy balance is also obvious, and the maximum magnitude of the changes to the sensible heat flux and latent heat flux is nearly 30%. To a certain extent, geographical factors(including longitude, latitude, and altitude) and climate factors(including temperature and precipitation) affect the surface energy fluxes. However, the surface net radiation is more closely related to latitude and altitude, sensible heat flux is more closely related to the monsoon rainfall and latitude, and latent heat flux and soil heat flux are more closely related to the monsoon rainfall.

Estimation of Turbulent Sensible Heat and Momentum Fluxes over a Heterogeneous Urban Area Using a Large Aperture Scintillometer

The accurate determination of surface-layer turbulent fluxes over urban areas is critical to understanding urban boundary layer （UBL） evolution. In this study, a remote-sensing technique using a large aperture scintillometer （LAS） was investigated to estimate surface-layer turbulent fluxes over a highly heterogeneous urban area. The LAS system, with an optical path length of 2.1 km, was deployed in an urban area characterized by a complicated land-use mix （residential houses, water body, bare ground, etc.）. The turbulent sensible heat （QH） and momentum fluxes （z） were estimated from the scintillation measurements obtained from the LAS system during the cold season. Three-dimensional LAS footprint modeling was introduced to identify the source areas （＂footprint＂） of the estimated turbulent fluxes. The analysis results showed that the LAS-derived turbulent fluxes for the highly heterogeneous urban area revealed reasonable temporal variation during daytime on clear days, in comparison to the land-surface process-resolving numerical modeling. A series of sensitivity tests indicated that the overall uncertainty in the LAS-derived daytime QH was within 20%-30% in terms of the influence of input parameters and the non- dimensional similarity function for the temperature structure function parameter, while the estimation errors in z were less sensitive to the factors of influence, except aerodynamic roughness length. The 3D LAS footprint modeling characterized the source areas of the LAS-derived turbulent fluxes in the heterogeneous urban area, revealing that the representative spatial scales of the LAS system deployed with the 2.1 km optical path distance ranged from 0.2 to 2 km2 （a ＂micro-a scale＂）, depending on local meteorological conditions.

Modeling diurnal variation of ground thermal radiance images using energy balance model and endmember composing technique

Modeling and analyzing dynamic changes of land thermal radiance scenes play an important role in thermal remote sensing. In this paper, the diurnal variation of ground surface thermal scene is mainly discussed. Firstly, based on the land surface energy balance equation, the diurnal variation of land surface temperatures （LSTs） over bare land covers were simulated by an analyt- ical thermal model with second harmonic terms, and the diurnal LST variation of vegetation canopy was simulated using the Cupid model. Secondly, normalized difference vegetation index （NDVI）, normalized difference water index （NDWI）, and ratio resident-area index （RRI） were used to evaluate the endmember abundance of four land cover types including vegetation, bare soil, impervious and water area, which were calculated from IKONOS visible and near infrared （VNIR） bands. Finally, the thermal radiance scenes at various times and view angles were modeled based on the linear-energy-mixing hypothesis. The re- suits showed that the simulated daily LST variations for vegetated and bare surfaces are correlated with the measured values with a maximum standard deviation of 2.7℃, that land thermal radiant textures with high-resolution are restored from the lin- ear-energy-mixing method, and that the information abundance of the scene are related to the distribution of land cover, the imaging time, and the view angle.

Inter-comparisons of thermodynamic sea-ice modeling results using various parameterizations of radiative flux

Radiative fluxes are of primary importance in the energy and mass balance of the sea-ice cover. Various parameterizations of the radiative fluxes are studied in a thermodynamic sea-ice model. Model outputs of the surface radiative and heat fluxes and mass balance are compared with observations. The contribution of short-wave radiation is limited to a long part of winter. Therefore, simple schemes are often sufficient. Errors in estimations of the short-wave radiation are due mainly to cloud effects and occasionally to multi-reflection between surface and ice crystals in the air. The long-wave radiation plays an important role in the ice surface heat and mass balance during most part of a winter. The effect of clouds on the accuracy of the simple radiative schemes is critical, which needs further attention. In general, the accuracy of an ice model depends on that of the radiative fluxes.

Analysis of the Surface Energy Closure for a Site in the Gobi Desert in Northwest China

The heat storage terms over an ideal（non-vegetated） horizontal desert surface may be very important and easily neglected in surface energy balance studies.In this paper,based on a field experiment over the Gobi Desert in the middle part of the Hexi Corridor in Northwest China（39 05 N,100 16 E;1457-m elevation）,we studied the energy budget closure and evaluated the contribution of the heat storage terms to the closure of the surface energy balance.There were imbalances of 8% and 15% in summer and winter,respectively,if the heat storage terms were not taken into account.For both seasons,a nearly perfect result of the surface energy closure（99%） was obtained by inclusion of the estimates of heat storage terms.The soil heat storage term improved the surface energy imbalance by about 6% in summer and 13% in winter,and the air enthalpy storage term improved it by about 0.6% in summer and 1% in winter,while the contribution of the atmospheric moisture changes could be ignored.

Comparative Study on Methods for Computing Soil Heat Storage and Energy Balance in Arid and Semi-Arid Areas

Observations collected in the Badan Jaran desert hinterland and edge during 19-23 August 2009 and in the Jinta Oasis during 12-16 June 2005 are used to assess three methods for calculating the heat storage of the5-20-cm soil layer.The methods evaluated include the harmonic method,the conduction-convection method,and the temperature integral method.Soil heat storage calculated using the harmonic method provides the closest match with measured values.The conduction-convection method underestimates nighttime soil heat storage.The temperature integral method best captures fluctuations in soil heat storage on sub-diurnal timescales,but overestimates the amplitude and peak values of the diurnal cycle.The relative performance of each method varies with the underlying land surface.The land surface energy balance is evaluated using observations of soil heat flux at 5-cm depth and estimates of ground heat flux adjusted to account for soil heat storage.The energy balance closure rate increases and energy balance is improved when the ground heat flux is adjusted to account for soil heat storage.The results achieved using the harmonic and temperature integral methods are superior to those achieved using the conduction-convection method.

Enhanced modeling of latent heat flux from urban surfaces in the Noah/single-layer urban canopy coupled model

The numerical modeling of the impacts of urban buildings in mesoscale meteorological models has gradually improved in recent years. Correctly representing the latent heat flux from urban surfaces is a key issue in urban land-atmosphere coupling studies but is a common weakness in current urban canopy models. Using the surface energy balance data at a height of 140 m from a 325 m meteorological tower in Beijing, we conducted a 1-year continuous off-line simulation by using a coupled land surface model and a single-layer urban canopy model and found that this model has a relatively large systematic error for simulated latent heat flux. To improve the numerical method for modeling latent heat flux from urban surfaces, we combined observational analysis and urban land surface model to derive an oasis effect coefficient for urban green areas; to develop a temporal variation formula for water availability in urban impervious surfaces; and to specify a diurnal profile and the maximum values of anthropogenic latent heat release for four seasons. These results are directly incorporated into the urban land surface model to improve model performance. In addition, this method serves as a reference for studies in other urban areas.

Estimating evapotranspiration from terrestrial groundwater-dependent ecosystems using Landsat images

Understanding and mitigating against the impact of groundwater extraction on groundwater-dependent ecosystems(GDE)requires information of evapotranspiration(ET)of these ecosystems.In this pilot study,we tested two remotesensing methods,Surface Energy Balance Algorithms for Land(SEBAL)and Vegetation Index/Temperature Trapezoid(VITT),for ET estimation from terrestrial GDEs.Multi-temporal Landsat Thematic Mapper(TM)and Enhanced Thematic Mapper(ETM
)images were used to derive vegetation indices and land surface temperatures for ET estimation.Radiative transfer model was used for atmospheric correction of the Landsat images.Field measurements were used to validate the remote sensing estimation of VI and surface temperature.Both methods have been implemented in a geographic information system(GIS)using automated scripts and ancillary GIS data for quality control process.Comparison of predicted ET by SEBAL to VITT model indicates relatively good agreement(R2
0.90)and promise for use in groundwater management.The average ET from woodland GDEs within the zone of influence of the pumping stations is in general lower than similar woodlands outside of the pumping area,particularly in summer seasons which demonstrates that the pumping regime has an impact on those GDEs.The study also demonstrates that even a simple physical ET model can provide useful information for groundwater management,and more broadly other applications in hydrologic modelling and digital earth studies.

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%.