Mapping regional turbulent heat fluxes via variational assimilation of land surface temperature data from polar orbiting satellites

摘要

Estimation of turbulent heat fluxes by assimilating sequences of land surface temperature (LST) measurements into variational data assimilation (VDA) frameworks has been the subject of several studies. The VDA approaches estimate turbulent heat fluxes by minimizing the difference between LST observations and estimations from the heat diffusion equation. The VDA methods have been tested only with high temporal resolution LST observations (e.g., from geostationary satellites) when applied at regional scales. Geostationary satellites can capture the diurnal cycle of LST, but they have a relatively low spatial resolution and mainly focus on low latitudes. To overcome these shortcomings, this study assimilates high spatial resolution LST data from polar orbiting satellites (e.g., Moderate Resolution Imaging Spectroradiometer, MODIS) into the combined-source (CS) and dual-source (DS) VDA schemes. An expression is developed to obtain an a priori evaporative fraction (EF) estimate from leaf area index (LAI) or apparent thermal inertia (ATI). The a priori EF estimate is used as an initial guess in the VDA approach. The results indicate that the VDA method is able to find the optimal value of EF by assimilating the low-temporal resolution MODIS LST data. The predicted turbulent heat fluxes from VDA are compared with the measurements from the large-aperture scintillometer at three sites (Arou, Daman, and Sidaoqiao) in the Heihe River Basin (located in northwest China). The findings indicate that the CS and DS VDA models perform well in various hydrological and vegetative conditions. The three-site-average root mean square errors (RMSEs) of sensible and latent heat fluxes estimates from the CS scheme are 37.44 W m(-2) and 94.30 W m(-2), respectively. The DS model reduces the above mentioned RMSEs by 19.82% and 21.37%, respectively. Overall, the results show that using the a priori EF estimate from the proposed expression in the VDA approach eliminates the need for the high reso