scholarly journals A continuum model for meltwater flow through compacting snow

2017 ◽  
Vol 11 (6) ◽  
pp. 2799-2813 ◽  
Author(s):  
Colin R. Meyer ◽  
Ian J. Hewitt

Abstract. Meltwater is produced on the surface of glaciers and ice sheets when the seasonal energy forcing warms the snow to its melting temperature. This meltwater percolates into the snow and subsequently runs off laterally in streams, is stored as liquid water, or refreezes, thus warming the subsurface through the release of latent heat. We present a continuum model for the percolation process that includes heat conduction, meltwater percolation and refreezing, as well as mechanical compaction. The model is forced by surface mass and energy balances, and the percolation process is described using Darcy's law, allowing for both partially and fully saturated pore space. Water is allowed to run off from the surface if the snow is fully saturated. The model outputs include the temperature, density, and water-content profiles and the surface runoff and water storage. We compare the propagation of freezing fronts that occur in the model to observations from the Greenland Ice Sheet. We show that the model applies to both accumulation and ablation areas and allows for a transition between the two as the surface energy forcing varies. The largest average firn temperatures occur at intermediate values of the surface forcing when perennial water storage is predicted.

2017 ◽  
Author(s):  
Colin R. Meyer ◽  
Ian J. Hewitt

Abstract. Meltwater is produced on the surface of glaciers and ice sheets when the seasonal surface energy forcing warms the snow to its melting temperature. This meltwater can run off the surface in streams or percolate through the porous snow and refreeze, which warms the subsurface through the release of latent heat. We model the percolation process from first principles using a continuum model that includes heat conduction, meltwater percolation and refreezing, as well as mechanical compaction. The model is forced by surface mass and energy balances. When the surface temperature reaches the melting point, we compute the amount of meltwater produced and allow it to percolate through the snow according to Darcy's law, or to run off the surface if the snow is already saturated. The model outputs the temperature, density, and water content profiles as well as the surface runoff and water storage. We compare the propagation of freezing fronts that occur in the model to observations from the Greenland ice sheet. The model applies to both accumulation and ablation areas and allows for a transition between the two as the surface energy forcing varies. The largest firn temperatures occur at intermediate values of the surface forcing when perennial water storage is predicted.


2016 ◽  
Vol 10 (5) ◽  
pp. 1933-1946 ◽  
Author(s):  
Michiel R. van den Broeke ◽  
Ellyn M. Enderlin ◽  
Ian M. Howat ◽  
Peter Kuipers Munneke ◽  
Brice P. Y. Noël ◽  
...  

Abstract. We assess the recent contribution of the Greenland ice sheet (GrIS) to sea level change. We use the mass budget method, which quantifies ice sheet mass balance (MB) as the difference between surface mass balance (SMB) and solid ice discharge across the grounding line (D). A comparison with independent gravity change observations from GRACE shows good agreement for the overlapping period 2002–2015, giving confidence in the partitioning of recent GrIS mass changes. The estimated 1995 value of D and the 1958–1995 average value of SMB are similar at 411 and 418 Gt yr−1, respectively, suggesting that ice flow in the mid-1990s was well adjusted to the average annual mass input, reminiscent of an ice sheet in approximate balance. Starting in the early to mid-1990s, SMB decreased while D increased, leading to quasi-persistent negative MB. About 60 % of the associated mass loss since 1991 is caused by changes in SMB and the remainder by D. The decrease in SMB is fully driven by an increase in surface melt and subsequent meltwater runoff, which is slightly compensated by a small ( <  3 %) increase in snowfall. The excess runoff originates from low-lying ( <  2000 m a.s.l.) parts of the ice sheet; higher up, increased refreezing prevents runoff of meltwater from occurring, at the expense of increased firn temperatures and depleted pore space. With a 1991–2015 average annual mass loss of  ∼  0.47 ± 0.23 mm sea level equivalent (SLE) and a peak contribution of 1.2 mm SLE in 2012, the GrIS has recently become a major source of global mean sea level rise.


2016 ◽  
Author(s):  
Michiel van den Broeke ◽  
Ellyn Enderlin ◽  
Ian Howat ◽  
Peter Kuipers Munneke ◽  
Brice Noël ◽  
...  

