An Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation of the ICESat-2 Mission

The Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission has been selected by NASA as a Decadal Survey mission, to be launched in 2016. Mission objectives are to measure land ice elevation, sea ice freeboard/ thickness and changes in these variables and to collect measurements over vegetation that will facilitate determination of canopy height, with an accuracy that will allow prediction of future environmental changes and estimation of sea-level rise. The importance of the ICESat-2 project in estimation of biomass and carbon levels has increased substantially, following the recent cancellation of all other planned NASA missions with vegetation-surveying lidars. Two innovative components will characterize the ICESat-2 lidar: (1) Collection of elevation data by a multi-beam system and (2) application of micropulse lidar (photon counting) technology. A micropulse photon-counting altimeter yields clouds of discrete points, which result from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of returned points to reflectors of interest including canopy and ground in forested areas. The objective of this paper is to derive and validate an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2-type data. Data are based on airborne observations with a Sigma Space micropulse lidar and vary with respect to signal strength, noise levels, photon sampling options and other properties. A mathematical algorithm is developed, using spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors and geostatistical classification parameters and hyperparameters. Validation shows that the algorithm works very well and that ground and canopy elevation, and hence canopy height, can be expected to be observable with a high accuracy during the ICESat-2 mission. A result relevant for instrument design is that even the two weaker beam classes considered can be expected to yield useful results for vegetation measurements (93.01-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp9) and 72.85% - 98.68% for 0.48 msp (msp4)). Resampling options affect results more than noise levels. The algorithm derived here is generally applicable for analysis of micropulse lidar altimeter data collected over forested areas as well as other surfaces, including land ice, sea ice and land surfaces

Assessment of NASA Airborne Laser Altimetry Data Using Ground-Based GPS Data Near Summit Station, Greenland

A series of NASA airborne lidars have been used in support of satellite laser altimetry missions. These airborne laser altimeters have been deployed for satellite instrument development, for spaceborne data validation, and to bridge the data gap between satellite missions. We used data from ground-based Global Positioning System (GPS) surveys of an 11 km long track near Summit Station, Greenland, to assess the surface–elevation bias and measurement precision of three airborne laser altimeters including the Airborne Topographic Mapper (ATM), the Land, Vegetation, and Ice Sensor (LVIS), and the Multiple Altimeter Beam Experimental Lidar (MABEL). Ground-based GPS data from the monthly ground-based traverses, which commenced in 2006, allowed for the assessment of nine airborne lidar surveys associated with ATM and LVIS between 2007 and 2016. Surface–elevation biases for these altimeters – over the flat, ice-sheet interior – are less than 0.12 m, while assessments of measurement precision are 0.09 m or better. Ground-based GPS positions determined both with and without differential post-processing techniques provided internally consistent solutions. Results from the analyses of ground-based and airborne data provide validation strategy guidance for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) elevation and elevation-change data products

Analytical Modeling and Performance Assessment of Micropulse Photon-counting Lidar System

The melting of polar ice sheets and evidence of global warming continue to remain prominent research interests among scientists. To better understand global volumetric change of ice sheets, NASA intends to launch Ice, Cloud and land Elevation Satellite-2 (ICESat-2) in 2017. ICESat-2 employs a high frequency photon-counting laser altimeter, which will provide significantly greater spatial sampling. However, the combined effects of sub-beam complex surfaces, as well as system effects on returning photon distribution have not been systematically studied. To better understand the effects of various system attributes and to help improve the theory behind lidar sensing of complex surfaces, an analytical model using a first principles 3-D Monte Carlo approach is developed to predict system performance. Based on the latest ICESat-2 design, this analytical model simulates photons which propagate from the laser transmitter to the scene, and reflected to the detector model. A radiometric model is also applied in the synthetic scene. Such an approach allows the study of surface elevation retrieval accuracy for landscapes, as well as surface reflectivities. It was found that ICESat-2 will have a higher precision on a smoother surface, and a surface with smaller diffuse albedo will on average result in smaller bias. Furthermore, an adaptive density-based algorithm is developed to detect the surface returns without any geometrical knowledge. This proposed approach is implemented using the aforementioned simulated data set, as well as airborne laser altimeter measurement. Qualitative and quantitative results are presented to show that smaller laser footprint, smoother surface, and lower noise rate will improve accuracy of ground height estimation. Meanwhile, reasonable detection accuracy can also be achieved in estimating both ground and canopy returns for data generated using Digital Imaging and Remote Sensing Image Generation (DIRSIG) model. This proposed approach was found to be generally applicable for surface and canopy finding from photon-counting laser altimeter data

Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation for the ICESat-2 Mission

