The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) launched on September 15th, 2018 and this mission offers an extraordinary opportunity to contribute to an assessment of forest resources at multiple spatial scales. This study served to develop a methodology for utilizing ICESat-2 data over vegetated areas. The specific objectives were to: (1) derive a simulated ICESat-2 photon-counting lidar (PCL) vegetation product using airborne lidar data, (2) examine the use of simulated PCL metrics for modelling aboveground biomass (AGB) along ICESat-2 profiles using a simulated ICESat-2 PCL vegetation product and reference AGB estimated from airborne lidar data, and (3) estimate forest canopy cover using simulated PCL canopy product data and airborne lidar-derived canopy cover. Using existing airborne lidar data for Sam Houston National Forest (SHNF) in Texas and known ICESat-2 track locations, PCL simulations were carried out. Three scenarios were analyzed in this study; 1) simulated data without the addition of noise, 2) processed simulated data for daytime, and 3) nighttime scenarios. Segments measuring 100 m along the proposed ICESat-2 tracks were used to extract simulated PCL metrics from each of the three data scenarios and spatially coincident, reference airborne lidar-estimated AGB and airborne lidar canopy cover estimates. Linear regression models were then developed with a subset of the simulated PCL segments to estimate AGB and canopy cover and their performance assessed using separate testing sets. AGB model testing with the simulated dataset without noise, nighttime and daytime scenarios resulted in R-2 values of 0.79, 0.79 and 0.63 respectively, with root mean square error (RMSE) values of 19.16 Mg/ha, 19.23 Mg/ha, and 25.35 Mg/ha. Predictive models for canopy cover (4.6 m) achieved R-2 values of 0.93, 0.75 and 0.63 and RMSE values of 6.36%, 12.33% and 15.01% for the simulated dataset without noise, nighttime and daytime scenarios respectively. Findings from this study suggest the potential of ICESat-2 for estimating AGB, given the photon detection rate and background noise anticipated by ICESat-2 under temperate forest settings, and especially in the data format provided by standard ICESat-2 data vegetation products.

NASA’s Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) was launched in September, 2018. The satellite carries a single instrument, ATLAS (Advanced Topographic Laser Altimeter System), a green wavelength, photon-counting lidar, enabling global measurement and monitoring of elevation with a primary focus on the cryosphere. Although bathymetric mapping was not one of the design goals for ATLAS, pre-launch work by our research team showed the potential to map bathymetry with ICESat-2, using data from MABEL (Multiple Altimeter Beam Experimental Lidar), NASA’s high-altitude airborne ATLAS emulator, and adapting the laser-radar equation for ATLAS specific parameters. However, many of the sensor variables were only approximations, which limited a full assessment of the bathymetric mapping capabilities of ICESat-2 during pre-launch studies. Following the successful launch, preliminary analyses of the geolocated photon returns have been conducted for a number of coastal sites, revealing several salient examples of seafloor detection in water depths of up to ~40 m. The geolocated seafloor photon returns cannot be taken as bathymetric measurements, however, since the algorithm used to generate them is not designed to account for the refraction that occurs at the air–water interface or the corresponding change in the speed of light in the water column. This paper presents the first early on-orbit validation of ICESat-2 bathymetry and quantification of the bathymetric mapping performance of ATLAS using data acquired over St. Thomas, U.S. Virgin Islands. A refraction correction, developed and tested in this work, is applied, after which the ICESat-2 bathymetry is compared against high-accuracy airborne topo-bathymetric lidar reference data collected by the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA). The results show agreement to within 0.43—0.60 m root mean square error (RMSE) over 1 m grid resolution for these early on-orbit data. Refraction-corrected bottom return photons are then inspected for four coastal locations around the globe in relation to Visible Infrared Imaging Radiometer Suite (VIIRS) Kd(490) data to empirically determine the maximum depth mapping capability of ATLAS as a function of water clarity. It is demonstrated that ATLAS has a maximum depth mapping capability of nearly 1 Secchi in depth for water depths up to 38 m and Kd(490) in the range of 0.05–0.12 m−1. Collectively, these results indicate the great potential for bathymetric mapping with ICESat-2, offering a promising new tool to assist in filling the global void in nearshore bathymetry.

