Estimating aboveground and organ biomass of plant canopies across the entire season of rice growth with terrestrial laser scanning

Abstract Non-destructive and accurate estimation of crop biomass is crucial for the quantitative diagnosis of growth status and timely prediction of grain yield. As an active remote sensing technique, terrestrial laser scanning (TLS) has become increasingly available in crop monitoring for its advantages in recording structural properties. Some researchers have attempted to use TLS data in the estimation of crop aboveground biomass, but only for part of the growing season. Previous studies rarely investigated the estimation of biomass for individual organs, such as the panicles in rice canopies, which led to the poor understanding of TLS technology in monitoring biomass partitioning among organs. The objective of this study was to investigate the potential of TLS in estimating the biomass for individual organs and aboveground biomass of rice and to examine the feasibility of developing universal models for the entire growing season. The field plots experiments were conducted in 2017 and 2018 and involved different nitrogen (N) rates, planting techniques and rice varieties. Three regression approaches, stepwise multiple linear regression (SMLR), random forest regression (RF) and linear mixed-effects (LME) modeling, were evaluated in estimating biomass with extensive TLS and biomass data collected at multiple phenological stages of rice growth across the entire season. The models were calibrated with the 2017 dataset and validated independently with the 2018 dataset. The results demonstrated that growth stage in LME modeling was selected as the most significant random effect on rice growth among the three candidates, which were rice variety, growth stage and planting technique. The LME models grouped by growth stage exhibited higher validation accuracies for all biomass variables over the entire season to varying degrees than SMLR models and RF models. The most pronounced improvement with a LME model was obtained for panicle biomass, with an increase of 0.74 in R2 (LME: R2 = 0.90, SMLR: R2 = 0.16) and a decrease of 1.15 t/ha in RMSE (LME: RMSE =0.79 t/ha, SMLR: RMSE =2.94 t/ha). Compared to SMLR and RF, LME modeling yielded similar estimation accuracies of aboveground biomass for pre-heading stages, but significantly higher accuracies for post-heading stages (LME: R2 = 0.63, RMSE =2.27 t/ha; SMLR: R2 = 0.42, RMSE =2.42 t/ha; RF: R2 = 0.57, RMSE =2.80 t/ha). These findings implied that SMLR was only suitable for the estimation of biomass at pre-heading stages and LME modeling performed remarkably well across all growth stages, especially for post-heading. The results suggest coupling TLS with LME modeling is a promising approach to monitoring rice biomass at post-heading stages at high accuracy and to overcoming the saturation of canopy reflectance signals encountered in optical remote sensing. It also has great potential in the monitoring of other crops in cloud-cover conditions and the instantaneous prediction of grain yield any time before harvest.

[1]  Sébastien Bauwens,et al.  Forest Inventory with Terrestrial LiDAR: A Comparison of Static and Hand-Held Mobile Laser Scanning , 2016 .

[2]  K. Keightley,et al.  Original paper: 3D volumetric modeling of grapevine biomass using Tripod LiDAR , 2010 .

[3]  Lloyd T. Wilson,et al.  Degree-days: An aid in crop and pest management , 1983 .

[4]  Qiang Cao,et al.  Multitemporal crop surface models: accurate plant height measurement and biomass estimation with terrestrial laser scanning in paddy rice , 2014 .

[5]  Nguyen Hieu Trung,et al.  A comparative analysis of multitemporal MODIS EVI and NDVI data for large-scale rice yield estimation , 2014 .

[6]  C. L. Todoroki,et al.  Linear mixed-effects models for estimating biomass and fuel loads in shrublands , 2010 .

[7]  S. Popescu,et al.  Satellite lidar vs. small footprint airborne lidar: Comparing the accuracy of aboveground biomass estimates and forest structure metrics at footprint level , 2011 .

[8]  M. Maziah,et al.  Understanding the shoot apical meristem regulation: A study of the phytohormones, auxin and cytokinin, in rice , 2015, Mechanisms of Development.

[9]  T.C.E. Cheng,et al.  Detection of wheat height using optimized multi-scan mode of LiDAR during the entire growth stages , 2019, Comput. Electron. Agric..

[10]  Hamid Hamraz,et al.  Forest understory trees can be segmented accurately within sufficiently dense airborne laser scanning point clouds , 2017, Scientific Reports.

[11]  Thuy Le Toan,et al.  Dependence of radar backscatter on coniferous forest biomass , 1992, IEEE Trans. Geosci. Remote. Sens..

