- Original Paper
- Published:
Remote estimation of canopy leaf area index and chlorophyll content in Moso bamboo (Phyllostachys edulis (Carrière) J. Houz.) forest using MODIS reflectance data
Annals of Forest Science volume 75, Article number: 33 (2018)
Abstract
Key message
We estimated the leaf area index (LAI) and canopy chlorophyll content (CC) of Moso bamboo forest by using statistical models based on MODIS data and field measurements. Results showed that the statistical model driven by MODIS data has the potential to accurately estimate LAI and CC, while the structure of the calibration models varied between on- and off-years because of the different leaf change and bamboo shoot production characteristics between these types of years.
Context
LAI and CC (gram per square meter of ground area) are important parameters for determining carbon exchange between Moso bamboo forest (Phyllostachys edulis (Carrière) J. Houz.) and the atmosphere.
Aims
This study evaluated the ability of a statistical model driven by MODIS data to accurately estimate the LAI and CC in Moso bamboo forest, and differences in the LAI and CC between on-years (years with great shoot production) and off-years (years with less shoot production) were analyzed.
Methods
The LAI and CC measurements were collected in Anji County, Zhejiang Province, China. Indicators of LAI and CC were calculated from MODIS data. Then, a regression analysis was used to build relationships between the LAI and CC and various indicators on the basis of leaf change and bamboo shoot production characteristics of Moso bamboo forest.
Results
LAI and CC were accurately estimated by using the regression analysis driven by MODIS-derived indicators with a relative root mean squared error (RMSEr) of 9.04 and 13.1%, respectively. The structure of the calibration models varied between on- and off-years. Long-term time series analysis from 2000 to 2015 showed that LAI and CC differed largely between on- and off-years.
Conclusion
This study demonstrates that LAI and CC of Moso bamboo forest can be estimated accurately by using a statistical model driven by MODIS-derived indicators, but attention should be paid to differences in the calibration models between on- and off-years.
1 Introduction
Leaf area index (LAI) and canopy chlorophyll content (CC, gram per square meter of ground area) are important indicators of the structural and biochemical characteristics of the canopy (Houborg et al. 2007). LAI, which is defined as half the total leaf area per unit ground surface area, is crucial in determining evapotranspiration and CO2 exchange between the biosphere and the atmosphere (Running et al. 1989; Sellers 1997; Houborg et al. 2007; Maire et al. 2011; Olsoy et al. 2016). Changes in CC are related to vegetation senescence, photosynthetic capacity due to photosynthetically active radiation absorbed by CC, and effects of stress on destruction of CC (Datt 1999). Both parameters are important inputs in ecosystem models designed to quantify the spatial and temporal distributions of carbon, water, and energy fluxes (Gitelson et al. 2006; Wu et al. 2009).
LAI is routinely available from satellite-based Earth observation instruments, such as the Moderate Resolution Imaging Spectroradiometer (MODIS). The MODIS LAI estimates have been used widely and successfully in ecosystem models. However, many previous validation studies on MODIS-derived LAI have suggested that they contain considerable errors (Fang and Liang 2005; Leuning et al. 2005; Cohen et al. 2006; De Kauwe et al. 2011; Fang et al. 2012; Xu et al. 2016). Compared with in situ observations, MODIS LAI was often overestimated for broadleaf forest, needle-leaf forest, savannah, and woody vegetation (Fang and Liang 2005; Heinsch et al. 2006; Fang et al. 2012; Xu et al. 2016), but underestimates for cropland and bamboo forest (Meyers and Hollinger 2004; Li et al. 2014; Xu et al. 2016). For example, the Collection 4 MODIS LAI product for broadleaf forest and evergreen needle-leaf forest overestimates by around 2.0–3.0 and 0.9 m2/m2, respectively (Heinsch et al. 2006). Results from De Kauwe et al. (2011) suggest that the Collection 5 MODIS LAI underestimates the upper range of in situ LAI measurements. The relatively large uncertainties (± 1.0 m2/m2) of the Collection 5 MODIS LAI product mean that it cannot satisfy the accuracy requirements of many systems used to observe global climate (± 0.5 m2/m2) (Fang et al. 2012). The application of the Collection 5 MODIS LAI product to Moso bamboo (Phyllostachys edulis (Carrière) J. Houz.) forest resulted in underestimates with a larger degree of error than that for deciduous broadleaf forest and evergreen needle-leaf forest (Xu et al. 2016). Furthermore, errors in the LAI will propagate into CC results because LAI is a driver of CC.
Alternative and effective methods for the retrieval of LAI and CC data in Moso bamboo forest are necessary for the improvement of the MODIS LAI product. There are two common methods for estimating LAI and CC based on remote sensing techniques, namely, physically based models and empirical–statistical models (Coops et al. 2003; Houborg et al. 2007, 2009; Maire et al. 2011; Zhang et al. 2012). Physically based models can accurately explain the relationships between biophysical variables and canopy reflectance because they explicitly describe the transfer and interaction of radiation inside the canopy and external factors (Jacquemoud et al. 2000; Houborg et al. 2007). The main limitations of physically based models are the expenses associated with the computational requirements and the non-uniqueness of the solution (Jacquemoud et al. 2000; Houborg et al. 2007). Empirical–statistical models are commonly used for the remote estimation of LAI and CC, as they are based on the strong relationships between LAI and CC and remote sensing reflectance data. Their main advantage is the ease with which they can provide estimations of large-scale biophysical variables, i.e., they there are simple and computationally efficient. The main problem with empirical–statistical models is their lack of generality (Houborg et al. 2007). The structure of empirical–statistical models depends on site-, time-, and species-specific factors (Colombo et al. 2003; Houborg et al. 2007) that require representative field data for their calibration.
Bamboos are distributed naturally in tropical, subtropical, and temperate regions of all the continents (Grattani et al. 2008). The entire area of bamboo forest around the world is about 31.5 million ha and accounts for roughly 0.8% of the global forest area (FAO 2010). In China, the bamboo area is more than 6 million ha, 70% of which is covered by Moso bamboo forest (Song et al. 2016). Moso bamboo has a large carbon storage capacity (Chen et al. 2009; Komatsu et al. 2010, 2012; Yen and Lee 2011), and thus it is important to develop accurate estimates of LAI and CC for these forests despite the inaccuracy in the MODIS estimates (Xu et al. 2016). Precise determinations of the LAI and CC of Moso bamboo forest is needed for accurate estimations of the transpiration and carbon fluxes of Moso bamboo forests.
A special characteristic of Moso bamboo forest is its ability to produce many bamboo shoots through an underground rhizome system. Unlike other vegetation types, a newly emerged bamboo shoot can reach its maximum height within 40 days (Yen 2016). According to the number of bamboo shoots, Moso bamboo forest can be classified as being in either an on-year or off-year state, i.e., years with great bamboo shoot production are known as on-years and years with less bamboo shoot production are known as off-years (Li et al. 1998; Zhou et al. 2011). The number of bamboo shoots in an on-year is considerably greater than that in an off-year during mid-March to mid-May. Moso bamboo forest commonly exists as an uneven-aged forest with one to 5-year-old culms in stands subjected to intensive management. Most of the culms in on-years have uneven-year ages (e.g., 1, 3, 5), while those in off-years have even-year ages (e.g., 2, 4). The leaf life span of Moso bamboo is 2 years, except for first-year shoots (Li et al. 1998). Leaves on culms of even-year ages are 1-year-old leaves, whereas those on culms of uneven-year ages are 2-year-old leaves (Kleinhenz et al. 2001). Therefore, most of the leaves in on-years are 2-year-old leaves and those in off-years are 1-year-old leaves. Rates of photosynthesis between 1-year-old leaves and 2-year-old leaves are significantly different (Kleinhenz et al. 2001). Two-year-old leaves fall in spring and are replaced by 1-year-old leaves (Li et al. 1998). Hence, the seasonal dynamics of LAI and CC differ between on- and off-years. Besides, the frequency distribution of culms among the age classes and the aboveground carbon storage were found to be significantly different between Moso bamboo stands that were subjected to two different management strategies, namely, intensive management and extensive management (Yen 2015), which resulted in differences in the LAI and CC between the different stands. Therefore, significantly different leaf change and bamboo shoot production characteristics between on- and off-years and different management strategies may cause challenges in estimating the LAI and CC of Moso bamboo forest based on MODIS reflectance data.
In this study, relationships between four vegetation indices from the Terra MODIS 8-day composite reflectance product (MOD09A1) and in situ LAI and CC observations were analyzed. Then, empirical–statistical models were calibrated with a regression method and validated by a leave-one-out cross-validation method. LAI and CC of Moso bamboo forest in the study site from 2000 to 2015 were predicted and compared with the MODIS LAI product (MOD15A2), and corrected MODIS LAI data were produced by using the “upper-envelope” smoothing method (Gu et al. 2006; Xu et al. 2016). Finally, the differences in LAI and CC of Moso bamboo forest between on- and off-years were analyzed. The objectives were (1) to test the ability of simple empirical–statistical models to estimate LAI and CC of Moso bamboo forest accurately by using MODIS-derived indicators, (2) to discover differences in the structure of the calibration models between on- and off-years, and (3) to evaluate the uncertainties in the MODIS LAI estimates of Moso bamboo forest.
2 Material and methods
2.1 Study area
The study area is located in Anji County in Zhejiang Province, China (Fig. 1). Here, bamboo forest encompasses an area of 757 km2, and it accounts for 56.47% of the forested area; of this, 79.30% is Moso bamboo. The canopy height of the Moso bamboo forest is approximately 11 m. The average annual precipitation is 1100–1900 mm, and the average annual temperature is 12.2–15.6 °C. Results from 55 sample plots showed that the mean diameter at breast height for the Moso bamboo in this study area is 9.3 cm (the range is between 6.9 and 11.3 cm) (Xu et al. 2013). The culm density is approximately 3235 culms per hectare (Xu et al. 2013). The average aboveground biomass of Moso bamboo in this study area was found to be 44.2 Mg/ha based on a point kriging interpolation method (Du et al. 2010). An eddy covariance (EC) tower (30.476° N, 119.673° E) was built in 2010 to measure the carbon fluxes of the Moso bamboo forest. The 1000 m × 1000 m area around the EC tower site comprises mainly Moso bamboo (86.1%) with small proportions of mixed forest (1.9%), cropland (8.6%), and buildings (3.4%). The EC system and its measurements are described in detail by Xu et al. (2013). In the study area, 2011 and 2015 were on-years and 2014 was off-years (Xu et al. 2013).
2.2 Canopy LAI and chlorophyll content
Twenty-five sample plots (Fig. 1) were selected within a 600 m × 600 m square around the EC tower to calculate the canopy LAI and leaf chlorophyll content (LCC). The size of each sample plot was 5 m × 5 m. We measured LCC and LAI once per month from April 2014 to December 2015. Furthermore, LAI data from April 2011 to July 2011 were also available.
A photograph of the canopy was taken at the center of each sample plot with a vegetation canopy analyzing system (WinSCANOPY2009a; Regent Instruments, Ste-Foy, Quebec, Canada). The canopy LAI was then calculated by using an LAI (2000G)-Log CI correction algorithm based on the photograph (Mao et al. 2017), which was taken as the measured LAI at the corresponding sample plot. The WinSCANOPY results processed by using the WinSCANOPY software with corrections for leaf clumping were close to the direct-harvest values for a tropical rain forest (Olivas et al. 2013).
Because most of the leaves on bamboo culms in on-years (except for new bamboo culms) or off-years are of the same age (Li et al. 1998), five leaves were randomly selected from each of three culms per sample plot and each of the 375 leaves was measured three times with a chlorophyll meter (CCM-200 plus; Opti-Sciences, Boston, USA). The average LCC was taken as representative ground truth data for the sample plot. A calibration equation was used to transform the values measured by the CCM-200 to the sum of the chlorophyll a and b per square meter of leaf (mg/cm2) (Eq. (1), Fig. S1, Gu 2013). Then, CC was defined as the total amount of chlorophyll per unit land area (g/m2), which was calculated by using product of the LAI and LCC (Gitelson et al. 2005).
where LCC is the sum of chlorophyll a and b per square meter of leaf (mg/cm2), and LCCCCM-200 is the value measured when using the CCM-200.
In reference to Wu et al. (2008), the canopy LAI and CC collected from the 25 sample plots were averaged to represent ground truth data for the single 500-m-resolution pixel covering the 25 sample plots. Then, the relationships between the mean LAI and CC obtained from the 25 plots and the MODIS-derived indicators from the pixel with central coordinates closest to the location of the EC tower site were analyzed to build the statistical model for estimating LAI and CC. The variability in LAI and CC across sample plots was not taken into account in the statistical model.
2.3 MODIS reflectance data and vegetation indices
Surface reflectance data were acquired from the Terra MODIS 8-day composite reflectance product (MOD09A1) at a 500-m resolution and were downloaded directly from the website of the Oak Ridge National Laboratory Distributed Active Archive Center. The MOD09A1 product includes the following seven bands: band 1 (red, 620–670 nm), band 2 (near infrared, NIR, 841–876 nm), band 3 (blue, 459–479 nm), band 4 (green, 545–565 nm), band 5 (mid-infrared, mid-IR, 1230–1250 nm), band 6 (mid-IR, 1628–1652 nm), and band 7 (mid-IR, 2105–2155 nm). There are uncertainties in MODIS reflectance data as a result of atmospheric effects, clouds, and snow. Outliers of band 3 reflectance above 0.05 were removed. Then, other band reflectance data were defined as bad if the blue band 3 reflectance data were bad.
Based on MOD09A1 reflectance data, four important vegetation indices were calculated. The Normalized Difference Vegetation Index (NDVI) is a good indicator for the estimation of canopy LAI because it combines the NIR band 2 and red band 1 reflectance data (Houborg et al. 2007). A non-linear transformation of NDVI, called the Wide Dynamic Range Vegetation Index (WDRVI), was also tested because WDRVI is linearly related to green LAI data (Gitelson et al. 2007; Guindin-Garcia et al. 2012). The simple ratio (SR) index was calculated by dividing the reflectance of band 2 by the reflectance of band 1, and it is highly sensitive to the LAI and CC (Chen and Cihlar 1996; Zhang et al. 2008). The Gitelson Green Index (GI) was also used because it is sensitive to changes in CC (Gitelson et al. 2005; Wu et al. 2009). The four indices can be calculated as follows:
where RNIR, Rred, and Rgreen are reflectance data in the NIR band 2, red band 1, and green band 4, respectively. α is equal to 0.1 (Guindin-Garcia et al. 2012).
2.4 Statistical analysis
The day of year (DOY) of the ground-measured biophysical parameters was not the same as the DOY of the MODIS 8-day reflectance data. Therefore, to ensure good correspondence between the ground-measured biophysical parameters and the MODIS 8-day reflectance data, the ground-measured biophysical parameters measured closest to the MODIS 8-day reflectance data in terms of DOY were paired and used for the relationship analysis. The relationships between LAI and CC in Moso bamboo forest and MODIS 8-day reflectance data were analyzed by using the correlation analysis method of the SPSS13.0 statistical software. Linear or non-linear regression models were built to estimate the LAI and CC in Moso bamboo forest. The data sets used in this study are presented in Table 1.
A data set from 3 years (2011, 2014, and 2015) was used for developing the regression models for the estimations of LAI and CC. A leave-one-out cross-validation method was used for the validation of the prediction accuracy of the regression models. The goodness of fit of the regression models was measured by calculating the coefficient of determination (R2), root mean squared error (RMSE), and relative RMSE (RMSEr) for the calibration and validation sets. Eight-day synthesized Collection 5 MODIS LAI (MOD15A2) data with a 1-km spatial resolution were compared with the LAI estimates. The time series of the remote sensing data and the biophysical parameters shown in Table 1 were used in this study.
Data availability statement
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
3 Results
3.1 Relationship between the LAI and CC
Relationships between CC and LAI differed between on-years (2011 and 2015) and off-years (Fig. 2). The relationship between LAI and CC in the off-year (R2 = 0.62) was weak relative to that during on-years (R2 = 0.91). CC during DOY 97 to DOY 137 in the off-year was obviously smaller than that in the on-years (Fig. 3a), but differences in the LAI in the same period between the off- and on-years were relatively small (Fig. 3b). Thus, the difference in CC was caused by the different change trends of the LCC between the off- and on-years during the period of DOY 97 to DOY 153 (Fig. 3c). The differences in CC between off- and on-years were also well represented by the visible red band 1 (Fig. 4). The reflectance of visible red band 1 in the off-year during DOY 65 to DOY 129 was greater than the reflectance in the on-years, which indirectly implies that there was a smaller CC in the off-year than that in the on-years. However, the reflectance of visible red band 1 in the off-year decreased earlier during DOY 129 to DOY 153 when CC in the off-year was still smaller than that in the on-years. Therefore, separation of sample points into on- and off-years is necessary for analyses of the relationships between LAI and CC and the vegetation indices.
3.2 Correlations of vegetation indices with the LAI and CC
The results in Figs. 5 and 6 show the relationships between the vegetation indices and the LAI and CC. The WDRVI, SR, and NDVI displayed significant linear relationships with the LAI and CC in on-years, with R2 values of over 0.7 (p < 0.01), thus implying that they are good indicators for the estimation of LAI in on-years (Figs. 5 and 6). The relationships between the vegetation indices and the LAI and CC in the off-year tended to be exponential. The slopes of the relationships between vegetation indices and CC were very different between on- and off-years, so for the same vegetation indices, CC in on- and off-years were different. Different relationships of LAI and CC to vegetation indices between on- and off-years were caused by the different change trends of the LAI and CC between the on- and off-year (Fig. 2). The WDRVI and SR (p < 0.01) were slightly accurate than the NDVI (p < 0.05) for estimating the LAI in the off-year. The GI had the weakest relationship with the LAI (Fig. 5d), whereas it had the strongest relationship with CC in the off-year (Fig. 6d).
3.3 Calibration and validation of the LAI and CC estimation models
Because the same capacity was exhibited by the WDRVI (α = 0.1) and SR for explaining the LAI and CC (Fig. 5), we chose the WDRVI (α = 0.1) to build a linear regression model for the estimation of the LAI in on-years and an exponential model for estimating the LAI in off-years by using the 21 sample data points collected in 2011, 2014, and 2015. Additionally, the WDRVI (α = 0.1) was used for the estimation of CC in on-years and the GI was used for the estimation of CC in off-years. The algorithm equations used for the estimations of LAI and CC were as follows:
Equation (6) yielded the following accuracy values for the LAI estimation: RMSE = 0.28 m2/m2 and RMSEr = 7.40% in on-years and RMSE = 0.24 m2/m2 and RMSEr = 6.04% in off-years. Equation (7) yielded the following accuracy values for the CC estimation: RMSE = 0.15 g/m2 and RMSEr = 7.42% in on-years and RMSE = 0.25 g/m2 and RMSEr = 14.29% in off-years. A leave-one-out cross-validation method was used to test the robustness and prediction accuracy of the equations. The values of R2, RMSE, and RMSEr are given in Table 2. The high R2 and low RMSEr values between the predicted data and measured data show the good predictive ability of the models for both on- and off-years. There were no obvious differences in the predictive ability between on- and off-years for the estimation of the LAI. However, the predictive ability for CC in off-years was slightly poorer than that in on-years. The models were able to estimate the LAI with accuracy values amounting to an RMSE = 0.34 m2/m2 and RMSEr = 9.04% and to estimate the CC with accuracy values amounting to an RMSE = 0.25 g/m2 and RMSEr = 13.09%. These validation results confirmed that accurate LAI and CC estimates of Moso bamboo forest can be obtained based on simple relationships between the LAI and CC and vegetation indices from MODIS reflectance data.
3.4 Comparison of the LAI estimates with the MODIS LAI product
The relationship between the validation results for the LAI and ground-measured LAI was close to the 1:1 line, thus implying that there is no problem with underestimation and overestimation (Fig. 7a). Compared with ground-measured LAI, the LAI of Moso bamboo forest from the MODIS LAI product (MOD15A2) had serious problems with underestimation and significant errors were encountered (RMSEr = 58.76%) (Fig. 7). The relationship between the measured LAI and the MODIS LAI corrected by using the “upper-envelope” smoothing method showed improvements, but the corrected MODIS LAI remained underestimated, i.e., RMSEr = 37.90% (Fig. 7c). Therefore, the LAI of Moso bamboo forest derived from the simple regression model (Eq. (6)) is more accurate than both the MODIS LAI product (MOD15A2) and the corrected MODIS LAI.
Variations in the 8-day LAI average from 2000 to 2015 showed that LAI estimates from this study ranged between 1.9 and 5.9 m2/m2 and MODIS LAI estimates of Moso bamboo forest ranged between 0.7 and 6.5 m2/m2 (Fig. 8a). Although the MODIS LAI time series can also capture the growing phases typical of Moso bamboo forest, the data tended to increase early, display a relatively short steady period during June–August, and then decrease rapidly after reaching the maximum value compared with the LAI estimates from this study (Fig. 8a). The LAI estimates from this study were significantly greater than the MODIS LAI products in the winter season (p < 0.01), thus implying that there is a large amount of underestimation in the MODIS LAI estimates of Moso bamboo forest. The relative difference (RD) between the LAI estimates from this study and the MODIS LAI estimates was very large (i.e., > 30%) for the periods extending from DOY 1 to DOY 129 (about May 10) and from DOY 257 (about September 15) to DOY 365 (Fig. 8b).
3.5 Differences in the LAI and CC between on- and off-years
Variations in the 8-day average LAI from 2000 to 2015 showed that the LAI in off-years was slightly smaller than that in on-years during DOY 1 to DOY 121 (end of April), whereas the LAI in off-years was about 12% (p < 0.01) larger than that in on-years during DOY 121 (beginning of May) to DOY 241 (end of August) (Fig. 9a). The LAI values remained stable in the on-years and only slightly decreased from DOY 1 to DOY 113 (mid-April). From DOY 113 (mid-April) to DOY 169 (mid-June), the LAI increased rapidly from 2.4 to 4.7 m2/m2 in the off-years because of the spread of new leaves, and the LAI increased rapidly from 2.6 to 4.1 m2/m2 in the on-years because of the spread of leaves of the bamboo shoots (Fig. 9a). The increase in LAI due to the flushing of new leaves on old bamboo in off-years was greater than that due to the flushing of leaves of the bamboo shoots in the on-years. After the new leaves and leaves of the bamboo shoots stopped spreading in mid-June, the LAI values reached their highest level and remained stable during DOY 169 (mid-June) to DOY 241 (end of August). The LAI decreased at the beginning of September for both the on- and off-years. Variations in the 8-day average of CC in the on- and off-years (Fig. 9b) showed that CC in the off-years was about 50% (p < 0.01) smaller than that in the on-years during DOY 1 to DOY 113 (mid-April). The CC in the off-years increased steeply after DOY 113 (mid-April), which is when new leaves were spreading, and it was comparable to the CC in the on-years after DOY 169 (mid-June), which is when the growth of new leaves had finished.
4 Discussion
4.1 Effect of various factors on the accuracy of LAI and CC estimations
The relationships of canopy reflectance to LAI and CC were likely confounded by other factors such as the understory component, canopy architecture, and sun/view angle, which are known to contribute to the reduction of accuracy in empirical models (Zhang et al. 2008). The main factor focused on here was the different leaf change and bamboo shoot production characteristics of Moso bamboo forest between on- and off-years, which caused the relationships between vegetation indices and CC were very different between on- and off-years. Furthermore, the factor caused the differences in the calibration model results between on- and off-years. For example, the relationships between vegetation indices and LAI and CC in the on-years were linear, whereas those relationships were exponential in the off-year that was evaluated in detail (Figs. 5 and 6). The reflectance of visible red band 1 decreased earlier in the off-year during DOY 129 to DOY 153 which also caused the structure of the calibration models varied between on- and off-years. The NDVI increased earlier than the CC in Amazon forests was also reported (Hilker et al. 2017). The GI showed a higher sensitivity than the other vegetation indices to the CC in the off-year. The reason for this is that reflectance at 550 nm (central wavelength of the green band) might be highly sensitive to CC in some types of leaves, e.g., uniformly young leaves (Datt 1999).
The retrieval of LAI data from remote sensing imagery can be hindered by the effect of the understory component on the relationships between the LAI and remote sensing reflectance data (De Kauwe et al. 2011; Qi et al. 2014). Leaf change in off-year will increase the effect of the understory component on the relationships between biophysical parameters and vegetation indices and bring about greater uncertainties in the models for estimating the LAI and CC in off-years compared to on-years (Table 2).
In this study, it is hard to assess whether the obtained calibration models are robust because the sample data sets used for the model calibration was limited in spatial variability. Five pixels randomly and spatially distributed in Anji County were selected. The LAI for the five pixels were calculated using the obtained calibration models (Fig. S2). Good relationships between the estimated LAI and MODIS LAI estimates were found and the MODIS LAI estimates were smaller than the estimated LAI in most of the pixels, which indirectly showed the ability of the obtained calibration models in explaining spatial variation of LAI. Many in situ measurements are required to assess accuracy of the obtained calibration models in estimating spatial variation of LAI in the future.
4.2 Comparison of the LAI estimates with the MODIS LAI product
The LAI estimates of Moso bamboo forest from this study were more accurate than the MODIS LAI product. Our findings suggest that the LAI of Moso bamboo forest was underestimated by the MODIS LAI product. The main causes of unrealistic MODIS LAI values can be attributed to cloud contamination in the surface reflectance input data, the MODIS LAI algorithm itself, and misclassifications (De Kauwe et al. 2011; Fang et al. 2013). In the latest MODIS LAI algorithm, land surfaces are classified into eight biomes, which is an insufficient number to describe actual vegetation types. For example, Moso bamboo forest is not included in the MODIS land cover type product (MOD12Q1); it is classified as mixed forest instead (Xu et al. 2016). There are obvious differences between mixed forest and Moso bamboo forest, which can be attributed to the special growth characteristics of Moso bamboo. The rapid growth from shoots to bamboo within 40 days (Yen 2016) and the changing leaf period are accompanied by changes in the LAI and CC. Therefore, the probable reason for the uncertainties in the MODIS LAI of Moso bamboo forest is that some parameters of mixed forest, such as the soil and leaf optical properties, are not suitable for estimating the LAI of Moso bamboo forest based on the lookup table method (Xu et al. 2016).
The proposed empirical models in this study have a greater ability than the MODIS LAI algorithm to capture the dynamic changes of the LAI of Moso bamboo forest. The MODIS LAI of Moso bamboo forest tends to rise early and to decrease quickly after a short constant period, which is a trend that also has been found in the MODIS LAI of deciduous broadleaf forest (Wang et al. 2005). The MODIS LAI product is used as a key variable for the MODIS gross primary productivity/net primary productivity (GPP/NPP) product, and therefore, errors in LAI will lead to uncertainty in the MODIS GPP/NPP product (Heinsch et al. 2006; De Kauwe et al. 2011; Qi et al. 2014). An underestimation of the MODIS GPP product for Moso bamboo forest compared with observations from the flux tower has been confirmed (Xu et al. 2013), and it might be strongly related to the underestimation of the MODIS LAI product. The LAI estimates from this study had higher accuracies and agreed better with the temporal and dynamic trends of observations than the MODIS LAI, thus implying that they are superior for capturing the important and special phenological phases of Moso bamboo forest and can be used for estimating the GPP of Moso bamboo forest.
4.3 Differences in LAI and CC between on- and off-years
There were obvious differences in LAI and CC between on-years and off-years, especially for CC data from DOY 1 to DOY 169. Previous studies have shown that the LAI increases rapidly as new leaf flushing starts in May and ends in June, and thereafter, it remains stable from July onward (Grattani et al. 2008; Jin et al. 2010). During the leaf flushing period, leaf area increases considerably from 270 to 835 mm2 (Jin et al. 2010). The mechanisms for new leaf flushing in the on- and off-years are different. New leaf flushing in the on-years is due to growth of the bamboo shoots, whereas that in the off-years results from 2-year-old leaf defoliation and other changes. Thus, changes in LAI from DOY 113 (mid-April) to DOY 241 (end of August) between the on- and off-years were different (Fig. 9a).
Leaves of old bamboo turn yellow after the flushing of new shoots in March during on-years because the growth of bamboo shoots consumes considerable quantities of carbohydrates and nutrients (Qiu 1984; Hu 2011; Song et al. 2016). Furthermore, the LCC in on-years will remain stable throughout the whole year (see Fig. 3c) as the increases in CC due to leaves of new bamboo becoming dark green compensate for the decreases in CC due to leaves of old bamboo becoming yellow. The LCC in off-years is relatively small and stable during DOY 1 to DOY 137 (mid-April) when the leaf changing process occurs (Grattani et al. 2008). It increases rapidly from DOY 137 (mid-April) to DOY 169 (mid-June) as the leaf growth process occurs and remains stable from DOY 169 (mid-June) to DOY 241 (end of August) as leaf color turns dark green (Grattani et al. 2008; Song et al. 2016). The leaf color turns yellow in November because of leaf senescence, which results in obvious decreases in the LAI and CC in both on-years and off-years.
5 Conclusions
Previous studies have revealed underestimates in the MODIS LAI and GPP products of Moso bamboo forest. In order to obtain estimates of biophysical parameters that are more accurate, simplified empirical–statistical models were developed based on an analysis of the relationships between vegetation indices and in situ LAI and CC observations. The calibration models for estimating LAI and CC between on- and off-years were different because of different leaf change and bamboo shoot production characteristics. This implies that it is necessary to build separate models for on- and off-years, which will improve the potential for retrieving accurate LAI and CC estimates. A relationship analysis between vegetation indices and LAI and CC showed that the WDRVI was a good indicator for estimations of LAI in on- and off-years and CC in off-years. The GI was the vegetation index most sensitive to CC in off-years. The models presented in this study produced better estimates of LAI than the MODIS LAI product (MOD15 A2). It was identified that the MODIS LAI time series has a problem with rising early, having a short constant period, and decreasing soon after reaching the maximum value. Therefore, caution must be observed when analyzing the phenological characteristics of Moso bamboo forest directly with the MODIS LAI product. The different variations in LAI and CC between on- and off-years were first discussed based on the LAI and CC estimates from this study. This study then presented simple methods to estimate the LAI and CC accurately from MODIS reflectance data based on the knowledge of complicated leaf change and bamboo shoot production characteristics of Moso bamboo forest between on- and off-years, and it was concluded that the statistical model driven by MODIS data had the potential to accurately estimate the LAI and CC. This study also identified superior biophysical parameters for both assessing the growth and health status of Moso bamboo forest and for estimating the amount of carbon fixation within Moso bamboo forest.
References
Chen JM, Cihlar J (1996) Retrieving leaf area index of boreal conifer forests using Landsat TM images. Remote Sens Environ 55:153–162. https://doi.org/10.1016/0034-4257(95)00195-6
Chen XG, Zhang XQ, Zhang YP, Booth T, He XH (2009) Changes of carbon stocks in bamboo stands in China during 100 years. For Ecol Manag 258:1489–1496. https://doi.org/10.1016/j.foreco.2009.06.051
Cohen WB, Maiersperger TK, Turner DP, Ritts WD, Pflugmacher D, Kennedy RE, Kirschbaum A, Running SW, Costa M, Gower ST (2006) MODIS land cover and LAI collection 4 product quality across nine sites in the western hemisphere. IEEE Trans Geosci Rem Sens 44:1843–1857. https://doi.org/10.1109/TGRS.2006.876026
Colombo R, Bellingeri D, Fasolini D, Marino CM (2003) Retrieval of leaf area index in different vegetation types using high resolution satellite data. Remote Sens Environ 86:120–131. https://doi.org/10.1016/S0034-4257(03)00094-4
Coops NC, Stone C, Culvenor DS, Chisholm LA, Merton RN (2003) Chlorophyll content in eucalypt vegetation at the leaf and canopy scales as derived from high resolution spectral data. Tree Physiol 23:23–31
Datt B (1999) A new reflectance index for remote sensing of chlorophyll content in higher plants: tests using eucalyptus leaves. J Plant Physiol 154:30–36. https://doi.org/10.1016/S0176-1617(99)80314-9
De Kauwe MG, Disney MI, Quaife T, Lewis P, Williams M (2011) An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest. Remote Sens Environ 115:767–780. https://doi.org/10.1016/j.rse.2010.11.004
Du H, Zhou G, Fan W, Ge H, Xu X, Shi Y, Fan W (2010) Spatial heterogeneity and carbon contribution of aboveground biomass of moso bamboo by using geostatistical theory. Plant Ecol 207:131–139. https://doi.org/10.1007/s11258-009-9659-3
Fang HL, Liang SL (2005) A hybrid inversion method for mapping leaf area index from MODIS data: experiments and application to broadleaf and needleleaf canopies. Remote Sens Environ 94:405–424. https://doi.org/10.1016/j.rse.2004.11.001
Fang HL, Li WJ, Myneni RB (2013) The impact of potential land cover misclassification on MODIS leaf area index (LAI) estimation: a statistical perspective. Remote Sens 5:830–844. https://doi.org/10.3390/rs5020830
Fang HL, Wei SS, Liang SL (2012) Validation of MODIS and CYCLOPES LAI products using global field measurement data. Remote Sens Environ 119:43–54. https://doi.org/10.1016/j.rse.2011.12.006
FAO (2010) Global forest resources assessment 2010: main report food and agriculture organization of the United Nations, Rome. https://doi.org/10.1016/s0264-8377(03)00003-6
Gitelson AA, Viña A, Ciganda V, Rundquist DC, Arkebauer TJ (2005) Remote estimation of canopy chlorophyll content in crops. Geophys Res Lett 32:L08403. https://doi.org/10.1029/2005GL022688
Gitelson AA, Viña A, Verma SB, Rundquist DC, Arkebauer TJ, Keydan G, Leavitt B, Ciganda V, Burba GG, Suyker AE (2006) Relationship between gross primary production and chlorophyll content in crops: implications for the synoptic monitoring of vegetation productivity. J Geophys Res 111:D08S11. https://doi.org/10.1029/2005JD006017
Gitelson AA, Wardlow BD, Keydan GP, Leavitt B (2007) An evaluation of MODIS 250-m data for green LAI estimation in crops. Geophys Res Lett 34:L20403. https://doi.org/10.1029/2007GL031620
Grattani L, Crescente MF, Varone L, Fabrini ED (2008) Growth pattern and photosynthetic activity of different bamboo species growing in the Botanical Garden of Rome. Flora 203:77–84. https://doi.org/10.1016/j.flora.2007.11.002
Gu CY (2013) Satellite-based retrieval of canopy parameters of Moso bamboo forest with PROSAIL radiative transfer model. Dissertation, Zhejiang A&F University [in Chinese]
Gu YX, Bélair S, Mahfouf J-F, Deblonde G (2006) Optimal interpolation analysis of leaf area index using MODIS data. Remote Sens Environ 104:283–296. https://doi.org/10.1016/j.rse.2006.04.021
Guindin-Garcia N, Gitelson AA, Arkebauer TJ, Shanahan J, Weiss A (2012) An evaluation of MODIS 8- and 16-day composite products for monitoring maize green leaf area index. Agric For Meteorol 161:15–25. https://doi.org/10.1016/j.agrformet.2012.03.012
Heinsch FA, Zhao M, Running SW, Kimball JS, Nemani RR, Davis KJ, Bolstad PV, Cook BD, Desai AR, Ricciuto DM, Law BE, Oechel WC, Kwon H, Luo HY, Wofsy SC, Dunn AL, Munger JW, Baldocchi DD, Xu LK, Hollinger DY, Richardson AD, Stoy PC, Siqueira MBS, Monson RK, Burns SP, Flanagan LB (2006) Evaluation of remote sensing based terrestrial productivity from MODIS using regional tower eddy flux network observations. IEEE Trans Geosci Rem Sens 44:1908–1925. https://doi.org/10.1109/TGRS.2005.853936
Hilker T, Galvão LS, Aragão LEOC, Moura YMD, Amaral CHD, Lyapustin AI, Wu J, Albert LP, Ferreira MJ, Anderson LO, dos Santos VAHF, Prohaska N, Tribuzy E, Ceron JVB, Saleska SR, Wang Y, de Carvalho Goncalves JF, de Oliveira Junior RC, Rodrigues JVFC, Garcia MN (2017) Vegetation chlorophyll estimates in the amazon from multi-angle MODIS observations and canopy reflectance model. Int J Appl Earth Obs 58:278–287. https://doi.org/10.1016/j.jag.2017.01.014
Houborg R, Anderson M, Daughtry C (2009) Utility of an image-based canopy reflectance modeling tool for remote estimation of LAI and leaf chlorophyll content at the field scale. Remote Sens Environ 113:259–274. https://doi.org/10.1016/j.rse.2008.09.014
Houborg R, Soegaard H, Boegh E (2007) Combining vegetation index and model inversion methods for the extraction of key vegetation biophysical parameters using Terra and Aqua MODIS reflectance data. Remote Sens Environ 106:39–58. https://doi.org/10.1016/j.rse.2006.07.016
Hu YB (2011) Effects of fertilizations on leaf characteristics and photosynthesis character during the shoot period to deciduous period of Phyllostarchys pubescens. Dissertation, Zhejiang A&F University [in Chinese]
Jacquemoud S, Bacour C, Poilvé H, Frangi J-P (2000) Comparison of four radiative transfer models to simulate plant canopies reflectance: direct and inverse mode. Remote Sens Environ 74:471–481. https://doi.org/10.1016/S0034-4257(00)00139-5
Jin XC, Jin AW, Song YD, Lou JF, Mei SM (2010) Canopy development and photosynthesis with fertilization during leaf expansion period of Phyllostarchys pubescens. J Zhejiang For Coll 27:57–62 [in Chinese]
Kleinhenz V, Midmore D J (2001) Aspects of bamboo agronomy. Adv agron 74:99–153. https://doi.org/10.1016/S0065-2113(01)74032-1
Komatsu H, Onozawa Y, Kume T, Tsuruta K, Kumagai T, Shinohara Y, Otsuki K (2010) Stand-scale transpiration estimates in a Moso bamboo forest: II comparison with coniferous forests. For Ecol Manag 260:1295–1302. https://doi.org/10.1016/j.foreco.2010.06.040
Komatsu H, Onozawa Y, Kume T, Tsuruta K, Shinohara Y, Otsuki K (2012) Canopy conductance for a Moso bamboo (Phyllostachys pubescens) forest in western Japan. Agric For Meteorol 156:111–120. https://doi.org/10.1016/j.agrformet.2012.01.004
Leuning R, Cleugh HA, Zeglin SJ, Hughes D (2005) Carbon and water fluxes over a temperate Eucalyptus forest and tropical wet/dry savanna in Australia: measurements and comparison with MODIS remote sensing estimates. Agric For Meteorol 129:153–173. https://doi.org/10.1016/j.agrformet.2004.12.004
Li R, Werger MJA, During HJ, Zhong ZC (1998) Biennial variation in production of new shoots in groves of the giant bamboo Phyllostachys pubescens in Sichuan. China Plant Eco 135:103–112. https://doi.org/10.1023/A:1009761428401
Li XJ, Xiao ZQ, Wang JD, Qu Y, Jin HA (2014) Dual ensemble Kalman filter assimilation method for estimating time series LAI. J Remote Sens 18:27–44
Maire GL, Marsden C, Verhoef W, Ponzoni FJ, Seen DL, Bégué A, Stape J-L, Nouvellon Y (2011) Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of eucalyptus, plantations. Remote Sens Environ 115:586–599. https://doi.org/10.1016/j.rse.2010.10.004
Mao F, Du H, Zhou G, Li X, Xu X, Li P, Sun S (2017) Coupled LAI assimilation and BEPS model for analyzing the spatiotemporal pattern and heterogeneity of carbon fluxes of the bamboo forest in Zhejiang Province, China. Agric For Meteorol 242:96–108. https://doi.org/10.1016/j.agrformet.2017.03.022
Meyers TP, Hollinger SE (2004) An assessment of storage terms in the surface energy balance of maize and soybean. Agric For Meteorol 125:105–115. https://doi.org/10.1016/j.agrformet.2004.03.001
Olivas PC, Oberbauer SF, Clark DB, Clark DA, Ryan MG, O’Brien JJ, Ordoñez H (2013) Comparison of direct and indirect methods for assessing leaf area index across a tropical rain forest landscape. Agric For Meteorol 177:110–116. https://doi.org/10.1016/j.agrformet.2013.04.010
Olsoy PJ, Mitchell JJ, Levia DF, Clark PE, Glenn NF (2016) Estimation of big sagebrush leaf area index with terrestrial laser scanning. Ecol Indic 61:815–821. https://doi.org/10.1016/j.ecolind.2015.10.034
Qi Y, Li F, Liu Z, Jin G (2014) Impact of understorey on overstorey leaf area index estimation from optical remote sensing in five forest types in northeastern China. Agric For Meteorol s198–199:72–80. https://doi.org/10.1016/j.agrformet.2014.08.001
Qiu FG (1984) The on-year and off-year of Phyllostarchys pubescens forests and their control. J Bamboo Res 3:62–69 [in Chinese]
Running SW, Nemani RR, Peterson DL, Band LE, Potts DF, Pierce LL, Spanner MA (1989) Mapping regional forest evapotranspiration and photosynthesis by coupling satellite data with ecosystem simulation. Ecology 70:1090–1101. https://doi.org/10.2307/1941378
Sellers PJ (1997) Modeling the exchange of energy, water, and carbon between continents and atmosphere. Science 275:602–609
Song X, Peng C, Zhou G, Gu H, Li Q, Zhang C (2016) Dynamic allocation and transfer of non-structural carbohydrates, a possible mechanism for the explosive growth of Moso bamboo (Phyllostachys heterocycla). Sci Rep 6:25908. https://doi.org/10.1038/srep25908
Wang Q, Tenhunen J, Dinh NQ, Reichstein M, Otieno D, Granier A, Pilegarrd K (2005) Evaluation of seasonal variation of MODIS derived leaf area index at two European deciduous broadleaf forest sites. Remote Sens Environ 96:475–484. https://doi.org/10.1016/j.rse.2005.04.003
Wu C, Zheng N, Quan T, Huang W (2008) Estimating chlorophyll content from hyperspectral vegetation indices: modeling and validation. Agric For Meteorol 148:1230–1241. https://doi.org/10.1016/j.agrformet.2008.03.005
Wu C, Zheng N, Quan T, Huang W, Rivard B, Feng J (2009) Remote estimation of gross primary production in wheat using chlorophyll-related vegetation indices. Agric For Meteorol 149:1015–1021. https://doi.org/10.1016/j.agrformet.2008.12.007
Xu X, Du H, Zhou G, Li P (2016) Method for improvement of MODIS leaf area index products based on pixel-to-pixel correlations. Eur J Remote Sens 49:57–72. https://doi.org/10.5721/EuJRS20164904
Xu X, Zhou G, Liu S, Du H, Mo L, Shi Y, Jiang H, Zhou Y, Liu E (2013) Implications of ice storm damages on the water and carbon cycle of bamboo forests in southeastern China. Agric For Meteorol 177:35–45. https://doi.org/10.1016/j.agrformet.2013.04.005
Yen TM (2015) Comparing aboveground structure and aboveground carbon storage of an age series of Moso bamboo forests subjected to different management strategies. J For Res 20:1–8. https://doi.org/10.1007/s10310-014-0455-0
Yen TM (2016) Culm height development, biomass accumulation and carbon storage in an initial growth stage for a fast-growing Moso bamboo (Phyllostachy pubescens). Bot Stud 57:1–9. https://doi.org/10.1186/s40529-016-0126-x
Yen TM, Lee JS (2011) Comparing aboveground carbon sequestration between Moso bamboo (Phyllostachys heterocycla) and China fir (Cunninghamia lanceolata) forests based on the allometric model. Forest Ecol Manag 261:995–1002. https://doi.org/10.1016/j.foreco.2010.12.015
Zhang Y, Chen JM, Miller JR, Noland TL (2008) Leaf chlorophyll content retrieval from airborne hyperspectral remote sensing imagery. Remote Sens Environ 112:3234–3247. https://doi.org/10.1016/j.rse.2008.04.005
Zhang Y, Qu Y, Wang J, Liang S, Liu Y (2012) Estimating leaf area index from modis and surface meteorological data using a dynamic Bayesian network. Remote Sens Environ 127:30–43. https://doi.org/10.1016/j.rse.2012.08.015
Zhou B, Li Z, Wang X, Cao Y, An Y, Deng Z, Letu G, Wang G, Gu L (2011) Impact of the 2008 ice storm on Moso bamboo plantations in Southeast China. J Geophys Res 116:G00H06. https://doi.org/10.1029/2009JG001234
Funding
This study was supported by the National Natural Science Foundation [31500520, 31370637, and 31670644] and Natural Science Foundation of Zhejiang Province [LQ15C160003].
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Handling Editor: Erwin Dreyer
Contribution of the co-authors
Xiaojun Xu, Huaqiang Du, and Guomo Zhou: conception and design, fieldwork, analysis and interpretation, statistical analysis, and writing the article. Fangjie Mao, Xuejian Li, Dien Zhu, Yangguang Li, and Lu Cui: fieldwork, analysis, and interpretation.
Rights and permissions
About this article
Cite this article
Xu, X., Du, H., Zhou, G. et al. Remote estimation of canopy leaf area index and chlorophyll content in Moso bamboo (Phyllostachys edulis (Carrière) J. Houz.) forest using MODIS reflectance data. Annals of Forest Science 75, 33 (2018). https://doi.org/10.1007/s13595-018-0721-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s13595-018-0721-y