Shezhou Luo

Local leaf size diversity (LLSD) is an essential functional indicator of plant biodiversity; however, massive challenges are encountered when quantifying it and decoding its global ecological patterns. To address this limitation, the present study defined a quantitative indicator of LLSD, termed coefficient of variation index (CVI), for the leaf sizes, regardless of plant species, collected in each sampling site. Then, we innovatively derived a set of global CVI values from a published dataset, which was obtained through a meta-analysis of global leaf area samples and their related climate factors. Our macroecological analyses indicate that the CVI values vary across continents and fluctuate with latitude. The global CVI values are predominantly influenced by the mean temperature of the coldest month during the growing season in the negative correlation mode. When two leading climate drivers are considered, the global CVI values are primarily influenced by the mean temperature during growing season and the mean annual sum precipitation. Overall, all of these contributions are pioneering in their implications for characterizing the global distribution and ecological patterns of LLSD and advancing the cutting-edge research domain of leaf functional biodiversity to a new quantitative stage.

Point cloud filtering is a crucial step in most airborne light detection and ranging (LiDAR) applications. Many filtering algorithms have been proposed, but the filtering effect has some limitations in complex environments. To improve the filtering effect in complex terrain, a multilevel adaptive filter (MAF) combining morphological reconstruction and thin plate spline (TPS) interpolation is proposed. The digital elevation model (DEM) generated in each iteration is used as the marker image for morphological reconstruction to extract ground pixels, and an adaptive residual threshold is achieved by using terrain gradient as a compensation. The benchmark dataset provided by the International Society for Photogrammetry and Remote Sensing (ISPRS) and another LiDAR dataset in northwestern China were used to evaluate the filtering performance of MAF. For the ISPRS benchmark dataset, MAF obtained the lowest average total error (3.72%) and highest average kappa coefficient (87.16%) compared with eight classic filtering algorithms. For the dataset in northwestern China, the DEM generated from the filtering result of MAF obtained higher accuracy than the filtering result of TerraScan. Overall, the MAF achieved promising results without considering the selection of filtering window, which may enhance the robustness and applicability of the algorithm in different environments.

Measuring forest aboveground biomass (AGB) at local to regional scales is critical to understanding their role in regional and global carbon cycles. The Three-North Shelterbelt Forest Program (TNSFP) is the largest ecological restoration project in the world, and has been ongoing for over 40 years. In this study, we developed models to estimate the planted forest aboveground biomass (PF_AGB) for Yulin, a typical area in the project. Surface reflectances in the study area from 1978 to 2013 were obtained from Landsat series images, and integrated forest z-scores were constructed to measure afforestation and the stand age of planted forest. Normalized difference vegetation index (NDVI) was combined with stand age to develop an initial model to estimate PF_AGB. We then developed additional models that added environment variables to our initial model, including climatic factors (average temperature, total precipitation, and total sunshine duration) and a topography factor (slope). The model which combined the total precipitation and slope greatly improved the accuracy of PF_AGB estimation compared to the initial model, indicating that the environmental variables related to water distribution indirectly affected the growth of the planted forest and the resulting AGB. Afforestation in the study area occurred mainly in the early 1980s and early 21st century, and the PF_AGB in 2003 was 2.3 times than that of 1998, since the fourth term TNSFP started in 2000. The PF_AGB in 2013 was about 3.33 times of that in 2003 because many young trees matured. The leave-one-out cross-validation (LOOCV) approach showed that our estimated PF_AGB had a significant correlation with field-measured data (correlation coefficient (r) = 0.89, p < 0.001, root mean square error (RMSE) = 6.79 t/ha). Our studies provided a method to estimate long time series PF_AGB using satellite repetitive measures, particularly for arid or semi-arid areas.

Optical methods require model inversion to infer plant area index (PAI) and woody area index (WAI) of leaf-on and leaf-off forest canopy from gap fraction or radiation attenuation measurements. Several inversion models have been developed previously, however, a thorough comparison of those inversion models in obtaining the PAI and WAI of leaf-on and leaf-off forest canopy has not been conducted so far. In the present study, an explicit 3D forest scene series with different PAI, WAI, phenological periods, stand density, tree species composition, plant functional types, canopy element clumping index, and woody component clumping index was generated using 50 detailed 3D tree models. The explicit 3D forest scene series was then used to assess the performance of seven commonly used inversion models to estimate the PAI and WAI of the leaf-on and leaf-off forest canopy. The PAI and WAI estimated from the seven inversion models and simulated digital hemispherical photography images were compared with the true PAI and WAI of leaf-on and leaf-off forest scenes. Factors that contributed to the differences between the estimates of the seven inversion models were analyzed. Results show that both the factors of inversion model, canopy element and woody component projection functions, canopy element and woody component estimation algorithms, and segment size are contributed to the differences between the PAI and WAI estimated from the seven inversion models. There is no universally valid combination of inversion model, needle-to-shoot area ratio, canopy element and woody component clumping index estimation algorithm, and segment size that can accurately measure the PAI and WAI of all leaf-on and leaf-off forest canopies. The performance of the combinations of inversion model, needle-to-shoot area ratio, canopy element and woody component clumping index estimation algorithm, and segment size to estimate the PAI and WAI of leaf-on and leaf-off forest canopies is the function of the inversion model as well as the canopy element and woody component clumping index estimation algorithm, segment size, PAI, WAI, tree species composition, and plant functional types. The impact of canopy element and woody component projection function measurements on the PAI and WAI estimation of the leaf-on and leaf-off forest canopy can be reduced to a low level (<4%) by adopting appropriate inversion models.

