The treatment of surface melt, runoff, and the snow-firn-ice transition in ice-sheet models (ISMs) is becoming increasingly important, as mobile liquid on Greenland and Antarctic flanks increases due to climate warming in the next century and beyond. Simple Positive Degree Day (PDD)-based box models used in some ISMs crudely capture liquid storage and refreezing, but need to be extended to include vertical structure through the whole firn-ice column, as in some regional climate models (RCMs). This is a necessary prelude to modeling the flow of mobile meltwater in channel-river-moulin systems, and routing to the base and/or margins of the ice sheet. More detailed column models of snow and firn exist, that include compaction, grain size, and other processes. Some focus on dry-snow zones, and have fine vertical resolution spanning the entire firn column with Lagrangian tracking of annual snow layers (e.g., FirnMICE: Lundin et al., J. Glac., 2017). However, they are mostly too computationally expensive for ISM applications, and are not designed for ablation zones with meltwater and bare ice in summer. More general models are used in some RCMs that include similar physics but with fewer layers, and are applicable both to accumulation and ablation zones. Here we formulate a new snow-firn model, similar to those in RCMs, for use within an ice-sheet model. A limited number of vertical layers is used (∼10), with Lagrangian tracking of layers, grain size evolution, compaction, ice lenses, liquid melting, storage, percolation and runoff. Surface melting is computed from linearized net atmospheric energy fluxes, not from PDDs. The model is tested using the FirnMICE experiments, and using gridded RACMO2 modern climate input over Greenland, seeking to balance model performance with computational efficiency.
more »
« less
A cold laboratory hyperspectral imaging system to map grain size and ice layer distributions in firn cores
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
- PAR ID:
- 10522336
- Publisher / Repository:
- The Cryosphere
- Date Published:
- Journal Name:
- The Cryosphere
- Volume:
- 18
- Issue:
- 4
- ISSN:
- 1994-0424
- Page Range / eLocation ID:
- 1925 to 1946
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. The net rate of snow accumulation b is predicted to increase over large areas of the Antarctic and Greenland ice sheets as the climate warms. Models disagree on how this will affect the thickness of the firn layer – the relatively low-density upper layer of the ice sheets that influences altimetric observations of ice sheet mass change and palaeo-climate reconstructions from ice cores. Here we examine how b influences firn compaction and porosity in a simplified model that accounts for mass conservation, dry firn compaction, grain-size evolution, and the impact of grain size on firn compaction. Treating b as a boundary condition and employing an Eulerian reference frame helps to untangle the factors controlling the b dependence of firn thickness. We present numerical simulations using the model, as well as simplified steady-state approximations to the full model, to demonstrate how the downward advection of porosity and grain size are both affected by b but have opposing impacts on firn thickness. The net result is that firn thickness increases with b and that the strength of this dependence increases with increasing surface grain size. We also quantify the circumstances under which porosity advection and grain-size advection balance exactly, which counterintuitively renders steady-state firn thickness independent of b. These findings are qualitatively independent of the stress-dependence of firn compaction and whether the thickness of the ice sheet is increasing, decreasing, or steady. They do depend on the grain-size dependence of firn compaction. Firn models usually ignore grain-size evolution, but we highlight the complex effect it can have on firn thickness when included in a simplified model. This work motivates future efforts to better observationally constrain the rheological effect of grain size in firn.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
-
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
-
Abstract The expansion of refrozen ice slabs in Greenland's firn may enhance meltwater runoff and increase surface mass loss. However, the impermeability of ice slabs and the pathways for meltwater export from these regions remain poorly characterized. Here, we present ice‐penetrating radar observations of extensive meltwater infiltration and refreezing beneath ice slabs in Northwest Greenland. We show that these buried ice complexes form where supraglacial streams or lakes drain through surface crevasses into relict firn beneath the ice slabs. This suggests that the firn can continue to buffer mass loss from surface meltwater runoff and limit meltwater delivery to the ice sheet bed even after ice slabs have formed. Therefore, a significant time lag may exist between the initial formation of ice slabs and the onset of complete surface runoff and seasonal meltwater drainage to the subglacial system in interior regions of the ice sheet.more » « less