Simultaneous inversion of multiple land surface parameters from MODIS optical-thermal observations

In this paper, we propose a unified algorithm for simultaneously retrieving six land surface parameters - Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), land surface albedo, Land Surface Emissivity (LSE), Land Surface Temperature (LST), and Upwelling Longwave radiation (LWUP) by exploiting MODIS visible-to-TIR observations. We incorporate a unified physical radiative transfer model into a data assimilation framework. The MODIS visible-to-TIR time series datasets include the daily surface reflectance product and MIR-to-TIR surface radiance, which are atmospherically corrected from the MODIS data using the Moderate Resolution Transmittance program (MODTRAN, ver. 5.0). LAI was first estimated using a data assimilation method that combines MODIS daily reflectance data and a LAI phenology model, and then the LAI was input to the unified radiative transfer model to simulate spectral surface reflectance and surface emissivity for calculating surface broadband albedo and emissivity, and FAPAR. LST was estimated from the MIR-TIR surface radiance data and the simulated emissivity, using an iterative optimization procedure. Lastly, LWUP was estimated using the LST and surface emissivity. The retrieved six parameters were extensively validated across six representative sites with different biome types, and compared with MODIS, GLASS, and GlobAlbedo land surface products. The results demonstrate that the unified inversion algorithm can retrieve temporally complete and physically consistent land surface parameters, and provides more accurate estimates of surface albedo, LST, and LWUP than existing products, with R2 values of 0.93 and 0.62, RMSE of 0.029 and 0.037, and BIAS values of 0.016 and 0.012 for the retrieved and MODIS albedo products, respectively, compared with field albedo measurements; R2 values of 0.95 and 0.93, RMSE of 2.7 and 4.2 K, and BIAS values of −0.6 and −2.7 K for the retrieved and MODIS LST products, respectively, compared with field LST measurements; and R2 values of 0.93 and 0.94, RMSE of 18.2 and 22.8 W/m2, and BIAS values of −2.7 and −14.6 W/m2 for the retrieved and MODIS LWUP products, respectively, compared with field LWUP measurements.