Canopy height serves as a good indicator of forest carbon content. Remote sensing-based direct estimations of canopy height are usually based on Light Detection and Ranging (LiDAR) or Synthetic Aperture Radar (SAR) interferometric data. LiDAR data is scarcely available for the Indian tropics, while Interferometric SAR data from commercial satellites are costly. High temporal decorrelation makes freely available Sentinel-1 interferometric data mostly unsuitable for tropical forests. Alternatively, other remote sensing and biophysical parameters have shown good correlation with forest canopy height. The study objective was to establish and validate a methodology by which forest canopy height can be estimated from SAR and optical remote sensing data using machine learning models i.e., Random Forest (RF) and Symbolic Regression (SR). Here, we analysed the potential of Sentinel-1 interferometric coherence and Sentinel-2 biophysical parameters to propose a new method for estimating canopy height in the study site of the Bhitarkanika wildlife sanctuary, which has mangrove forests. The results showed that interferometric coherence, and biophysical variables (Leaf Area Index (LAI) and Fraction of Vegetation Cover (FVC)) have reasonable correlation with canopy height. The RF model showed a Root Mean Squared Error (RMSE) of 1.57 m and R2 value of 0.60 between observed and predicted canopy heights; whereas, the SR model through genetic programming demonstrated better RMSE and R2 values of 1.48 and 0.62 m, respectively. The SR also established an interpretable model, which is not possible via any other machine learning algorithms. The FVC was found to be an essential variable for predicting forest canopy height. The canopy height maps correlated with ICESat-2 estimated canopy height, albeit modestly. The study demonstrated the effectiveness of Sentinel series data and the machine learning models in predicting canopy height. Therefore, in the absence of commercial and rare data sources, the methodology demonstrated here offers a plausible alternative for forest canopy height estimation.

The leaf area index (LAI) is a key indicator of the status of forest ecosystems that is important for understanding global carbon and water cycles as well as terrestrial surface energy balances and the impacts of climate change. Machine learning (ML) methods offer promising ways of generating spatially explicit LAI data covering large regions based on optical images. However, there have been few efforts to analyze the LAI in heterogeneous subtropical forests with complex terrain by fusing high-resolution multi-sensor data from the Sentinel-1 Synthetic Aperture Radar (SAR), Sentinel-2 Multi Spectral Instrument (MSI), and Advanced Land Observing Satellite-1 digital elevation model (DEM). Here, forest LAI mapping was performed by integrating the MSI, SAR, and DEM data using a stacking learning (SL) approach that incorporates distinct predictions from a set of optimized individual ML algorithms. The method’s performance was evaluated by comparison to field forest LAI measurements acquired in Xingguo and Gandong of subtropical China. The results showed that the addition of the SAR and DEM images using the SL model compared to the inputs of only optical images reduced the mean absolute error (MAE) and root mean square error (RMSE) by 26% and 18%, respectively, in Xingguo, and by 12% and 8%, respectively, in Gandong. Furthermore, the combination of all images had the best prediction performance. SL was found to be more robust and accurate than conventional individual ML models, while the MAE and RMSE were decreased by 71% and 64%, respectively, in Xingguo, and by 68% and 59%, respectively, in Gandong. Therefore, the SL model using the three-source data combination produced satisfied prediction accuracy with the coefficients of determination (R2), MAE, and RMSE of 0.96, 0.17, and 0.28, respectively, in Xingguo and 0.94, 0.30, and 0.47, respectively, in Gandong. This study revealed the potential of the SL algorithm for retrieving the forest LAI using multi-sensor data in areas with complex terrain.

