How Universal Is the Relationship between Remotely Sensed Vegetation Indices and Crop Leaf Area Index? A Global Assessment

Leaf Area Index (LAI) is a key variable that bridges remote sensing observations to the quantification of agroecosystem processes. In this study, we assessed the universality of the relationships between crop LAI and remotely sensed Vegetation Indices (VIs). We first compiled a global dataset of 1459 in situ quality-controlled crop LAI measurements and collected Landsat satellite images to derive five different VIs including Simple Ratio (SR), Normalized Difference Vegetation Index (NDVI), two versions of the Enhanced Vegetation Index (EVI and EVI2), and Green Chlorophyll Index (CIGreen). Based on this dataset, we developed global LAI-VI relationships for each crop type and VI using symbolic regression and Theil-Sen (TS) robust estimator. Results suggest that the global LAI-VI relationships are statistically significant, crop-specific, and mostly non-linear. These relationships explain more than half of the total variance in ground LAI observations (R2 >0.5), and provide LAI estimates with RMSE below 1.2 m2/m2. Among the five VIs, EVI/EVI2 are the most effective, and the crop-specific LAI-EVI and LAI-EVI2 relationships constructed by TS, are robust when tested by three independent validation datasets of varied spatial scales. While the heterogeneity of agricultural landscapes leads to a diverse set of local LAI-VI relationships, the relationships provided here represent global universality on an average basis, allowing the generation of large-scale spatial-explicit LAI maps. This study contributes to the operationalization of large-area crop modeling and, by extension, has relevance to both fundamental and applied agroecosystem research.

[1]  Weiwei Shao,et al.  Improving leaf area index simulation of IBIS model and its effect on water carbon and energy—A case study in Changbai Mountain broadleaved forest of China , 2015 .

[2]  Hiroshi Matsuyama,et al.  Improving the estimation of leaf area index by using remotely sensed NDVI with BRDF signatures , 2010 .

[3]  A. J. Stern,et al.  Application of MODIS derived parameters for regional crop yield assessment , 2005 .

[4]  T. Carlson,et al.  On the relation between NDVI, fractional vegetation cover, and leaf area index , 1997 .

[5]  F. Baret,et al.  Potentials and limits of vegetation indices for LAI and APAR assessment , 1991 .

[6]  B. Markham,et al.  Summary of Current Radiometric Calibration Coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI Sensors , 2009 .

[7]  J. Hill,et al.  Use of coupled canopy structure dynamic and radiative transfer models to estimate biophysical canopy characteristics , 2005 .

[8]  Yann Kerr,et al.  The AACES field experiments: SMOS calibration and validation across the Murrumbidgee River catchment , 2012 .

[9]  N. Goel,et al.  Influences of canopy architecture on relationships between various vegetation indices and LAI and Fpar: A computer simulation , 1994 .

[10]  Hod Lipson,et al.  Distilling Free-Form Natural Laws from Experimental Data , 2009, Science.

[11]  David J. Harding,et al.  Amazon forests maintain consistent canopy structure and greenness during the dry season , 2014, Nature.

[12]  A. Cutini,et al.  Estimation of leaf area index with the Li-Cor LAI 2000 in deciduous forests , 1998 .

[13]  Samuel C. Zipper,et al.  Untangling the effects of shallow groundwater and soil texture as drivers of subfield‐scale yield variability , 2015 .

[14]  N. C. Strugnell,et al.  First operational BRDF, albedo nadir reflectance products from MODIS , 2002 .

[15]  K. Moffett,et al.  Remote Sens , 2015 .

[16]  Gerald J. Hahn,et al.  Applied Regression Analysis (2nd Ed.) , 2012 .

[17]  A. Viña,et al.  Comparison of different vegetation indices for the remote assessment of green leaf area index of crops , 2011 .

[18]  Yuri A. Gritz,et al.  Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. , 2003, Journal of plant physiology.

[19]  N. Breda.  Ground-based measurements of leaf area index: a review of methods, instruments and current controversies. , 2003, Journal of experimental botany.