NASA's Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission is a decadal survey mission (2016 launch). The mission objectives are to measure land ice elevation, sea ice freeboard, and changes in these variables, as well as to collect measurements over vegetation to facilitate canopy height determination. Two innovative components will characterize the ICESat-2 lidar: 1) collection of elevation data by a multibeam system and 2) application of micropulse lidar (photon-counting) technology. A photon-counting altimeter yields clouds of discrete points, resulting from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of the returned points to reflectors of interest. The objective of this paper is to derive an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2 data, based on airborne observations with a Sigma Space micropulse lidar. The mathematical algorithm uses spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors, and geostatistical classification parameters and hyperparameters. Validation shows that ground and canopy elevation, and hence canopy height, can be expected to be observable with high accuracy by ICESat-2 for all expected beam energies considered for instrument design (93.01%-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp), and 72.85%-98.68% for 0.48 msp). The algorithm derived here is generally applicable for elevation determination from photoncounting lidar altimeter data collected over forested areas, land ice, sea ice, and land surfaces, as well as for cloud detection

Radiometric Assessment of ICESat-2 over Vegetated Surfaces

The ice, cloud, and land elevation satellite-2 (ICESat-2) is providing global elevation measurements to the science community. ICESat-2 measures the height of the Earth’s surface using a photon counting laser altimeter, ATLAS (advanced topographic laser altimetry system). As a photon counting system, the number of reflected photons per shot, or radiometry, is a function primarily of the transmitted laser energy, solar elevation, surface reflectance, and atmospheric scattering and attenuation. In this paper, we explore the relationship between detected scattering and attenuation in the atmosphere against the observed radiometry for three general forest types, as well as the radiometry as a function of day versus night. Through this analysis, we found that ATLAS strong beam radiometry exceeds the pre-launch design cases for boreal and tropical forests but underestimates the predicted radiometry over temperate forests by approximately half a photon. The weak beams, in contrast, exceed all pre-launch conditions by a factor of two to six over all forest types. We also observe that the signal radiometry from day acquisitions is lower than night acquisitions by 10% and 40% for the strong and weak beams, respectively. This research also found that the detection ratio between each beam-pair was lower than the predicted 4:1 values. This research also presents the concept of ICESat-2 radiometric profiles; these profiles provide a path for calculating vegetation structure. The results from this study are intended to be informative and perhaps serve as a benchmark for filtering or analysis of the ATL08 data products over vegetated surfaces

Advances in Waveform and Photon Counting Lidar Processing for Forest Vegetation Applications

Full waveform (FW) and photon counting LiDAR (PCL) data have garnered greater attention due to increasing data availability, a wealth of information they contain and promising prospects for large scale vegetation mapping. However, many factors such as complex processing steps and scarce non-proprietary tools preclude extensive and practical uses of these data for vegetation characterization. Therefore, the overall goal of this study is to develop algorithms to process FW and PCL data and to explore their potential in real-world applications. Study I explored classical waveform decomposition methods such as the Gaussian decomposition, Richardson–Lucy (RL) deconvolution and a newly introduced optimized Gold deconvolution to process FW LiDAR data. Results demonstrated the advantages of the deconvolution and decomposition method, and the three approaches generated satisfactory results, while the best performances varied when different criteria were used. Built upon Study I, Study II applied the Bayesian non-linear modeling concepts for waveform decomposition and quantified the propagation of error and uncertainty along the processing steps. The performance evaluation and uncertainty analysis at the parameter, derived point cloud and surface model levels showed that the Bayesian decomposition could enhance the credibility of decomposition results in a probabilistic sense to capture the true error of estimates and trace the uncertainty propagation along the processing steps. In study III, we exploited FW LiDAR data to classify tree species through integrating machine learning methods (the Random forests (RF) and Conditional inference forests (CF)) and Bayesian inference method. Results of classification accuracy highlighted that the Bayesian method was a superior alternative to machine learning methods, and rendered users with more confidence for interpreting and applying classification results to real-world tasks such as forest inventory. Study IV focused on developing a framework to derive terrain elevation and vegetation canopy height from test-bed sensor data and to pre-validate the capacity of the upcoming Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission. The methodology developed in this study illustrates plausible ways of processing the data that are structurally similar to expected ICESat-2 data and holds the potential to be a benchmark for further method adjustment once genuine ICESat-2 are available

Comparison of 16-Channel Laser Photoreceivers for Topographic Mapping

