The objective of this paper is to demonstrate a new capability to detect faint stratospheric aerosols in atmospheric lidar data from NASA’s Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations (CALIPSO) Mission (12 June 2006 30 June 2023) at high resolution. Faint aerosols from wildfires, distant volcanic eruptions, and Asian and Saharan dust storms play important roles in the Earth-Atmosphere system, as they provide essential air quality and pollution, airline safety, climate radiative forcing, monitoring of fires and volcanic eruptions, environmental safety, and weather forecasting. Tenuous aerosols are recorded in the CALIOP atmospheric lidar data from the CALIPSO mission, but especially extremely faint stratospheric aerosol layers often escape detection and classification partly or entirely in the current CALIOP data analysis scheme. To solve this problem, we introduce a new algorithm for detection and height determination of atmospheric layers, the Density-Dimension Algorithm for CALIOP data analysis (CALIOP-DDA). The DDA is an auto-adaptive algorithm that builds on concepts from artificial intelligence and spatial statistics. Core steps are calculation of a density field and application of a threshold function for signal-noise separation. Stratospheric aerosol detection is aided by the tropopause split concept. The CALIOP-DDA facilitates detection of extremely faint stratospheric aerosols from various sources, including distant wildfires and volcanic eruptions, in nighttime and day-time CALIOP data, even in presence of complex types of other cloud and aerosol layers across a large range of optical thicknesses, while retaining the full 334m along-track, 30m height resolution, without creating false positives. CALIOPDDA results are evaluated by comparison to layer heights derived from airborne validation experiments, conducted using the Cloud Physics Lidar (CPL). In conclusion, the CALIOP-DDA holds promise as the algorithmic basis for a future improved, high-resolution CALIPSO data product.

The objective is to demonstrate that spatially variable ice-surface roughness is an important, and so-far overlooked, component of melt processes in the Greenland Ice Sheet (GrIS), linking the analysis of spatial ice-surface roughness (SISR) derived from satellite laser altimeter data with a new regional energy balance model. Specific results are: (1) SISR can be calculated from ICESat Geoscience Laser Altimeter System (GLAS) data for the GrIS. (2) Seasonal and interannual changes of SISR are reflected in satellite-altimetry-based roughness maps, demonstrated for GLAS data. (3) Current regional climate models largely utilize roughness as a constant value in space and time. Here, we develop a new Regional Energy Balance Model (REBM) that is sensitive to SISR, using the same physical principles as the Integrated Forecasting System (IFS), constrained by climate fields of the Regional Atmospheric Climate Model (RACMO2). (4) A control study is carried out to ensure that REBM works correctly. (5) Melt energy of the GrIS calculated by REBM using SISR results in much higher melt values than predicted by regional climate models (RACMO2, IFS), which underestimate melting. (6) Application of the approach using REBM to the GrIS, driven by seasonally and interannually variable SISR, highlights regional and temporal differences in sensible heat flux and thus melt differences. In summary, SISR explains some of the discrepancy between observed and modeled melting in the GrIS. This study serves as a proof of concept for an approach that establishes mathematical and physical concepts for linking satellite altimeter measurements, SISR and energy balance modeling.

A document by Camden Opfer. Click on the document to view its contents.

The objectives of this paper are to investigate the tradeoffs between a physically constrained neural network and a deep, convolutional neural network and to design a combined ML approach (“VarioCNN”). Our solution is provided in the framework of a cyberinfrastructure that includes a newly designed ML software, GEOCLASS-image, modern high-resolution satellite image datasets (Maxar WorldView data) and instructions/descriptions that may facilitate solving similar spatial classification problems. Combining the advantages of the physically-driven connectionist-geostatistical classification method with those of an efficient CNN, VarioCNN, provides a means for rapid and efficient extraction of complex geophysical information from submeter resolution satellite imagery. A retraining loop overcomes the difficulties of creating a labeled training data set.Computational analyses and developments are centered on a specific, but generalizable, geophysical problem: The classification of crevasse types that form during the surge of a glacier system. A surge is a glacial catastrophe, an acceleration of a glacier to typically 100-200 times its normal velocity, which for a marine-terminating glacier leads to sudden and substantial mass transfer from the cryosphere to the oceans, contributing significantly to sea-level-rise. The sudden and rapid acceleration characteristic of a surge results in formation of crevasses, whose spatial characteristics provide informants on the ice-dynamic processes that occur during the surge. GEOCLASS-image is applied to study the current (2016-2024) surge in the Negribreen Glacier System, Svalbard. The geophysical result is a description of the structural evolution and expansion of the surge, based on crevasse types that capture ice deformation in 6 simplified classes.

