Estimation of crop gross primary production (GPP): I. impact of MODIS observation footprint and impact of vegetation BRDF characteristics

a b s t r a c t Accurate estimation of gross primary production (GPP) is essential for carbon cycle and climate change studies. Three AmeriFlux crop sites of maize and soybean were selected for this study. Two of the sites were irrigated and the other one was rainfed. The normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the green band chlorophyll index (CIgreen), and the green band wide dynamic range vegetation index (WDRVIgreen) were computed from the moderate resolution imaging spectroradiometer (MODIS) surface reflectance data. We examined the impacts of the MODIS observation footprint and the vegetation bidirectional reflectance distribution function (BRDF) on crop daily GPP esti- mation with the four spectral vegetation indices (VIs - NDVI, EVI, WDRVIgreen and CIgreen) where GPP was predicted with two linear models, with and without offset: GPP = a × VI × PAR and GPP = a × VI × PAR + b. Model performance was evaluated with coefficient of determination (R 2 ), root mean square error (RMSE), and coefficient of variation (CV). The MODIS data were filtered into four categories and four experiments were conducted to assess the impacts. The first experiment included all observations. The second experi- ment only included observations with view zenith angle (VZA) ≤ 35 ◦ to constrain growth of the footprint size,which achieved a better grid cell match with the agricultural fields. The third experiment included only forward scatter observations with VZA ≤ 35 ◦ . The fourth experiment included only backscatter obser- vations with VZA ≤ 35◦. Overall, the EVI yielded the most consistently strong relationships to daily GPP under all examined conditions. The model GPP = a × VI × PAR + b had better performance than the model GPP = a × VI × PAR, and the offset was significant for most cases. Better performance was obtained for the irrigated field than its counterpart rainfed field. Comparison of experiment 2 vs. experiment 1 was used to examine the observation footprint impact whereas comparison of experiment 4 vs. experiment 3 was used to examine the BRDF impact. Changes in R2, RMSE,CV and changes in model coefficients "a" and "b" (experiment 2 vs. experiment 1; and experiment 4 vs. experiment 3) were indicators of the impacts. The second experiment produced better performance than the first experiment, increasing R2 (↑0.13) and

[1]  C. Tucker,et al.  Multi-angle implementation of atmospheric correction for MODIS (MAIAC): 3. Atmospheric correction , 2012 .