[20]  J. Goudriaan,et al.  Monitoring rice reflectance at field level for estimating biomass and LAI , 1998 .

[21]  C. Justice,et al.  Development of vegetation and soil indices for MODIS-EOS , 1994 .

[22]  S. Liang,et al.  Validation of MODIS and CYCLOPES LAI products using global field measurement data , 2012 .

[23]  Edward J. Kim,et al.  The NAFE'05/CoSMOS Data Set: Toward SMOS Soil Moisture Retrieval, Downscaling, and Assimilation , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[24]  D. Deering.  Rangeland reflectance characteristics measured by aircraft and spacecraft sensors , 1979 .

[25]  Ghassem R. Asrar,et al.  Estimates of leaf area index from spectral reflectance of wheat under different cultural practices and solar angle , 1985 .

[26]  M. Guérif,et al.  Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications , 2005 .

[27]  Rasmus Houborg,et al.  Retrieving Leaf Area Index From Landsat Using MODIS LAI Products and Field Measurements , 2014, IEEE Geoscience and Remote Sensing Letters.

[28]  J. Tenhunen,et al.  On the relationship of NDVI with leaf area index in a deciduous forest site , 2005 .

[29]  A. Kuusk.  A Markov chain model of canopy reflectance , 1995 .

[30]  John Fox,et al.  Applied Regression Analysis and Generalized Linear Models , 2008 .

[31]  S. T. Gower,et al.  Direct and Indirect Estimation of Leaf Area Index, fAPAR, and Net Primary Production of Terrestrial Ecosystems , 1999 .

[32]  C. Tucker,et al.  Vegetation dynamics and rainfall sensitivity of the Amazon , 2014, Proceedings of the National Academy of Sciences.

[33]  A. Huete,et al.  A comparison of vegetation indices over a global set of TM images for EOS-MODIS , 1997 .

[34]  J. Fox,et al.  Applied Regression Analysis and Generalized Linear Models , 2008 .

[35]  J. Hicke,et al.  Global synthesis of leaf area index observations: implications for ecological and remote sensing studies , 2003 .

[36]  Joshua Gray,et al.  Mapping leaf area index using spatial, spectral, and temporal information from multiple sensors , 2011 .

[37]  Jörg Kaduk,et al.  Comparing global models of terrestrial net primary productivity (NPP): importance of vegetation structure on seasonal NPP estimates , 1999 .

[38]  Alessandro Anav,et al.  Global Data Sets of Vegetation Leaf Area Index (LAI)3g and Fraction of Photosynthetically Active Radiation (FPAR)3g Derived from Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI3g) for the Period 1981 to 2011 , 2013, Remote. Sens..

[39]  Pete Smith,et al.  The net biome production of full crop rotations in Europe , 2010 .

[40]  Roshanak Darvishzadeh,et al.  Inversion of a Radiative Transfer Model for Estimation of Rice Canopy Chlorophyll Content Using a Lookup-Table Approach , 2012, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[41]  Brigitte Leblon,et al.  Non-destructive estimation of potato leaf chlorophyll from canopy hyperspectral reflectance using the inverted PROSAIL model , 2007, International Journal of Applied Earth Observation and Geoinformation.

[42]  Anatoly A. Gitelson,et al.  An evaluation of MODIS 8- and 16-day composite products for monitoring maize green leaf area index , 2012 .

[43]  S. Liang,et al.  Real-time retrieval of Leaf Area Index from MODIS time series data , 2011 .

[44]  Adrian V. Rocha,et al.  Advantages of a two band EVI calculated from solar and photosynthetically active radiation fluxes , 2009 .

[45]  J. Cihlar,et al.  Radiometric normalization of multitemporal high-resolution satellite images with quality control for land cover change detection , 2002 .

[46]  A. Gitelson,et al.  Relationships between gross primary production, green LAI, and canopy chlorophyll content in maize: Implications for remote sensing of primary production , 2014 .

[47]  S. Liang,et al.  Retrieving leaf area index using a genetic algorithm with a canopy radiative transfer model , 2003 .

