NASA's ICESat-2 mission was launched in 2018, carrying a photon-counting laser altimeter, with a primary objective of measuring height changes across Earth's surface. ICESat-2 has provided measurements of ice surface height between 88º N and S, repeated four times per year, with high vertical accuracy and along-track spatial resolution. Its accuracy and coverage has enabled near-complete recovery of height changes across the ice sheets, capturing subtle changes in the interior, and rapid changes along the dynamic margins with steep slopes and the floating peripheral ice shelves. The ICESat-2 Science Team has developed a suite of algorithms that produce along-track and gridded land ice height products at various levels of processing, all freely available at the National Snow and Ice Data Center. Here, we describe three higher-level land-ice data products derived from ATL06 and their underlying algorithms: along-track height change (ATL11), digital elevation model (ATL14) and gridded surface height change (ATL15). We demonstrate the suitability of each data product for studying different ice sheet regions. We then show height changes for Greenland and Antarctica from ATL15 during the first 6 years of the ICESat-2 mission (October 2018-December 2024), illustrating how ICESat-2 measurements can distinguish the multi-year trends from seasonal fluctuations.

The process of laser light reflecting from surfaces made of scattering materials that do not strongly absorb at the wavelength of the laser can involve reflections from hundreds or thousands of individual grains, which can introduce delays in the time between light entering and leaving the surface. These time-of-flight biases depend on the grain size and density of the medium, and thus they can result in spatially and temporally varying surface height biases estimated from laser altimeters, such as NASA's ICESat-2 (Ice Cloud, and land Elevation Satellite-2) mission. Modeling suggests that ICESat-2 might experience a bias difference as large as 0.1–0.2 m between coarse-grained melting snow and fine-grained wintertime snow (Smith et al., 2018), which exceeds the mission's requirement to measure seasonal height differences to an accuracy better than 0.1 m (Markus et al., 2017). In this study, we investigate these biases using a model of subsurface scattering, laser altimetry measurements from NASA's ATM (Airborne Topographic Mapper) system, and grain size estimates based on optical imagery of the ice sheet. We demonstrate that distortions in the shapes of waveforms measured using ATM are related to the optical grain size of the surface estimated using optical reflectance measurements and show that they can be used to estimate an effective grain radius for the surface. Using this effective grain radius as a proxy for the severity of subsurface scattering, we use our model with grain size estimates from optical imagery to simulate corrections for biases in ICESat-2 data due to subsurface scattering and demonstrate that, on the basis of large-scale averages, the corrections calculated based on the satellite optical imagery match the biases in the data. This work demonstrates that waveform-based altimetry data can measure the optical properties of granular surfaces and that corrections based on optical grain size estimates can correct for subsurface-scattering biases in ICESat-2 data.