[2]  Elizabeth M. Middleton,et al.  Variations of Foliage Chlorophyll fAPAR and Foliage Non-Chlorophyll fAPAR (fAPAR$_{\rm chl}$, fAPAR$_{\rm non\mathchar "702D chl}$) at the Harvard Forest , 2013, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[3]  Alan H. Strahler,et al.  The Moderate Resolution Imaging Spectroradiometer (MODIS): land remote sensing for global change research , 1998, IEEE Trans. Geosci. Remote. Sens..

[4]  Clemens Beckstein,et al.  Characterization of ecosystem responses to climatic controls using artificial neural networks , 2010 .

[5]  Christopher B. Field,et al.  Substrate limitations for heterotrophs: Implications for models that estimate the seasonal cycle of atmospheric CO2 , 1996 .

[6]  J. Randerson,et al.  Global net primary production: Combining ecology and remote sensing , 1995 .

[7]  T. Arkebauer,et al.  Gross primary production and ecosystem respiration of irrigated maize and irrigated soybean during a growing season , 2005 .

[8]  Alfredo Huete,et al.  Terrestrial carbon sinks in the Brazilian Amazon and Cerrado region predicted from MODIS satellite data and ecosystem modeling , 2009 .

[9]  R. Fensholt,et al.  Evaluation of MODIS LAI, fAPAR and the relation between fAPAR and NDVI in a semi-arid environment using in situ measurements , 2004 .

[10]  S. Goward,et al.  Evaluation of the NOAA/NASA Pathfinder AVHRR Land Data Set for global primary production , 1995 .

[11]  W. Oechel,et al.  On the use of MODIS EVI to assess gross primary productivity of North American ecosystems , 2006 .

[12]  Bo-Cai Gao,et al.  Using EO-1 Hyperion to Simulate HyspIRI Products for a Coniferous Forest: The Fraction of PAR Absorbed by Chlorophyll $(\hbox{fAPAR}_{\rm chl})$ and Leaf Water Content (LWC) , 2012, IEEE Transactions on Geoscience and Remote Sensing.

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

[14]  A-Xing Zhu,et al.  Developing a continental-scale measure of gross primary production by combining MODIS and AmeriFlux data through Support Vector Machine approach , 2007 .

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

[16]  David P. Roy,et al.  MODIS land data storage, gridding, and compositing methodology: Level 2 grid , 1998, IEEE Trans. Geosci. Remote. Sens..

[17]  Maosheng Zhao,et al.  A Continuous Satellite-Derived Measure of Global Terrestrial Primary Production , 2004 .

[18]  Thomas Hilker,et al.  Dynamics of spectral bio-indicators and their correlations with light use efficiency using directional observations at a Douglas-fir forest , 2009 .

[19]  Jing M. Chen,et al.  Predicting gross primary production from the enhanced vegetation index and photosynthetically active radiation: Evaluation and calibration , 2011 .

[20]  S. Wofsy,et al.  Midday values of gross CO2 flux and light use efficiency during satellite overpasses can be used to directly estimate eight-day mean flux , 2005 .

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

[22]  Wenjiang Huang,et al.  Remote estimation of gross primary production in wheat using chlorophyll-related vegetation indices , 2009 .

[23]  Anatoly A. Gitelson,et al.  Remote estimation of gross primary productivity in crops using MODIS 250m data , 2013 .

[24]  C. Justice,et al.  A Revised Land Surface Parameterization (SiB2) for Atmospheric GCMS. Part II: The Generation of Global Fields of Terrestrial Biophysical Parameters from Satellite Data , 1996 .

[25]  W. Oechel,et al.  A new model of gross primary productivity for North American ecosystems based solely on the enhanced vegetation index and land surface temperature from MODIS , 2008 .

[26]  A. Viña,et al.  Relationship between gross primary production and chlorophyll content in crops: Implications for the synoptic monitoring of vegetation productivity , 2006 .

[27]  Christopher M. Gough,et al.  Remote sensing of canopy light use efficiency in temperate and boreal forests of North America using MODIS imagery , 2012 .

[28]  S. Goward,et al.  Vegetation canopy PAR absorptance and the normalized difference vegetation index - An assessment using the SAIL model , 1992 .

[29]  Youngwook Kim,et al.  Spectral compatibility of vegetation indices across sensors: band decomposition analysis with Hyperion data , 2010 .

[30]  J. Randerson,et al.  Terrestrial ecosystem production: A process model based on global satellite and surface data , 1993 .

[31]  Nicholas C. Coops,et al.  Improving predictions of forest growth using the 3-PGS model with observations made by remote sensing , 2010 .

[32]  A. Dalcher,et al.  A Simple Biosphere Model (SIB) for Use within General Circulation Models , 1986 .

[33]  Lawrence A. Corp,et al.  Utilizing in situ directional hyperspectral measurements to validate bio-indicator simulations for a corn crop canopy , 2010, Ecol. Informatics.

[34]  Anatoly A. Gitelson,et al.  Application of chlorophyll-related vegetation indices for remote estimation of maize productivity , 2011 .

[35]  Andrew E. Suyker,et al.  Synoptic Monitoring of Gross Primary Productivity of Maize Using Landsat Data , 2008, IEEE Geoscience and Remote Sensing Letters.

[36]  S. Verma,et al.  Modeling gross primary production of irrigated and rain-fed maize using MODIS imagery and CO2 flux tower data , 2011 .

[37]  D. P. Turner,et al.  Scaling net primary production to a MODIS footprint in support of Earth observing system product validation , 2004 .

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

[39]  S. Running,et al.  MODIS Leaf Area Index (LAI) And Fraction Of Photosynthetically Active Radiation Absorbed By Vegetation (FPAR) Product , 1999 .

[40]  D. Randall,et al.  A Revised Land Surface Parameterization (SiB2) for Atmospheric GCMS. Part I: Model Formulation , 1996 .

[41]  Anatoly A. Gitelson,et al.  Remote estimation of gross primary productivity in soybean and maize based on total crop chlorophyll content , 2012 .

[42]  I. Simmonds,et al.  Modeling δ18O in tropical precipitation and the surface ocean for present‐day climate , 2006 .

[43]  Thomas Hilker,et al.  Linking foliage spectral responses to canopy-level ecosystem photosynthetic light-use efficiency at a Douglas-fir forest in Canada , 2009 .

[44]  R. Myneni,et al.  On the relationship between FAPAR and NDVI , 1994 .

[45]  C. Tucker.  Red and photographic infrared linear combinations for monitoring vegetation , 1979 .

[46]  T. A. Black,et al.  Separating physiologically and directionally induced changes in PRI using BRDF models , 2008 .

[47]  Andrew E. Suyker,et al.  Estimating daily gross primary production of maize based only on MODIS WDRVI and shortwave radiation data , 2011 .

[48]  Ramakrishna R. Nemani,et al.  Evaluation of remote sensing based terrestrial productivity from MODIS using regional tower eddy flux network observations , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[49]  Maosheng Zhao,et al.  Drought-Induced Reduction in Global Terrestrial Net Primary Production from 2000 Through 2009 , 2010, Science.

[50]  Nicholas C. Coops,et al.  Characterizing spatial representativeness of flux tower eddy-covariance measurements across the Canadian Carbon Program Network using remote sensing and footprint analysis , 2012 .

[51]  A. Arneth,et al.  Phenology and gross primary production of two dominant savanna woodland ecosystems in Southern Africa , 2013 .

[52]  Yujie Wang,et al.  An automatic cloud mask algorithm based on time series of MODIS measurements , 2008 .

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

[54]  Yujie Wang,et al.  Multiangle implementation of atmospheric correction (MAIAC): 1. Radiative transfer basis and look-up tables , 2011 .

[55]  A. Bondeau,et al.  Comparing global models of terrestrial net primary productivity (NPP): analysis of differences in light absorption and light‐use efficiency , 1999 .

[56]  C. Justice,et al.  A revised land surface parameterization (SiB2) for GCMs. Part III: The greening of the Colorado State University general circulation model , 1996 .

[57]  Thomas Hilker,et al.  Tracking plant physiological properties from multi-angular tower-based remote sensing , 2011, Oecologia.

[58]  W. Oechel,et al.  A continuous measure of gross primary production for the conterminous United States derived from MODIS and AmeriFlux data , 2010, Remote Sensing of Environment.

[59]  S. Goward,et al.  Vegetation canopy PAR absorptance and NDVI: An assessment for ten tree species with the SAIL model , 1997 .

[60]  E. Vermote,et al.  Seasonal and inter-annual variation in view angle effects on MODIS vegetation indices at three forest sites , 2011 .

[61]  N. Coops,et al.  Simulating gross primary production across a chronosequence of coastal Douglas-fir forest stands with a production efficiency model. , 2010 .

[62]  W. Cohen,et al.  Scaling Gross Primary Production (GPP) over boreal and deciduous forest landscapes in support of MODIS GPP product validation , 2003 .

[63]  E. Davidson,et al.  Satellite-based modeling of gross primary production in an evergreen needleleaf forest , 2004 .

[64]  D. Roberts,et al.  On intra-annual EVI variability in the dry season of tropical forest A case study with MODIS and hyperspectral data , 2011 .

[65]  A. Gitelson,et al.  Remote estimation of crop gross primary production with Landsat data , 2012 .

[66]  Alfredo Huete,et al.  An empirical investigation of cross-sensor relationships of NDVI and red/near-infrared reflectance using EO-1 Hyperion data , 2006 .

[67]  Andrew E. Suyker,et al.  REMOTE ESTIMATION OF GROSS PRIMARY PRODUCTION IN MAIZE , 2011 .

[68]  A. Huete,et al.  Overview of the radiometric and biophysical performance of the MODIS vegetation indices , 2002 .

[69]  Markus Reichstein,et al.  Improving canopy processes in the Community Land Model version 4 (CLM4) using global flux fields empirically inferred from FLUXNET data , 2011 .

[70]  Yujie Wang,et al.  Multiangle implementation of atmospheric correction (MAIAC): 2. Aerosol algorithm , 2011 .