Estimation of spatially distributed surface energy fluxes using remotely‐sensed data for agricultural fields

Land surface energy fluxes are required in many environmental studies, including hydrology, agronomy and meteorology. Surface energy balance models simulate microscale energy exchange processes between the ground surface and the atmospheric layer near ground level. Spatial variability of energy fluxes limits point measurements to be used for larger areas. Remote sensing provides the basis for spatial mapping of energy fluxes. Remote‐sensing‐based surface energy flux‐mapping was conducted using seven Landsat images from 1997 to 2002 at four contiguous crop fields located in Polk County, northwestern Minnesota. Spatially distributed surface energy fluxes were estimated and mapped at 30 m pixel level from Landsat Thematic Mapper and Enhanced Thematic Mapper images and weather information. Net radiation was determined using the surface energy balance algorithm for land (SEBAL) procedure. Applying the two‐source energy balance (TSEB) model, the surface temperature and the latent and sensible heat fluxes were partitioned into vegetation and soil components and estimated at the pixel level. Yield data for wheat and soybean from 1997 to 2002 were mapped and compared with latent heat (evapotranspiration) for four of the fields at pixel level. The spatial distribution and the relation of latent heat flux and Bowen ratio (ratio of sensible heat to latent heat) to crop yield were studied. The root‐mean‐square error and the mean absolute percentage of error between the observed and predicted energy fluxes were between 7 and 22 W m⁻² and 12 and 24% respectively. Results show that latent heat flux and Bowen ratio were correlated (positive and negative) to the yield data. Wheat and soybean yields were predicted using latent heat flux with mean R² = 0·67 and 0·70 respectively, average residual means of −4·2 bushels/acre and 0·11 bushels/acre respectively, and average residual standard deviations of 16·2 bushels/acre and 16·6 bushels/acre respectively (1 bushel/acre ≈ 0·087 m³ ha⁻¹). The flux estimation procedure from the SEBAL‐TSEB model was useful and applicable to agricultural fields. Copyright © 2005 John Wiley & Sons, Ltd.     

Application of remote sensing to forest hydrology

The paper reviews the application of remote sensing to forest hydrology. After discussing the general advantages and disadvantages of satellite remote sensing, the estimation of precipitation, changes in soil moisture, runoff, radiation components, sensible heat flux, latent heat flux, soil heat flux, changes in energy storage in biomass, primary production and monitoring the extent, type and density of forests are reviewed. Finally, the paper looks forward to future developments and concludes that these are likely to come from the use of multitemporal data, combined analysis of different types of remotely sensed data and of remotely sensed and ground data, improved image analysis techniques and combining satellite data with models.     

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.     

A simplified climate model with combined atmospheric-hydrological processes

Abstract

Global climate change can be reproduced in detail by using three-dimensional general circulation models (GCMs). However, such complex models require super-computers and extensive hours of computational time for a single attempt at reproducing long term climate change. An alternative approach is to make simplifying assumptions that retain the essential physics for the desired simulation. Energy balance and Radiative-convective models are examples of such models. The model in this study follows the simplified approach using physics-based climate processes as well as interactions between atmospheric and hydrological processes. The vertically and latitudinally averaged mean temperature and mean water vapour content between 30°N-50°N latitudes are considered as atmospheric state variables while soil and sea temperatures and water storage amount are considered for describing the behaviour of the hydrological system. Temperatures in both the atmosphere and ground are calculated by a thermal energy equation that considers the physically-based processes of shortwave radiation, longwave radiation, sensible heat flux, and latent heat flux. Precipitation and evaporation processes transport moisture between the atmosphere and ground. In this study, the radiation parameterization of the simplified climate model is tested in the investigation of the various effects of global warming due to doubling and quadrupling of CO₂. Changes of temperature, soil water content, evaporation rate and precipitation rate are investigated by numerical experiments. The simplified climate model provides acceptable simulation of climate change and holds promise for practical investigations such as the interactions of physical processes in the evolution of drought phenomena.     

A remote sensing model of CO2 flux for wheat and studying of regional distribution

A model for calculating CO₂ flux in the wheat field and an algorithm for estimating CO₂ flux in the mejonal scale were presented using the remote sensing data and supplementary micpo-met~orological data. First of all a-longertenn measurement wae carried out during winter wheat growing period in Yucheng Experimental Station udng the spectmradiometer system, the thermal infrared radiometer system, the Bowen-ratio device as well as the eddy-correlation device. Two kinds of issues concerning remote sensing and CO₂ flux can be obtained. Based on the obeervations a remote sensing model was estabilished. Then when the NOAA-AVHRR passed over the experimental area simultaneous measurements were carried out with the satellites. A regional distribution image for CO₂ flux over wheat canopy in North China (500×500 km²) was made using the supplementary ground data and NOAA-AVHRR remote sensing data which was calibrated by the synchronous observation. The sources and sinks for CO₂ fluxes in the region can be seen obviously. Project supported by the National Natural Science Foundation of China (Grant Nos. 49671058, 49890330)     

