As climate warms and the transition from a perennial to a seasonal Arctic sea-ice cover is imminent, understanding melt ponding is central to understanding changes in the new Arctic. National Aeronautics and Space Administration (NASA)’s Ice, Cloud and land Elevation Satellite (ICESat-2) has the capacity to provide measurements and monitoring of the onset of melt in the Arctic and on melt progression. Yet ponds are currently not identified on the ICESat-2 standard sea-ice products, in which only a single surface is determined. The objective of this article is to introduce a mathematical algorithm that facilitates automated detection of melt ponds in the ICESat-2 Advanced Topographic Laser Altimeter System (ATLAS) data, retrieval of two surface heights, pond surface and bottom, and measurements of depth and width of melt ponds. With ATLAS, ICESat-2 carries the first spaceborne multibeam micropulse photon-counting laser altimeter system, operating at 532-nm frequency. ATLAS data are recorded as clouds of discrete photon points. The Density-Dimension Algorithm for bifurcating sea-ice reflectors (DDA-bifurcate-seaice) is an autoadaptive algorithm that solves the problem of pond detection near the 0.7-m nominal along-track spacing of ATLAS data, utilizing the radial basis function for calculation of a density field and a threshold function that automatically adapts to changes in the background, apparent surface reflectance, and some instrument effects. The DDA-bifurcate-seaice is applied to large ICESat-2 datasets from the 2019 and 2020 melt seasons in the multiyear Arctic sea-ice region. Results are evaluated by comparison with those from a manually forced algorithm.
more »
« less
Airborne Validation of ICESat-2 ATLAS Data over Crevassed Surfaces and Other Complex Glacial Environments: Results from Experiments of Laser Altimeter and Kinematic GPS Data Collection from a Helicopter over a Surging Arctic Glacier (Negribreen, Svalbard)
The topic of this paper is the airborne evaluation of ICESat-2 Advanced Topographic Laser Altimeter System (ATLAS) measurement capabilities and surface-height-determination over crevassed glacial terrain, with a focus on the geodetical accuracy of geophysical data collected from a helicopter. To obtain surface heights over crevassed and otherwise complex ice surface, ICESat-2 data are analyzed using the density-dimension algorithm for ice surfaces (DDA-ice), which yields surface heights at the nominal 0.7 m along-track spacing of ATLAS data. As the result of an ongoing surge, Negribreen, Svalbard, provided an ideal situation for the validation objectives in 2018 and 2019, because many different crevasse types and morphologically complex ice surfaces existed in close proximity. Airborne geophysical data, including laser altimeter data (profilometer data at 905 nm frequency), differential Global Positioning System (GPS), Inertial Measurement Unit (IMU) data, on-board-time-lapse imagery and photographs, were collected during two campaigns in summers of 2018 and 2019. Airborne experiment setup, geodetical correction and data processing steps are described here. To date, there is relatively little knowledge of the geodetical accuracy that can be obtained from kinematic data collection from a helicopter. Our study finds that (1) Kinematic GPS data collection with correction in post-processing yields higher accuracies than Real-Time-Kinematic (RTK) data collection. (2) Processing of only the rover data using the Natural Resources Canada Spatial Reference System Precise Point Positioning (CSRS-PPP) software is sufficiently accurate for the sub-satellite validation purpose. (3) Distances between ICESat-2 ground tracks and airborne ground tracks were generally better than 25 m, while distance between predicted and actual ICESat-2 ground track was on the order of 9 m, which allows direct comparison of ice-surface heights and spatial statistical characteristics of crevasses from the satellite and airborne measurements. (4) The Lasertech Universal Laser System (ULS), operated at up to 300 m above ground level, yields full return frequency (400 Hz) and 0.06–0.08 m on-ice along-track spacing of height measurements. (5) Cross-over differences of airborne laser altimeter data are −0.172 ± 2.564 m along straight paths, which implies a precision of approximately 2.6 m for ICESat-2 validation experiments in crevassed terrain. (6) In summary, the comparatively light-weight experiment setup of a suite of small survey equipment mounted on a Eurocopter (Helicopter AS-350) and kinematic GPS data analyzed in post-processing using CSRS-PPP leads to high accuracy repeats of the ICESat-2 tracks. The technical results (1)–(6) indicate that direct comparison of ice-surface heights and crevasse depths from the ICESat-2 and airborne laser altimeter data is warranted. Numerical evaluation of height comparisons utilizes spatial surface roughness measures. The final result of the validation is that ICESat-2 ATLAS data, analyzed with the DDA-ice, facilitate surface-height determination over crevassed terrain, in good agreement with airborne data, including spatial characteristics, such as surface roughness, crevasse spacing and depth, which are key informants on the deformation and dynamics of a glacier during surge.
more »
« less
- PAR ID:
- 10346536
- Date Published:
- Journal Name:
- Remote Sensing
- Volume:
- 14
- Issue:
- 5
- ISSN:
- 2072-4292
- Page Range / eLocation ID:
- 1185
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The ATLAS sensor onboard the ICESat-2 satellite is a photon-counting lidar (PCL) with a primary mission to map Earth's ice sheets. A secondary goal of the mission is to provide vegetation and terrain elevations, which are essential for calculating the planet's biomass carbon reserves. A drawback of ATLAS is that the sensor does not provide reliable terrain height estimates in dense, high-closure forests because only a few photons reach the ground through the canopy and return to the detector. This low penetration translates into lower accuracy for the resultant terrain model. Tropical forest measurements with ATLAS have an additional problem estimating top of canopy because of frequent atmospheric phenomena such as fog and low clouds that can be misinterpreted as top of the canopy. To alleviate these issues, we propose using a ConvPoint neural network for 3D point clouds and high-density airborne lidar as training data to classify vegetation and terrain returns from ATLAS. The semantic segmentation network provides excellent results and could be used in parallel with the current ATL08 noise filtering algorithms, especially in areas with dense vegetation. We use high-density airborne lidar data acquired along ICESat-2 transects in Central American forests as a ground reference for training the neural network to distinguish between noise photons and photons lying between the terrain and the top of the canopy. Each photon event receives a label (noise or signal) in the test phase, providing automated noise-filtering of the ATL03 data. The terrain and top of canopy elevations are subsequently aggregated in 100 m segments using a series of iterative smoothing filters. We demonstrate improved estimates for both terrain and top of canopy elevations compared to the ATL08 100 m segment estimates. The neural network (NN) noise filtering reliably eliminated outlier top of canopy estimates caused by low clouds, and aggregated root mean square error (RMSE) decreased from 7.7 m for ATL08 to 3.7 m for NN prediction (18 test profiles aggregated). For terrain elevations, RMSE decreased from 5.2 m for ATL08 to 3.3 m for the NN prediction, compared to airborne lidar reference profiles.ICESat-2LidarPoint cloudNoise filteringmore » « less
-
The unique measurement capabilities of ICESat-2 allow high spatiotemporal resolution of complex ice-dynamic processes that occur during a surge. Detailed and precise mapping of height changes on surge glaciers has previously escaped observations from space due to limited resolution of space-borne altimeter data and the surface characteristics of glaciers during surge such as heavy crevassing. This makes geophysical interpretation of deformation and assessment of mass transfer difficult. In this paper, we present an approach that facilitates analysis of the evolution of geophysical processes during a surge, including height changes, crevassing, mass transfer and roughness evolution. We utilize all data from 2 years of ICESat-2 observations collected during the mature phase of the Negribreen Glacier System (NGS) surge in 2019 and 2020. The progression of the NGS surge has resulted in large-scale elevation changes and wide-spread crevassing making it an ideal case study to demonstrate ICESat-2 measurement capabilities, which are maximized when coupled with the Density Dimension Algorithm for Ice (DDA-ice). Results show the expansion of the surge in upper Negribreen which demonstrates the unique ability of ICESat-2/DDA-ice to measure a rapidly changing surge glacier and provide the best estimates for cryospheric changes and their contributions to sea-level rise.more » « less
-
Accurate measurements of terrain elevation are crucial for many ecological applications. In this study, we sought to assess new global three-dimensional Earth observation data acquired by the spaceborne Light Detection and Ranging (LiDAR) missions Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) and Global Ecosystem Dynamics Investigation (GEDI). For this, we examined the “ATLAS/ICESat-2 L3A Land and Vegetation Height”, version 5 (20 × 14 m and 100 × 14 m segments) and the “GEDI Level 2A Footprint Elevation and Height Metrics”, version 2 (25 m circle). We conducted our analysis across four land cover classes (bare soil, herbaceous, forest, savanna), and six forest types (temperate broad-leaved, temperate needle-leaved, temperate mixed, tropical upland, tropical floodplain, and tropical secondary forest). For assessment of terrain elevation estimates from spaceborne LiDAR data we used high resolution airborne data. Our results indicate that both LiDAR missions provide accurate terrain elevation estimates across different land cover classes and forest types with mean error less than 1 m, except in tropical forests. However, using a GEDI algorithm with a lower signal end threshold (e.g., algorithm 5) can improve the accuracy of terrain elevation estimates for tropical upland forests. Specific environmental parameters (terrain slope, canopy height and canopy cover) and sensor parameters (GEDI degrade flags, terrain estimation algorithm; ICESat-2 number of terrain photons, terrain uncertainty) can be applied to improve the accuracy of ICESat-2 and GEDI-based terrain estimates. Although the goodness-of-fit statistics from the two spaceborne LiDARs are not directly comparable since they possess different footprint sizes (100 × 14 m segment or 20 × 14 m segment vs. 25 m circle), we observed similar trends on the impact of terrain slope, canopy cover and canopy height for both sensors. Terrain slope strongly impacts the accuracy of both ICESat-2 and GEDI terrain elevation estimates for both forested and non-forested areas. In the case of GEDI the impact of slope is, however, partly caused by horizontal geolocation error. Moreover, dense canopies (i.e., canopy cover higher than 90%) affect the accuracy of spaceborne LiDAR terrain estimates, while canopy height does not, when considering samples over flat terrains. Our analysis of the accuracy and precision of current versions of spaceborne LiDAR products for different vegetation types and environmental conditions provides insights on parameter selection and estimated uncertainty to inform users of these key global datasets.more » « less
-
Canopy heights and vertical profiles were analyzed for 12 airborne lidar tracks acquired over forests of the mid-Juruá region, Brazil. Canopy height models were classified at 1m resolution as floodplain, terrace, hillslope, or interfluvial flat; floodplains were further separated according to Horton-Strahler (HS) stream order. RH95 canopy heights, and vertical profiles at 1m intervals, were aggregated to 30m scale and compared with Copernicus DEM heights, using a DEM transform, the Relative Terrain Height (RTH). Median canopy height ranged from 15.4 m for the Juruá floodplain to 25.5 m for hillslopes; maximum canopy heights varied from 37.4 m to 60.0 m. A strong correlation between RTH and median canopy height (r = 0.75) was found for the Juruá floodplain tracks. Vertical profiles of Juruá floodplain tracks showed that the height above ground of maximum returns increased monotonically with RTH height. Our results clearly show the influence of floodplain topography on forest canopy structure.more » « less