Abstract. We assess the recent contribution of the Greenland ice sheet (GrIS) to sea level change. We use the mass budget method, which quantifies ice sheet mass balance (MB) as the difference between surface mass balance (SMB) and solid ice discharge across the grounding line (D). A comparison with independent gravity change observations from GRACE shows good agreement for the overlapping period 2002–2015, giving confidence in the partitioning of recent GrIS mass changes. The estimated 1995 value of D and the 1958–1995 average value of SMB are similar at 411 and 418 Gt yr-1, respectively, suggesting that ice flow in the mid-nineties was well adjusted to the average annual mass input, reminiscent of an ice sheet in approximate balance. Starting in the early to mid-1990's, SMB decreased while D increased, leading to quasi-persistent negative MB. About 60 % of the associated mass loss since 1991 is caused by changes in SMB and the remainder by D. The decrease in SMB is fully driven by an increase in surface melt and subsequent meltwater runoff, which is slightly compensated by a small (< 3 %) increase in snowfall. The excess runoff originates from low-lying (< 2000 m a.s.l.) parts of the ice sheet; higher up, increased refreezing prevents runoff of meltwater to occur, at the expense of increased firn temperatures and depleted pore space. With a 1991–2015 average annual mass loss of ~ 0.47 ± 0.23 mm sea level equivalent (SLE) and a peak contribution of 1.2 mm SLE in 2012, the GrIS has recently become a major source of global mean sea level rise.


2012 ◽  
Vol 6 (4) ◽  
pp. 3101-3147 ◽  
Author(s):  
X. Fettweis ◽  
B. Franco ◽  
M. Tedesco ◽  
J. H. van Angelen ◽  
J. T. M. Lenaerts ◽  
...  

Abstract. We report future projections of Surface Mass Balance (SMB) over the Greenland ice sheet (GrIS) obtained with the regional climate model MAR, forced by the outputs of three CMIP5 General Circulation Models (GCMs) when considering two different warming scenarios (RCP 4.5 and RCP 8.5). The GCMs selected in this study have been chosen according to their ability to simulate the current climate over Greenland. Our results indicate that in a warmer climate (i) the mass gained due to increased precipitation over GrIS does not compensate the mass lost through increased run-off; (ii) the surface melt increases non-linearly with rising temperatures due to the positive feedback between surface albedo and melt, associated with the expansion of bare ice zones which, in addition, decreases the ice sheet refreezing capacity; (iii) most of the precipitation is expected to fall as rainfall in summer, which further increases surface melt; (iv) no considerable change is expected on the length of the melting season, since heavier winter snowfall dampens the melt increase at the end of spring; (v) the increase of meltwater run-off versus temperature anomalies is dependent of the GCM-forced MAR ability to simulate the current climate; (vi) the MAR-simulated SMB changes can be approximated using the annual accumulated snowfall and summer 600 hPa temperature increase simulated by the forcing GCMs. In view of the large range in the CMIP5 future projections for the same future scenario, the GCM-based SMB approximations allow us to estimate what future projections are most likely within the CMIP5 multi-model ensemble. In 2100, the ensemble mean projects a sea level rise, resulting from a GrIS SMB decrease, estimated to be +4 &amp;pm; 2 cm and +9 &amp;pm; 4 cm for the RCP 4.5 and RCP 8.5 scenarios, respectively. The GrIS SMB should remain positive with respect to RCP 4.5 scenario and becomes negative around 2070 in the case of the RCP 8.5 scenario since a global warming >+3 °C is needed. However, these future projections do not consider the positive melt-elevation feedback because the ice sheet topography is fixed through the whole simulation. In this regard, the MAR simulations suggest a cumulative ice sheet thinning by 2100 of ~100–200 m in the ablation zone. This highlights the importance of coupling climate models to an ice sheet model to consider the future response of both surface processes and ice-dynamic changes, and their mutual feedbacks to rising temperatures.


2021 ◽  
Author(s):  
Max Brils ◽  
Peter Kuipers Munneke ◽  
Willem Jan van de Berg ◽  
Achim Heilig ◽  
Baptiste Vandercrux ◽  
...  

&lt;p&gt;Recent studies indicate that a declining surface mass balance will dominate the Greenland Ice Sheet&amp;#8217;s (GrIS) contribution to 21&lt;sup&gt;st&lt;/sup&gt; century sea level rise. It is therefore crucial to understand the liquid water balance of the ice sheet and its response to increasing temperatures and surface melt if we want to accurately predict future sea level rise. The ice sheet firn layer covers ~90% of the GrIS and provides pore space for storage and refreezing of meltwater. Because of this, the firn layer can retain up to ~45% of the surface meltwater and thus act as an efficient buffer to ice sheet mass loss. However, in a warming climate this buffer capacity of the firn layer is expected to decrease, amplifying meltwater runoff and sea-level rise. Dedicated firn models are used to understand how firn layers evolve and affect runoff. Additionally, firn models are used to estimate the changing thickness of the firn layer, which is necessary in altimetry to convert surface height change into ice sheet mass loss.&lt;/p&gt;&lt;p&gt;Here, we present the latest version of our firn model IMAU-FDM. With respect to the previous version, changes have been made to the handling of the freshly fallen snow, the densification rate of the firn and the conduction of heat. These changes lead to an improved representation of firn density and temperature. The results have been thoroughly validated using an extensive dataset of density and temperature measurements that we have compiled covering 126 different locations on the GrIS. Meltwater behaviour in the model is validated with upward-looking GPR measurements at Dye-2. Lastly, we present an in-depth look at the evolution firn characteristics at some typical locations in Greenland.&lt;/p&gt;&lt;p&gt;Dedicated, stand-alone firn models offer various benefits to using a regional climate model with an embedded firn model. Firstly, the vertical resolution for buried snow and ice layers can be larger, improving accuracy. Secondly, a stand-alone firn model allows for spinning up the model to a more accurate equilibrium state. And thirdly, a stand-alone model is more cost- and time-effective to use. Firn models are increasingly capable of simulating the firn layer, but areas with large amounts of melt still pose the greatest challenge.&lt;/p&gt;