[12]  Nicanor Z. Saliendra,et al.  Bowen ratio and closed chamber carbon dioxide flux measurements over sagebrush steppe vegetation , 2001 .

[13]  C. Bennington,et al.  Use and Misuse of Mixed Model Analysis of Variance in Ecological Studies , 1994 .

[14]  Qi Chen Retrieving vegetation height of forests and woodlands over mountainous areas in the Pacific Coast region using satellite laser altimetry , 2010 .

[15]  G. Bareth,et al.  Multi-temporal Crop Surface Models derived from terrestrial laser scanning for accurate plant height measurement and biomass estimation of barley , 2014 .

[16]  Colin Bradley,et al.  Reverse Engineering Physical Models Employing Wrap-Around B-Spline Surfaces and Quadrics , 1996 .

[17]  W. K. Lauenroth,et al.  Small rainfall events: An ecological role in semiarid regions , 1982, Oecologia.

[18]  Stefan Sommer,et al.  Addressing the complexity in non-linear evolution of vegetation phenological change with time-series of remote sensing images , 2013 .

[19]  In-Sue Oh,et al.  Fixed- versus random-effects models in meta-analysis: model properties and an empirical comparison of differences in results. , 2009, The British journal of mathematical and statistical psychology.

[20]  Horst Bischof,et al.  On-line Random Forests , 2009, 2009 IEEE 12th International Conference on Computer Vision Workshops, ICCV Workshops.

[21]  Weixing Cao,et al.  Development of critical nitrogen dilution curve of Japonica rice in Yangtze River Reaches , 2013 .

[22]  Y. Miao,et al.  Precise plant height monitoring and biomass estimation with Terrestrial Laser Scanning in paddy rice , 2013 .

[23]  A. K. Milne,et al.  The potential of synthetic aperture radar (SAR) for quantifying the biomass of Australia's woodlands. , 2000 .

[24]  T. Reeves,et al.  The Cereal of the World's Poor Takes Center Stage , 2002, Science.

[25]  J. Eitel,et al.  Estimating aboveground biomass and leaf area of low-stature Arctic shrubs with terrestrial LiDAR , 2015 .

[26]  J. Schjoerring,et al.  Reflectance measurement of canopy biomass and nitrogen status in wheat crops using normalized difference vegetation indices and partial least squares regression , 2003 .

[27]  Ryotaro Kamimura,et al.  A Stepwise AIC Method for Variable Selection in Linear Regression , 2007 .

[28]  Nora Tilly,et al.  Fusion of Plant Height and Vegetation Indices for the Estimation of Barley Biomass , 2015, Remote. Sens..

[29]  Dong Li,et al.  Spectroscopic Estimation of Biomass in Canopy Components of Paddy Rice Using Dry Matter and Chlorophyll Indices , 2017, Remote. Sens..

[30]  Toshihiro Sakamoto,et al.  Refined shape model fitting methods for detecting various types of phenological information on major U.S. crops , 2018 .

[31]  K. Omasa,et al.  Estimation of vertical plant area density profiles in a rice canopy at different growth stages by high-resolution portable scanning lidar with a lightweight mirror , 2012 .

[32]  C. Kao,et al.  Senescence of Rice Leaves : VII. PROLINE ACCUMULATION IN SENESCING EXCISED LEAVES. , 1982, Plant physiology.

[33]  V. Carey,et al.  Mixed-Effects Models in S and S-Plus , 2001 .

[34]  M. Jeuffroy,et al.  Diagnosis tool for plant and crop N status in vegetative stage Theory and practices for crop N management , 2008 .

[35]  Marguerite Madden,et al.  A linear mixed-effects model of biomass and volume of trees using Landsat ETM+ images , 2007 .

[36]  D. Long,et al.  Mapping straw yield using on-combine light detection and ranging (lidar) , 2013 .

[37]  Hongzhi An,et al.  Fast stepwise procedures of selection of variables by using AIC and BIC criteria , 1989 .

[38]  Guijun Yang,et al.  Winter wheat biomass estimation based on spectral indices, band depth analysis and partial least squares regression using hyperspectral measurements , 2014 .

[39]  J. Parfitt,et al.  Food waste within food supply chains: quantification and potential for change to 2050 , 2010, Philosophical Transactions of the Royal Society B: Biological Sciences.

