skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Observing the evolution of summer melt on multiyear sea ice with ICESat-2 and Sentinel-2
We investigate sea ice conditions during the 2020 melt season, when warm air temperature anomalies in spring led to early melt onset, an extended melt season, and the second-lowest September minimum Arctic ice extent observed. We focus on the region of the most persistent ice cover and examine melt pond depth retrieved from Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) using two distinct algorithms in concert with a time series of melt pond fraction and ice concentration derived from Sentinel-2 imagery to obtain insights about the melting ice surface in three dimensions. We find the melt pond fraction derived from Sentinel-2 in the study region increased rapidly in June, with the mean melt pond fraction peaking at 16 % ± 6 % on 24 June 2020, followed by a slow decrease to 8 % ± 6 % by 3 July, and remained below 10 % for the remainder of the season through 15 September. Sea ice concentration was consistently high (>95 %) at the beginning of the melt season until 4 July, and as floes disintegrated, it decreased to a minimum of 70 % on 30 July and then became more variable, ranging from 75 % to 90 % for the remainder of the melt season. Pond depth increased steadily from a median depth of 0.40 m ± 0.17 m in early June and peaked at 0.97 m ± 0.51 m on 16 July, even as melt pond fraction had already started to decrease. Our results demonstrate that by combining high-resolution passive and active remote sensing we now have the ability to track evolving melt conditions and observe changes in the sea ice cover throughout the summer season.  more » « less
Award ID(s):
1942356 2325430 1835256
PAR ID:
10473351
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Copernicus Publications
Date Published:
Journal Name:
The Cryosphere
Volume:
17
Issue:
9
ISSN:
1994-0424
Page Range / eLocation ID:
3695 to 3719
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Melt ponds forming on Arctic sea ice in summer significantly reduce the surface albedo and impact the heat and mass balance of the sea ice. Therefore, their areal coverage, which can undergo rapid change, is crucial to monitor. We present a revised method to extract melt pond fraction (MPF) from Sentinel‐2 satellite imagery, which is evaluated by MPF products from higher‐resolution satellite and helicopter‐borne imagery. The analysis of melt pond evolution during the MOSAiC campaign in summer 2020, shows a split of the Central Observatory (CO) into a level ice and a highly deformed ice part, the latter of which exhibits exceptional early melt pond formation compared to the vicinity. Average CO MPFs are 17% before and 23% after the major drainage. Arctic‐wide analysis of MPF for years 2017–2021 shows a consistent seasonal cycle in all regions and years. 
    more » « less
  2. Sea ice growth and decay are critical processes in the Arctic climate system, but comprehensive observations are very sparse. We analyzed data from 23 sea ice mass balance buoys (IMBs) deployed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020 to investigate the seasonality and timing of sea ice thermodynamic mass balance in the Arctic Transpolar Drift. The data reveal four stages of the ice season: (I) onset of ice basal freezing, mid-October to November; (II) rapid ice growth, December–March; (III) slow ice growth, April–May; and (IV) melting, June onward. Ice basal growth ranged from 0.64 to 1.38 m at a rate of 0.004–0.006 m d–1, depending mainly on initial ice thickness. Compared to a buoy deployed close to the MOSAiC setup site in September 2012, total ice growth was about twice as high, due to the relatively thin initial ice thickness at the MOSAiC sites. Ice growth from the top, caused by surface flooding and subsequent snow-ice formation, was observed at two sites and likely linked to dynamic processes. Snow reached a maximum depth of 0.25 ± 0.08 m by May 2, 2020, and had melted completely by June 25, 2020. The relatively early onset of ice basal melt on June 7 (±10 d), 2019, can be partly attributed to the unusually rapid advection of the MOSAiC floes towards Fram Strait. The oceanic heat flux, calculated based on the heat balance at the ice bottom, was 2.8 ± 1.1 W m–2 in December–April, and increased gradually from May onward, reaching 10.0 ± 2.6 W m–2 by mid-June 2020. Subsequently, under-ice melt ponds formed at most sites in connection with increasing ice permeability. Our analysis provides crucial information on the Arctic sea ice mass balance for future studies related to MOSAiC and beyond. 
    more » « less
  3. This dataset contains data used in the paper: Revisiting the Last Ice Area Projections from a High-Resolution Global Earth System Model - Fol et al (2025). Results are organized in excel files or numpy arrays with the dataset name, variable and ensemble member (for simulations) in the name of the file. See below for more information on what variables are included in the files and their structure.  CESM_HR :  CESM_HR_SIAFluxes - per ensemble member: Timeseries of monthly SIA flux per gate. CESM_HR_fluxesCAADiv.npy, _fluxesQEIDiv.npy, _fuxesQEIDivMeltSeason.npy : Timeseries of annual divergence over the Queen Elizabeth Islands and the Southern Canadian Arctic Archipelago derived from monthly SIA fluxes at the entry and exit gates. CESM_HR_ThicknessDistribution.xlsx : Thickness distribution for the LIA-N, QEI, and CAA-S computed from the simulated thickness distribution (aicen001, aicen002, aicen003, aicen004, aicen005). CESM_HR_tendencies - per ensemble member- per region (LIA-N, QEI, CAA-S). Timeseries of melt season integrated thermodynamic, dynamic (advection and ridging terms) sea ice area loss. CESH_HR_sitPanArctic - per ensemble member: Timeseries of pan-Arctic mean may sea ice thickness. CESM_HR_sieSept and CESM_HR_sieMarch - per ensemble member- per region (LIA-N, QEI, CAA-S) and pan-Arctic : Timeseries of March or September sea ice extent, sea ice area. CESM_HR_sic - per ensemble member- per region (LIA-N, QEI, CAA-S): Timeseries of mean sea ice concentration for grid cells having more than 15% of SIC (no open water).  CESM_HR_meltSeason - per ensemble member- per region (LIA-N, QEI, CAA-S): Timeseries of annual freeze and melt onset dates allowing the definition of the melt season based on the thermodynamic sea ice area tendency crossing 0.  CESM_HR_mean_mayThickness.npy and CESM_HR_meanseptconc.npy: Results for map of the mean september sea ice concentration and may ice thickness for 1981-2000, 2001-2020, 2021-2040 and 2041-2060. CESM_LR :  CESM_LR_sieSept - per region (LIA-N, QEI, CAA-S) and pan-Arctic: Timeseries of September sea ice extent and sea ice area. CESM_LR_sitPanArctic.xlsx: Timeseries of pan-Arctic mean May sea ice thickness. CESM_LR_tendencies - per region (LIA-N, QEI, CAA-S). Timeseries of melt season integrated thermodynamic, dynamic (advection and ridging terms) sea ice area loss. CESM_LR_meltSeason - per region (LIA-N, QEI, CAA-S): Timeseries of annual freeze and melt onset dates allowing the definition of the melt season based on the thermodynamic sea ice area tendency crossing 0. CESM2_LE:  CESM2_LE_CAA_sept, _LIAN_sept, _QEI_sept, panArctic: Mean September sea ice extent and sea ice area per region. There is one excel tab per ensemble member in each file. CESM2_LE_pan_Arctic_hi_may: Pan-Arctic mean May sea ice thickness. There is one excel tab per ensemble member in each file. PIOMAS:  PIOMAS_panArctic_hi.xlsx: Timeseries of mean may sea ice thickness. PIOMAS_mean_1981_2000_mean_mayThickness.npy : Results for map of the mean May ice thickness for 1981-2000 and 2001-2020. Observations: CIS_marchSept_1982_1990_sie - per region (QEI and CAA-S) :  Timeseries of March and September mean sea ice extent and area. CIS ice charts do not fully cover the LIA-N.  NSDICCDR_1979_2023_sia - per region (LIA-N, QEI, CAA-S) and pan-Arctic : Timeseries of March and September mean sea ice extent and area. NSDICCDR_1981_2000_mean_septConc.npy: Results for map of the mean September sea ice concentration for 1981-2000 and 2001-2020. NSDICCDR_1981_2000_LIAN_sic.xlsx: Timeseries of monthly mean sea ice concentration in the LIA-N. NSDICCDR_1981_2000_QEI_sic.xlsx: Timeseries of monthly mean sea ice concentration in the QEI. These results are derived from the following datasets:  Ensemble members 1 and 3 of simulations from the high-resolution Community Earth System Model version 1.3 (CESM1.3-HR) produced for the International Laboratory for High-Resolution Earth System Prediction (iHESP) by the Qingdao National Laboratory for Marine Science and Technology (QNLM), Texas A&M University (TAMU), and the U.S. National Center for Atmospheric Research (NCAR). The lower resolution simulation is also used (CESM1.3-LR). (Chang et al., 2020; Zhang et al., 2020). The 100-ensemble members Community Earth System Model (version 2) Large Ensemble (CESM2-LE) (Danabasoglu et al., 2020). Satellite-derived monthly mean SIA fluxes through entry and exit gates of the CAA and Nares Strait (Howell et a., 2019; 2021; 2023; 2024; Smedsrud et al., 2017; Kwok, 2006). The National Snow and Ice Data Center (NSIDC) Climate Data Record (CDR) (version 4) sea ice concentration, stored on a 25 x 25 km polar stereographic grid centered on the North Pole from 1979 to 2023 (Meier et al., 2021). The gridded version of the regional Canadian Ice Service (CIS) Digital Archive ice charts from the Eastern and Western Arctic regions (Tivy et al., 2011). The Pan-Arctic Ice Ocean Modeling and Assimilation System (PIOMAS) assimilated sea ice concentration and derived ice thickness distribution estimates in the Arctic from 1978 to 2022 (Zhang et al., 2000).   
    more » « less
  4. 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
  5. As part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC), the HELiX uncrewed aircraft system (UAS) was deployed over the sea ice in the central Arctic Ocean during summer 2020. Albedo measurements were obtained with stabilized pyranometers, and melt pond fraction was calculated from orthomosaic imagery from a surface-imaging multispectral camera. This study analyzed HELiX flight data to provide insights on the temporal and spatial evolution of albedo and melt pond fraction of the MOSAiC floe during the melt season as it drifted south through Fram Strait. The surface albedo distributions showed peak values changing from high albedo (0.55–0.6) to lower values (0.3) as the season advanced. Inspired by methods developed for satellite data, an algorithm was established to retrieve melt pond fraction from the orthomosaic images. We demonstrate that the near-surface observations of melt pond fraction were highly dependent on sample area, offering insight into the influence of subgrid scale features and spatial heterogeneity in satellite observations. Vertical observations conducted with the HELiX were used to quantify the influence of melt pond scales on observed surface albedo as a function of sensor footprint. These scaling results were used to link surface-based measurements collected during MOSAiC to broader-scale satellite data to investigate the influence of surface features on observed albedo. Albedo values blend underlying features within the sensor footprint, as determined by the melt pond size and concentration. This study framed the downscaling (upscaling) problem related to the airborne (surface) observations of surface albedo across a variety of spatial scales. 
    more » « less