Precise measurements of Arctic sea ice mass balance are necessary to understand the rapidly changing sea ice cover and its representation in climate models. During the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, we made repeat point measurements of snow and ice thickness on primarily level first- and second-year ice (FYI, SYI) using ablation stakes and ice thickness gauges. This technique enabled us to distinguish surface and bottom (basal) melt and characterize the importance of oceanic versus atmospheric forcing. We also evaluated the time series of ice growth and melt in the context of other MOSAiC observations and historical mass balance observations from the Surface Heat Budget of the Arctic (SHEBA) campaign and the North Pole Environmental Observatory (NPEO). Despite similar freezing degree days, average ice growth at MOSAiC was greater on FYI (1.67 m) and SYI (1.23 m) than at SHEBA (1.45 m, 0.53 m), due in part to initially thinner ice and snow conditions on MOSAiC. Our estimates of effective snow thermal conductivity, which agree with SHEBA results and other MOSAiC observations, are unlikely to explain the difference. On MOSAiC, FYI grew more and faster than SYI, demonstrating a feedback loop that acts to increase ice production after multi-year ice loss. Surface melt on MOSAiC (mean of 0.50 m) was greater than at NPEO (0.18 m), with considerable spatial variability that correlated with surface albedo variability. Basal melt was relatively small (mean of 0.12 m), and higher than NPEO observations (0.07 m). Finally, we present observations showing that false bottoms reduced basal melt rates in some FYI cases, in agreement with other observations at MOSAiC. These detailed mass balance observations will allow further investigation into connections between the carefully observed surface energy budget, ocean heat fluxes, sea ice, and ecosystem at MOSAiC and during other campaigns.

Some of the largest climatic changes in the Arctic have been observed in Alaska and the surrounding marginal seas. Near-surface air temperature (T2m), precipitation (P), snowfall, and sea ice changes have been previously documented, often in disparate studies. Here, we provide an updated, long-term trend analysis (1957–2021; n = 65 years) of such parameters in ERA5, NOAA U.S. Climate Gridded Dataset (NClimGrid), NOAA National Centers for Environmental Information (NCEI) Alaska climate division, and composite sea ice products preceding the upcoming Fifth National Climate Assessment (NCA5) and other near-future climate reports. In the past half century, annual T2m has broadly increased across Alaska, and during winter, spring, and autumn on the North Slope and North Panhandle (T2m > 0.50°C decade^{-1). Precipitation has also increased across climate divisions and appears strongly interrelated with temperature–sea ice feedbacks on the North Slope, specifically with increased (decreased) open water (sea ice extent). Snowfall equivalent (SFE) has decreased in autumn and spring, perhaps aligned with a regime transition of snow to rain, while winter SFE has broadly increased across the state. Sea ice decline and melt-season lengthening also have a pronounced signal around Alaska, with the largest trends in these parameters found in the Beaufort Sea. Alaska’s climatic changes are also placed in context against regional and contiguous U.S. air temperature trends and show ∼50% greater warming in Alaska relative to the lower-48 states. Alaska T2m increases also exceed those of any contiguous U.S. subregion, positioning Alaska at the forefront of U.S. climate warming.}

The microstructure of the uppermost portions of a melting Arctic sea ice cover has a disproportionately large influence on how incident sunlight is reflected and absorbed in the ice/ocean system. The surface scattering layer (SSL) effectively backscatters solar radiation and keeps the surface albedo of melting ice relatively high compared to ice with the SSL manually removed. Measurements of albedo provide information on how incoming shortwave radiation is partitioned by the SSL and have been pivotal to improving climate model parameterizations. However, the relationship between the physical and optical properties of the SSL is still poorly constrained. Until now, radiative transfer models have been the only way to infer the microstructure of the SSL. During the MOSAiC expedition of 2019–2020, we took samples and, for the first time, directly measured the microstructure of the SSL on bare sea ice using X-ray micro-computed tomography. We show that the SSL has a highly anisotropic, coarse, and porous structure, with a small optical diameter and density at the surface, increasing with depth. As the melting surface ablates, the SSL regenerates, maintaining some aspects of its microstructure throughout the melt season. We used the microstructure measurements with a radiative transfer model to improve our understanding of the relationship between physical properties and optical properties at 850 nm wavelength. When the microstructure is used as model input, we see a 10–15% overestimation of the reflectance at 850 nm. This comparison suggests that either a) spatial variability at the meter scale is important for the two in situ optical measurements and therefore a larger sample size is needed to represent the microstructure or b) future work should investigate either i) using a ray-tracing approach instead of explicitly solving the radiative transfer equation or ii) using a more appropriate radiative transfer model.