[48]  J. Dungan,et al.  Generating global Leaf Area Index from Landsat: Algorithm formulation and demonstration , 2012 .

[49]  A. Huete,et al.  Development of a two-band enhanced vegetation index without a blue band , 2008 .

[50]  R. Fernandes,et al.  Parametric (modified least squares) and non-parametric (Theil–Sen) linear regressions for predicting biophysical parameters in the presence of measurement errors , 2005 .

[51]  A. Gitelson.  Wide Dynamic Range Vegetation Index for remote quantification of biophysical characteristics of vegetation. , 2004, Journal of plant physiology.

[52]  P. J. Curran,et al.  Multispectral remote sensing for the estimation of green leaf area index , 1983, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[53]  Frédéric Baret,et al.  GEOV1: LAI and FAPAR essential climate variables and FCOVER global time series capitalizing over existing products. Part1: Principles of development and production , 2013 .

[54]  Crystal B. Schaaf,et al.  Assessment of the status of the development of standards for the Terrestrial Essential Climate Variables , 2008 .

[55]  R. Dickinson,et al.  Evaluation of the Utility of Satellite-Based Vegetation Leaf Area Index Data for Climate Simulations , 2001 .

[56]  Didier Tanré,et al.  Atmospherically resistant vegetation index (ARVI) for EOS-MODIS , 1992, IEEE Trans. Geosci. Remote. Sens..

[57]  Jihua Wang,et al.  A comparative analysis of spectral vegetation indices to estimate crop leaf area index , 2013, Intell. Autom. Soft Comput..

[58]  J. C. Price,et al.  Leaf area index estimation from visible and near-infrared reflectance data , 1995 .

[59]  J. Chen.  Optically-based methods for measuring seasonal variation of leaf area index in boreal conifer stands , 1996 .

[60]  Flávio Jorge Ponzoni,et al.  Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations , 2012, Remote. Sens..

[61]  Karin S. Fassnacht,et al.  Relationships between leaf area index and Landsat TM spectral vegetation indices across three temperate zone sites , 1999 .

[62]  S. Weisberg,et al.  Diagnostics for heteroscedasticity in regression , 1983 .

[63]  Kenneth A. Sudduth,et al.  Comparing LAI Estimates of Corn and Soybean from Vegetation Indices of Multi-resolution Satellite Images , 2012 .

[64]  R. Casa,et al.  LAI retrieval from multiangular image classification and inversion of a ray tracing model , 2005 .

[65]  A. Viña,et al.  Remote estimation of leaf area index and green leaf biomass in maize canopies , 2003 .

[66]  R. Confalonieri,et al.  Evaluation of LAI-2000 for leaf area index monitoring in paddy rice , 2006 .

[67]  J. Harlan,et al.  Spectral estimation of Green leaf area index of oats , 1985 .

[68]  Hongliang Fang,et al.  Seasonal variation of leaf area index (LAI) over paddy rice fields in NE China: Intercomparison of destructive sampling, LAI-2200, digital hemispherical photography (DHP), and AccuPAR methods , 2014 .

[69]  Richard F. Gunst,et al.  Applied Regression Analysis , 1999, Technometrics.

[70]  Andrew E. Suyker,et al.  Annual carbon dioxide exchange in irrigated and rainfed maize-based agroecosystems , 2005 .

[71]  D. Watson.  Comparative Physiological Studies on the Growth of Field Crops: I. Variation in Net Assimilation Rate and Leaf Area between Species and Varieties, and within and between Years , 1947 .

[72]  Ying Gao,et al.  The Soil Moisture Active Passive Experiments (SMAPEx): Toward Soil Moisture Retrieval From the SMAP Mission , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[73]  Jan G. P. W. Clevers,et al.  Experimental Sentinel-2 LAI estimation using parametric, non-parametric and physical retrieval methods - A comparison , 2015 .

[74]  Clement Atzberger,et al.  Estimation of Leaf Area Index Using DEIMOS-1 Data: Application and Transferability of a Semi-Empirical Relationship between two Agricultural Areas , 2013, Remote. Sens..