Leaf area index (LAI) is a critical vegetation structural variable and is essential in the feedback of vegetation to the climate system. The advancement of the global Earth Observation has enabled the development of global LAI products and boosted global Earth system modeling studies. This overview provides a comprehensive analysis of LAI field measurements and remote sensing estimation methods, the product validation methods and product uncertainties, and the application of LAI in global studies. First, the paper clarifies some definitions related to LAI and introduces methods to determine LAI from field measurements and remote sensing observations. After introducing some major global LAI products, progresses made in temporal compositing and prospects for future LAI estimation are analyzed. Subsequently, the overview discusses various LAI product validation schemes, uncertainties in global moderate resolution LAI products, and high resolution reference data. Finally, applications of LAI in global vegetation change, land surface modeling, and agricultural studies are presented. It is recommended that (1) continued efforts are taken to advance LAI estimation algorithms and provide high temporal and spatial resolution products from current and forthcoming missions; (2) further validation studies be conducted to address the inadequacy of current validation studies, especially for underrepresented regions and seasons; and (3) new research frontiers, such as machine learning algorithms, light detection and ranging technology, and unmanned aerial vehicles be pursued to broaden the production and application of LAI.

ICESat-2（Ice, Cloud, and land Elevation Satellite-2）是美国NASA（National Aeronautics and Space Administration）在2018年发射的激光测高卫星,其上搭载的激光测高系统ATLAS（Advanced Topographic Laser Altimeter System）采用微脉冲多波束光子计数激光雷达系统,因其低能耗、高探测灵敏度、高重复频率的特性极大改善了沿轨采样密度,但也使获取的数据中包含大量的噪声,如何有效实现光子点云去噪分类成为后续应用的关键,也是当前研究的热点和难点,为此本文提出一种基于卷积神经网络的光子点云去噪和分类算法。首先将光子点云按照沿轨和高程方向划分格网,去除明显的噪声光子,并将每个粗信号光子点栅格化为影像;然后基于少量样本构建的卷积神经网络分类模型实现光子点云精去噪和分类;最后利用机载激光雷达数据进行验证,并与ATL08产品的去噪分类结果进行对比。结果表明,对于裸地和森林区域,卷积神经网络算法均能有效去除噪声光子,特别对于森林区域,可同时实现去噪和分类;其中,裸地区域地表计算的R²和RMSE分别为1.0和0.72 m,森林区域地表和树冠计算的R²分别为1.0和0.70, RMSE分别为1.11 m和4.99 m。本文利用深度学习算法实现光子点云去噪分类,在裸地和森林区域均取得了较好的结果,为后续光子点云数据处理提供了参考。

The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission employs a micro-pulse photon-counting LiDAR system for mapping and monitoring the biomass and carbon of terrestrial ecosystems over large areas. In preparation for ICESat-2 data processing and applications, this paper aimed to develop and validate an effective algorithm for better estimating ground elevation and vegetation height from photon-counting LiDAR data. Our new proposed algorithm consists of three key steps. Firstly, the noise photons were filtered out using a noise removal algorithm based on localized statistical analysis. Secondly, we classified the signal photons into canopy photons and ground photons by conducting a series of operations, including elevation frequency histogram building, empirical mode decomposition (EMD), and progressive densification. At the same time, we also identified the top of canopy (TOC) photons from canopy photons by percentile statistics method. Thereafter, the ground and TOC surfaces were generated from ground photons and TOC photons by cubic spline interpolation, respectively. Finally, the ground elevation and vegetation height were estimated by retrieved ground and TOC surfaces. The results indicate that the noise removal algorithm is effective in identifying background noise and preserving signal photons. The retrieved ground elevation is more accurate than the retrieved vegetation height, and the results of nighttime data are better than those of the corresponding daytime data. Specifically, the root-mean-square error (RMSE) values of ground elevation estimates range from 2.25 to 6.45 m for daytime data and 2.03 to 6.03 m for nighttime data. The RMSE values of vegetation height estimates range from 4.63 to 8.92 m for daytime data and 4.55 to 8.65 m for nighttime data. Our algorithm performs better than the previous algorithms in estimating ground elevation and vegetation height due to lower RMSE values. Additionally, the results also illuminate that the photon classification algorithm effectively reduces the negative effects of slope and vegetation coverage. Overall, our paper provides an effective solution for estimating ground elevation and vegetation height from micro-pulse photon-counting LiDAR data.

