Estimation of Daily Actual Evapotranspiration using Microwave and Optical Vegetation Indices for Clear and Cloudy Sky Conditions

Abstract

Evapotranspiration (ET) is a significant hydrological process. It can be studied and estimated using remote sensing based methods at multiple spatial and temporal scales. Most commonly and widely used remote sensing based methods to estimate actual evapotranspiration (AET) are a) methods based on energy balance equations, b) vegetation coefficient based method and c) contextual methods. These three methods require reflectance and land surface temperature (LST) data measured at optical and thermal portion of the electromagnetic spectrum. However, these data are available only for clear sky conditions and fail to be retrieved under overcast conditions creating gaps in the data, which result in discontinuous of AET product. Moreover, energy balance equation based methods and evaporative fraction (EF) based contextual methods are difficult to apply over overcast conditions. In this context, vegetation coefficient based (Tasumi et al., 2005; Allen et al., 2005) and microwave remote sensing based methods can be applied under cloudy sky conditions (Sun et al., 2012), since microwave radiations can penetrate through clouds, but these data are available at coarse resolution. In the vegetation coefficient method temporal upscaling can be avoided. Therefore in this research vegetation coefficient based method is employed over Cauvery basin to estimate daily AET for clear and cloudy sky conditions. Required critical variables for this method such as reference evapotranspiration (ETo) and vegetation coefficients are obtained using LST and optical vegetation indices for all sky conditions. In this study, all sky conditions refer to both clear and cloudy sky conditions.

Most important variable for estimation of ETo using radiation and temperature based models is air temperature (Ta). In this study, for better accuracy of Ta, two satellite based approaches namely, Temperature Vegetation Index (TVX) and Advance Statistical

Approaches (ASA) were evaluated. In the TVX approach, in addition to traditional Normalized Difference Vegetation Index (NDVI), other vegetation indices such as Enhanced Vegetation Index (EVI) and Global Vegetation Moisture Index (GVMI) were also examined. In case of ASA, bootstrap technique was used to generate calibration and validation samples and Levenberg Marquardt algorithm was used to find the solution of the models. The better of the Ta results obtained out of these two approaches were employed in the ETo models and are referred as Ta based ETo models. Instead of Ta, processed LST data obtained directly from the satellite (Aqua/Moderate Resolution Imaging Spectroradiometer (MODIS)) was applied in the ETo models and these are referred as LST based ETo models. These Ta and LST based Hargreaves-Samani (H-S), Makkink (Makk) and Penman Monteith Temperature (PMT) models were evaluated by comparing with the FAO56 PM model. Additionally, simple LST based equation (SLBE) proposed by Rivas et al. (2004) was also examined. Required solar radiation (Rs) data for ETo estimation was obtained from Kalpana1/VHRR satellite data. Results implied that, Ta based PMT model performed better than the Ta based H-S, Makk and SLBE with less RMSE, MAPE and MBE values for all land cover classes and for various climatic regions for clear sky conditions. LST based H-S, PMT, Makk and Ta based Makkink advection models predominantly overestimated ETo for the study region. In the case of TVX approach, to estimate maximum Ta (Tmax), GVMI performed better than NDVI and EVI. Nevertheless, TVX approach poorly estimated Tmax in comparison with statistical approach. ASA performed better for both Tmax and minimum Ta. This study demonstrates the applicability of satellite based Ta and ETo models by considering very few variables for clear sky conditions.

Spatially distributed vegetation coefficients (Kv) data with high temporal resolution is another important variable in vegetation coefficient method for daily AET estimation and also it is in demand for crop condition assessment, irrigation scheduling, etc. But available Kv models application hinders because of two main reasons i.e 1) Spectral reflectance based Kv accounts only for transpiration factor but not evaporation, which fails to account for total AET. 2) Required optical spectral reflectances are available only during clear sky conditions, which creates gaps in the Kv data. Hence there is a necessity of a model which accounts for both transpiration and evaporation factors and also gap filling method, which produces accurate continuous quantification of Kv values. Therefore, different combinations of EVI, GVMI and temperature vegetation dryness index (TVDI) have been employed in linear and non linear regression techniques to obtain best model. This best Kv model had been compared with Guershman et al. (2009) Kv model. To fill the gaps in the data, initially, temporal fitting of Kv values have been examined using Savitsky-Goley (SG) filter for three years of data (2012 to 2014), but this fails when sufficient high quality Kv values were unavailable. In this regard, three gap filling techniques namely regression, Artificial Neural Networks (ANNs) and interpolation techniques have been analyzed. Microwave polarization difference index (MPDI) has been employed in ANN technique to estimate Kv values under cloudy sky conditions. The results revealed that the combination of GVMI and TVDI using linear regression technique performed better than other combinations and also yielded better results than Guershman et al. (2009) Kv model. Furthermore, the results indicated that SG filter can be used for temporal fitting and for filling the gaps, regression technique can be used as it performed better than other techniques for Berambadi station.