1998 ◽  
Vol 09 (08) ◽  
pp. 1535-1543 ◽  
Author(s):  
Somalee Datta ◽  
Sidney Redner

The phenomenon of clogging in depth filtration is investigated, in which a dirty fluid is "cleaned" by the trapping of dirt particles within the pore space during flow through a porous medium. This gives rise to a self-generated gradient percolation process which exhibits a power law distribution for the density of trapped particles at downstream distance x from the input, both in idealized and lattice networks. Implications for efficient filter design are also mentioned.


1991 ◽  
Vol 37 (126) ◽  
pp. 296-300
Author(s):  
R. D Moore

AbstractEnergy exchange between the atmosphere and a melting glacier surface is mediated by the presence of a water layer. Under conditions of rapid melt and/or heavy rainfall, the possibility exists that a supraglacial run-off layer can advect sensible heat and influence the spatial variations of melt. The potential magnitude of such advection was investigated by numerically solving differential equations expressing the mass and energy balances of a two-dimensional run-off layer. Solutions were obtained for conditions typical of rainfall events, in which the potential for supraglacial heat advection should be maximal. The solutions indicate that advection cannot influence macro-scale melt patterns and surface morphology, except perhaps under heavy rainfall and/or rapid melt conditions, but can possibly cause micro-scale variations in ice melt. One-dimensional energy-balance models, which have normally been applied over glacier surfaces, should remain valid for most conditions.


2019 ◽  
Vol 66 (255) ◽  
pp. 61-73
Author(s):  
David W. Ashmore ◽  
Douglas W. F. Mair ◽  
David O. Burgess

AbstractThe retention of meltwater in the accumulation area of the Greenland ice sheet and other Arctic ice masses buffers their contribution to sea level change. However, sustained warming also results in impermeable ice layers or ‘ice slabs’ that seal the underlying pore space. Here, we use a 1-D, physically based, high-resolution model to simulate the surface mass balance (SMB), percolation, refreezing, ice layer formation and runoff from across the high-elevation area of Devon Ice Cap, Canada, from 2001 to 2016. We vary the thickness of the ‘impermeable’ ice layer at which underlying firn becomes inaccessible to meltwater. Thick near-surface ice layers are established by an initial deep percolation, the formation of decimetre ice layers and the infilling of interleaving pore space. The cumulative SMB increases by 48% by varying impermeable layer thickness between 0.01 and 5 m. Within this range we identify narrower range (0.25–1 m) that can simulate both the temporal variability in SMB and the observed near-surface density structure. Across this range, cumulative SMB variation is limited to 6% and 45–49% of mass retention takes place within the annually replenished snowpack. Our results indicate cooler summers after intense mid-2000s warming have led to a partial replenishment of pore space.


2018 ◽  
Author(s):  
Stefan R. M. Ligtenberg ◽  
Peter Kuipers Munneke ◽  
Brice P. Y. Noël ◽  
Michiel R. van den Broeke

Abstract. By providing pore space for storage or refreezing of meltwater, the Greenland ice sheet firn layer strongly modulates runoff. Correctly representing the firn layer is therefore crucial for Greenland (surface) mass balance studies. Here, we present an improved simulation of the Greenland firn layer with the firn model IMAU-FDM forced by the latest output of the regional climate model RACMO2, version 2.3p2. In the percolation zone, much improved agreement with firn density and temperature observations is found. A full simulation of Greenland firn at high temporal (10 days) and spatial (11 km) resolution is available for the period 1960–2016.


2013 ◽  
Vol 7 (2) ◽  
pp. 1833-1870 ◽  
Author(s):  
S. MacDonell ◽  
C. Kinnard ◽  
T. Mölg ◽  
L. Nicholson ◽  
J. Abermann

Abstract. Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s−1), high shortwave radiation receipt (mean annual 295 W m−2) and low precipitation rates (mean annual 45 mm w.e.). Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.


Sign in / Sign up

Export Citation Format

Share Document