Forest biomass is an important parameter for quantifying and understanding biological and physical processes on the Earth’s surface. Rapid, reliable, and objective estimations of forest biomass are essential to terrestrial ecosystem research. The Geoscience Laser Altimeter System (GLAS) produced substantial scientific data for detecting the vegetation structure at the footprint level. This study combined GLAS data with MODIS/BRDF (Bidirectional Reflectance Distribution Function) and ASTER GDEM data to estimate forest aboveground biomass (AGB) in Xishuangbanna, Yunnan Province, China. The GLAS waveform characteristic parameters were extracted using the wavelet method. The ASTER DEM was used to compute the terrain index for reducing the topographic influence on the GLAS canopy height estimation. A neural network method was applied to assimilate the MODIS BRDF data with the canopy heights for estimating continuous forest heights. Forest leaf area indices (LAIs) were derived from Landsat TM imagery. A series of biomass estimation models were developed and validated using regression analyses between field-estimated biomass, canopy height, and LAI. The GLAS-derived canopy heights in Xishuangbanna correlated well with the field-estimated AGB (R² = 0.61, RMSE = 52.79 Mg/ha). Combining the GLAS estimated canopy heights and LAI yielded a stronger correlation with the field-estimated AGB (R² = 0.73, RMSE = 38.20 Mg/ha), which indicates that the accuracy of the estimated biomass in complex terrains can be improved significantly by integrating GLAS and optical remote sensing data.

Accurate land cover classification information is a critical variable for many applications. This study presents a method to classify land cover using the fusion data of airborne discrete return LiDAR (Light Detection and Ranging) and CASI (Compact Airborne Spectrographic Imager) hyperspectral data. Four LiDAR-derived images (DTM, DSM, nDSM, and intensity) and CASI data (48 bands) with 1 m spatial resolution were spatially resampled to 2, 4, 8, 10, 20 and 30 m resolutions using the nearest neighbor resampling method. These data were thereafter fused using the layer stacking and principal components analysis (PCA) methods. Land cover was classified by commonly used supervised classifications in remote sensing images, i.e., the support vector machine (SVM) and maximum likelihood (MLC) classifiers. Each classifier was applied to four types of datasets (at seven different spatial resolutions): (1) the layer stacking fusion data; (2) the PCA fusion data; (3) the LiDAR data alone; and (4) the CASI data alone. In this study, the land cover category was classified into seven classes, i.e., buildings, road, water bodies, forests, grassland, cropland and barren land. A total of 56 classification results were produced, and the classification accuracies were assessed and compared. The results show that the classification accuracies produced from two fused datasets were higher than that of the single LiDAR and CASI data at all seven spatial resolutions. Moreover, we find that the layer stacking method produced higher overall classification accuracies than the PCA fusion method using both the SVM and MLC classifiers. The highest classification accuracy obtained (OA = 97.8%, kappa = 0.964) using the SVM classifier on the layer stacking fusion data at 1 m spatial resolution. Compared with the best classification results of the CASI and LiDAR data alone, the overall classification accuracies improved by 9.1% and 19.6%, respectively. Our findings also demonstrated that the SVM classifier generally performed better than the MLC when classifying multisource data; however, none of the classifiers consistently produced higher accuracies at all spatial resolutions.

The Leaf Area Index (LAI) is one of the most important parameters, which controls biological and physical processes associated with vegetation on the Earth's surface, such as photosynthesis, respiration, transpiration, carbon and nutrient cycles, and rainfall interception. Therefore, rapid, reliable and objective estimations of LAI are essential. In this study, we used a simplified approach to compute laser penetration index (LPI) from airborne discrete‐return LiDAR data, and LPI was computed based on corrected echo intensity for the first time. Using the variable of LPI, we built an LAI estimation model based on Beer‐Lambert law. This approach was applied to a forest area in Dayakou, Gansu Province of China. The accuracy of the corrected intensity‐derived LAI estimation model was compared with that of uncorrected intensity‐derived and echo counts‐derived model. The results show that the corrected echoes intensity can improve the accuracy of LAI estimation. To assess validity and generalization of the model, we validated the optimum model using the Leave‐One‐Out Cross‐Validation (LOOCV) procedure, and show that the model has no overfitting and is more general. Finally, we examined the accuracy of predicted LAIs with 16 field‐measured LAIs which were not involved in the modeling process and indicated that LAI can be estimated with a high accuracy in mountains area by corrected echo intensity. The comparison indicates the LiDAR‐derived LAI (R²=0.825, RMSE=0.165) is much higher than that of the LAI from Landsat TM images (R²=0.605, RMSE=0.257). It means that airborne LiDAR data can be used to obtain high‐accuracy LAI estimation and provide reliable data for ecological environment research.