Land Surface Temperature (LST) with high spatiotemporal resolution is required in the estimation of ETo to obtain AET. MODIS is one of the most commonly used sensors owing to its high spatial and temporal availability over the globe, but is incapable of providing LST data under cloudy conditions, resulting in gaps in the data. In contrast, microwave measurements have a capability to penetrate under clouds. The current study proposes a methodology by exploring this property to predict high spatiotemporal resolution LST under cloudy conditions during daytime and night time without employing in-situ LST measurements. To achieve this, ANN based models were employed for different land cover classes, utilizing MPDI at finer resolution with ancillary data. MPDI was derived using resampled (from 0.250 to 1 km) brightness temperatures (Tb) at 36.5 GHz channel of dual polarization from Advance Microwave Scanning Radiometer (AMSR)-Earth Observing System and AMSR2 sensors. The proposed methodology was quantitatively evaluated through three performance measures namely correlation coefficient (r), Nash Sutcliffe Efficiency (NSE) and Root Mean Square Error (RMSE). Results revealed that during daytime, AMSR-E(AMSR2) derived LST under clear sky conditions corresponds well with MODIS LST resulting in values of r ranging from 0.76(0.78) to 0.90(0.96), RMSE from 1.76(1.86) K to 4.34(4.00) K and NSE from

0.58(0.61) to 0.81(0.90) for different land cover classes. For night time, r values ranged from 0.76(0.56) to 0.87(0.90), RMSE from 1.71(1.70) K to 2.43(2.12) K and NSE from 0.43 (0.28) to 0.80(0.81) for different land cover classes. RMSE values found between predicted LST and MODIS LST during daytime under clear sky conditions were within acceptable limits. Under cloudy conditions, results of microwave derived LST were evaluated with Ta which indicated that the approach performed well with RMSE values lesser than the results obtained under clear sky conditions for land cover classes for both day and nighttimes. These predicted LSTs can be applied for the estimation of soil

moisture in hydrological studies, in climate studies, ecology, urban climate and environmental studies, etc.

AET was estimated for all sky conditions using vegetation coefficient method. Essential parameter ETo under cloudy conditions was estimated using LST and Ta based PMT and H-S models and required solar radiation (Rs) in these two models estimated using equation proposed by Samani (2000). In this equation it was found that the differences between LSTmax or Tmax and LSTmin or Tmin could able to capture the variations due to cloudy sky conditions and hence can be used for estimating ETo under cloudy sky conditions. Results revealed that the estimated Rs correlated well with observed Rs for Berambadi station under cloudy conditions for the year 2013. PMT based ETo values were corresponded with observed ETo under cloudy sky condition. The difference between LST and Ta was less during cloudy conditions, therefore LST or Ta can be used as the only input in temperature based PMT model to estimate ETo. AET estimated correlated well with the observed AET values for clear and cloudy sky conditions. In addition, AET estimated using vegetation coefficient method was compared with two source energy balance (TSEB) method developed by Nishida et al. (2003) under clear sky conditions. It was found that the improved vegetation coefficient method performed better than the TSEB method for Berambadi station.

Other microwave vegetation indices such as Microwave Vegetation Indices (MVIs) and Emissivity Difference Vegetation Index (EDVI) are available in literature. Therefore in this study, MVIs are used to predict LST under cloudy conditions using proposed methodology to check whether the MVIs could yield better LST values. Results showed that MPDI performed better than MVIs to predict LST under cloudy sky conditions. Furthermore, MPDI obtained using dual polarizations of 37 GHz channel Tb has advantage

of having fine spatial resolution compared to MVIs, as it requires Tb of 19 GHz in addition to Tb of 37 GHz channel which is of coarse resolution and therefore uncertainties resulting from re-sampling technique can be minimized.

x