A document by Thomas Trantow. Click on the document to view its contents.

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. 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 reported on the ICESat-2 standard sea-ice products because of the low resolution of the products, in which only a single surface is determined. The objective of this paper is to introduce a mathematical algorithm that facilitates automated detection of melt ponds in ICESat-2 ATLAS data, retrieval of two surface heights, pond surface and bottom, and measurements of depth and width of melt ponds. With the Advanced Topographic Laser Altimeter System (ATLAS), ICESat-2 carries the first space-borne multi-beam micro-pulse 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 auto-adaptive algorithm that solves the problem of pond detection near  the 0.7m nominal alongtrack resolution of ATLAS data, utilizing the radial basis function for calculation of a density field and a threshold function that automatically adapts to changes in background, apparent surface reflectance and some instrument effects. The DDA-bifurcate-seaice is applied to large ICESat-2 data sets from the 2019 and 2020 melt seasons in the multi-year Arctic sea-ice region. Results are evaluated by comparison to those from a manually forced algorithm.

Observations reveal end of summer Arctic sea ice extent is declining at an accelerating rate. Model projections underestimate this decline and continue to have a broad spread in forecasted September sea ice extent. This suggests some important summer processes, such as melt pond formation and evolution, may not be properly represented in current models. Melt ponds form on the sea ice surface as snow melts, and pools in low lying areas on the sea ice surface. The evolution of the ponds depends on snow depth, ice thickness, and surface conditions. Melt water may spread across a level surface, or be confined to depressions between sea ice ridges. Ponds decrease the albedo of the surface and enhance the positive ice albedo feedback, accelerating further melt. Until recently, Arctic-wide observations of individual melt ponds were not available. ICESat-2, a photon counting laser altimeter launched in 2018, provides high resolution detail of sea ice and snow topography due to its unique combination of a small footprint (~12 m) and high-resolution along-track sampling (0.7 m). The green laser (532 nm) is able to penetrate water, enabling melt pond depth measurements. We have developed methods to track the melt pond surface and bathymetry in ICESat-2 data to determine melt pond depth. We also track melt pond evolution through application of a sea ice classification algorithm to 10 m resolution Sentinel-2 imagery. The combination of these two datasets allows for an evolving, three-dimensional view of the melting sea ice surface. We focus on the evolution of summer melt on multiyear ice in the Central Arctic north of Greenland and Canada in 2020. Our findings are put in context of existing literature on melt pond depth, volume, and evolution. We also discuss our results in relation to the melt pond fraction north of the Fram Strait, where we expect different ice conditions in the vicinity of the 2020 MOSAiC field studies. Observational data products comprising melt pond fraction and pond depth are being developed for public distribution. These products may be of interest to those studying under-ice light and biology, as well as modelers who are interested in understanding the evolution of melt pond parameters for model initialization and validation.

The Advanced Topographic Laser Altimeter System (ATLAS) was launched aboard the Ice Cloud and land-Elevation Satellite-2 (ICESat-2) satellite in September 2018. ATLAS is a single wavelength (532 nm) lidar system designed to acquire high resolution measurements of the earth’s surface while also obtaining atmospheric backscatter from molecules, clouds, and aerosols. Because ATLAS is optimized for altimetry, the atmospheric data acquired is unique in many respects and requires non-standard analysis techniques. For example, the high repetition rate laser limits the vertical extent of the profiles to just 14 km and causes atmospheric scattering from above 15 km to be added to the scattering in the lower 0 -14 km profile. In addition, the limited vertical range of the acquired profiles renders it difficult to compute the magnitude of the solar background and hinders the application of standard calibration techniques. Despite these limitations, techniques have been developed to successfully produce data products that have value to the atmospheric community for cloud and aerosol research and are currently available at the National Snow and Ice Data Center (NSIDC). In this paper we describe the ICESat-2 atmospheric channel and the methods used to process the ATLAS raw photon count data to obtain calibrated backscatter and higher level products such as layer heights and type, blowing snow, column optical depth and apparent surface reflectance.