Numerical simulation and data assimilation of the water-energy cycle over semiarid northeastern China

The default fractional vegetation cover and terrain height were replaced by the estimated fractional vegetation cover, which was calculated by the Normalized Difference Vegetation Index(NDVI) of Earth Observing System Moderate-Resolution Imaging Spectroradiometer(EOS-MODIS) and the Digital Elevation Model of the Shuttle Radar Topography Mission(SRTM) system. The near-surface meteorological elements over northeastern China were assimilated into the three-dimensional variational data assimilation system(3DVar) module in the Weather Research and Forecasting(WRF) model. The structure and daily variations of air temperature, humidity, wind and energy fields over northeastern China were simulated using the WRF model. Four groups of numerical experiments were performed, and the simulation results were analyzed of latent heat flux, sensible heat flux, and their relationships with changes in the surface energy flux due to soil moisture and precipitation over different surfaces. The simulations were compared with observations of the stations Tongyu, Naiman, Jinzhou, and Miyun from June to August, 2009. The results showed that the WRF model achieves high-quality simulations of the diurnal characteristics of the surface layer temperature, wind direction, net radiation, sensible heat flux, and latent heat flux over semiarid northeastern China in the summer. The simulated near-surface temperature, relative humidity, and wind speed were improved in the data assimilation case(Case 2) compared with control case(Case 1). The simulated sensible heat fluxes and surface heat fluxes were improved by the land surface parameterization case(Case 3) and the combined case(Case 4). The simulated temporal variations in soil moisture over the northeastern arid areas agree well with observations in Case 4, but the simulated precipitation should be improved in the WRF model. This study could improve the land surface parameters by utilizing remote sensing data and could further improve atmospheric elements with a data assimilation system. This work provides an effective attempt at combining multi-source data with different spatial and temporal scales into numerical simulations. The assimilation datasets generated by this work can be applied to research on climate change and environmental monitoring of arid lands, as well as research on the formation and stability of climate over semiarid areas.     

Energy balance above a boreal coniferous forest: a difference in turbulent fluxes between snow‐covered and snow‐free canopies

To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day⁻¹, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.     

Using a remote sensing field experiment to investigate flux-footprint relations and flux sampling distributions for tower and aircraft-based observations

During a remote sensing field experiment conducted in the Southern Great Plains in 1997 (SGP97), tower and aircraft-based flux observations were collected over one of the main study sites in central Oklahoma. This is an agricultural region and contains primarily grassland/pasture and winter wheat, which was recently harvested leaving a significant number of fields either as wheat stubble or plowed bare soil. Multi-spectral data obtained by aircraft provided high-resolution (30 m) spatially-distributed vegetation cover and surface temperature information over the study area. The spatial variations in these surface states strongly affect the partitioning of surface fluxes between sensible and latent heat. These data, together with coarser resolution (5 km) satellite data, are used in a remote sensing-based energy balance modeling system that disaggregates flux estimates from 5 km to 30-m resolution. The resulting high-resolution flux maps provide a means for evaluating whether tower and aircraft-based flux measurements sample a full range in flux conditions for this landscape. In addition, this remote sensing-based modeling system can be used to investigate the influence of variability in these key surface states on tower and aircraft measurements through flux-footprint modeling. Under the light wind and unstable conditions that existed during the observations, highest correlation between aircraft and modeled estimated heat and water vapor fluxes were obtained using different flux-footprint estimates. More specifically, the source area for heat was estimated to be much closer to the aircraft flight line than for water vapor.     

Heat discharge from summit crater of Mihara volcano, Ooshima, Japan

Using an infrared radiation thermometer, measurements by remote sensing were carried out for surface temperatures of the bottom-floor and wall of the central pit and its surrounding crater bottom of Mihara volcano, Ooshima, Japan. Isotherms for the wall of the central pit are compared with the thermal and other surface manifestations. Heat discharge from the summit crater of Mihara volcano is estimated on the basis of various meteorological and other techniques, and for convenience of computation it is divided into the following three categories; latent and sensible heat transferred by ascending volcanic gas; heat emitted from the bottom-floor and wall of the central pit; and heat transported conductively and convectively from the underground heat source through the western, southern and eastern parts of the crater bottom. The computed heat discharge is estimated to be 1.4 x 10⁷ cal s⁻¹, which is about one order of magnitude less than that of the minor eruption on 28 February 1974. The methodology developed here for estimation of the heat discharge from a volcano can also be applied to other volcanoes.     

