Abstract. The Greenland and Antarctic ice sheets are covered in a layer of porous firn. Knowledge of firn structure improves our understanding of ice sheet mass balance, supra- and englacial hydrology, and ice core paleoclimate records. While macroscale firn properties, such as firn density, are relatively easy to measure in the field or lab, more intensive measurements of microstructural properties are necessary to reduce uncertainty in remote sensing observations of mass balance, model meltwater infiltration, and constrain ice age – gas age differences in ice cores. Additionally, as the duration and extent of surface melting increases, refreezing meltwater will greatly alter firn structure. Field observations of firn grain size and ice layer stratigraphy are required to test and validate physical models that simulate the ice-sheet-wide evolution of the firn layer. However, visually measuring grain size and ice layer distributions is tedious, is time-consuming, and can be subjective depending on the method. Here we demonstrate a method to systematically map firn core grain size and ice layer stratigraphy using a near-infrared hyperspectral imager (NIR-HSI; 900–1700 nm). We scanned 14 firn cores spanning ∼ 1000 km across western Greenland’s percolation zone with the NIR-HSI mounted on a linear translation stage in a cold laboratory. We leverage the relationship between effective grain size, a measure of NIR light absorption by firn grains, and NIR reflectance to produce high-resolution (0.4 mm) maps of effective grain size and ice layer stratigraphy. We show the NIR-HSI reproduces visually identified ice layer stratigraphy and infiltration ice content across all cores. Effective grain sizes change synchronously with traditionally measured grain radii with depth, although effective grains in each core are 1.5× larger on average, which is largely related to the differences in measurement techniques. To demonstrate the utility of the firn stratigraphic maps produced by the NIR-HSI, we track the 2012 melt event across the transect and assess its impact on deep firn structure by quantifying changes to infiltration ice content and grain size. These results indicate that NIR-HSI firn core analysis is a robust technique that can document deep and long-lasting changes to the firn column from meltwater percolation while quickly and accurately providing detailed firn stratigraphy datasets necessary for firn research applications.
more »
« less
In situ effective snow grain size mapping using a compact hyperspectral imager
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
- Award ID(s):
- 2012091
- PAR ID:
- 10253899
- Date Published:
- Journal Name:
- Journal of Glaciology
- Volume:
- 67
- Issue:
- 261
- ISSN:
- 0022-1430
- Page Range / eLocation ID:
- 49 to 57
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Snow albedo, a measure of the amount of solar radiation that is reflected at the snow surface, plays a critical role in Earth’s climate and in regional hydrology because it is a primary driver of snowmelt timing. Satellite multi-spectral remote sensing provides a multi-decade record of land surface reflectance, from which snow albedo can be retrieved. However, this observational record is challenging to assess because discretein situobservations are not well suited for validation of snow properties at the spatial resolution of satellites (tens to hundreds of meters). For example, snow grain size, a primary driver of snow albedo, can vary at the sub-meter scale driven by changes in aspect, elevation, and vegetation. Here, we present a new uncrewed aerial vehicle hyperspectral imaging (UAV-HSI) method for mapping snow surface properties at high resolution (20 cm). A Resonon near-infrared HSI was flown on a DJI Matrice 600 Pro over the meadow encompassing Swamp Angel Study Plot in Senator Beck Basin, Colorado. Using a radiative transfer forward modeling approach, effective snow grain size and albedo maps were produced from measured surface reflectance. Coincident ground observations were used for validation; relative to retrievals from a field spectrometer the mean grain size difference was 2 μm, with an RMSE of 12 μm, and the mean broadband albedo was within 1% of that measured near the center of the flight area. Even though the snow surface was visually homogenous, the maps showed spatial variability and coherent patterns in the freshly fallen snow. To demonstrate the potential for UAV-HSI to be used to improve validation of satellite retrievals, the high-resolution maps were used to assess grain size and albedo retrievals, and subpixel variability, across 17 Landsat 9 OLI pixels from a satellite overpass with similar conditions two days following the flight. Although Landsat 9 did not capture the same range of values and spatial variability as the UAV-HSI, on average the comparison showed good agreement, with a mean grain size difference of 9 μm and the same broadband albedo (86%).more » « less
-
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
-
Hyperspectral imaging (HSI) technology has been increasingly used in Earth and planetary sciences. This imaging technique has been successfully tested on ice cores using VNIR (visible and near-infrared, 380-1000 nm) (Garzonio et al., 2018) and near-infrared (900 - 1700 nm) (McDowell et al, 2023) line-scan cameras. Results show that HSI data greatly expand ice core line-scan imaging capabilities, previously used with gray or RGB cameras (see summary in Dey et al., 2023). Combinations of selected HSI bands from the hyperspectral data cube improve feature detection in ice core stratigraphy, and map distribution of volcanic material, dust, air bubbles, fractures, and ice crystals in ice cores. Captured spectral information provides unique fingerprints for specific materials present in ice cores. This method helps to guide ice core sampling because it provides non-destructive, rapid visualization of microstructural properties, layering, bubble contents, increases in dust, or presence of tephra material. Precise identification of these atmospheric components is important for understanding past climate drivers reconstructed from ice cores. As part of the COLDEX project (Brook et al., this meeting) we adapted the SPECIM SisuSCS HSI system for ice core imaging. The ice core scanning system is housed inside the ca. -20ºC main NSF ICF freezer, and externally computer-controlled. The operator monitors scanning operations and communicates with personnel inside of the freezer via radio. The system is equipped with a SPECIM FX10 camera that measures up to 224 bands in the VNIR range. We modified the ice core holder tray and installed a heated enclosure for the camera. The system uses SCHOTT DCR III Fiber Optic light sources with an OSL2BIR bulb from Thorlabs. IR filters are removed to extend the light spectral range beyond the 700 nm limit without heating the ice core surface during rapid (<5 minutes) scanning of an entire meter-long section. Emitted light enters ice at a 45º angle from two top and two bottom light sources. To calibrate absolute reflectance we use three Spectralon panels with 100, 50 and 20% reflectance values with every scan as well as several secondary reflective standards and USAF targets for geometric corrections. We are developing Python-based open source data processing routines and currently comparing HSI data with existing ice core physical and chemical measurements. The goal is to fully integrate the ice core HSI system with ice core processing at the NSF ICF. Dey et al., 2023. Application of Visual Stratigraphy from Line-Scan Images to Constrain Chronology and Melt Features of a Firn Core from Coastal Antarctica. Journal of Glaciology 69(273): 179–90. https://doi.org/10.1017/jog.2022.59.Garzonio et al., 2018. A Novel Hyperspectral System for High Resolution Imaging of Ice Cores: Application to Light-Absorbing Impurities and Ice Structure. Cold Regions Science and Technology 155: 47–57. https://doi.org/10.1016/j.coldregions.2018.07.005.McDowell et al., 2023. A Cold Laboratory Hyperspectral Imaging System to Map Grain Size and Ice Layer Distributions in Firn Cores. Preprint. Ice sheets/Instrumentation. https://doi.org/10.5194/egusphere-2023-2351.more » « less
-
Abstract Understanding soil organic carbon (SOC) response to global change has been hindered by an inability to map SOC at horizon scales relevant to coupled hydrologic and biogeochemical processes. Standard SOC measurements rely on homogenized samples taken from distinct depth intervals. Such sampling prevents an examination of fine‐scale SOC distribution within a soil horizon. Visible near‐infrared hyperspectral imaging (HSI) has been applied to intact monoliths and split cores surfaces to overcome this limitation. However, the roughness of these surfaces can influence HSI spectra by scattering reflected light in different directions posing challenges to fine‐scale SOC mapping. Here, we examine the influence of prescribed surface orientation on reflected spectra, develop a method for correcting topographic effects, and calibrate a partial least squares regression (PLSR) model for SOC prediction. Two empirical models that account for surface slope, aspect, and wavelength and two theoretical models that account for the geometry of the spectrometer were compared using 681 homogenized soil samples from across the United States that were packed into sample wells and presented to the spectrometer at 91 orientations. The empirical approach outperformed the more complex geometric models in correcting spectra taken at non‐flat configurations. Topographically corrected spectra reduced bias and error in SOC predicted by PLSR, particularly at slope angles greater than 30°. Our approach clears the way for investigating the spatial distributions of multiple soil properties on rough intact soil samples.more » « less