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: Mapping liquid water content in snow at the millimeter scale: an intercomparison of mixed-phase optical property models using hyperspectral imaging and in situ measurements
Abstract. It is well understood that the distribution and quantity of liquid water in snow is relevant for snow hydrology and avalanche forecasting, yet detecting and quantifying liquid water in snow remains a challenge from themicro- to the macro-scale. Using near-infrared (NIR) spectral reflectancemeasurements, previous case studies have demonstrated the capability toretrieve surface liquid water content (LWC) of wet snow by leveraging shifts in the complex refractive index between ice and water. However, different models to represent mixed-phase optical properties have been proposed, including (1) internally mixed ice and water spheres, (2) internally mixed water-coated ice spheres, and (3) externally mixed interstitial ice and water spheres. Here, from within a controlled laboratory environment, we determined the optimal mixed-phase optical property model for simulating wet snow reflectance using a combination of NIR hyperspectral imaging, radiative transfer simulations (Discrete Ordinate Radiative Transfer model, DISORT), and an independent dielectric LWC measurement (SLF Snow Sensor). Maps of LWC were produced by finding the lowest residual between measured reflectance and simulated reflectance in spectral libraries, generated for each model with varying LWC and grain size, and assessed against the in situ LWC sensor. Our results show that the externally mixed model performed the best, retrieving LWC with an uncertainty of ∼1 %, while the simultaneously retrieved grain size better represented wet snow relative to the established scaled band area method. Furthermore, the LWC retrieval method was demonstrated in the field by imaging a snowpit sidewall during melt conditions and mapping LWC distribution in unprecedented detail, allowing for visualization of pooling water and flow features.  more » « less
Award ID(s):
2012091
PAR ID:
10332390
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
16
Issue:
1
ISSN:
1994-0424
Page Range / eLocation ID:
43 to 59
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Effective snow grain radius ( r e ) is mapped at high resolution using near-infrared hyperspectral imaging (NIR-HSI). The NIR-HSI method can be used to quantify r e spatial variability, change in r e due to metamorphism, and visualize water percolation in the snowpack. Results are presented for three different laboratory-prepared snow samples (homogeneous, ice lens, fine grains over coarse grains), the sidewalls of which were imaged before and after melt induced by a solar lamp. The spectral reflectance in each ~3 mm pixel was inverted for r e using the scaled band area of the ice absorption feature centered at 1030 nm, producing r e maps consisting of 54 740 pixels. All snow samples exhibited grain coarsening post-melt as the result of wet snow metamorphism, which is quantified by the change in r e distributions from pre- and post-melt images. The NIR-HSI method was compared to r e retrievals from a field spectrometer and X-ray computed microtomography (micro-CT), resulting in the spectrometer having the same mean r e and micro-CT having 23.9% higher mean r e than the hyperspectral imager. As compact hyperspectral imagers become more widely available, this method may be a valuable tool for assessing r e spatial variability and snow metamorphism in field and laboratory settings. 
    more » « less
  2. Abstract. Surface melting on the Antarctic Ice Sheet has been monitored by satellite microwave radiometry for over 40 years. Despite this long perspective, our understanding of the microwave emission from wet snow is still limited, preventing the full exploitation of these observations to study supraglacial hydrology. Using the Snow Microwave Radiative Transfer (SMRT) model, this study investigatesthe sensitivity of microwave brightness temperature to snow liquid water content at frequencies from 1.4 to 37 GHz. We first determine the snowpack properties for eight selected coastal sites byretrieving profiles of density, grain size and ice layers from microwave observations when the snowpack is dry during wintertime. Second, a series of brightness temperature simulations is run with added water. The results show that (i) a small quantity of liquid water (≈0.5 kg m−2) can be detected, but the actual quantity cannot be retrieved out of the full range of possible water quantities; (ii) the detection of a buried wet layer is possible up to a maximum depth of 1 to 6 m depending on the frequency (6–37 GHz) and on the snow properties (grain size, density) at each site; (iii) surface ponds and water-saturated areas may prevent melt detection, but the current coverage of these waterbodies in the large satellite field of view is presently too small in Antarctica to have noticeable effects; and (iv) at 1.4 GHz, while the simulations are less reliable, we found a weaker sensitivity to liquid water and the maximal depth of detection is relatively shallow (<10 m) compared to the typical radiation penetration depth in dry firn (≈1000 m) at this low frequency. These numerical results pave the way for the development of improved multi-frequency algorithms to detect melt intensity and the depth of liquid water below the surface in the Antarctic snowpack. 
    more » « less
  3. 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
  4. Abstract Mixed‐phase clouds contribute to substantial uncertainties in global climate models due to their complex microphysical properties. Former model evaluations almost exclusively rely on satellite observations to assess cloud phase distributions globally. This study investigated mixed‐phase cloud properties using near global‐scale in situ observation data sets from 14 flight campaigns in combination with collocated output from a global climate model. The Southern Hemisphere (SH) shows significantly higher occurrence frequencies and higher mass fractions of supercooled liquid water than Northern Hemisphere (NH) based on observations at 0.2 and 100 km horizontal scales. Such hemispheric asymmetry is not captured by the model. The model also consistently overestimates liquid water content (LWC) in all cloud phases but shows ice water content (IWC) biases that vary with phase. Key processes contributing to model biases in phase partition can be identified through the combination of evaluation of phase frequency, liquid mass fraction, LWC and IWC. 
    more » « less
  5. Cloud phase and relative humidity (RH) distributions at −67° to 0°C over the Southern Ocean during austral summer are compared between in situ airborne observations and global climate simulations. A scale-aware comparison is conducted using horizontally averaged observations from 0.1 to 50 km. Cloud phase frequencies, RH distributions, and liquid mass fraction are found to be less affected by horizontal resolutions than liquid and ice water content (LWC and IWC, respectively), liquid and ice number concentrations (Ncliqand Ncice, respectively), and ice supersaturation (ISS) frequency. At −10° to 0°C, observations show 27%–34% and 17%–37% of liquid and mixed phases, while simulations show 60%–70% and 3%–4%, respectively. Simulations overestimate (underestimate) LWC and Ncliqin liquid (mixed) phase, overestimate Ncicein mixed phase, underestimate IWC in ice and mixed phases, and underestimate (overestimate) liquid mass fraction below (above) −5°C, indicating that observational constraints are needed for different cloud phases. RH frequently occurs at liquid saturation in liquid and mixed phases for all datasets, yet the observed RH in ice phase can deviate from liquid saturation by up to 20%–40% at −20° to 0°C, indicating that the model assumption of liquid saturation for coexisting ice and liquid is inaccurate for low liquid mass fractions (<0.1). Simulations lack RH variability for partial cloud fractions (0.1–0.9) and underestimate (overestimate) ISS frequency for cloud fraction <0.1 (≥0.6), implying that improving RH subgrid-scale parameterizations may be a viable path to account for small-scale processes that affect RH and cloud phase heterogeneities. Two sets of simulations (nudged and free-running) show very similar results (except for ISS frequency) regardless of sample sizes, corroborating the statistical robustness of the model–observation comparisons. 
    more » « less