Impurities of glacier ice : accumulation, transport and albedo
Abstract
The Greenland ice sheet is the largest ice mass in the Northern Hemisphere and
has experienced accelerating mass loss in recent decades. An increase in surface
melt is the major cause for the loss of mass. Albedo is a major control of surface
melt and has decreased over the entire Greenland ice sheet, especially at lower
elevations near the margin where glacier ice is exposed for part of each year. The
albedos of ice and snow are lowered by dark impurities such as mineral dust and
black carbon, a by-product of combustion. These impurities are a heterogeneous
group of highly absorbent particles which reach the ice sheet via atmospheric
transport from distant and local sources. Parts of these impurities are buried in the
accumulation zone and become part of the moving ice. After a significant transport
time, up to tens of thousands of years, these englacial impurities re-emerge in
the ablation zone near the margin. These re-emerging impurities, together with
those directly deposited, accumulate on the ice surface. Once located on the ice
surface, these impurities darken the ice surface for several years and lower the
albedo, causing more melt.
Current sea-level projections rarely include albedo as a dynamic model component.
In models which consider albedo, snow albedo is often treated with
sophisticated methods, while ice albedo is still treated as a constant. Albedo
lowering is a major cause of the current mass loss of the Greenland ice sheet,
and the role of ice albedo will increase under a warmer climate. Therefore, this
study presents a model with a dynamic component of ice albedo.
The model framework includes the effect of impurities on the mass balance and
ice sheet geometry. The framework consists of an ice dynamics component which
is linked to a module of englacial impurity transport. This component feeds into
the impurity accumulation module, which deals with both snow and ice. An albedo module which accounts for mineral dust and black carbon accumulation uses this
output to derive daily albedo values. A simplified surface energy balance model
is used to derive the surface melt rate and surface mass balance, which is then
fed back into the ice dynamics component.
This model framework is used to investigate the role of melt-out and impurity
accumulation on the melt of the Greenland ice sheet. For that purpose, simplified
geometry and different temperature pathways are used to simulate the evolution
of the ice sheet over 1000 years. Due to the feedback between melt-out, ice
albedo and impurity accumulation, the role of impurities is disproportionally larger
in warmer scenarios. In the warmest scenario (RCP8.5), a conservative estimate
for the additional mass loss due to impurity accumulation in the year 3000 was
7%.
Melt-out of dust is the largest source of impurities on the ice surface. Darkening
is not always dominated by dust, however, due to the high absorption of black
carbon. The amount of impurities from melt-out depends on the englacial impurity
concentration and surface melt. The englacial impurity concentration, in turn,
relies on the computed age of the ice and a time series of impurity concentration.
Therefore, the accuracy of the transport scheme, which provides the age of the
ice, is crucial for overall accuracy.
A semi-Lagrangian transport scheme of second-order accuracy was implemented
in the 3D ice sheet model SICOPOLIS (SImulation COde for POLythermal
Ice Sheets). The model was applied to the ice sheets of Greenland and Antarctica.
Artificial ice cores of 18O, a proxy for surface temperature, were compared to ice
core data. The results of the second-order scheme were identical to the results of
the first-order scheme in the ice sheet interior. The results deviated substantially,
however, in the outer regions near the margin.
The results emphasise that the role of ice albedo and impurities for the surface
mass balance will be even greater under warmer conditions. Furthermore, the
presented model framework is not limited to Greenland but can also be adapted
for valley glaciers and paleo-ice sheets. Iskappen på Grønland er den største ismassen på den nordlige halvkule og massetapet
her har vært akselererende de siste tiårene. En økning i overflatesmeltingen
er den viktigste årsaken til tap av masse. Albedoen, som er hovedfaktoren for
overflatesmeltingen, har avtatt over hele iskappen på Grønland, og særlig i lavere
høyde nær iskanten hvor breisen er eksponert for et visst tidsrom årlig. Albedoen
for is og snø er redusert pga mørkere partikler og urenheter, som mineralstøv
og black carbon (sot), et bi-produkt av forbrenning. Disse urenhetene er en heterogen
gruppe av svært absorberende partikler som transporteres til iskappen
via atmosfæren fra lokale og fjerntliggende kilder. Deler av disse urenhetene
opptas i akkumulasjonssonen og blir del av den bevegelige isen. Etter lang tids
transport, opptil titusener av år, vil disse englasiale urenhetene komme frem i
ablasjonssonen nær brekanten. Både englasialt transporterte urenhetene og
dirkete avsatte urenheter akkumuleres på isoverflaten. Når urenhetene først er
avsatt på isoverflaten vil de bidra til en mørkere isoverflate med redusert albedo,
og derved bidra til økt smelting.
Nåværende havnivå prognoser inkluderer sjelden albedo som en dynamisk
modellkomponent. I modeller som inkluderer albedo blir snøalbedoen ofte justert
med sofistikerte metoder, mens isalbedoen blir holdt konstant. Albedo er hoveddriveren
av det pågående massetapet for Grønlandsisen, og isalbedoens rolle vil
forsterkes under et varmere klima. Av denne grunn presenteres det i dette studie
en modell med en dynamisk komponent for isalbedo.
Dette modellrammeverket blir brukt til å undersøke hvilken rolle utsmelting og
akkumulasjon av urenheter spiller på avsmeltingen av iskappen på Grønland.
En forenklet geometri og forskjellige temperaturtraseer blir brukt til å simulere
utviklingen av iskappen over 1000 år. På grunn av tilbakekoblingsmekanismene mellom utsmelting, isalbedo, og akkumulasjon av urenheter er rollen av urenheter
uproporsjonalt større i varmere scenarier. I det varmeste scenariet (RCP8.5) er et
konservativt estimat for massetap på grunn av akkumulasjon av urenheter i året
3000 på 7 %.
Utsmelting av støv er den største kilden til urenheter på isoverflaten. Sverting av
overflaten er allikevel ikke alltid dominert av støv, pga. den høye absorbsjonen av
black carbon. Mengden av urenheter pga. utsmelting avhenger av den englasiale
urenhetskonsentrasjonen og overflatesmelting. Den englasiale konsentrasjonen
avhenger igjen av den beregnede alderen av isen og av en tidsserie av urenhetskonsentrasjoner.
Derfor er nøyaktigheten av transportskjema, som gir alderen
på isen, essensiell for den generelle nøyaktigheten.
Et semi-Lagrangian transport skjema av annen ordens nøyaktighet ble implementert
i den tre dimensjaonale iskappemodellen SICOPOLIS (SImulation COde
for POLythermal Ice Sheets). Modellen ble anvendt for iskappene på Grønland og
Antarktis. Genererte iskjerne data av 18O, en proxy for overflate temperatur, ble
sammenlignet med iskjerne data. Resultatet av annen ordens skjemaet var identisk
til resultatet til første ordens skjemaet i det indre av iskappen, men resultatene
avvek vesentlig i de ytre regionene nær iskanten.
Resultatene understreker rollen isalbedo og urenheter har for overfalte massebalansen,
en rolle som også vil bli enda større ved varmere forhold. Videre er det
presenterte modellrammeverket ikke begrenset til Grønland men det kan også
tilpasses dalbreer og paleo-iskapper.