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: SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice
Abstract. Accurate modeling of cryospheric surface albedo is essential for ourunderstanding of climate change as snow and ice surfaces regulate the globalradiative budget and sea-level through their albedo and massbalance. Although significant progress has been made using physicalprinciples to represent the dynamic albedo of snow, models of glacier icealbedo tend to be heavily parameterized and not explicitly connected withphysical properties that govern albedo, such as the number and size of airbubbles, specific surface area (SSA), presence of abiotic and biotic lightabsorbing constituents (LACs), and characteristics of any overlyingsnow. Here, we introduce SNICAR-ADv4, an extension of the multi-layertwo-stream delta-Eddington radiative transfer model with theadding–doubling solver that has been previously applied to represent snowand sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnelreflectance and transmittance between overlying snow and higher-densityglacier ice, scattering by air bubbles of varying sizes, and numerous typesof LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadbandalbedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03for bare and bubble-free ice under direct light. Our results indicate thatrepresenting ice with a density of 650 kg m−3 as snow with norefractive Fresnel layer, as done previously, generally overestimates theBBA by an average of 0.058. However, because mostnaturally occurring ice surfaces are roughened “white ice”, we recommendmodeling a thin snow layer over bare ice simulations. We find optimalagreement with measurements by representing cryospheric media with densitiesless than 650 kg m−3 as snow and larger-density media as bubbly icewith a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedoimpacts from LACs with changing ice SSA, with peak impact per unit mass ofLACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACsactually increase the albedo. SNICAR-ADv4 represents smooth transitionsbetween snow, firn, and ice surfaces and accurately reproduces measuredspectral albedos of a variety of glacier surfaces. This work paves the wayfor adapting SNICAR-ADv4 to be used in land ice model components of Earthsystem models.  more » « less
Award ID(s):
1712695
PAR ID:
10382458
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
16
Issue:
4
ISSN:
1994-0424
Page Range / eLocation ID:
1197 to 1220
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. The Snow, Ice, and Aerosol Radiative (SNICAR) model has been used in various capacities over the last 15 years to model the spectral albedo of snow with light-absorbing constituents (LACs). Recent studies have extended the model to include an adding-doubling two-stream solver and representations of non-spherical ice particles; carbon dioxide snow; snow algae; and new types of mineral dust, volcanic ash, and brown carbon. New options also exist for ice refractive indices and solar-zenith-angle-dependent surface spectral irradiances used to derive broadband albedo. The model spectral range was also extended deeper into the ultraviolet for studies of extraterrestrial and high-altitude cryospheric surfaces. Until now, however, these improvements and capabilities have not been merged into a unified code base. Here, we document the formulation and evaluation of the publicly available SNICAR-ADv3 source code, web-based model, and accompanying library of constituent optical properties. The use of non-spherical ice grains, which scatter less strongly into the forward direction, reduces the simulated albedo perturbations from LACs by ∼9 %–31 %, depending on which of the three available non-spherical shapes are applied. The model compares very well against measurements of snow albedo from seven studies, though key properties affecting snow albedo are not fully constrained with measurements, including ice effective grain size of the top sub-millimeter of the snowpack, mixing state of LACs with respect to ice grains, and site-specific LAC optical properties. The new default ice refractive indices produce extremely high pure snow albedo (>0.99) in the blue and ultraviolet part of the spectrum, with such values only measured in Antarctica so far. More work is needed particularly in the representation of snow algae, including experimental verification of how different pigment expressions and algal cell concentrations affect snow albedo. Representations and measurements of the influence of liquid water on spectral snow albedo are also needed. 
    more » « less
  2. 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. 
    more » « less
  3. Abstract Snow is the most reflective natural surface on Earth. Since fresh snow on bare sea ice increases the surface albedo, the impact of summer snow accumulation can have a negative radiative forcing effect, which would inhibit sea ice surface melt and potentially slow sea‐ice loss. However, it is not well known how often, where, and when summer snowfall events occur on Arctic sea ice. In this study, we used in situ and model snow depth data paired with surface albedo and atmospheric conditions from satellite retrievals to characterize summer snow accumulation on Arctic sea ice from 2003 to 2017. We found that, across the Arctic, ∼2 snow accumulation events occurred on initially snow‐free conditions each year. The average snow depth and albedo increases were ∼2 cm and 0.08, respectively. 16.5% of the snow accumulation events were optically thick (>3 cm deep) and lasted 2.9 days longer than the average snow accumulation event (3.4 days). Based on a simple, multiple scattering radiative transfer model, we estimated a −0.086 ± 0.020 W m−2change in the annual average top‐of‐the‐atmosphere radiative forcing for summer snowfall events in 2003–2017. The following work provides new information on the frequency, distribution, and duration of observed snow accumulation events over Arctic sea ice in summer. Such results may be particularly useful in understanding the impacts of ephemeral summer weather on surface albedo and their propagating effects on the radiative forcing over Arctic sea ice, as well as assessing climate model simulations of summer atmosphere‐ice processes. 
    more » « less
  4. Abstract Feedbacks between ice melt, glacier flow and ocean circulation can rapidly accelerate ice loss at tidewater glaciers and alter projections of sea-level rise. At the core of these projections is a model for ice melt that neglects the fact that glacier ice contains pressurized bubbles of air due to its formation from compressed snow. Current model estimates can underpredict glacier melt at termini outside the region influenced by the subglacial discharge plume by a factor of 10–100 compared with observations. Here we use laboratory-scale experiments and theoretical arguments to show that the bursting of pressurized bubbles from glacier ice could be a source of this discrepancy. These bubbles eject air into the seawater, delivering additional buoyancy and impulses of turbulent kinetic energy to the boundary layer, accelerating ice melt. We show that real glacier ice melts 2.25 times faster than clear bubble-free ice when driven by natural convection in a laboratory setting. We extend these results to the geophysical scale to show how bubble dynamics contribute to ice melt from tidewater glaciers. Consequently, these results could increase the accuracy of modelled predictions of ice loss to better constrain sea-level rise projections globally. 
    more » « less
  5. 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. 
    more » « less