Evapotranspiration estimation considering anthropogenic heat based on remote sensing in urban area

Urbanization influences hydrologic cycle significantly on local, regional even global scale. With urbanization the water resources demand for dense population sharpened, thus it is a great challenge to ensure water supply for some metropolises such as Beijing. Urban area is traditionally considered as the area with lower evapotranspiration (ET) on account of the impervious surface and the lower wind speed. For most remote sensing models, the ET, defined as latent heat in energy budget, is estimated as the difference between net radiation and sensible heat. The sensible heat is generally higher in urban area due to the high surface temperature caused by heat island, therefore the latent heat (i.e. the ET) in urban area is lower than that in other region. We estimated water consumption from 2003 to 2012 in Beijing based on water balance method and found that the annual mean ET in urban area was about 654 mm. However, using Surface Energy Balance System (SEBS) model, the annual mean ET in urban area was only 348 mm. We attributed this inconsistence to the impact of anthropogenic heat and quantified this impact on the basis of the night-light maps. Therefore, a new model SEBS-Urban, coupling SEBS model and anthropogenic heat was developed to estimate the ET in urban area. The ET in urban area of Beijing estimated by SEBS-Urban showed a good agreement with the ET from water balance method. The findings from this study highlighted that anthropogenic heat should be included in the surface energy budget for a highly urbanized area.     

Mapping of H2S fluxes from the ground using copper passive samplers: An application study at the Zolforata di Pomezia degassing area (Alban Hills,Central Italy)

A new method for measuring H₂S mass flux from the ground, based on the digital analysis of the interference colours produced by the sulphidation of copper passive samplers (CPS), is proposed and discussed in this article. CPS sulphidation has a wide range of linear responses to H₂S doses and can be used together the accumulation chamber method to estimate gas fluxes from natural degassing areas. These are often characterized by the presence of vent centred degassing areas (VCDAs), which are recognizable from the absence or rarefaction of vegetation due to high acid gas concentrations in the soil pores and in the air at ground level. A reference emission curve, accounting for the advective and diffusive components of the flux, can be modelled and used to estimate the total H₂S mass released from each VCDA. The application of this method can be supported by remote sensing analysis that helps identify VCDAs in the field in perivolcanic H₂S degassing areas.As an illustrative application, H₂S gas fluxes from the ground were measured in spring 2007 at the Zolforata di Pomezia degassing area (ZPDA, Alban Hills, Central Italy) using an accumulation chamber internally equipped with CPS. H₂S peak fluxes were measured over the vents after remote sensing assisted identification of the VCDAs. Further measurements were carried out in two ponds and one artificial channel bordering the study area. The total atmospheric flux released at the ZPDA, estimated to be about 1207.6 kg day^? 1, was calculated as the summation of the fluxes from all the H₂S sources, the background flux being negligible.     

Seasonal variations between sampling and classical mean turbulent heat flux estimates in the eastern North Atlantic

The two commonly used statistical measures of the air-sea heat flux, the sampling and classical means, have been compared using hourly reports over a 7-year-period from a weather ship stationed in the NE Atlantic. The sampling mean is the average over all flux estimates in a given period, where individual flux estimates are determined from ship reports of meteorological variables using the well-known bulk formulae. The classical mean is the flux derived by substituting period-averaged values for each of the meteorological variables into the bulk formula (where the averaging period employed is the same as that over which the fluxes are to be determined). Monthly sampling and classical means are calculated for the latent and sensible heat fluxes. The monthly classical mean latent heat flux is found to overestimate the sampling mean by an amount which increases from 1–2 W m⁻² in summer to 7 W m⁻² in winter, on average, over the 7-year-period. In a given winter month, the excess may be as great as 15 W m⁻², which represents about 10% of the latent heat flux. For the sensible heat flux, any seasonal variation between the two means is of the order of 1 W m⁻² and is not significant compared to the interannual variation. The discrepancy between the two means for the latent heat flux is shown to arise primarily from a negative correlation between the wind speed and sea-air humidity difference, the effects of which are implicitly included in the sampling method but not in the classical. The influence of the dominant weather conditions on the sign and magnitude of this correlation are explored, and the large negative values that it takes in winter are found to depend on the typical track of the mid-latitude depressions with respect to the position sampled. In conclusion, it is suggested that sampling means should be employed where possible in future climatological studies.     