[75]  H. Dolman,et al.  Terrestrial essential climate variables for climate change assessment, mitigation and adaptation (GTOS 52) , 2008 .

[76]  C. Bacour,et al.  Comparison of four radiative transfer models to simulate plant canopies reflectance: direct and inverse mode. , 2000 .

[77]  J. Nagy,et al.  Growth and yield of cotton in response to a free-air carbon dioxide enrichment (FACE) environment , 1994 .

[78]  Keith F. Lewin,et al.  Growth and yield of cotton exposed to free‐air CO2 enrichment (Face) , 1992 .

[79]  E. Pattey,et al.  Assessment of vegetation indices for regional crop green LAI estimation from Landsat images over multiple growing seasons , 2012 .

[80]  A. Viña,et al.  Green leaf area index estimation in maize and soybean: Combining vegetation indices to achieve maximal sensitivity , 2012 .

[81]  Yuri Knyazikhin,et al.  Retrieval of canopy biophysical variables from bidirectional reflectance Using prior information to solve the ill-posed inverse problem , 2003 .

[82]  C. Tucker,et al.  Climate-Driven Increases in Global Terrestrial Net Primary Production from 1982 to 1999 , 2003, Science.

[83]  G. Hoogenboom,et al.  Integration of MODIS LAI and vegetation index products with the CSM–CERES–Maize model for corn yield estimation , 2011 .

[84]  Luiz E. O. C. Aragão,et al.  Can MODIS EVI monitor ecosystem productivity in the Amazon rainforest? , 2014 .

[85]  C. Jordan.  Derivation of leaf-area index from quality of light on the forest floor , 1969 .

[86]  S. Running,et al.  Synergistic algorithm for estimating vegetation canopy leaf area index and fraction of absorbed photosynthetically active , 1998 .

[87]  A. Huete.  A soil-adjusted vegetation index (SAVI) , 1988 .

[88]  N. H. Brogea,et al.  Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density , 2022 .

[89]  James W. Jones,et al.  The DSSAT cropping system model , 2003 .

[90]  Ranga B. Myneni,et al.  Estimation of global leaf area index and absorbed par using radiative transfer models , 1997, IEEE Trans. Geosci. Remote. Sens..

[91]  John R. Miller,et al.  Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture , 2004 .

[92]  F. Baret,et al.  GEOV1: LAI, FAPAR essential climate variables and FCOVER global time series capitalizing over existing products. Part 2: Validation and intercomparison with reference products , 2013 .

[93]  D. Watson,et al.  Comparative physiological studies on the growth of field crops. III. The effect of infection with beet yellows and beet mosaic viruses on the growth and yield of the sugar-beet root crop , 1953 .

[94]  D. Roy,et al.  Continental-scale Validation of MODIS-based and LEDAPS Landsat ETM+ Atmospheric Correction Methods , 2012 .

[95]  Alessandro Anav,et al.  Evaluation of Land Surface Models in Reproducing Satellite Derived Leaf Area Index over the High-Latitude Northern Hemisphere. Part II: Earth System Models , 2013, Remote. Sens..

[96]  Alfonso Fernández-Manso,et al.  Spectral unmixing , 2012 .

[97]  V. Capelozzi,et al.  The minimum sum of absolute errors regression: a robust alternative to the least squares regression. , 1999, Statistics in medicine.

[98]  A. Viña,et al.  Remote estimation of canopy chlorophyll content in crops , 2005 .

[99]  Yann Kerr,et al.  Leaf area index estimates using remotely sensed data and BRDF models in a semiarid region. , 2000 .

[100]  James Hansen,et al.  Assimilation of remotely sensed soil moisture and vegetation with a crop simulation model for maize yield prediction , 2013 .

[101]  Michael H. Young,et al.  Monitoring Vegetation Phenological Cycles in Two Different Semi-Arid Environmental Settings Using a Ground-Based NDVI System: A Potential Approach to Improve Satellite Data Interpretation , 2010, Remote. Sens..