Topographic mapping lidar instruments must be able to detect extremely weak laser return signals from high altitudes including orbital distance. The signals have a wide dynamic range caused by the variability in atmospheric transmission and surface reflectance under a fast moving spacecraft. Ideally, lidar detectors should be able to detect laser signal return pulses at the single photon level and produce linear output for multiple photon events. Silicon avalanche photodiode (APO) detectors have been used in most space lidar receivers to date. Their sensitivity is typically hundreds of photons per pulse, and is limited by the quantum efficiency, APO gain noise, dark current, and preamplifier noise. NASA is pursuing three approaches for a 16-channel laser photoreceiver for use on the next generation direct-detection airborne and spacebome lidars. We present our measurement results and a comparison of their performance

The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2): Science Requirements, Concept, and Implementation

The Ice, Cloud, and land Elevation Satellite (ICESat) mission used laser altimetry measurements to determine changes in elevations of glaciers and ice sheets, as well as sea ice thickness distribution. These measurements have provided important information on the response of the cryosphere (Earths frozen surfaces) to changes in atmosphere and ocean condition. ICESat operated from 2003-2009 and provided repeat altimetry measurements not only to the cryosphere scientific community but also to the ocean, terrestrial and atmospheric scientific communities. The conclusive assessment of significant ongoing rapid changes in the Earths ice cover, in part supported by ICESat observations, has strengthened the need for sustained, high accuracy, repeat observations similar to what was provided by the ICESat mission. Following recommendations from the National Research Council for an ICESat follow-on mission, the ICESat-2 mission is now under development for planned launch in 2018. The primary scientific aims of the ICESat-2 mission are to continue measurements of sea ice freeboard and ice sheet elevation to determine their changes at scales from outlet glaciers to the entire ice sheet, and from 10s of meters to the entire polar oceans for sea ice freeboard. ICESat carried a single beam profiling laser altimeter that produced approximately 70 m diameter footprints on the surface of the Earth at approximately 150 m along-track intervals. In contrast, ICESat-2 will operate with three pairs of beams, each pair separated by about 3 km across-track with a pair spacing of 90 m. Each of the beams will have a nominal 17 m diameter footprint with an along-track sampling interval of 0.7 m. The differences in the ICESat-2 measurement concept are a result of overcoming some limitations associated with the approach used in the ICESat mission. The beam pair configuration of ICESat-2 allows for the determination of local cross-track slope, a significant factor in measuring elevation change for the outlet glaciers surrounding the Greenland and Antarctica coasts. The multiple beam pairs also provide improved spatial coverage. The dense spatial sampling eliminates along-track measurement gaps, and the small footprint diameter is especially useful for sea surface height measurements in the often narrow leads needed for sea ice freeboard and ice thickness retrievals. The ICESat-2 instrumentation concept uses a low energy 532 nm (green) laser in conjunction with single-photon sensitive detectors to measure range. Combining ICESat-2 data with altimetry data collected since the start of the ICESat mission in 2003, such as Operation IceBridge and ESAs CryoSat-2, will yield a 15+ year record of changes in ice sheet elevation and sea ice thickness. ICESat-2 will also provide information of mountain glacier and ice cap elevations changes, land and vegetation heights, inland water elevations, sea surface heights, and cloud layering and optical thickness

Potential of Forest Parameter Estimation Using Metrics from Photon Counting LiDAR Data in Howland Research Forest

ICESat-2 is the new generation of NASA’s ICESat (Ice, Cloud and land Elevation Satellite) mission launched in September 2018. We investigate the potential of forest parameter estimation using metrics from photon counting LiDAR data, using an integrated dataset including photon counting LiDAR data from SIMPL (the Slope Imaging Multi-polarization Photon-counting LiDAR), airborne small footprint LiDAR data from G-LiHT and a stem map in Howland Research Forest, USA. First, we propose a noise filtering method based on a local outlier factor (LOF) with elliptical search area to separate the ground and canopy surfaces from noise photons. Next, a co-registration technique based on moving profiling is applied between SIMPL and G-LiHT data to correct geolocation error. Then, we calculate height metrics from both SIMPL and G-LiHT. Finally, we investigate the relationship between the two sets of metrics, using a stem map from field measurement to validate the results. Results of the ground and canopy surface extraction show that our methods can detect the potential signal photons effectively from a quite high noise rate environment in relatively rough terrain. In addition, results from co-registration between SIMPL and G-LiHT data indicate that the moving profiling technique to correct the geolocation error between these two datasets achieves favorable results from both visual and statistical indicators validated by the stem map. Tree height retrieval using SIMPL showed error of less than 3 m. We find good consistency between the metrics derived from the photon counting LiDAR from SIMPL and airborne small footprint LiDAR from G-LiHT, especially for those metrics related to the mean tree height and forest fraction cover, with mean R 2 value of 0.54 and 0.6 respectively. The quantitative analyses and validation with field measurements prove that these metrics can describe the relevant forest parameters and contribute to possible operational products from ICESat-2