The inherent biases of different long-term gridded surface soil moisture (SSM) products, unconstrained by the in situ observations, implies different spatio-temporal patterns. In this study, the Random Forest (RF) model was trained to predict SSM from relevant land surface feature variables (i.e., land surface temperature, vegetation indices, soil texture, and geographical information) and precipitation, based on the in situ soil moisture data of the International Soil Moisture Network (ISMN.). The results of the RF model show an RMSE of 0.05 m3 m−3 and a correlation coefficient of 0.9. The calculated impurity-based feature importance indicates that the Antecedent Precipitation Index affects most of the predicted soil moisture. The geographical coordinates also significantly influence the prediction (i.e., RMSE was reduced to 0.03 m3 m−3 after considering geographical coordinates), followed by land surface temperature, vegetation indices, and soil texture. The spatio-temporal pattern of RF predicted SSM was compared with the European Space Agency Climate Change Initiative (ESA-CCI) soil moisture product, using both time-longitude and latitude diagrams. The results indicate that the RF SSM captures the spatial distribution and the daily, seasonal, and annual variabilities globally.

Wetland vegetation is an important component of wetland ecosystems and plays a crucial role in the ecological functions of wetland environments. Accurate distribution mapping and dynamic change monitoring of vegetation are essential for wetland conservation and restoration. The development of unoccupied aerial vehicles (UAVs) provides an efficient and economic platform for wetland vegetation classification. In this study, we evaluated the feasibility of RGB imagery obtained from the DJI Mavic Pro for wetland vegetation classification at the species level, with a specific application to Honghu, which is listed as a wetland of international importance. A total of ten object-based image analysis (OBIA) scenarios were designed to assess the contribution of five machine learning algorithms to the classification accuracy, including Bayes, K-nearest neighbor (KNN), support vector machine (SVM), decision tree (DT), and random forest (RF), multi-feature combinations and feature selection implemented by the recursive feature elimination algorithm (RFE). The overall accuracy and kappa coefficient were compared to determine the optimal classification method. The main results are as follows: (1) RF showed the best performance among the five machine learning algorithms, with an overall accuracy of 89.76% and kappa coefficient of 0.88 when using 53 features (including spectral features (RGB bands), height information, vegetation indices, texture features, and geometric features) for wetland vegetation classification. (2) The RF model constructed by only spectral features showed poor classification results, with an overall accuracy of 73.66% and kappa coefficient of 0.70. By adding height information, VIs, texture features, and geometric features to construct the RF model layer by layer, the overall accuracy was improved by 8.78%, 3.41%, 2.93%, and 0.98%, respectively, demonstrating the importance of multi-feature combinations. (3) The contribution of different types of features to the RF model was not equal, and the height information was the most important for wetland vegetation classification, followed by the vegetation indices. (4) The RFE algorithm effectively reduced the number of original features from 53 to 36, generating an optimal feature subset for wetland vegetation classification. The RF based on the feature selection result of RFE (RF-RFE) had the best performance in ten scenarios, and provided an overall accuracy of 90.73%, which was 0.97% higher than the RF without feature selection. The results illustrate that the combination of UAV-based RGB imagery and the OBIA approach provides a straightforward, yet powerful, approach for high-precision wetland vegetation classification at the species level, in spite of limited spectral information. Compared with satellite data or UAVs equipped with other types of sensors, UAVs with RGB cameras are more cost efficient and convenient for wetland vegetation monitoring and mapping.

