An official website of the United States government
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
Crary Ice Rise region dh/dt estimates
v2.0 of this dataset includes: All surface-elevation change (dh/dt) data from ICESat, CryoSat-2, and ICESat-2 altimetry missions necessary to reproduce figures and analysis from Verboncoeur et al. (2024) ('*dhdt_smb'); a file containing x,y positions of the ad-hoc reference tracks formed around ICESat ground tracks ('xy_is_masked.csv'); a folder containing delineated boundaries used in analysis ('SHAPES.zip'); folders containing raw subsetted ICESat data ('IS_data.zip') and a folder containing CryoCloud scripts for downloading ICESat-2 data ('IS2_processing_cryocloud.zip')
more »
« less
- Award ID(s):
- 2049302
- PAR ID:
- 10596051
- Publisher / Repository:
- Zenodo
- Date Published:
- Format(s):
- Medium: X
- Right(s):
- Creative Commons Attribution 4.0 International
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
NASA’s ICESat-2 has been providing sea ice freeboard measurements across the polar regions since October 2018. In spite of the outstanding spatial resolution and precision of ICESat-2, the spatial sparsity of the data can be a critical issue for sea ice monitoring. This study employs a geostatistical approach (i.e., ordinary kriging) to characterize the spatial autocorrelation of the ICESat-2 freeboard measurements (ATL10) to estimate weekly freeboard variations in 2019 for the entire Ross Sea area, including where ICESat-2 tracks are not directly available. Three variogram models (exponential, Gaussian, and spherical) are compared in this study. According to the cross-validation results, the kriging-estimated freeboards show correlation coefficients of 0.56–0.57, root mean square error (RMSE) of ~0.12 m, and mean absolute error (MAE) of ~0.07 m with the actual ATL10 freeboard measurements. In addition, the estimated errors of the kriging interpolation are low in autumn and high in winter to spring, and low in southern regions and high in northern regions of the Ross Sea. The effective ranges of the variograms are 5–10 km and the results from the three variogram models do not show significant differences with each other. The southwest (SW) sector of the Ross Sea shows low and consistent freeboard over the entire year because of the frequent opening of wide polynya areas generating new ice in this sector. However, the southeast (SE) sector shows large variations in freeboard, which demonstrates the advection of thick multiyear ice from the Amundsen Sea into the Ross Sea. Thus, this kriging-based interpolation of ICESat-2 freeboard can be used in the future to estimate accurate sea ice production over the Ross Sea by incorporating other remote sensing data.more » « 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
-
The raw data for the associated manuscript is organized here into three categories: 1) relating to the measurement and analysis of the native recluse spiders loop junctions, 2) raw images found in the figures throughout the manuscript, and 3) relating to the experiments testing the effect that junction angle has on the strength of two intersecting tapes. It is recommended to browse the data files in Tree mode, which will make the files appear in folders reflecting this organization. 1) Loxosceles Loop Junction Images and Analysis The folder titled, SEM Raw Images, has all of the scanning electron microscopy (SEM) images taken of the native recluse loop junctions. Some images are close-ups of individual junctions and others take a broader perspective (macro) of many loop junctions in series. Where possible several close-up images of the individual junctions are accompanied with a macro image. These images were imported into ImageJ where the junction angle was measured. The measurements for all 41 loop junctions observed are in the folder titled, Raw Data Files in the file titled, Loxosceles Loop Junction Angle Measurements.txt. The folder titled, Raw Data Files contains, in addition to the angle measurements, the raw data for analyzing the strength of individual loop junctions. The data is in native MATLAB data format. These datasets include the complete tensile data and the cross-sectional area data for each spiders silk. The MATLAB code titled, Figure_2A_2B_code, processes the raw tensile data from the natural recluse spiders loop junctions. This data is plotted as two representative curves in Figure 2A and as a complete set as a histogram in Figure 2B. The MATLAB code titled, Figure_7_code, processes and plots the loop junction data found in, Loxosceles Loop Junction Angle Measurements.txt and executed the model of a random set of recluse loops. This code can be executed to generate Figure 7. The folder titled, Raw Data Files, must be open in MATLAB to run this code! This code uses the MATLAB function, areacalculation, to calculate the junction area for a given junction angle. 2) Raw Images This folder is organized by the respective figure in the manuscript where each image can be found. Additional metadata for each image can be found accompanying each image. 3) Tensile Data and Analysis This folder contains all of the raw tensile data for all tape-tape junction experiments conducted. All of the tensile data is in the folder titled, Raw Data Test Files. Within this folder is a .txt file for each sample tested. The file names are critical to the figure codes working properly because they contain the information for the junction angle and iterations. The file names are in the format year-month-day_trialnumber_junctionangle.txt. Also in the Raw Data Test Files folder are two functions used within some of the figure codes: fbfill and areacalculation. These functions will be used in the figure codes to properly analyze the data. To generate any figure using the MATLAB code in this folder, first open the code in MATLAB. Then within MATLAB, open the folder Raw Data Test Files. Only with this folder open in MATLAB will the code be able to find the correct raw data .txt files. The rest of the contents of this folder are MATLAB codes for specific figures in the manuscript. The only exception to this is the code titled, surfaceenergy_code, which is executed to calculate the phenomenological surface energy for the tapes used in these experiments.more » « less