Given the importance that melt ponds have on the energy balance of summer sea ice, there have been several efforts to develop pan-Arctic datasets using satellite data. Here we intercompare three melt pond data sets that rely on multi-frequency optical satellite data. Early in the melt season, the three data sets have similar spatial patterns in melt pond fraction, but this agreement weakens as the melt season progresses despite relatively high interannual correlations in pond fractions between the data products. Most of the data sets do not exhibit trends towards increased melt pond fractions from 2002 to 2011 despite overall Arctic warming and earlier melt onset. Further comparisons are made against higher resolution optical data to assess relative accuracy. These comparisons reveal the challenges in retrieving melt ponds from coarse resolution satellite data, and the need to better discriminate between leads, small open water areas and melt ponds. Finally, we assess melt pond data sets as a function of ice type and how well they correlate with surface albedo. As expected, melt pond fractions are negatively correlated with surface albedo, though the strength of the correlation varies across products and regions. Overall, first-year ice has larger melt pond fractions than multi-year ice.

The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 x 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

The Arctic is warming faster than anywhere on Earth, and with these warming temperatures, there is likely to be more precipitation falling as rain. This precipitation phase change will have profound impacts on the hydrologic cycle, energy balance, and snow and sea ice mass budgets. Here, we examine the number of rainfall days in the Arctic from three reanalysis: ERA-Interim, ERA5 and MERRA-2 over 1980-2016. We show that the number of rainfall days has increased over this period, predominantly in the autumn and in the North Atlantic and Peripheral Seas, and the length of the rain season has increased in all reanalyses. This is positively correlated to the number of days with above freezing air temperatures and a lengthening of the warm season. ERA-Interim produces significantly more rainfall days than other reanalyses and CloudSat observations, as well as significantly more rainfall when temperatures are below freezing. Investigation into the cloud microphysics schemes revealed that the scheme employed by ERA-Interim allowed for mixed-phase clouds to form rain at temperatures below freezing following a temperature-dependent phase partitioning function between 250K and 273K. This simple diagnostic treatment erroneously overestimates rain at temperatures below 273K and produces unrealistic rainfall compared to ERA5 and MERRA-2. This work highlights the importance of having accurate physics and improving microphysical schemes in models for simulating precipitation in the Arctic and the caution that is warranted for interpreting reanalysis trends.

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.

Some of the largest climatic changes in the Arctic have been observed in Alaska and surrounding marginal seas. Near-surface air temperature (T2m), precipitation (P), snowfall, and sea ice changes have been previously documented, often in disparate studies. Here we provide an updated, long-term trend analysis (1957–2021; n=65 years) of such parameters in ERA5, NOAA NClimGrid, NOAA NCEI Alaska climate division, and composite sea ice products preceding the upcoming Fifth National Climate Assessment (NCA5) and other near-future climate reports. In the past half century, annual T2m has broadly increased across Alaska, and during winter, spring, and autumn on the North Slope and North Panhandle (T2m>0.50°C decade^-1). P has also increased across climate divisions, and appears strongly interrelated with temperature-sea ice feedbacks on the North Slope, specifically with increased (decreased) open water (sea ice extent). Snowfall equivalent (SFE) has decreased in autumn and spring, perhaps aligned with a regime transition of snow to rain, while winter SFE has broadly increased across the state. Sea ice decline and melt season lengthening also have a pronounced signal around Alaska, with the largest trends in these parameters found in the Beaufort Sea. Alaska’s climatic changes are also placed in context against regional and contiguous U.S. air temperature trends, and show ~50% greater warming in Alaska relative to the lower-48 states. Alaska T2m increases also exceed those of any contiguous U.S. sub-region, positioning Alaska at the forefront of U.S. climate warming.

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.

Our understanding of Arctic sea ice and its wide-ranging influence is deeply rooted in observation. Advancing technologies have profoundly improved our ability to observe Arctic sea ice, document its processes and properties, and describe atmosphere-ice-ocean interactions with unprecedented detail. Yet, our progress toward better understanding the Arctic sea ice system is mired by the stark disparities between observations that tend to be siloed by method, scientific discipline, and application. This article presents a review and philosophical design for observing sea ice and accelerating our understanding of the Arctic sea ice system. We give a brief history of Arctic sea ice observations and showcase the 2018 melt season within the context of five observational themes: spatial heterogeneity, temporal variability, cross-disciplinary science, scalability, and retrieval uncertainty. We synthesize buoy data, optical imagery, satellite retrievals, and airborne measurements to demonstrate how disparate data sets can be woven together to transcend issues of observational scale. The results show that there are limitations to interpreting any single data set alone. However, many of these limitations can be surmounted by combining observations that cross spatial and temporal scales. We conclude the article with pathways toward enhanced coordination across observational platforms in order to: (1) optimize the scientific, operational, and community return on observational investments, and (2) facilitate a richer understanding of Arctic sea ice and its role in the climate system.