[102]  Michael E. Schaepman,et al.  A review on reflective remote sensing and data assimilation techniques for enhanced agroecosystem modeling , 2007, Int. J. Appl. Earth Obs. Geoinformation.

[103]  Zhiqiang Xiao,et al.  Reprocessing the MODIS Leaf Area Index products for land surface and climate modelling , 2011 .

[104]  Frédéric Baret,et al.  Review of methods for in situ leaf area index determination Part I. Theories, sensors and hemispherical photography , 2004 .

[105]  Subhash C. NarulaI,et al.  The Minimum Sum of Absolute Errors Regression: A State of the Art Survey , 1982 .

[106]  Y. Knyazikhin,et al.  Validation and intercomparison of global Leaf Area Index products derived from remote sensing data , 2008 .

[107]  Clement Atzberger,et al.  Comparative analysis of different retrieval methods for mapping grassland leaf area index using airborne imaging spectroscopy , 2015, Int. J. Appl. Earth Obs. Geoinformation.

[108]  F. Baret,et al.  Review of methods for in situ leaf area index (LAI) determination: Part II. Estimation of LAI, errors and sampling , 2004 .

[109]  R. Colombo,et al.  Retrieval of leaf area index in different vegetation types using high resolution satellite data , 2003 .

[110]  T. O. Kvålseth.  Cautionary Note about R 2 , 1985 .

[111]  W. Cohen,et al.  An improved strategy for regression of biophysical variables and Landsat ETM+ data. , 2003 .

[112]  Nadine Gobron,et al.  Optical remote sensing of vegetation: Modeling, caveats, and algorithms , 1995 .

[113]  Wolfram Mauser,et al.  Optimal Exploitation of the Sentinel-2 Spectral Capabilities for Crop Leaf Area Index Mapping , 2012, Remote. Sens..

[114]  A. Gitelson,et al.  Estimating green LAI in four crops: Potential of determining optimal spectral bands for a universal algorithm , 2014 .

[115]  Jingfeng Huang,et al.  Evaluation of MODIS surface reflectance products for wheat leaf area index (LAI) retrieval , 2008 .

[116]  P. Sen.  Estimates of the Regression Coefficient Based on Kendall's Tau , 1968 .

[117]  W. Verhoef.  Light scattering by leaf layers with application to canopy reflectance modelling: The SAIL model , 1984 .

[118]  S. Running,et al.  Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data , 2002 .

[119]  G. Foody,et al.  Sub-pixel land cover composition estimation using a linear mixture model and fuzzy membership functions , 1994 .

[120]  D. Diner,et al.  Estimation of vegetation canopy leaf area index and fraction of absorbed photosynthetically active radiation from atmosphere‐corrected MISR data , 1998 .

[121]  M. Weissa,et al.  Review of methods for in situ leaf area index ( LAI ) determination Part II . Estimation of LAI , errors and sampling , 2003 .

[122]  Nadine Gobron,et al.  On Line Validation Exercise (OLIVE): A Web Based Service for the Validation of Medium Resolution Land Products. Application to FAPAR Products , 2014, Remote. Sens..

[123]  O. Hagolle,et al.  LAI, fAPAR and fCover CYCLOPES global products derived from VEGETATION: Part 1: Principles of the algorithm , 2007 .

[124]  Edward J. Kim,et al.  The NAFE'06 data set: towards soil moisture retrieval at intermediate resolution , 2008 .

[125]  Samuel C. Zipper,et al.  Using evapotranspiration to assess drought sensitivity on a subfield scale with HRMET, a high resolution surface energy balance model , 2014 .

[126]  H. Theil.  A Rank-Invariant Method of Linear and Polynomial Regression Analysis , 1992 .

[127]  C. Justice,et al.  Atmospheric correction of visible to middle-infrared EOS-MODIS data over land surfaces: Background, operational algorithm and validation , 1997 .

[128]  Robert E. Wolfe,et al.  A Landsat surface reflectance dataset for North America, 1990-2000 , 2006, IEEE Geoscience and Remote Sensing Letters.