A remote sensing model for monitoring soil evaporation based on differential thermal inertia and its validation

To develop geosciences quantification and multi-dimensional researches will be an inevitable trend in the 21st century. The interaction between the land surface and the atmosphere not only serves as an important component in geosciences quantification, bu…     

基于双温度差分双源模型的遥感方法反演四川汶川地区的潜热通量

研究表明,在很多地震前都出现了潜热通量异常。文中基于Modis遥感数据,结合研究区(四川省汶川地区)的高光谱数据及地面观测的气象数据,利用一种简化的双层模型(双温度差分双源模型)反演显热通量,并通过地表能量平衡方程计算研究区的潜热通量。还探讨了该模型的特点,并对龙门山断裂带通过地区进行了潜热通量的遥感反演,发现在2008年5月12日汶川MS8.0地震前48h内,研究区内的震中地区、龙门山断裂周边地区潜热通量值出现了明显的高值,5月12日中午12时7分沿龙门山断裂带潜热通量值达到300～400W/m2,可能是汶川8.0级地震前的异常反应     

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.     

Remote sensing parameterization of the processes of energy and water cycle over desertification areas

In order to understand the processes of land surface-atmosphere interaction over de-sertification area, it is indispensable to utilize of satellite remote sensing. Two scenes of LandsatTM were used to produce a set of maps of surface reflectance, MSAVI, vegetation coverage, sur-face temperature, net radiation, soil heat flux, sensible heat flux and latent heat flux. Statisticalanalysis based on these maps revealed some quantitative significant land surface characteristics.Future developments of the method are also discussed.     

A land surface model incorporated with soil freeze/thaw and its application in GAME/Tibet

Land surface process is of great importance in global climate change, moisture and heat exchange in the interface of the earth and atmosphere, human impacts on the environment and eco- system, etc. Soil freeze/thaw plays an important role in cold land surface processes. In this work the diurnal freeze/thaw effects on energy partition in the context of GAME/Tibet are studied. A sophisti- cated land surface model is developed, the particular aspect of which is its physical consideration of soil freeze/thaw and vapor flux. The simultaneous water and heat transfer soil sub-model not only reflects the water flow from unfrozen zone to frozen fringe in freezing/thawing soil, but also demon- strates the change of moisture and temperature field induced by vapor flux from high temperature zone to low temperature zone, which makes the model applicable for various circumstances. The modified Picard numerical method is employed to help with the water balance and convergence of the numerical scheme. Finally, the model is applied to analyze the diurnal energy and water cycle char- acteristics over the Tibetan Plateau using the Game/Tibet datasets observed in May and July of 1998. Heat and energy transfer simulation shows that: (i) There exists a negative feedback mechanism between soil freeze/thaw and soil temperature/ground heat flux; (ii) during freezing period all three heat fluxes do not vary apparently, in spite of the fact that the negative soil temperature is higher than that not considering soil freeze; (iii) during thawing period, ground heat flux increases, and sensible heat flux decreases, but latent heat flux does not change much; and (iv) during freezing period, soil temperature decreases, though ground heat flux increases.     

Remote sensing parameterization of the processes of energy and water cycle over desertification areas

In order to understand the processes of land surface-atmosphere interaction over desertification area, it is indispensable to utilize of satellite remote sensing. Two scenes of Landsat TM were used to produce a set of maps of surface reflectance, MSAVI, vegetation coverage, surface temperature, net radiation, soil heat flux, sensible heat flux and latent heat flux. Statistical analysis based on these maps revealed some quantitative significant land surface characteristics. Future developments of the method are also discussed.     

Water and heat fluxes above a lowland dipterocarp forest in Peninsular Malaysia

We measured the fluxes of sensible and latent heat between a low‐land dipterocarp forest in Peninsular Malaysia and the atmosphere. No clear seasonal or interannual changes in latent heat flux were found from 2003 to 2005, while sensible heat flux sometimes fluctuated depending on the fluctuation of incoming radiation between wet and dry seasons. The evapotranspiration rates averaged for the period between 2003 and 2005 were 2·77 and 3·61 mm day^?1 using eddy covariance data without and with an energy balance correction, respectively. Average precipitation was 4·74 mm day^?1. Midday surface conductance decreased with an increasing atmospheric water vapour pressure deficit and thus restricted the excess water loss on sunny days in the dry season. However, the relationship between the surface conductance and vapour pressure deficit did not significantly decline with an increase in volumetric soil water content even during a period of extremely low rainfall. Copyright © 2009 John Wiley & Sons, Ltd.