Combining FAO Model and Vegetation Indices to Estimate Crop Coefficient Using Principle Component Analysis

آب و خاک

Article 2, Volume 31, Issue 5 - Serial Number 55, November 2018, Pages 1290-1301   

Document Type: Research Article

DOI: 10.22067/jsw.v31i5.63990

Author

Laleh Parviz^*

Azarbaijan Shahid Madani University

Abstract

Introduction: The globally growing demand for water has shown the need for its efficient and judicial utilization, and particularly in agriculture being single largest consumer of water. Crop evapotranspiration represents crop water demand and governed by weather and crop conditions and most of the current water demand models are non-spatial models, they use point data. Global scale satellite images can solve these problems. According to the high performance of satellite indices, it is necessary to estimate crop coefficient using combination of reflectance and thermal bands. The aim of this research was to estimate the effective crop coefficient of potato using vegetation indices and principle component analysis.
Materials and Methods: Principle component analysis (PCA) was used for effective crop coefficient estimation. Modeling of associations between vegetation indices and crop coefficient were conducted using artificial neural network. In the present study, NDVI, RI, EVI, SAVI, MSAVI, NVSWI, TVX, TVI, mNDVI and mTVI were the used as vegetation indices. PCA is designed to transform the original variables into new and uncorrelated variables (axes), namely the principal components, which are linear combinations of the original variables. The new axes lie along the directions of maximum variance. PCA provides an objective procedure of finding indices and information on the most meaningful parameters, which describes a whole data set affording data reduction with minimum loss of original information. Artificial neural networks are a computational model which is based on a large collection of simple neural units, loosely analogous to the observed behavior of a biological brain's axons. RMSE, MAE and MARE were the statistics used for investigating the performance of crop coefficient of vegetation indices with FAO crop coefficient.
Results and Discussion: Eleven MODIS vegetation indices are derived in the period of 2013 to 2016 for potato over the limited area between Tabriz and Bostanabad. The last year was considered as the validation period. According to the FAO-56 paper, the lengths of initial stage, crop development stage, mid-season stage, late season stage were considered to be 25, 30, 45, 30 days, respectively. The vegetation indices were derived using MODIS sensor with 2×2 pixels. The PCA showed that with increasing the number of components, the eigenvalues decreased. The analysis indicated that the three first components accounted for the 85.45 % of the total variance of data and the eigenvalues of them were greater than 1, the three first components were thus selected. NDVI, RI, TVI, MSAVI and NVSWI in the first component, mNDVI in the second component and LST in the third component had the highest coefficients. NDVI in the first component with high coefficient indicted the importance of index in the crop coefficient determination. The coefficients of SAVI and MSAVI were higher than NDVI. From the three investigations on the kind of principle component, the first investigation led to a 55.75 % decrease in RMSE relative to the second and third investigations. The first and second components together had less error rather than third component. The average of MAE for first, second and third investigations was, respectively, 0.17, 0.22 and 0.2. Therefore, component with exact values of particular vectors resulted in a reduced error. The sensitivity of artificial neural network led to an increase in the simulation accuracy (for example the RMSE decreased from epoch 100 to 50 was 48.27%). Crop coefficient estimation using vegetation indices of principle component analysis was underestimated about 1% in the validation period. Overestimation and underestimation were found in the initial and crop development stages, respectively.
Conclusions: The quantities of statistics showed the improvement of artificial neural network performance with combination of vegetation indices using principle component analysis. The vegetation indices with reflectance bands performed well. The combination of thermal and reflectance bands enhanced the vegetation indices efficiency. In addition to NDVI index for crop coefficient estimation, improvement of indices according to the study area condition increased the indices performance. The kind of mathematical equations of indices can increase the indices performance which using the same bands with different equations have different results. The selected component of principle component analysis has important role in increasing the simulation accuracy. The error reduction of simulated crop coefficients can increase the precision of irrigation consumption and agricultural planning which the principle component analysis has more important role.

Keywords

Artificial neural network; Components; SAVI; Water demand

References