Arctic cyclones are a fundamental component of Arctic climate, influencing atmospheric heat and moisture transport into the region and surface energy, moisture, and momentum fluxes. Arctic cyclones can also drive changes in sea ice and energize ocean waves. Here we investigate a record low sea level pressure (SLP) Arctic cyclone which formed in East Greenland and tracked NE over the Barents and Kara seas between 21 and 27 January 2022. At its peak intensity on 24 January, the cyclone reached an estimated depth of 932.2 mb at 79.5°N 20°E. North of 70°N, this is the lowest SLP in the ERA-5 reanalysis over 1979 to present. The cyclone resulted in a record (over the period 1979–2022) weekly loss of regional sea ice area and surface wind speeds, and generated ocean waves exceeding 8 m that impinged on sea ice in the Barents sea, observed via satellite altimetry as large waves-in-sea ice up to 2 m in amplitude more than 100 km into the ice pack. Surface heat fluxes were strongly impacted by the cyclone, with record atmosphere-to-surface turbulent fluxes. However, the direct atmospheric thermodynamic impact on sea ice loss was modest, and the record sea ice changes were likely mainly driven by dynamical and/or ocean processes. While the storm was well predicted up to 8 days in advance, subsequent changes in sea ice cover were not, likely due to biases in the forecasts' sea ice initial conditions and missing physics in the forecast model such as wave-sea ice interaction.

The magnitude, spectral composition, and variability of the Arctic sea ice surface albedo are key to understanding and numerically simulating Earth’s shortwave energy budget. Spectral and broadband albedos of Arctic sea ice were spatially and temporally sampled by on-ice observers along individual survey lines throughout the sunlit season (April–September, 2020) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. The seasonal evolution of albedo for the MOSAiC year was constructed from spatially averaged broadband albedo values for each line. Specific locations were identified as representative of individual ice surface types, including accumulated dry snow, melting snow, bare and melting ice, melting and refreezing ponded ice, and sediment-laden ice. The area-averaged seasonal progression of total albedo recorded during MOSAiC showed remarkable similarity to that recorded 22 years prior on multiyear sea ice during the Surface Heat Budget of the Arctic Ocean (SHEBA) expedition. In accord with these and other previous field efforts, the spectral albedo of relatively thick, snow-free, melting sea ice shows invariance across location, decade, and ice type. In particular, the albedo of snow-free, melting seasonal ice was indistinguishable from that of snow-free, melting second-year ice, suggesting that the highly scattering surface layer that forms on sea ice during the summer is robust and stabilizing. In contrast, the albedo of ponded ice was observed to be highly variable at visible wavelengths. Notable temporal changes in albedo were documented during melt and freeze onset, formation and deepening of melt ponds, and during melt evolution of sediment-laden ice. While model simulations show considerable agreement with the observed seasonal albedo progression, disparities suggest the need to improve how the albedo of both ponded ice and thin, melting ice are simulated.

During the Arctic melt season, relatively fresh meltwater layers can accumulate under sea ice as a result of snow and ice melt, far from terrestrial freshwater inputs. Such under-ice meltwater layers, sometimes referred to as under-ice melt ponds, have been suggested to play a role in the summer sea ice mass balance both by isolating the sea ice from saltier water below, and by driving formation of 'false bottoms' below the sea ice. Such layers form at the interface of the fresher under-ice layer and the colder, saltier seawater below. During the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) expedition in the Central Arctic, we observed the presence of under-ice meltwater layers and false bottoms throughout July 2020 at primarily first-year ice locations. Here, we examine the distribution, prevalence, and drivers of under-ice ponds and the resulting false bottoms during this period. The average thickness of observed false bottoms and freshwater equivalent of under-ice meltwater layers was 0.08 m, with false bottom ice comprised of 74–87% FYI melt and 13–26% snow melt. Additionally, we explore these results using a 1D model to understand the role of dynamic influences on decoupling the ice from the seawater below. The model comparison suggests that the ice-ocean friction velocity was likely exceptionally low, with implications for air-ice-ocean momentum transfer. Overall, the prevalence of false bottoms was similar to or higher than noted during other observational campaigns, indicating that these features may in fact be common in the Arctic during the melt season. These results have implications for the broader ice-ocean system, as under-ice meltwater layers and false bottoms provide a source of ice growth during the melt season, potentially reduce fluxes between the ice and the ocean, isolate sea ice primary producers from pelagic nutrient sources, and may alter light transmission to the ocean below.

