An official website of the United States government
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
Grain-size evolution controls the accumulation dependence of modelled firn thickness
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
- Award ID(s):
- 1935438
- PAR ID:
- 10394597
- Date Published:
- Journal Name:
- The Cryosphere
- Volume:
- 16
- Issue:
- 9
- ISSN:
- 1994-0424
- Page Range / eLocation ID:
- 3413 to 3430
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Reconstructions of paleoatmospheric H2using polar firn air and ice cores would lead to a better understanding of the H2biogeochemical cycle and how it is influenced by climate change and human activity. In this study, the permeability, diffusivity, and solubility of H2are determined experimentally in ice Ih at temperatures relevant to polar ice sheets (199–253 K). The experimental data are used in conjunction with simplified diffusion models to assess the implications for: (a) Diffusion of H2from pressurized closed bubbles to open pores in polar firn, (b) diffusive smoothing of H2gradients in the ice sheet, and (c) post‐coring diffusive losses of H2from ice core samples. The results indicate that diffusive equilibrium between open and closed pores is likely achieved in the firn lock‐in zone. Diffusive smoothing of atmospheric variations is significant and should be accounted for in atmospheric reconstructions on millennial time scales. Diffusive losses from a bubbly ice sample are sufficiently slow that samples may be meaningfully analyzed for H2after storage on the order of a year. These results suggest that the mobility of H2in ice should not preclude the reconstruction of paleoatmospheric H2from firn air and ice cores.more » « less
-
Abstract. Viscous flow in ice is often described by the Glen flow law – anon-Newtonian, power-law relationship between stress and strain rate with astress exponent n ∼ 3. The Glen law is attributed tograin-size-insensitive dislocation creep; however, laboratory and fieldstudies demonstrate that deformation in ice can be strongly dependent ongrain size. This has led to the hypothesis that at sufficiently lowstresses, ice flow is controlled by grain boundary sliding, which explicitly incorporates the grain size dependence of ice rheology. Experimental studiesfind that neither dislocation creep (n ∼ 4) nor grain boundarysliding (n ∼ 1.8) have stress exponents that match the value ofn ∼ 3 in the Glen law. Thus, although the Glen law provides anapproximate description of ice flow in glaciers and ice sheets, itsfunctional form is not explained by a single deformation mechanism. Here weseek to understand the origin of the n ∼ 3 dependence of theGlen law by using the “wattmeter” to model grain size evolution in ice.The wattmeter posits that grain size is controlled by a balance between themechanical work required for grain growth and dynamic grain size reduction.Using the wattmeter, we calculate grain size evolution in two end-membercases: (1) a 1-D shear zone and (2) as a function of depth within anice sheet. Calculated grain sizes match both laboratory data and ice coreobservations for the interior of ice sheets. Finally, we show thatvariations in grain size with deformation conditions result in an effectivestress exponent intermediate between grain boundary sliding and dislocationcreep, which is consistent with a value of n = 3 ± 0.5 over the rangeof strain rates found in most natural systems.more » « less
-
null (Ed.)Climate change is extensively affecting ice sheets resulting in accelerating mass loss in recent decades. Assessment of this reduction and its causes is required to project future ice mass loss. Annual snow accumulation is an important component of the surface mass balance of ice sheets. While in situ snow accumulation measurements are temporally and spatially limited due to their high cost, airborne radar sounders can achieve ice sheet wide coverage by capturing and tracking annual snow layers in the radar images or echograms. In this paper, we use deep learning to uniquely identify the position of each annual snow layer in the Snow Radar echograms taken across different regions over the Greenland ice sheet. We train with more than 15,000 images generated from radar echograms and estimate the thickness of each snow layer within a mean absolute error of 0.54 to 7.28 pixels, depending on dataset. A highly precise snow layer thickness can help improve weather models and, thus, support glaciological studies. Such a well-trained deep learning model can be used with ever-growing datasets to aid in the accurate assessment of snow accumulation on the dynamically changing ice sheets.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5194/tc-16-3413-2022
