Assimilating leaf area index of three typical types of subtropical forest in China from MODIS time series data based on the integrated ensemble Kalman filter and PROSAIL model

Abstract Subtropical forest ecosystems play essential roles in the global carbon cycle and in carbon sequestration functions, which challenge the traditional understanding of the main functional areas of carbon sequestration in the temperate forests of Europe and America. The leaf area index (LAI) is an important biological parameter in the spatiotemporal simulation of the carbon cycle, and it has considerable significance in carbon cycle research. Dynamic retrieval based on remote sensing data is an important method with which to obtain large-scale high-accuracy assessments of LAI. This study developed an algorithm for assimilating LAI dynamics based on an integrated ensemble Kalman filter using MODIS LAI data, MODIS reflectance data, and canopy reflectance data modeled by PROSAIL, for three typical types of subtropical forest (Moso bamboo forest, Lei bamboo forest, and evergreen and deciduous broadleaf forest) in China during 2014–2015. There were some errors of assimilation in winter, because of the bad data quality of the MODIS product. Overall, the assimilated LAI well matched the observed LAI, with R2 of 0.82, 0.93, and 0.87, RMSE of 0.73, 0.49, and 0.42, and aBIAS of 0.50, 0.23, and 0.03 for Moso bamboo forest, Lei bamboo forest, and evergreen and deciduous broadleaf forest, respectively. The algorithm greatly decreased the uncertainty of the MODIS LAI in the growing season and it improved the accuracy of the MODIS LAI. The advantage of the algorithm is its use of biophysical parameters (e.g., measured LAI) in the LAI assimilation, which makes it possible to assimilate long-term MODIS LAI time series data, and to provide high-accuracy LAI data for the study of carbon cycle characteristics in subtropical forest ecosystems.

[1]  W. Verhoef,et al.  Coupled soil–leaf-canopy and atmosphere radiative transfer modeling to simulate hyperspectral multi-angular surface reflectance and TOA radiance data , 2007 .

[2]  W. Verhoef,et al.  PROSPECT+SAIL models: A review of use for vegetation characterization , 2009 .

[3]  S. Jacquemoud Inversion of the PROSPECT + SAIL Canopy Reflectance Model from AVIRIS Equivalent Spectra: Theoretical Study , 1993 .

[4]  Improvement of MODIS LAI Product in China , 2008 .

[5]  F. Baret,et al.  PROSPECT: A model of leaf optical properties spectra , 1990 .

[6]  R. B. Jackson,et al.  A Large and Persistent Carbon Sink in the World’s Forests , 2011, Science.

[7]  Liu Jiyuan,et al.  Multi-scale observation and cross-scale mechanistic modeling on terrestrial ecosystem carbon cycle , 2005 .

[8]  Chunlin Huang,et al.  A simplified data assimilation method for reconstructing time-series MODIS NDVI data , 2009 .

[9]  Guo-mo Zhou,et al.  [Retrieval of leaf area index of moso bamboo forest with Landsat Thematic Mapper image based on PROSAIL canopy radiative transfer model]. , 2013, Ying yong sheng tai xue bao = The journal of applied ecology.

[10]  F. Woodward,et al.  Terrestrial Gross Carbon Dioxide Uptake: Global Distribution and Covariation with Climate , 2010, Science.

[11]  Lufeng Mo,et al.  Current and potential carbon stocks in Moso bamboo forests in China. , 2015, Journal of environmental management.

[12]  Yanxia Zhao,et al.  Assimilating remote sensing information with crop model using Ensemble Kalman Filter for improving LAI monitoring and yield estimation , 2013 .

[13]  Roberta E. Martin,et al.  PROSPECT-4 and 5: Advances in the leaf optical properties model separating photosynthetic pigments , 2008 .

[14]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

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

[16]  Rasmus Fensholt,et al.  MODIS leaf area index products: from validation to algorithm improvement , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[17]  Guoyi Zhou,et al.  Estimates of soil respiration and net primary production of three forests at different succession stages in South China , 2006 .

[18]  Jinling Song,et al.  Improvement of spatially continuous forest LAI retrieval by integration of discrete airborne LiDAR and remote sensing multi-angle optical data , 2014 .

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

[20]  A. Kuusk,et al.  A reflectance model for the homogeneous plant canopy and its inversion , 1989 .

[21]  Guo-shuai Zhao,et al.  [Vegetation net primary productivity in Northeast China in 2000-2008: simulation and seasonal change]. , 2011, Ying yong sheng tai xue bao = The journal of applied ecology.

[22]  J. Chen,et al.  Retrieving Leaf Area Index of Boreal Conifer Forests Using Landsat TM Images , 1996 .

[23]  S. Leblanc,et al.  Derivation and validation of Canada-wide coarse-resolution leaf area index maps using high-resolution satellite imagery and ground measurements , 2002 .

[24]  J. Mahfouf,et al.  Optimal interpolation analysis of leaf area index using MODIS data , 2006 .

[25]  Zhou Guo-mo Dynamic change of Phyllostachys edulis forest canopy parameters and their relationships with photosynthetic active radiation in the bamboo shooting growth phase , 2012 .

[26]  Sandra A. Brown Measuring carbon in forests: current status and future challenges. , 2002, Environmental pollution.

[27]  Jing Chen,et al.  Estimating seasonal variations of leaf area index using litterfall collection and optical methods in four mixed evergreen–deciduous forests , 2015 .

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

[29]  T. Meng Intraspecific and Interspecific competition analysis of community dominant plant populations based on Voronoi diagram , 2007 .

[30]  Philip Lewis,et al.  Assimilating canopy reflectance data into an ecosystem model with an Ensemble Kalman Filter , 2008 .

[31]  G. Evensen Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[32]  Philippe Ciais,et al.  The carbon balance of terrestrial ecosystems in China , 2009, Nature.

[33]  Weimin Ju,et al.  Spatial and temporal variations of forest LAI in China during 2000–2010 , 2012 .

[34]  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.

[35]  Xiaojun Xu,et al.  Spatiotemporal heterogeneity of Moso bamboo aboveground carbon storage with Landsat Thematic Mapper images: a case study from Anji County, China , 2013 .

[36]  Soroosh Sorooshian,et al.  Dual state-parameter estimation of hydrological models using ensemble Kalman filter , 2005 .

[37]  Dennis McLaughlin,et al.  An integrated approach to hydrologic data assimilation: interpolation, smoothing, and filtering , 2002 .

[38]  Guirui Yu,et al.  An old-growth subtropical Asian evergreen forest as a large carbon sink , 2011 .

[39]  Shuguang Liu,et al.  Old-Growth Forests Can Accumulate Carbon in Soils , 2006, Science.

[40]  S. Plummer,et al.  Perspectives on combining ecological process models and remotely sensed data , 2000 .