1- Adamala S., Rajwade Y.A., and Krishna Reddy Y.V. 2016. Estimation of wheat crop evapotranspiration using NDVI vegetation index. Journal of Applied and Natural Science, 8(1):159-166.

2- Agapiou A., Hadjimitsis D.G., and Alexakis D.D. 2012. Evaluation of broadband and narrowband vegetation indices for the identification of archaeological crop marks. Remote Sensing, 4:3892-3919.

3- Amirmoradi K., and Bahmani O. 2014. Estimation the daily soil temperate using neural artificial network. Journal of Soil Research (Soil Science and Water), 28(3):543-556. (In Persian)

4- El-Shirbeny M.A., Ali A-E M., Badr M.A., and Bauomy E.M. 2014. Assessment of wheat crop coefficient using remote sensing techniques. World Research Journal of Agricultural Sciences, 1(2):012-016.

5- Er-Raki S., Chehbouni A., Guemouria N., Duchemin B., Ezzahar J., and Hadria R. 2007. Combining FAO-56 model and ground-based remote sensing to estimate water consumptions of wheat crops in a semi-arid region. Agriculture Water Management, 87(1):41-54.

6- Gholami M.A., Mirlatifi S.M., Mohammadi K., and Alimohammadi A. 2010. Estimating crop coefficient and actual evapotranspiration of wheat by remote sensing: A case study, gorgan rood command area, Golestan, Iran. Iranian Journal of lrrigation and Drainage, 2(4):222-231. (In Persian with English abstract)

7- Gontia N.K., and Tiwari K.N. 2010. Estimation of crop coefficient and evapotranspiration of wheat (Triticum aestivum) in an irrigation command using remote sensing and GIS. Water Resources Managemnet, 24:1399-1414.

8- Jiang Z., Huete A.R., Didan K., and Miura T. 2008. Development of a two-band enhanced vegetation index without a blue band. Remote Sensing of Environment, 112: 3833–3845.

9- Kamble B., Kilic A., and Hubbard K. 2013. Estimating crop coefficients using remote sensing-based vegetation index. Remote Sensing, 5:1588-1602.

10- Karamouz M., and Araghinejad S. 2006. Advanced hydrology. 2nd ed. Amirkabir University of Technology, Tehran. (In Persian)

11- Mohammadi H., Ahmadi A., Moradi F., Abbasi A., Poustini K., Joudi M., and Fatehi F. 2011. Evaluation of critical traits for improving wheat yield under drought stress. Iranian Journal of Agricultural Science, 42(2):373-385. (In Persian)

12- Pimpale A.R., Rajankar P.B., Wadatkar S.B., and Ramteke I.K. 2014. Determination of spatial crop coefficient of chickpea using remote sensing and GIS. American International Journal of Research in Formal, Applied and Natural Sciences, 59-64.

13- Price J.C. 1984. Land surface temperature measurements from the split window channel of the NOAA 7 advanced very high resolution radiometer. Journal of Geophysical Research, 89(5):7231–7237.

14- Richard G.A., Luis S.P., Dirk R., and Martin S. 2006. Crop evapotranspiration (guidelines for computing crop water requirements). FAO Irrigation and Drainage, Paper No. 56.

15- Rostami A., and Raeini-Sarjanz M. 2016. Remotely sensed measurements of apple orchard actual evapotranspiration and plant coefficient using MODIS images and SEBAL algorithm (Case Study: Ahar plain, Iran). Journal of Agricultural and Meteorology, 4(1):32-43. (In Persian with English abstract)

16- Shahabfar A., Ghulam A., and Eitzinger J. 2012. Drought monitoring in Iran using the perpendicular drought indices. International Journal of Applied Earth Observation and Geoinformation, 18:119–127.

17- Sheikholeslami N., Ghahraman B., Mosaedi A., Davary K., and Mohejerpour M. 2014. Estimating reference evapotranspiration by using principal component analysis (PCA) and the development of a regression model (MLR-PCA) (Case Study: Mashhad Station). Journal of Water and Soil, 28(2):420-429. (In Persian with English abstract)