13.02.2013 Views

Download the entire proceedings as an Adobe PDF - Eastern Snow ...

Download the entire proceedings as an Adobe PDF - Eastern Snow ...

Download the entire proceedings as an Adobe PDF - Eastern Snow ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

247<br />

63 rd EASTERN SNOW CONFERENCE<br />

Newark, Delaware USA 2006<br />

Qu<strong>an</strong>tifying <strong>the</strong> Effect of Anisotropic Properties in <strong>Snow</strong><br />

for Modelling Meltwater Retention<br />

INTRODUCTION<br />

C. E. BØGGILD 1<br />

In cold snowpacks of <strong>the</strong> Arctic meltwater retention is a signific<strong>an</strong>t factor for timing <strong>an</strong>d<br />

magnitude of runoff. Despite being <strong>an</strong> import<strong>an</strong>t component in Arctic hydrology <strong>an</strong>d glacier m<strong>as</strong>s<br />

bal<strong>an</strong>ce <strong>the</strong> meltwater infiltration <strong>an</strong>d subsequent re-freezing in cold snow is ra<strong>the</strong>r little qu<strong>an</strong>tified<br />

in literature. The qu<strong>an</strong>tity of super-imposed ice (SI) h<strong>as</strong> previously been modeled in onedimensional<br />

vertical profiles. But, since <strong>the</strong>se approaches rely on few points <strong>the</strong> spatial<br />

distribution of SI still needs to be determined more precisely. The problem of moving interface<br />

<strong>as</strong>sociated with SI growth calls for high precision in numerical models, which again greatly<br />

enh<strong>an</strong>ce <strong>the</strong> computational dem<strong>an</strong>d. This computational dem<strong>an</strong>d c<strong>an</strong> likely be <strong>the</strong> re<strong>as</strong>on for <strong>the</strong><br />

lack of SI treatment in most cryospheric models.<br />

Bøggild (submitted) proposes <strong>an</strong> <strong>an</strong>alytical solution to <strong>the</strong> ice warming, from which <strong>the</strong><br />

temporal SI-qu<strong>an</strong>tity c<strong>an</strong> be derived directly without comprehensive numerical modeling.<br />

However, this solution is only valid for iso<strong>the</strong>rmal <strong>an</strong>d isotropic condition, which likely never<br />

occurs in nature. Here is presented <strong>the</strong> result from statistical fitting of <strong>the</strong>rmal gradients <strong>an</strong>d<br />

<strong>the</strong>rmal diffusivity to <strong>the</strong> resulting effect on meltwater refreezing by SI formation. The purpose is<br />

to extend <strong>the</strong> existing <strong>an</strong>alytical solution to <strong>an</strong>isotropic condition also. And, since <strong>the</strong> <strong>the</strong>ory<br />

behind <strong>an</strong> <strong>an</strong>alytical solution is only valid for isotropic condition <strong>the</strong> approach here h<strong>as</strong> to rely on<br />

statistical methods.<br />

METHODS<br />

Bøggild et al (2005) found that for iso<strong>the</strong>rmal <strong>an</strong>d isotropic condition <strong>the</strong> growth rate of SI is:<br />

H ( t)<br />

=<br />

2γ<br />

αt<br />

where H is SI, t is time α <strong>an</strong>d γ are const<strong>an</strong>ts. With a temperature gradient <strong>an</strong>d <strong>an</strong>isotropic<br />

conditions SI is modelled by temporal heat diffusion modelling followed by:<br />

dH<br />

dt<br />

[ ρsi<br />

− ρws(<br />

1−<br />

ω ) ]<br />

L<br />

= xm<br />

∂T<br />

ciρi<br />

∫ dx<br />

∂t<br />

0<br />

1<br />

The University Centre in Svalbard (UNIS), PO Box 156, N-9171 Longyearbyen, Norway<br />

e-mail: carl.egede.boggild@unis.no

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!