ICESat-2 Applications for Investigating Emerging Volcanoes

Submarine volcanism in shallow waters (<100 m), particularly in remote settings, is difficult to monitor quantitatively and, in the rare formation of islands, it is challenging to understand the rapid-paced erosion. However, these newly erupted volcanic islands become observable to airborne and/or satellite remote sensing instruments. NASA’s ICESat-2 satellite laser altimeter, combined with visible imagery (optical and microwave), provide a novel method of evaluating the elevation characteristics of newly emerged volcanoes and their subaerial eruption products. Niijima Fukutoku-Okanoba (NFO) is a submarine volcano 1300 km south of Tokyo (Ogasawara Archipelago of Japan) that periodically breaches the ocean surface to create new islands that are subsequently eroded. The recent eruption in August 2021 is a rare opportunity to investigate this island evolution using high-resolution satellite datasets with geodetic-quality ICESat-2 altimetry. Lansdat-8 and Planet imagery provide a qualitative analysis of the exposed volcanic deposits, while ICESat-2 products provide elevation profiles necessary to quantify the physical surface structures. This investigation determines an innovative application for ICESat-2 data in evaluating newly emerged islands and how the combination of satellite remote sensing (visible and lidar) to investigate these short-lived volcanic features can improve our understanding of the volcanic island system in ways not previously possible.

Exploring the Relationship between Forest Canopy Height and Canopy Density from Spaceborne LiDAR Observations

Forest structure is a useful proxy for carbon stocks, ecosystem function and species diversity, but it is not well characterised globally. However, Earth observing sensors, operating in various modes, can provide information on different components of forests enabling improved understanding of their structure and variations thereof. The Ice, Cloud and Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS), providing LiDAR footprints from 2003 to 2009 with close to global coverage, can be used to capture elements of forest structure. Here, we evaluate a simple allometric model that relates global forest canopy height (RH100) and canopy density measurements to explain spatial patterns of forest structural properties. The GLA14 data product (version 34) was applied across subdivisions of the World Wildlife Federation ecoregions and their statistical properties were investigated. The allometric model was found to correspond to the ICESat GLAS metrics (median mean squared error, MSE: 0.028; inter-quartile range of MSE: 0.022–0.035). The relationship between canopy height and density was found to vary across biomes, realms and ecoregions, with denser forest regions displaying a greater increase in canopy density values with canopy height, compared to sparser or temperate forests. Furthermore, the single parameter of the allometric model corresponded with the maximum canopy density and maximum height values across the globe. The combination of the single parameter of the allometric model, maximum canopy density and maximum canopy height values have potential application in frameworks that target the retrieval of above-ground biomass and can inform on both species and niche diversity, highlighting areas for conservation, and potentially enabling the characterisation of biophysical drivers of forest structure.

Light detection and ranging (LIDAR) laser altimeter for the Martian Moons Exploration (MMX) spacecraft

An altimeter is a critical instrument in planetary missions, for both safe operations and science activities. We present required specifications and link budget calculations for light detection and ranging (LIDAR) onboard the Martian Moons Exploration (MMX) spacecraft. During the mission phase, this LIDAR will continuously measure the distance between the spacecraft and its target. The time-series distance provides important diagnostic information for safe spacecraft operations and important information for geomorphological studies. Because MMX is a sample return mission, its LIDAR must accommodate physical disturbances on the Martian satellite surface. This resulted in changes to the optical system design. Graphical abstract

Supraglacial lake bathymetry automatically derived from ICESat-2 constraining lake depth estimates from multi-source satellite imagery