Accurate multidecadal radiative flux records are vital to understand Arctic amplification and constrain climate model uncertainties. Uncertainty in the NASA Clouds and the Earth's Radiant Energy System (CERES)-derived irradiances is larger over sea ice than any other surface type and comes from several sources. The year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic provides a rare opportunity to explore uncertainty in CERES-derived radiative fluxes. First, a systematic and statistically robust assessment of surface shortwave and longwave fluxes was conducted using in situ measurements from MOSAiC flux stations. The CERES Synoptic 1degree (SYN1deg) product overestimates the downwelling shortwave flux by +11.40 Wm^-2 and underestimates the upwelling shortwave flux by –15.70 Wm^-2 and downwelling longwave fluxes by –12.58 Wm^-2 at the surface during summer. In addition, large differences are found in the upwelling longwave flux when the surface approaches the melting point (approximately 0°C). The biases in downwelling shortwave and longwave fluxes suggest that the atmosphere represented in CERES is too optically thin. The large negative bias in upwelling shortwave flux can be attributed in large part to lower surface albedo (–0.15) in satellite footprint relative to surface sensors. Additionally, the results show that the spectral surface albedo used in SYN1deg overestimates albedo in visible and mid-infrared bands. A series of radiative transfer model perturbation experiments are performed to quantify the factors contributing to the differences. The CERES-MOSAiC broadband albedo differences (approximately 20 Wm^-2) explain a larger portion of the upwelling shortwave flux difference than the spectral albedo shape differences (approximately 3 Wm^-2). In addition, the differences between perturbation experiments using hourly and monthly MOSAiC surface albedo suggest that approximately 25% of the sea ice surface albedo variability is explained by factors not correlated with daily sea ice concentration variability. Biases in net shortwave and longwave flux can be reduced to less than half by adjusting both albedo and cloud inputs toward observed values. The results indicate that improvements in the surface albedo and cloud data would substantially reduce the uncertainty in the Arctic surface radiation budget derived from CERES data products.

Melt ponds on sea ice play an important role in the Arctic climate system. Their presence alters the partitioning of solar radiation: decreasing reflection, increasing absorption and transmission to the ice and ocean, and enhancing melt. The spatiotemporal properties of melt ponds thus modify ice albedo feedbacks and the mass balance of Arctic sea ice. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition presented a valuable opportunity to investigate the seasonal evolution of melt ponds through a rich array of atmosphere-ice-ocean measurements across spatial and temporal scales. In this study, we characterize the seasonal behavior and variability in the snow, surface scattering layer, and melt ponds from spring melt to autumn freeze-up using in situ surveys and auxiliary observations. We compare the results to satellite retrievals and output from two models: the Community Earth System Model (CESM2) and the Marginal Ice Zone Modeling and Assimilation System (MIZMAS). During the melt season, the maximum pond coverage and depth were 21% and 22 ± 13 cm, respectively, with distribution and depth corresponding to surface roughness and ice thickness. Compared to observations, both models overestimate melt pond coverage in summer, with maximum values of approximately 41% (MIZMAS) and 51% (CESM2). This overestimation has important implications for accurately simulating albedo feedbacks. During the observed freeze-up, weather events, including rain on snow, caused high-frequency variability in snow depth, while pond coverage and depth remained relatively constant until continuous freezing ensued. Both models accurately simulate the abrupt cessation of melt ponds during freeze-up, but the dates of freeze-up differ. MIZMAS accurately simulates the observed date of freeze-up, while CESM2 simulates freeze-up one-to-two weeks earlier. This work demonstrates areas that warrant future observation-model synthesis for improving the representation of sea-ice processes and properties, which can aid accurate simulations of albedo feedbacks in a warming climate.

We assess the influence of snow on sea ice in experiments using the Community Earth System Model version 2 for a preindustrial and a 2xCO2 climate state. In the preindustrial climate, we find that increasing simulated snow accumulation on sea ice results in thicker sea ice and a cooler climate in both hemispheres. The sea ice mass budget response differs fundamentally between the two hemispheres. In the Arctic, increasing snow results in a decrease in both congelation sea ice growth and surface sea ice melt due to the snow's impact on conductive heat transfer and albedo, respectively. These factors dominate in regions of perennial ice but have a smaller influence in seasonal ice areas. Overall, the mass budget changes lead to a reduced amplitude in the annual cycle of ice thickness. In the Antarctic, with increasing snow, ice growth increases due to snow–ice formation and is balanced by larger basal ice melt, which primarily occurs in regions of seasonal ice. In a warmer 2xCO2 climate, the Arctic sea ice sensitivity to snow depth is small and reduced relative to that of the preindustrial climate. In contrast, in the Antarctic, the sensitivity to snow on sea ice in the 2xCO2 climate is qualitatively similar to the sensitivity in the preindustrial climate. These results underscore the importance of accurately representing snow accumulation on sea ice in coupled Earth system models due to its impact on a number of competing processes and feedbacks that affect the melt and growth of sea ice.