[40]  W D Hutchison,et al.  Comparison of economic and plant phenology-based thresholds for management of Lygus lineolaris (Hemiptera: Miridae) in Minnesota strawberries. , 2003, Journal of economic entomology.

[41]  M. Cheung A model for integrating fixed-, random-, and mixed-effects meta-analyses into structural equation modeling. , 2008, Psychological methods.

[42]  Changying Li,et al.  In-field High Throughput Phenotyping and Cotton Plant Growth Analysis Using LiDAR , 2018, Front. Plant Sci..

[43]  Hui Feng,et al.  High-throughput volumetric reconstruction for 3D wheat plant architecture studies , 2016 .

[44]  Ning Zhang,et al.  Retrieval of three-dimensional tree canopy and shade using terrestrial laser scanning (TLS) data to analyze the cooling effect of vegetation , 2016 .

[45]  N. West,et al.  Measurements of Terrestrial Vegetation , 1989 .

[46]  De-Shuang Huang,et al.  Lidar signal denoising using least-squares support vector machine , 2005, IEEE Signal Processing Letters.

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

[48]  D. Ehlert,et al.  Laser rangefinder-based measuring of crop biomass under field conditions , 2009, Precision Agriculture.

[49]  Qi Chen,et al.  Modeling aboveground tree woody biomass using national-scale allometric methods and airborne lidar , 2015 .

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

[51]  Lalit Kumar,et al.  Estimation of Winter Wheat Biomass and Yield by Combining the AquaCrop Model and Field Hyperspectral Data , 2016, Remote. Sens..

[52]  Hengbiao Zheng,et al.  Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery , 2018, Precision Agriculture.

[53]  J. Eitel,et al.  LiDAR based biomass and crop nitrogen estimates for rapid, non-destructive assessment of wheat nitrogen status , 2014 .

[54]  Yan Li,et al.  Towards an operational system for regional-scale rice yield estimation using a time-series of Radarsat ScanSAR images , 2003 .

[55]  P. Thenkabail,et al.  Hyperspectral Vegetation Indices and Their Relationships with Agricultural Crop Characteristics , 2000 .

[56]  David Saah,et al.  Integration of airborne lidar and vegetation types derived from aerial photography for mapping aboveground live biomass , 2012 .

[57]  Ke Zhang,et al.  Estimation of Rice Growth Parameters Based on Linear Mixed-Effect Model Using Multispectral Images from Fixed-Wing Unmanned Aerial Vehicles , 2019, Remote. Sens..

[58]  R. B. Waite,et al.  The application of visual estimation procedures for monitoring pasture yield and composition in exclosures and small plots , 1994 .

[59]  G. Bareth,et al.  Terrestrial laser scanning for plant height measurement and biomass estimation of maize , 2014 .

[60]  K. Omasa,et al.  Estimating vertical plant area density profile and growth parameters of a wheat canopy at different growth stages using three-dimensional portable lidar imaging , 2009 .

[61]  Jose A. Jiménez-Berni,et al.  High Throughput Determination of Plant Height, Ground Cover, and Above-Ground Biomass in Wheat with LiDAR , 2018, Front. Plant Sci..

[62]  Jose A. Jiménez-Berni,et al.  Proximal Remote Sensing Buggies and Potential Applications for Field-Based Phenotyping , 2014 .

[63]  Yutaka Sato,et al.  Formation, Maintenance and Function of the Shoot Apical Meristem in Rice , 2006, Plant Molecular Biology.

[64]  Teemu Hakala,et al.  Uncertainty in multispectral lidar signals caused by incidence angle effects , 2018, Interface Focus.

[65]  Jihua Wang,et al.  Assimilation of Two Variables Derived from Hyperspectral Data into the DSSAT-CERES Model for Grain Yield and Quality Estimation , 2015, Remote. Sens..

[66]  Tee-Ann Teo,et al.  Radiometric Block Adjustment for Multi-Strip Airborne Waveform Lidar Data , 2015, Remote. Sens..

[67]  Romeo M. Visperas,et al.  Grain and dry matter yields and partitioning of assimilates in japonica/indica hybrid rice , 2002 .

[68]  Xin-jie Wang,et al.  [Simulation of the branch biomass for Chinese fir plantation using the linear mixed effects model]. , 2015, Ying yong sheng tai xue bao = The journal of applied ecology.

[69]  J. Eitel,et al.  An automated method to quantify crop height and calibrate satellite-derived biomass using hypertemporal lidar , 2016 .