We introduce an algorithm (Watta) which automatically calculates supraglacial lake bathymetry and detects potential ice layers along tracks of the ICESat-2 (Ice, Cloud, and Land Elevation Satellite) laser altimeter. Watta uses photon heights estimated by the ICESat-2 ATL03 product and extracts supraglacial lake surface, bottom, and depth corrected for refraction and (sub-)surface ice cover in addition to producing surface heights at the native resolution of the ATL03 photon cloud. These measurements are used to constrain empirical estimates of lake depth from satellite imagery, which were thus far dependent on sparse sets of in situ measurements for calibration. Imagery sources include Landsat 8 Operational Land Imager (OLI), Sentinel-2, and high-resolution Planet Labs PlanetScope and SkySat data, used here for the first time to calculate supraglacial lake depths. The Watta algorithm was developed and tested using a set of 46 lakes near Sermeq Kujalleq (Jakobshavn) glacier in western Greenland, and we use multiple imagery sources (available for 45 of these lakes) to assess the use of the red vs. green band to extrapolate depths along a profile to full lake volumes. We use Watta-derived estimates in conjunction with high-resolution imagery from both satellite-based sources (tasked over the season) and nearly simultaneous Operation IceBridge CAMBOT (Continuous Airborne Mapping By Optical Translator) imagery (on a single airborne flight) for a focused study of the drainage of a single lake over the 2019 melt season. Our results suggest that the use of multiple imagery sources (both publicly available and commercial), in combination with altimetry-based depths, can move towards capturing the evolution of supraglacial hydrology at improved spatial and temporal scales.

A Method of Extracting High-Accuracy Elevation Control Points from ICESat-2 Altimetry Data

Due to its high ranging accuracy, spaceborne laser altimetry technology can improve the accuracy of satellite stereo mapping without ground control points. In the past, full-waveform ICE, CLOUD, and Land Elevation Satellite (ICESat) laser altimeter data have been used as one of the main data sources for global elevation control. As a second-generation satellite, ICESat-2 is equipped with an altimeter using photon counting mode. This can further improve the application capability for stereo mapping because of the six laser beams with high along-track repetition frequency, which can provide more detailed ground contour descriptions. Previous studies have addressed how to extract high-accuracy elevation control points from ICESat data. However, these methods cannot be directly applied to ICESat-2 data because of the different modes of the laser altimeters. Therefore, in this paper, we propose a method using comprehensive evaluation labels that can extract high-accuracy elevation control points that meet the different level elevation accuracy requirements for large scale mapping from the ICESat-2 land-vegetation along-track product. The method was verified using two airborne lidar data sets. In flat, hilly, and mountainous areas, by using our method to extract the terrain elevation, the root-mean-square error of elevation control points decrease from 1.249–2.094 m, 2.237–3.225 m, and 2.791–4.822 m to 0.262–0.429 m, 0.484–0.596 m, and 0.611–1.003 m, respectively. The results show that the extraction elevations meet the required accuracy for large scale mapping.

A Wind–Wave-Dependent Sea Spray Volume Flux Model Based on Field Experiments

Sea spray can contribute significantly to the exchanges of heat and momentum across the air–sea interface. However, while critical, sea spray physics are typically not included in operational atmospheric and oceanic models due to large uncertainties in their parameterizations. In large part, this is because of the scarcity of in-situ sea spray observations which prevent rigorous validation of existing sea spray models. Moreover, while sea spray is critically produced through the fundamental interactions between wind and waves, traditionally, sea spray models are parameterized in terms of wind properties only. In this study, we present novel in-situ observations of sea spray derived from a laser altimeter through the adoption of the Beer–Lambert law. Observations of sea spray cover a broad range of wind and wave properties and are used to develop a wind–wave-dependent sea spray volume flux model. Improved performance of the model is observed when wave properties are included, in contrast to a parameterization based on wind properties alone. The novel in-situ sea spray observations and the predictive model derived here are consistent with the classic spray model in both trend and magnitude. Our model and novel observations provide opportunities to improve the prediction of air–sea fluxes in operational weather forecasting models.

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.1918 $\pm$ 2.385~m along straight paths over generally crevassed terrain, which implies a precision of approximately 2.4~m for ICESat-2 validation experiments. (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. 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.

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.1918 $\pm$ 2.385~m along straight paths over generally crevassed terrain, which implies a precision of approximately 2.4~m for ICESat-2 validation experiments. (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. 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.