The albedo and transmittance of melting, first-year Arctic sea ice were measured during two cruises of the Impacts of Climate on the Eco-Systems and Chemistry of the Arctic Pacific Environment (ICESCAPE) project during the summers of 2010 and 2011. Spectral measurements were made for both bare and ponded ice types at a total of 19 ice stations in the Chukchi and Beaufort Seas. These data, along with irradiance profiles taken within boreholes, laboratory measurements of the optical properties of core samples, ice physical property observations, and radiative transfer model simulations are employed to describe representative optical properties for melting first-year Arctic sea ice. Ponded ice was found to transmit roughly 4.4 times more total energy into the ocean, relative to nearby bare ice. The ubiquitous surface-scattering layer and drained layer present on bare, melting sea ice are responsible for its relatively high albedo and relatively low transmittance. Light transmittance through ponded ice depends on the physical thickness of the ice and the magnitude of the scattering coefficient in the ice interior. Bare ice reflects nearly three-quarters of the incident sunlight, enhancing its resiliency to absorption by solar insolation. In contrast, ponded ice absorbs or transmits to the ocean more than three-quarters of the incident sunlight. Characterization of the heat balance of a summertime ice cover is largely dictated by its pond coverage, and light transmittance through ponded ice shows strong contrast between first-year and multiyear Arctic ice covers.

The seasonal evolution of melt ponds has been well documented on multiyear and landfast first-year sea ice, but is critically lacking on drifting, first-year sea ice, which is becoming increasingly prevalent in the Arctic. Using 1 m resolution panchromatic satellite imagery paired with airborne and in situ data, we evaluated melt pond evolution for an entire melt season on drifting first-year and multiyear sea ice near the 2011 Applied Physics Laboratory Ice Station (APLIS) site in the Beaufort and Chukchi seas. A new algorithm was developed to classify the imagery into sea ice, thin ice, melt pond, and open water classes on two contrasting ice types: first-year and multiyear sea ice. Surprisingly, melt ponds formed ~3 weeks earlier on multiyear ice. Both ice types had comparable mean snow depths, but multiyear ice had 0–5 cm deep snow covering ~37% of its surveyed area, which may have facilitated earlier melt due to its low surface albedo compared to thicker snow. Maximum pond fractions were 53 ± 3% and 38 ± 3% on first-year and multiyear ice, respectively. APLIS pond fractions were compared with those from the Surface Heat Budget of the Arctic Ocean (SHEBA) field campaign. APLIS exhibited earlier melt and double the maximum pond fraction, which was in part due to the greater presence of thin snow and first-year ice at APLIS. These results reveal considerable differences in pond formation between ice types, and underscore the importance of snow depth distributions in the timing and progression of melt pond formation.

Snow plays a key role in the growth and decay of Arctic sea ice. In winter, it insulates sea ice from cold air temperatures, slowing sea ice growth. From spring into summer, the albedo of snow determines how much insolation is absorbed by the sea ice and underlying ocean, impacting ice melt processes. Knowledge of the contemporary snow depth distribution is essential for estimating sea ice thickness and volume, and for understanding and modeling sea ice thermodynamics in the changing Arctic. This study assesses spring snow depth distribution on Arctic sea ice using airborne radar observations from Operation IceBridge for 2009–2013. Data were validated using coordinated in situ measurements taken in March 2012 during the BRomine, Ozone, and Mercury EXperiment (BROMEX) field campaign. We find a correlation of 0.59 and root-mean-square error of 5.8 cm between the airborne and in situ data. Using this relationship and IceBridge snow thickness products, we compared the recent results with data from the 1937, 1954–1991 Soviet drifting ice stations. The comparison shows thinning of the snow pack, from 35.1 ± 9.4 cm to 22.2 ± 1.9 cm in the western Arctic, and from 32.8 ±D 9.4 cm to 14.5 ± 1.9 cm in the Beaufort and Chukchi seas. These changes suggest a snow depth decline of 37 ± 29% in the western Arctic and 56 ± 33% in the Beaufort and Chukchi seas. Thinning is negatively correlated with the delayed onset of sea ice freeze-up during autumn.