National Aeronautics and Space Administration’s (NASA’s) Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) provides rich insights over the Earth’s surface through elevation data collected by its Advanced Topographic Laser Altimeter System (ATLAS) since its launch in September 2018. While this mission is primarily aimed at capturing ice measurements, ICESat-2 also provides data over vegetated areas, offering the capability to gain insights into ecosystem structure and the potential to contribute to the sustainable management of forests. This study involved an examination of the utility of ICESat-2 for estimating forest aboveground biomass (AGB). The objectives of this study were to: (1) investigate the use of canopy metrics for estimating AGB, using data extracted from an ICESat-2 transect over forests in south-east Texas; (2) compare the accuracy for estimating AGB using data from the strong beam and weak beam; and (3) upscale predicted AGB estimates using variables from Landsat multispectral imagery and land cover and canopy cover maps, to generate a 30 m spatial resolution AGB map. Methods previously developed with simulated ICESat-2 data over Sam Houston National Forest (SHNF) in southeast Texas were adapted using actual data from an adjacent ICESat-2 transect over similar vegetation conditions. Custom noise filtering and photon classification algorithms were applied to ICESat-2’s geolocated photon data (ATL03) for one beam pair, consisting of a strong and weak beam, and canopy height estimates were retrieved. Canopy height parameters were extracted from 100 m segments in the along-track direction for estimating AGB, using regression analysis. ICESat-2-derived AGB estimates were then extrapolated to develop a 30 m AGB map for the study area, using vegetation indices from Landsat 8 Operational Land Imager (OLI), National Land Cover Database (NLCD) landcover and canopy cover, with random forests (RF). The AGB estimation models used few canopy parameters and suggest the possibility for applying well-developed methods for modeling AGB with airborne light detection and ranging (lidar) data, using processed ICESat-2 data. The final regression model achieved a R2 and root mean square error (RMSE) value of 0.62 and 24.63 Mg/ha for estimating AGB and RF model evaluation with a separate test set yielded a R2 of 0.58 and RMSE of 23.89 Mg/ha. Findings provide an initial look at the ability of ICESat-2 to estimate AGB and serve as a basis for further upscaling efforts.

Spatially continuous estimates of forest aboveground biomass (AGB) are essential to supporting the sustainable management of forest ecosystems and providing invaluable information for quantifying and monitoring terrestrial carbon stocks. The launch of the Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) on September 15th, 2018 offers an unparalleled opportunity to assess AGB at large scales using along-track samples that will be provided during its three-year mission. The main goal of this study was to investigate deep learning (DL) neural networks for mapping AGB with ICESat-2, using simulated photon-counting lidar (PCL)-estimated AGB for daytime, nighttime, and no noise scenarios, Landsat imagery, canopy cover, and land cover maps. The study was carried out in Sam Houston National Forest located in south-east Texas, using a simulated PCL-estimated AGB along two years of planned ICESat-2 profiles. The primary tasks were to investigate and determine neural network architecture, examine the hyper-parameter settings, and subsequently generate wall-to-wall AGB maps. A first set of models were developed using vegetation indices calculated from single-date Landsat imagery, canopy cover, and land cover, and a second set of models were generated using metrics from one year of Landsat imagery with canopy cover and land cover maps. To compare the effectiveness of final models, comparisons with Random Forests (RF) models were made. The deep neural network (DNN) models achieved R2 values of 0.42, 0.49, and 0.50 for the daytime, nighttime, and no noise scenarios respectively. With the extended dataset containing metrics calculated from Landsat images acquired on different dates, substantial improvements in model performance for all data scenarios were noted. The R2 values increased to 0.64, 0.66, and 0.67 for the daytime, nighttime, and no noise scenarios. Comparisons with Random forest (RF) prediction models highlighted similar results, with the same R2 and root mean square error (RMSE) range (15–16 Mg/ha) for daytime and nighttime scenarios. Findings suggest that there is potential for mapping AGB using a combinatory approach with ICESat-2 and Landsat-derived products with DL.

The ancient topography of the Tibetan Plateau and its role in biotic evolution are still poorly understood, mostly due to a lack of fossil evidence. Our discovery of ∼47-Mya plant fossils from a present elevation of 4,850 m in central Tibet, diminishes, significantly, that lack of knowledge. The fossils represent a humid subtropical vegetation and some of the 70 different plant forms show affinity to Early-Middle Eocene floras in both North America and Europe. Using leaf architecture, we calculate that the forest grew at ∼1,500-m elevation within an east–west trending valley under a monsoonal climate. Our findings highlight the complexity of Tibet’s ancient landscape and emphasize the importance of Tibet in the history of global biodiversity.