The Ice, Cloud, and land Elevation Satellite‐2 (ICESat‐2) is a NASA Earth observing satellite mission that provides global elevation measurements using the Advanced Topographic Laser Altimetry System (ATLAS). ICESat‐2 was launched in September 2018 and completed its prime mission of 3 years of on‐orbit science data collection in December 2021. ICESat‐2, as the successor mission to ICESat (2003–2009) (Schutz et al., 2005, https://doi.org/10.1029/2005gl024009 ), was designed to provide global elevation measurements of Earth's surfaces. Changes in elevation, such as those over glaciers, ice sheets and sea ice, are some of the most critical observations for characterizing and understanding Earth's dynamic processes and the response to climate variability. The overarching scientific goals of ICESat‐2 are associated primarily with the cryosphere, but from a space‐based platform, the altimeter measurements serve a wide range of science disciplines. Prior to launch during the early mission development phase, the Level 1 Science Requirements were established, which at the time were some of the most stringent metrics created for space‐based altimetry. These requirements were the primary drivers of both the instrument technology development and the mission operational strategies. Here, we evaluate each of the science requirements using the science data collected over the prime mission timeline of 3 years. We conclude from our analyses that the mission has successfully met each of the Level 1 Science Requirements. Further, we evaluate the onboard consumables (fuel and laser energy) and demonstrate that the satellite's operational lifetime could potentially last an additional ∼10 years.

Antarctic ice shelves buttress the flow of the ice sheet but are vulnerable to increased basal melting from contact with a warming ocean and increased mass loss from calving due to changing flow patterns. Channels and similar features at the bases of ice shelves have been linked to enhanced basal melting and observed to intersect the grounding zone, where the greatest melt rates are often observed. The ice shelf of Thwaites Glacier is especially vulnerable to basal melt and grounding zone retreat because the glacier has a retrograde bed leading to a deep trough below the grounded ice sheet. We use digital surface models from 2010-2022 to investigate the evolution of its ice-shelf channels, grounding zone position, and the interactions between them. We find that the highest sustained rates of grounding zone retreat (up to 0.7 kmyr-1) are associated with high basal melt rates (up to similar to 250 myr-1) and are found where ice-shelf channels intersect the grounding zone, especially atop steep local retrograde slopes where subglacial channel discharge is expected. We find no areas with sustained grounding zone advance, although some secular retreat was distal from ice-shelf channels. Pinpointing other locations with similar risk factors could focus assessments of vulnerability to grounding zone retreat.

This investigation quantifies the dynamic geometries of Antarctic active subglacial lakes by leveraging satellite altimetry data and an image processing algorithm. Using high-resolution altimetry satellite data from 2010 to 2023 and a computer-vision-based contour finding algorithm to delineate contours of ice-surface, height-change anomalies indicative of lake activity, we generated evolving outlines of 131 previously identified Antarctic active subglacial lakes. Evolving outlines provided a novel time series of wetted area to quantify changes in the bed substrate available for geochemical and microbial processes. Evolving outlines revealed previously unseen transitory lake lobes that refine preexisting time series of height change and ice-volume displacement, from which we can infer subglacial water flux. The findings underscore the value of employing evolving subglacial lake outlines to record active subglacial lake geometric changes across Antarctica and suggest broader glaciological applications, such as the identification of new active subglacial lake candidates.

Abstract. Antarctic ice shelves buttress the flow of the ice sheet, tempering sea level rise, but they are vulnerable to basal melting from contact with the ocean, as well as mass loss due to fracture and calving. Melt channels and similar features at the bases of ice shelves have been linked to enhanced basal melting and observed to intersect the grounding line, where the greatest melt rates are often observed. The ice shelf of Thwaites Glacier is especially vulnerable to melt and subsequent retreat of the grounding line because the glacier has an inland–sloping bed leading to a deep trough below the grounded ice sheet. We use digital surface models from 2010–2022 to investigate the evolution of ice–shelf basal channels and a proxy for the grounding line on the Thwaites Glacier ice shelf. We find that the highest sustained rates of grounding–line retreat (up to 0.7 km a-1) are associated with high melt rates (up to ~250 m a-1) near the intersections of basal channels with the grounding zone, steep local retrograde slopes, and where subglacial channel discharge is expected. Detailed observations of basal channels collocated with regions of grounding–line retreat will further elucidate the complicated processes occurring at the ice–ocean interface and hopefully lead to more accurate estimates of current and future ice–shelf melting and evolution.

Differential penetration of green laser light into snow and ice has long been considered a possible cause of range and thus elevation bias in laser altimeters. Over snow, ice, and water, green photons can penetrate the surface and experience multiple scattering events in the subsurface volume before being scattered back to the surface and subsequently the instrument's detector, therefore biasing the range of the measurement. Newly formed sea ice adjacent to open-water leads provides an opportunity to identify differential penetration without the need for an absolute reference surface or dual-color lidar data. We use co-located, coincident high-resolution natural-color imagery and airborne lidar data to identify surface and ice types and evaluate elevation differences between those surfaces. The lidar data reveals that apparent elevations of thin ice and finger-rafted thin ice can be several tens of centimeters below the water surface of surrounding leads, but not over dry snow. These lower elevations coincide with broadening of the laser pulse, suggesting that subsurface volume scattering is causing the pulse broadening and elevation shift. To complement our analysis of pulse shapes and help interpret the physical mechanism behind the observed elevation biases, we match the waveform shapes with a model of scattering of light in snow and ice that predicts the shape of lidar waveforms reflecting from snow and ice surfaces based on the shape of the transmitted pulse, the surface roughness, and the optical scattering properties of the medium. We parameterize the scattering in our model based on the scattering length Lscat, the mean distance a photon travels between isotropic scattering events. The largest scattering lengths are found for thin ice that exhibits the largest negative elevation biases, where scattering lengths of several centimeters allow photons to build up considerable range biases over multiple scattering events, indicating that biased elevations exist in lower-level Airborne Topographic Mapper (ATM) data products. Preliminary analysis of ICESat-2 ATL10 data shows that a similar relationship between subsurface elevations (restored negative freeboard) and “pulse width” is present in ICESat-2 data over sea ice, suggesting that biased elevations caused by differential penetration likely also exist in lower-level ICESat-2 data products. The spatial correlation of observed differential penetration in ATM data with surface and ice type suggests that elevation biases could also have a seasonal component, increasing the challenge of applying a simple bias correction.

Abstract. Surface-mass-balance (SMB) and firn-densification (FD) models are widely used in altimetry studies as a tool to separate atmospheric-driven from ice-dynamics-driven ice-sheet mass changes and to partition observed volume changes into ice-mass changes and firn-air-content changes. Until now, SMB models have been principally validated based on comparison with ice core and weather station data or comparison with widely separated flight radar-survey flight lines. Firn-densification models have been primarily validated based on their ability to match net densification over decades, as recorded in firn cores, and the short-term time-dependent component of densification has rarely been evaluated at all. The advent of systematic ice-sheet-wide repeated ice-surface-height measurements from ICESat-2 (the Ice Cloud, and land Elevation Satellite, 2) allows us to measure the net surface-height change of the Greenland ice sheet at quarterly resolution and compare the measured surface-height differences directly with those predicted by three FD–SMB models: MARv3.5.11 (Modèle Atmosphérique Régional version 3.5.11) and GSFCv1.1 and GSFCv1.2 (the Goddard Space Flight Center FD–SMB models version 1.1 and 1.2). By segregating the data by season and elevation, and based on the timing and magnitude of modelled processes in areas where we expect minimal ice-dynamics-driven height changes, we investigate the models' accuracy in predicting atmospherically driven height changes. We find that while all three models do well in predicting the large seasonal changes in the low-elevation parts of the ice sheet where melt rates are highest, two of the models (MARv3.5.11 and GSFCv1.1) systematically overpredict, by around a factor of 2, the magnitude of height changes in the high-elevation parts of the ice sheet, particularly those associated with melt events. This overprediction seems to be associated with the melt sensitivity of the models in the high-elevation part of the ice sheet. The third model, GSFCv1.2, which has an updated high-elevation melt parameterization, avoids this overprediction.

It contains the code, terminus data, flowlines, and hybridDEMs used in the paper titled "Compounding seasonal variations in outlet glacier dynamics revealed by high-resolution observations".

Ice losses from the Greenland and Antarctic ice sheets have accelerated since the 1990s, accounting for a significant increase in the global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume, and in Earth's gravity field. Between 1992 and 2020, the ice sheets contributed 21.0±1.9 mm to global mean sea level, with the rate of mass loss rising from 105 Gt yr−1 between 1992 and 1996 to 372 Gt yr−1 between 2016 and 2020. In Greenland, the rate of mass loss is 169±9 Gt yr−1 between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86 Gt yr−1 in 2017 to 444 Gt yr−1 in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (82±9 Gt yr−1) and, to a lesser extent, from the Antarctic Peninsula (13±5 Gt yr−1). East Antarctica remains close to a state of balance, with a small gain of 3±15 Gt yr−1, but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at https://doi.org/10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021).

SlideRule is a software framework for on demand science data processing on Cloud computing infrastructure that involves several innovative open-source software packages working in concert. SlideRule development has been guided by a fundamental change in the way Earth science data is distributed by institutions and interacted with by data users. Historically, institutions pre-generate datasets which are archived and made available via data centers. Data users download parts of those datasets from the data centers and perform local processing and analysis. SlideRule models an alternative approach that is centered around services . Institutions make science data available via web services, and data users integrate those web services directly into their existing data analysis systems. This new approach is enabled by co-locating low-level instrument data with scalable compute resources, and permits efficient generation of both global standard products and custom user-defined derived datasets. The development of SlideRule was led by the University of Washington in collaboration with NASA Goddard and has been funded by NASA’s ICESat-2 mission 1 . SlideRule is designed to be extensible in anticipation of upcoming Earth System Observatory (ESO) Missions, and takes advantage of ICESat-2 being one of the first NASA missions to have archival data products hosted in commercial cloud object storage (AWS S3 in us-west-2). The University of Washington’s deployment of SlideRule runs in AWS us-west-2 and consists of a single autoscaling cluster of EC2 instances that have access to NASA’s Cumulus datasets in S3, and NASA’s Common Metadata Repository (CMR). Users can access the system via an HTTP API, as well as a Python client.

Conversion of altimetry-derived ice-sheet volume change to mass requires an understanding of the evolution of the combined ice and air content within the firn column. In the absence of suitable techniques to observe the changes to the firn column across the entirety of an ice sheet, the firn column processes are typically modeled. Here, we present new simulations of firn processes over the Greenland and Antarctic ice sheets (GrIS and AIS) using the Community Firn Model and atmospheric reanalysis variables for more than four decades. A data set of more than 250 measured depth–density profiles from both ice sheets provides the basis of the calibration of the dry-snow densification scheme. The resulting scheme results in a reduction in the rate of densification, relative to a commonly used semi-empirical model, through a decreased dependence on the accumulation rate, a proxy for overburden stress. The 1980–2020 modeled firn column runoff, when combined with atmospheric variables from MERRA-2, generates realistic mean integrated surface mass balance values for the Greenland (+390 Gt yr−1) and Antarctic (+2612 Gt yr−1) ice sheets when compared to published model-ensemble means. We find that seasonal volume changes associated with firn air content are on average approximately 2.5 times larger than those associated with mass fluxes from surface processes for the AIS and 1.5 times larger for the GrIS; however, when averaged over multiple years, ice and air-volume fluctuations within the firn column are of comparable magnitudes. Between 1996 and 2019, the Greenland Ice Sheet lost nearly 5 % of its firn air content, indicating a reduction in the total meltwater retention capability. Nearly all (94 %) of the meltwater produced over the Antarctic Ice Sheet is retained within the firn column through infiltration and refreezing.

The National Aeronautics and Space Administration (NASA)’s Operation IceBridge (OIB) was a 13‐year (2009–2021) airborne mission to survey land and sea ice across the Arctic, Antarctic, and Alaska. Here, we review OIB’s goals, instruments, campaigns, key scientific results, and implications for future investigations of the cryosphere. OIB’s primary goal was to use airborne laser altimetry to bridge the gap in fine‐resolution elevation measurements of ice from space between the conclusion of NASA’s Ice, Cloud, and land Elevation Satellite (ICESat; 2003–2009) and its follow‐on, ICESat‐2 (launched 2018). Additional scientific requirements were intended to contextualize observed elevation changes using a multisensor suite of radar sounders, gravimeters, magnetometers, and cameras. Using 15 different aircraft, OIB conducted 968 science flights, of which 42% were repeat surveys of land ice, 42% were surveys of previously unmapped terrain across the Greenland and Antarctic ice sheets, Arctic ice caps, and Alaskan glaciers, and 16% were surveys of sea ice. The combination of an expansive instrument suite and breadth of surveys enabled numerous fundamental advances in our understanding of the Earth’s cryosphere. For land ice, OIB dramatically improved knowledge of interannual outlet‐glacier variability, ice‐sheet, and outlet‐glacier thicknesses, snowfall rates on ice sheets, fjord and sub‐ice‐shelf bathymetry, and ice‐sheet hydrology. Unanticipated discoveries included a reliable method for constraining the thickness within difficult‐to‐sound incised troughs beneath ice sheets, the extent of the firn aquifer within the Greenland Ice Sheet, the vulnerability of many Greenland and Antarctic outlet glaciers to ocean‐driven melting at their grounding zones, and the dominance of surface‐melt‐driven mass loss of Alaskan glaciers. For sea ice, OIB significantly advanced our understanding of spatiotemporal variability in sea ice freeboard and its snow cover, especially through combined analysis of fine‐resolution altimetry, visible imagery, and snow radar measurements of the overlying snow thickness. Such analyses led to the unanticipated discovery of an interdecadal decrease in snow thickness on Arctic sea ice and numerous opportunities to validate sea ice freeboards from satellite radar altimetry. While many of its data sets have yet to be fully explored, OIB’s scientific legacy has already demonstrated the value of sustained investment in reliable airborne platforms, airborne instrument development, interagency and international collaboration, and open and rapid data access to advance our understanding of Earth’s remote polar regions and their role in the Earth system.

A series of traverses has been conducted for validation of the National Aeronautics and Space Administration Ice, Cloud, and land Elevation Satellite 2 (ICESat‐2) on the flat interior of the Antarctic ice sheet. Global Navigation Satellite System data collected on three separate 88S Traverses intersect 20% of the ICESat‐2 reference ground tracks and have precisions of better than ±7 cm and biases of less than ∼4 cm. Data from these traverses were used to assess heights from ICESat‐2, CryoSat‐2, and Airborne Topographic Mapper (ATM). ICESat‐2 heights have better than ±3.3 cm bias and better than ±7.2 cm precision. ATM heights have better than 9.3 cm bias and better than ±9.6 cm precision. CryoSat‐2 heights have −38.9 cm of bias and ±47.3 cm precision. These best case results are from the flat ice‐sheet interior but provide a characterization of the quality of satellite and airborne altimetry.

Ice core archives are well suited for reconstructing rapid past climate changes at high latitudes. Despite this, few records currently exist from coastal Greenlandic ice caps due to their remote nature, limiting our long-term understanding of past maritime and coastal climate variability across this rapidly changing Arctic region. Here, we reconstruct regionally representative glacier surface mass balance and climate variability over the last two thousand years (~169–2015 ce) using an ice core collected from the Nuussuaq Peninsula, west Greenland. We find indications of abrupt regional hydroclimate shifts, including an up to 20% decrease in average snow accumulation during the transition from the Medieval Warm Period (950–1250 ce) to Little Ice Age (1450–1850 ce), followed by a subsequent >40% accumulation increase from the early 18th to late 20th centuries ce. These coastal changes are substantially larger than those previously reported from interior Greenland records. Moreover, we show that the strong relationship observed today between Arctic temperature rise and coastal ice cap decay contrasts with that of the last millennium, during which periods of warming led to snowfall-driven glacial growth. Taken together with modern observations, the ice core evidence could indicate a recent reversal in the response of west Greenland ice caps to climate change. Coastal west Greenland ice caps fluctuated strongly compared to the interior in response to rapid Common Era changes in snow accumulation, according to modelling of proxy records developed from a Nuussuaq Peninsula ice core covering the last 2,000 years.

A system of subglacial lakes drained on Thwaites Glacier from 2012–2014. To improve coverage for subsequent drainage events, we extended the elevation and ice-velocity time series on Thwaites Glacier through austral winter 2019. These new observations document a second drainage cycle in 2017/18 and identified two new lake systems located in the western tributaries of Thwaites and Haynes glaciers. In situ and satellite velocity observations show temporary < 3 % speed fluctuations associated with lake drainages. In agreement with previous studies, these observations suggest that active subglacial hydrology has little influence on thinning and retreat of Thwaites Glacier on decadal to centennial timescales.

Abstract. Conversion of altimetry-derived ice-sheet volume change to mass requires an understanding of the evolution of the combined ice and air content within the firn column. In the absence of suitable techniques to observe the changes to the firn column across the entirety of an ice sheet, the firn column processes are typically modelled. Here, we present new 40-year simulations of firn processes over the Greenland and Antarctic Ice Sheets using the Community Firn Model and atmospheric reanalysis variables. A dataset of more than 250 measured depth-density profiles from both ice sheets provides the basis of the calibration of the dry-snow densification scheme. The resulting scheme results in a reduction in the rate of densification, relative to a commonly used semi-empirical model, through a decreased dependence on the accumulation rate, a proxy for overburden stress. The modelled firn column runoff, when combined with atmospheric variables from MERRA-2, generates realistic mean integrated surface mass balance values for the Greenland (+361 Gt yr−1) and Antarctic (+2623 Gt yr−1) ice sheets when compared to published model-ensemble means. We find that seasonal volume changes associated with firn air content are approximately 3 times larger than those associated with surface mass balance; however, when averaged over multiple years, ice and air-volume fluctuations within the firn column are of comparable magnitudes. Between 1996 and 2019, the Greenland Ice Sheet lost more than 5 % of its firn air content indicating a reduction in the total meltwater retention capability. Nearly all (>98 %) of the meltwater produced over the Antarctic Ice Sheet is retained within the firn column through infiltration and refreezing.

The speed of Greenland's fastest glacier, Jakobshavn Isbrae, has varied substantially since its speedup in the late 1990s. Here we present observations of surface velocity, mélange rigidity, and surface elevation to examine its behaviour over the last decade. Consistent with earlier results, we find a pronounced cycle of summer speedup and thinning followed by winter slowdown and thickening. There were extended periods of rigid mélange in the winters of 2016-17 and 2017-18, concurrent with terminus advances ~6 km farther than in the several winters prior. These terminus advances to shallower depths caused slowdowns, leading to substantial thickening, as has been noted elsewhere. The extended periods of rigid mélange coincide well with a period of cooler waters in Disko Bay. Thus, along with the relative timing of the seasonal slowdown, our results suggest that the ocean's dominant influence on Jakobshavn Isbrae is through its effect on winter mélange rigidity, rather than summer submarine melting. The elevation time series also reveals that in summers when the area upstream of the terminus approaches flotation, large surface depressions can form, which eventually become the detachment points for major calving events. It appears that as elevations near flotation, basal crevasses can form, which initiates a necking process that forms the depressions. The elevation data also show that steep cliffs often evolve into short floating extensions, rather than collapsing catastrophically due to brittle failure. Finally, summer 2019 speeds were slightly faster than the prior two summers, leaving it unclear whether the slowdown is ending.