During the Arctic summer season, snow atop the sea ice melts and pools into low-lying areas on the surface. These melt ponds reduce surface albedo and increase solar absorption in the Arctic Ocean. Throughout the summer, melt ponds grow, drain, and connect, through a complex drainage system. Current melt pond schemes in sea ice models, such as the level-ice scheme in the Los Alamos Sea Ice Model (CICE), rely on a linear relationship between pond depth and fraction to predict the evolution of pond growth as the snow and sea ice melt. Although the inclusion of melt ponds in models has been shown to improve forecasts of end-of-summer sea ice extent, observations of melt pond depth and fraction guiding these models are from SHEBA, a spatially-limited field campaign which occurred over 20 years ago. Until recently, melt ponds characteristics have been difficult to resolve from spaceborne platforms due to their small size (10s - 100s m in diameter), and indistinguishable radiometric similarity to open water. Here we show that new, high-resolution laser altimetry measurements from ICESat-2 (IS2), combined with coincident high-resolution satellite imagery, provides a three-dimensional view of the melting sea ice cover. IS2, launched in September 2018, has now observed two summer melt seasons in the Arctic. IS2 operates at 532 nm, a wavelength that penetrates low turbidity water, and can therefore be used to capture the bathymetry of shallow water features. Building on previous work, we demonstrate IS2’s ability to detect and measure melt ponds on multiyear sea ice. We validate the existence of melt ponds with high resolution (10 m) visible imagery from the Sentinel-2 (S2) MultiSpectral Instrument. We apply the “density dimension algorithm – bifurcate” (DDA-bifurcate), an auto-adaptive algorithm utilizing data aggregation with the ability to track two surfaces, as well as a second algorithm that tracks melt pond surface and bottom, to derive melt pond depth for dozens of melt ponds in 2019 and 2020. Applying a sea ice surface classification algorithm to S2 imagery, we are able to determine melt pond fraction. We compare our findings of coincident melt pond fraction and depth with the melt pond parameterization used in the level-ice scheme in CICE. We discuss our results in the context of the existing literature on pond depth and volume.

Glacial acceleration is the largest source of uncertainty in sea-level-rise assessment, according to the Intergovernmental Panel on Climate Change. Of the different types of glacial acceleration, surging is the least understood. In this paper, we demonstrate how a combination of automated algorithms dedicated to analysis of two entirely different observation types - satellite altimetry from NASA’s ICESat-2 and satellite imagery from Planet SkySat - can aid in advancing glaciology, utilizing state-of-the art remote sensing /Earth observation technology. NASA’s Ice, Cloud and land Elevation Satellite ICESat-2, launched 15~September~2018, carries the first space-borne multi-beam micro-pulse photon counting laser altimeter system, the Advanced Topographic Laser Altimeter System (ATLAS). ATLAS observations are collected in three pairs of weak and strong beams with 0.7m nominal along-track spacing (under clear-sky conditions). The recording of the observations as a photon-point cloud requires a dedicated algorithm for identification of signal photons and determination of surface heights. As a solution, we developed the density-dimension algorithm for ice surfaces, the DDA-ice. ATLAS data analyzed with the DDA-ice allow determination of heights over heavily crevassed ice surfaces, which are characteristics of accelerating glaciers. The study presented here builds on a special multi-component data set, obtained through synoptic observations of an Arctic glacier system during surge (Negribreen, Svalbard): Airborne altimeter and image data collected during our ICESat-2 validation campaign, and SkySat image data from a special acquisition collected as part of NASA’s Commercial Smallsat Data Acquisitions Pilot program. These are complemented by WorldView (Maxar) and ESA Sentinel-1 data. With a spatial resolution of 0.7-0.86m, SkySat data and WorldView lend themselves to automated classification of crevasse types. Altogether, we obtain a characterization in 3 dimensions that allows discrimination of ice-surface types from surging glaciers (Negribreen) and continuously fast-moving and accelerating glaciers (Jakobshavn Isbrae) based on morphological characteristics.