Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-05T01:06:35.376Z Has data issue: false hasContentIssue false

In situ quantification of supraglacial cryoconite morphodynamics using time-lapse imaging: an example from Svalbard

Published online by Cambridge University Press:  08 September 2017

Tristram D.L. Irvine-Fynn
Affiliation:
Department of Geography, University of Sheffield, Winter Street, Sheffield S10 2TN, UK Institute of Geography and Earth Sciences, Aberystwyth University, Penglais Campus, Aberystwyth SY23 3DB, UK E-mail: [email protected]
Jonathan W. Bridge
Affiliation:
Cell-Mineral Research Centre, Kroto Research Institute, University of Sheffield, Broad Lane, Sheffield S3 7HQ, UK
Andrew J. Hodson
Affiliation:
Department of Geography, University of Sheffield, Winter Street, Sheffield S10 2TN, UK
Rights & Permissions [Opens in a new window]

Abstract

There is growing recognition of the significance of biologically active supraglacial dust (cryoconite) for glacial mass balance and ecology. Nonetheless, the processes controlling the distribution, transport and fate of cryoconite particles in the glacial system remain somewhat poorly understood. Here, using a 216 hour time series of plot-scale (0.04 m2) images, we quantify the small-scale dynamics of cryoconite on Longyearbreen, Svalbard. We show significant fluctuations in the apparent cryoconite area and dispersion of cryoconite over the plot, within the 9 day period of observations. However, the net movement of cryoconite across the ice surface averaged only 5.3 mm d−1. High-resolution measurements of cryoconite granule motion showed constant, random motion but weak correlation with meteorological forcing factors and no directional trends for individual particle movement. The high-resolution time-series data suggest that there is no significant net transport of dispersed cryoconite material across glacier surfaces. The areal coverage and motion of particles within and between cryoconite holes appears to be a product of differential melting leading to changes in plot-scale microtopography, local meltwater flow dynamics and weather-dependent events. These subtle processes of cryoconite redistribution may be significant for supraglacial albedo and have bearing on the surface energy balance at the glacier scale.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2011

Introduction

The surface of Arctic glaciers typically exhibits a widespread yet spatially variable coverage of granular biologically active dust termed cryoconite. Cryoconite consists of significant quantities of organic material, algae and bacteria bound together in and around inorganic, mineral particles (Reference Takeuchi, Kohshima, Goto-Azuma and KoernerTakeuchi and others, 2001a, Reference Takeuchi, Nishiyama and Li2010; Reference HodsonHodson and others, 2008; Reference Langford, Hodson, Banwart and BøggildLangford and others, 2010). The presence of such material in the supraglacial landscape has a profound effect upon ice melt through its influence on surface albedo (e.g. Reference BøggildBøggild, 1997; Reference Takeuchi, Kohshima and SekoTakeuchi and others, 2001b; Reference TakeuchiTakeuchi, 2009; Reference Bøggild, Brandt, Brown and WarrenBøggild and others, 2010) and on nutrient cycling, particularly nitrogen and carbon (Reference Hodson, Mumford, Kohler and WynnHodson and others, 2005a, Reference Hodson2007, Reference Hodson2010; Reference Stibal, Tranter, Benning and RehákStibal and others, 2008; Reference Anesio, Hodson, Fritz, Psenner and SattlerAnesio and others, 2009, Reference Anesio2010). There is a growing literature examining cryoconite microbial ecology and its spatial variability (e.g. De Smet and Van Rompu, 1994; Reference Säwström, Mumford, Marshall, Hodson and Laybourn-ParrySäwström and others, 2002; Reference Christner, Kvito and ReeveChristner and others, 2003; Reference Mueller and PollardMueller and Pollard, 2004; Reference Porazinska, Fountain, Nylen, Tranter, Virginia and WallPorazinska and others, 2004; Reference Stibal, Šabacká and KaštovskáStibal and others, 2006; Reference Anesio, Mindl, Laybourn-Parry, Hodson and SattlerAnesio and others, 2007; Reference Foreman, Sattler, Mikucki, Porazinska and PriscuForeman and others, 2007; Reference Langford, Hodson, Banwart and BøggildLangford and others, 2010; Reference Segawa, Takeuchi, Ishida, Kanda and KohshimaSegawa and others, 2010; Reference EdwardsEdwards and others, 2011). Recent work has highlighted the intimate coupling between cryoconite and supraglacial hydrology both in Arctic and Antarctic settings (Reference Fountain, Tranter, Nylen, Lewis and MuellerFountain and others, 2004; Reference HodsonHodson and others, 2007; Reference Irvine-FynnIrvine-Fynn, 2008; Reference MacDonell and FitzsimonsMacDonell and Fitzsimons, 2008; Reference EdwardsEdwards and others, 2011). Glaciers are increasingly understood to represent an important source of labile carbon to connected hydrological and ecological systems (Reference HoodHood and others, 2009). Yet, despite this growing perception, the processes and conditions controlling cryoconite distribution, transport and fate in the supraglacial environment remain relatively poorly understood. The residence time of cryoconite on ice surfaces, lateral rates of movement of cryoconite towards ice margins or supraglacial streams, and the interactions of particles with each other and with ice surface microtopography are all, to date, unquantified.

Using imagery acquired from a camera mounted on an unmanned aerial vehicle (UAV), Reference HodsonHodson and others (2007) demonstrated a spatial gradient in the ratio of cryoconite form over Midtre Lovénbreen, a small valley glacier in Svalbard: distributed cryoconite forms were prevalent at lower elevations while isolated cryoconite holes dominated up-glacier areas. This contrast was inferred to be the result of hydrological mobilization: at lower elevations, greater melt rates and marginally steeper surface slopes result in frequent or periodic disruption and redistribution of cryoconite holes and their contents (Reference HodsonHodson and others, 2007, Reference Hodson2008). Exploring this in greater detail, Reference HodsonHodson and others (2010) and Reference Irvine-Fynn, Bridge and HodsonIrvine-Fynn and others (2010) used plot-scale image data collected from Longyearbreen, another valley glacier in Svalbard. They conjectured a net downslope movement of cryoconite interspersed with ‘relief inversions’; during sunny periods cryoconite melts into the ice due to preferential radiation absorption, creating depressions (Reference McIntyreMcIntyre, 1984; Reference Podgorny and GrenfellPodgorny and Grenfell, 1996; Reference Bøggild, Brandt, Brown and WarrenBøggild and others, 2010) and lessening the apparent cryoconite area relative to the ice surface plane. Conversely, ice melt through localized water flow, and during cloudy periods when dominated by turbulent energy fluxes, results in the melt-out of cryoconite (Reference McIntyreMcIntyre, 1984; Reference HodsonHodson and others, 2008), thereby reducing the contrasts in microtopography and increasing its apparent area.

Here, building upon the earlier work of Reference Irvine-Fynn, Bridge and HodsonIrvine-Fynn and others (2010) and Reference HodsonHodson and others (2010), we present a full exploration of a 9 day hourly time-series image set for a typical supraglacial plot. The datasets were acquired with commercially available digital cameras in an Arctic setting, where during summer months cryoconite holes tend not to display ice lids. The aim of this work is to quantify the nature of cryoconite mobility within the supraglacial environment, by conducting the first detailed plot-scale analysis of cryoconite granule distribution and mobility under multiple diurnal cycles.

Materials and Methods

Field site and cryoconite imaging

Over a 9 day period (21–29 July 2008; day of year (DOY) 203–211) a time-lapse sequence of 216 images was collected at a plot on Longyearbreen (78°1049 N, 15°3021′′ E). The image datasets were captured as uncompressed 7.08 Mpix (3072 × 2304) digital colour (RGB) JPEG files, using a Pentax Optio WP30 digital camera triggered at hourly intervals. The camera installation comprised a vertical aluminium pole drilled and frozen into the glacier ice with the camera attached to a crossbar; the arrangement was immobile although the ice surface ablation resulted in a change in the camera-to-ice distance of ∼0.4 m over the imaging period. The plot, averaging ∼0.12 m2, was considered to be representative of an ablating Arctic glacier surface, exhibiting an irregular, time-variable microtopography. The imaged plot sloped 2° down-glacier to the northeast (52° from north). A 47 mm filter paper was placed on the ice surface to act as a control for image scale; the filter paper was present throughout the image sequence. The camera was oriented such that the local surface slope was parallel to the vertical (Y) axis in images. The broader glaciological characteristics of the site are described by Reference Etzelmüller, Ödegård, Vatne, Mysterud, Tonning and SollidEtzelmüller and others (2000), while subsets of the data and further study-site details have been detailed previously (Reference HodsonHodson and others, 2010; Reference Irvine-Fynn, Bridge and HodsonIrvine-Fynn and others, 2010).

Concurrent meteorological data were collected at hourly intervals. At the camera location, time series of shortwave incident radiation, IR, and reflected radiation were recorded using horizontal, cross-calibrated up- and down-looking hemispheric pyranometers (Middleton and Skye Instruments, respectively). Supplementary meteorological data including air temperature, T a, were taken from the Gruvefjellet automatic weather station (AWS) located ∼4 km north of Longyearbreen but at an elevation similar to that of the imaging site.

Digital image processing

The purpose of image pre-processing was to provide a dataset normalized to a common spatial scale to account for changes in the position of the camera relative to the ablating ice surface. Analysis of this processed dataset was conducted with two objectives: (1) quantification of changes in the fractional areal coverage and distribution of cryoconite over the plot at hourly and diurnal timescales; and (2) quantification of individual particle motion vectors at hourly resolution over the course of each day. We used the public domain Java-based image processing software ImageJ (W.S. Rasband, http://imagej.nih.gov/ij/).

Pre-processing

Image pre-processing followed the method outlined by Reference Irvine-Fynn, Bridge and HodsonIrvine-Fynn and others (2010) for supraglacial cryoconite. Briefly, RGB colour images were split into red, green and blue colour channels, and blue channel data (as 8-bit, greyscale) retained for analysis since, in this colour space, red/brown cryoconite appears darker (lower pixel values) than surrounding blue, white and grey ice and shadows. The 216 images were grouped into 24-image ‘stacks’, each representing 1 day (midnight to midnight). Subsequent analyses were conducted on each 24-image stack separately.

Rescaling to normalize spatial scale

The spatial resolution of each image (pixels mm−1) was obtained by recording the apparent diameter (in pixels) of the 47 mm filter disc. The difference between this value and a target resolution of 10 pixels mm−1 was calculated and used as a scaling factor to normalize all images to the same size and spatial resolution. This step is essential for monitoring relative changes in position of cryoconite within the field-ofview over time. The final processed images showed a region of interest (ROI) centred on the centre of each original camera image, measuring ∼0.2 m × 0.2 m. This process does not account for changes in position of the field-of-view due to unconstrained camera ‘wobble’. To characterize this, we tracked and triangulated the positions of the five least-mobile cryoconite granules (see ‘Particle motion’ below) through each 24-image stack to obtain the apparent motion of a ‘virtual’ fixed reference point. The mean displacement of this point during any 24 hour period was only 2.04 ± 1.04 mm, with the exception of DOY 208 when a manual adjustment in camera position resulted in ∼36 mm displacement of the field of view. In this instance, this latter displacement vector was subtracted from subsequent particle position data.

Spatial moments

Regions of cryoconite were classified using pixel thresholding as described by Reference Irvine-Fynn, Bridge and HodsonIrvine-Fynn and others (2010). The total apparent cryoconite area (ACA; m2 cryoconite / m2 plot area) was evaluated as the fraction of pixels classified as cryoconite in each rescaled image. For each image, the centre of mass (COM) and XY covariance (a measure of dispersion of mass about the centroid) were calculated using built-in ImageJ functions: these two measures are the first and second spatial moments, respectively, and were used to quantify the net movement and distribution of cryoconite mass across the field of view.

Image differencing

To visualize changes in the spatial distribution of cryoconite across the plot over time, images obtained at 1200 h on DOY 203, 205, 207, 209 and 211 were processed to show areas of cryoconite loss and gain. Using the ‘Image arithmetic’ function in ImageJ, data for DOY 205 were subtracted, pixel-by-pixel, from data at DOY 203, resulting in an image in which negative pixel values indicated change from cryoconite (dark) to ice (bright), positive pixel values indicated change from ice to cryoconite, and zero values indicated regions of no substantial change. The procedure was repeated for successive image pairs at 48 hour intervals to DOY 211.

Particle motion

In each 24-image stack, the movement of a random selection of 50 individual cryoconite granules was tracked using the ImageJ plug-in MTrackJ (E. Meijering, http://www.imagescience.org/meijering/software/mtrackj/). The origin, hourly displacement (mm) and total displacement from the origin (mm) of each granule were recorded and averaged to obtain characteristic particle-motion data at hourly and diurnal timescales. From these data, motion path-length and -direction vectors could be derived.

Results and Analysis

Temporal variation in cryoconite areal coverage

Figure 1 illustrates the temporal change in ACA compared with air temperature, T a, and incident solar radiation, IR, the primary meteorological variables assumed to drive cryoconite mobility through ice melt. Note that ACA for a subset of these data was previously reported but used a subtly larger ROI within the measurement plot, so values of ACA in comparable images are slightly different (mean difference in ACA is 2.6 ± 1.3 m2 m−2 (cf. Reference Irvine-Fynn, Bridge and HodsonIrvine-Fynn and others, 2010).

Fig. 1. Hourly time series of ACA (black curve) including uncertainty (greyed area), incident radiation at the site (dotted curve) and air temperature at similar elevation (data from Gruvefjellet AWS: grey curve).

The ACA showed substantial fluctuations (maxima ∼0.195 and minima 0.058) with ∼4 day periodicity. Correlations between ACA and IR or T a were weak (r = −0.002 and r = −0.02, respectively), and covariance between IR and T a was only moderate (r = 0.4). Improvements in the correlation coefficients were found through cross-correlation: ACA lagged IR (r = 0.27, lag 14 hours) and T a (r = 0.37, lag 13 hours). Within the 9 day period of observations, ACA showed a statistically significant decrease (r 2 = 0.31, p < 0.01). However, given the relatively large fluctuation within this period, it is impossible to extrapolate a longer-term trend on this basis.

Spatial variation in cryoconite areal coverage

Figure 2 shows the change in the spatial distribution of cryoconite within the ROI for sequential 48 hour intervals processed using image differencing as described above. These data emphasize the substantial reworking of the cryoconite deposit within the plot, despite the fact that ACAs measured at the start and end of the observation period are only marginally different. Figure 3a plots the COM position over time and clearly indicates a trend in cryoconite movement at 57.1° from the downslope direction. Averaged over the observation period, the daily rate of COM motion is 5.3 mm, of which 2.9 mm is downslope (Y-direction) and 4.4 mm cross-slope (X-direction). The concurrent trend in the second spatial moment away from zero (Fig. 3b) indicates that cryoconite became more distributed or dispersed across the ROI, although like ACA this dataset showed strong fluctuations on hourly and diurnal timescales and the trend within the observation period cannot be extrapolated to longer timescales.

Fig. 2. Difference images highlighting spatial changes in cryoconite distribution at 48 hour intervals showing substantial reworking of the deposits. Note that downslope is parallel to the Y-axis, top to bottom.

Fig. 3. Plots of the change in (a) COM for cryoconite within the ∼0.04 m2 ROI, with each day highlighted with differing marker styles, and (b) the XY covariance indicating the degree of cryoconite dispersion over the ROI.

Motion of individual cryoconite granules

The motion tracks recorded for individual cryoconite granules were used to generate an hourly time series of mean particle motion rate (mm h−1). Plotted alongside the IR time series (Fig. 4), these data show an apparently diurnal signal. However, the correlation between the datasets, although significant, was weak (r = 0.25), as was the correlation with T a (r = 0.27; data not shown). Cross-correlation analysis suggested a lag time of 1 hour for mean granule movement against both IR and T a (r increased by 0.02 for both variables). Accumulated over each day, the net particle displacement of the 50 measured granules is shown in Figure 5a. Over the period of observations, this mean daily particle displacement was 8.6 mm, but the orientation of granule movement vectors showed no preferential direction, either downslope (180°) or in the direction of net cryoconite COM movement shown in Figure 3a (123°). However, the magnitude of mean daily particle displacement did, broadly, appear to correspond to the presumed driving variable, IR (Fig. 5b).

Fig. 4. Mean particle displacement as a function of time for sample of 50 random cryoconite granules; the coincident series of IR at the site is shown for comparison (grey dotted curve). Note that the velocity series is not continuous: data were examined at 24 hour intervals (n = 23 for 1 day), and images for which rescaling uncertainty occurred were omitted to avoid introduction of bias.

Fig. 5. (a) Compass plots illustrating net daily motion (magnitude and direction) for 50 cryoconite granules for DOY 203–211. Outer plot limits are 25 mm, with a 12.5 mm axis marked; downslope is oriented at 180°. (b) Daily sum of hourly-averaged IR data (bars) plotted alongside daily mean particle displacement (curve).

Cryoconite mobility and supraglacial energy balance

A full energy-balance model (EBM) estimating hourly energy fluxes to the ice surface was run following Reference Brock and ArnoldBrock and Arnold (2000). Site-specific details were used and adjustments made to the incident solar radiation algorithm to account for high-latitude solar geometry (after Reference Kreith and KreiderKreith and Kreider, 1978). The time series of albedo (derived from the two pyrometer records) was corrected for surface slope effects following Reference ArendtArendt (1999) and included within the model; a static value of 0.0253 mm was employed as a representative ice roughness (Reference Hodson, Kohler and BrinkhausHodson and others, 2005b). Particle motion compared with energy-balance components at hourly timescales revealed correlations |r| < 0.02 for sensible and latent heat fluxes, and an inverse relationship with long-wave radiative flux (r = −0.14); a stronger positive correlation between motion and shortwave radiation is apparent (r = 0.25; no lag was found with cross-correlation analysis), and regression analysis showed the slope of the relationship was significantly different from zero (p = 0.05). Figure 6 suggests that a number of outlying data points may conceal a more robust trend. These analyses, albeit with weak associations, support the assertion that, at the hourly timescale, shortwave radiation may provide a driving force for cryoconite granule redistribution.

Fig. 6. Scatter plot of hourly particle motion (from a sample of 50 granules) vs shortwave radiation flux predicted using an adjusted Reference Brock and ArnoldBrock and Arnold (2000) energy-balance model.

Crucially, comparison between the ACA time series and the plot-scale albedo measurements (Fig. 7) shows a significant relationship between the two variables over the 9 day period of observation (r 2 = 0.243, p < 0.05). This finding associates the process of relief inversions (the apparent melting in and out) and the fluctuations in cryoconite coverage and dispersion with ice surface albedo. The wider implication of this association is that the short-term variability in ACA may be a significant contributing factor in defining short-term variability in the energy balance over the snow-free glacier ice and thus provide some control on the variability in melt rates.

Fig. 7. Scatter plot of ACA vs albedo (the ratio of reflected/incident radiation) corrected for the effect of surface slope.

Discussion and Conclusions

We have used quantitative analysis of time-lapse digital image data to examine, in detail, a small (sub-metre) patch of dispersed cryoconite at unprecedented spatial and temporal resolution over a 9 day period. Cryoconite coverage varied substantially at frequencies of several days, which corresponds to typical synoptic weather pattern periodicity in western Svalbard (e.g. Reference Hanssen-Bauer, Solås and SteffensenHanssen-Bauer and others, 1990; Reference Hodson, Gurnell, Washington, Tranter, Clark and HagenHodson and others, 1998). Although, statistically, relationships were weak, the analyses suggest some degree of meteorological forcing of the ACA, with an inverse and 1 hour lagged relationship between energy receipt and ACA, suggestive of relief inversion processes at the diurnal scale. The fluctuation in ACA is coupled with a progressive but rather small lateral movement in the centre of mass for the cryoconite across the plot. Assuming a typical 70 day High Arctic melt season (Reference Hodson, Tranter and VatneHodson and others, 2000), this movement is insignificant at the scale of the glacier, being equivalent to 0.37 m a−1. Therefore, based on these preliminary data from a typical patch of dispersed cryoconite, and knowledge of distinctly linear relationships between cryoconite area and mass (Reference Cook, Hodson, Telling, Anesio, Irvine-Fynn and BellasCook and others, 2010), we find little evidence for mass transfer of significant quantities of cryoconite across the ice surface by ‘washing’ by meltwater flow over (or through) the so-called ‘weathering crust’ (Reference Müller and KeelerMüller and Keeler, 1969). This stability has significant implications for the cryoconite ecosystem, in that it increases the effective residence time of cryoconite on the glacier surface, particularly compared with meltwater fluxes, thereby enabling the formation of complex microbial communities, which are tolerant of the prevailing conditions (Reference EdwardsEdwards and others, 2011).

Conversely, the motion of individual cryoconite particles within the patch was pervasive, persistent and apparently random. Although the daily mean particle displacement was of the same order as the net cryoconite movement (8.6 and 5.3 mm d−1, respectively), the mean distance travelled by granules was significantly larger, at 18.3 mm d−1, emphasizing the tortuosity of granule motion. Inspection of individual motion tracks (Fig. 8) indicates that motion was dominated by two processes: random ‘shuffling’ of particles (leading to very small net displacement) and larger, rapid, directed displacements of particles associated with melting-induced flows, often between neighbouring holes. Critically, then, the major control on granule motion appears to be the microtopography of the weathering crust and the relative elevations of perched water levels generated and then released by progressive (or fluctuating) processes of relief inversion and local hydraulics.

Fig. 8. Cryoconite distribution at 0000 and 2300 h on DOY 207, showing initial positions (a) and motion tracks (b) of selected cryoconite granules. Of these, a subset (labelled 2:n and circled in black) show small, random ‘shuffling’ over the 23 hour period. Others (labelled 1:n) are associated with significant reworking of cryoconite hole morphology.

Importantly, the observations presented here relate to a patch of associated but distributed cryoconite holes, which contrasts in form with other distributions of cryoconite, such as in-stream and fully dispersed cryoconite (Reference HodsonHodson and others, 2007), deep, isolated holes, or the larger cryoconite deposits observed in the field as ‘cryoconite mantles’ (Reference HodsonHodson and others, 2008) or lobate forms. Additional data are required to characterize the morpho-dynamics of cryoconite in these forms, and over longer timescales, before more general estimates of cryoconite redistribution and transport phenomena at the glacier scale can be drawn. However, the methods we present here provide a rapid, reproducible protocol for collecting and analysing such datasets.

A key disadvantage of the current imaging method is the lack of quantitative depth perception from a single camera. Small-scale photogrammetry is under development to address this issue and enable more detailed quantitative investigation of the development of ice surface micro-topography and cryoconite coverage over the course of the melt season. Our preliminary analyses here indicate that this is a significant driver of cryoconite granule redistribution and hence of patch-scale ACA and, more significantly, glacier surface albedo fluctuations. Quantifying and understanding the small-scale, diurnal interactions between cryoconite granules and microtopography, using time-lapse image sequences, may thus critically enhance our ability to model the dynamic coupling between cryoconite, surface roughness, albedo and melt generation over larger areas, at seasonal timescales.

Acknowledgements

Fieldwork was funded by the Leverhulme Trust Research Fellowship RF/4/RFG/2007/0398 awarded to A.J.H. T.I.-F. and A.J.H. acknowledge UK Natural Environment Research Council (NERC) Standard Grant NE/G006253/1, and the RGS Peter Fleming Award (2009) enabled thoughtful discussions in the field. T.I.-F. also recognizes backing from the Climate Change Consortium of Wales (C3W) initiative, enabling completion of the manuscript. J.W.B. was supported in this work by S. Banwart and NERC grant NE/ C521001/1. The Gruvefjellet meteorological data were made available by researchers based at The University Centre in Svalbard (UNIS). A. Edwards is thanked for useful comments on the manuscript. Comments from two anonymous reviewers are appreciated.

References

Anesio, A.M., Mindl, B., Laybourn-Parry, J., Hodson, A.J. and Sattler, B.. 2007. Viral dynamics in cryoconite holes on a high Arctic glacier (Svalbard). J. Geophys. Res., 112(G4), G04S31. (10.1029/2006JG000350.)Google Scholar
Anesio, A.M., Hodson, A.J., Fritz, A., Psenner, R. and Sattler, B.. 2009. High microbial activity on glaciers: importance to the global carbon cycle. Global Change Biol., 15(4), 955960.Google Scholar
Anesio, A.M. and 6 others. 2010. Carbon fluxes through bacterial community on glacier surfaces. Ann. Glaciol., 51(56), 3240.CrossRefGoogle Scholar
Arendt, A. 1999. Approaches to modelling the surface albedo of a high Arctic glacier. Geogr. Ann., 81A(4), 477487.CrossRefGoogle Scholar
Bøggild, C.E. 1997. Different melt regimes indicated by surface albedo measurements at the Greenland ice sheet margin – application of TM image. EARSeL Adv. Remote Sens., 5, Yearbook 1997, 8288.Google Scholar
Bøggild, C.E., Brandt, R.E., Brown, K.J. and Warren, S.G.. 2010. The ablation zone in northeast Greenland: ice types, albedos and impurities. J. Glaciol., 56(195), 101113.Google Scholar
Brock, B.W. and Arnold, N.S.. 2000. A spreadsheet-based (Microsoft Excel) point surface energy balance model for glacier and snowmelt studies. Earth Surf. Process. Landf., 25(6), 649658.Google Scholar
Christner, B.C., Kvito, B.H. and Reeve, J.N.. 2003. Molecular identification of bacteria and eukarya inhabiting an Antarctic cryoconite hole. Extremophiles, 7(3), 177183.CrossRefGoogle ScholarPubMed
Cook, J., Hodson, A., Telling, J., Anesio, A., Irvine-Fynn, T. and Bellas, C.. 2010. The mass–area relationship within cryoconite holes and its implications for primary production. Ann. Glaciol., 51(56), 106110.CrossRefGoogle Scholar
De Smet, W.H. and van Rompu, E.A.. 1994. Rotifera and Tardigrada from some cryoconite holes on a Spitsbergen (Svalbard) glacier. Belg. J. Zool., 124(1), 2737.Google Scholar
Edwards, A. and 7 others. 2011. Possible interactions between bacterial diversity, microbial activity and supraglacial hydrology of cryoconite holes in Svalbard. ISME J., 5(1), 150160.Google Scholar
Etzelmüller, B., Ödegård, R.S., Vatne, G., Mysterud, R.S., Tonning, T. and Sollid, J.L.. 2000. Glacier characteristics and sediment transfer system of Longyearbreen and Larsbreen, western Spitsbergen. Nor. Geogr. Tidsskr., 54(4), 157168.Google Scholar
Foreman, C.M., Sattler, B., Mikucki, D.L., Porazinska, D.L. and Priscu, J.C.. 2007. Metabolic activity and diversity of cryoconites in the Taylor Valley, Antarctica. J. Geophys. Res., 112(G4), G04S32. (10.1029/2006JG000358.)Google Scholar
Fountain, A.G., Tranter, M., Nylen, T.H., Lewis, K.J. and Mueller, D.R.. 2004. Evolution of cryoconite holes and their contribution to meltwater runoff from glaciers in the McMurdo Dry Valleys, Antarctica. J. Glaciol., 50(168), 3545.Google Scholar
Hanssen-Bauer, I., Solås, M.K. and Steffensen, E.L.. 1990. The climate of Spitsbergen. Oslo, Norwegian Meteorological Institute. (DNMI Research Report 39/90.)Google Scholar
Hodson, A.J., Gurnell, A.M., Washington, R., Tranter, M., Clark, M.J. and Hagen, J.O.. 1998. Meteorological and runoff time-series characteristics in a small, high-Arctic glaciated basin, Svalbard. Hydrol. Process., 12(3), 509526.Google Scholar
Hodson, A.J., Tranter, M. and Vatne, G.. 2000. Contemporary rates of chemical denudation and atmospheric CO2 sequestration in glacier basins: an Arctic perspective. Earth Surf. Process. Landf., 25(13), 14471471.Google Scholar
Hodson, A.J., Mumford, P.N., Kohler, J. and Wynn, P.M.. 2005a. The High Arctic glacial ecosystem: new insights from nutrient budgets. Biogeochemistry, 72(2), 233256.Google Scholar
Hodson, A., Kohler, J. and Brinkhaus, M.. 2005b. Multi-year water and surface energy budget of a high-latitude polythermal glacier: evidence for overwinter storage in a dynamic subglacial reservoir. Ann. Glaciol., 42, 4246.Google Scholar
Hodson, A.J. and 10 others. 2007. A glacier respires: quantifying the distribution and respiration CO2 flux of cryoconite across an entire Arctic supraglacial ecosystem. J. Geophys. Res., 112(G4), G04S36. (10.1029/2007JG000452.)Google Scholar
Hodson, A. and 7 others. 2008. Glacial ecosystems. Ecol. Monogr., 78(1), 4167.CrossRefGoogle Scholar
Hodson, A. and 6 others. 2010. The structure, biological activity and biogeochemistry of cryoconite aggregates upon an Arctic valley glacier: Longyearbreen, Svalbard. J. Glaciol., 56(196), 349362.Google Scholar
Hood, E. and 6 others. 2009. Glaciers as a source of ancient and labile organic matter to the marine environment. Nature, 462(7276), 10441047.CrossRefGoogle Scholar
Irvine-Fynn, T.D.L. 2008. Modelling runoff from the maritime Arctic cryosphere: water storage and routing at Midtre Lovénbreen. (PhD thesis, University of Sheffield.)Google Scholar
Irvine-Fynn, T.D.L., Bridge, J.W. and Hodson, A.J.. 2010. Rapid quantification of cryoconite: granule geometry and in situ supraglacial extents, using examples from Svalbard and Greenland. J. Glaciol., 56(196), 297308.Google Scholar
Kreith, F. and Kreider, J.F.. 1978. Principles of solar engineering. New York, McGraw-Hill.Google Scholar
Langford, H., Hodson, A., Banwart, S. and Bøggild, C.. 2010. The microstructure and biogeochemistry of Arctic cryoconite granules. Ann. Glaciol., 51(56), 8794.Google Scholar
MacDonell, S. and Fitzsimons, S.. 2008. The formation and hydrological significance of cryoconite holes. Progr. Phys. Geogr., 32(6), 595610.CrossRefGoogle Scholar
McIntyre, N.F. 1984. Cryoconite hole thermodynamics. Can. J. Earth Sci., 21(2), 152156.Google Scholar
Mueller, D.R. and Pollard, W.H.. 2004. Gradient analysis of cryoconite ecosystems from two polar glaciers. Polar Biol., 27(2), 6674.Google Scholar
Müller, F. and Keeler, C.M.. 1969. Errors in short-term ablation measurements on melting ice surfaces. J. Glaciol., 8(52), 91105.Google Scholar
Podgorny, I.A. and Grenfell, T.C.. 1996. Absorption of solar energy in a cryoconite hole. Geophys. Res. Lett., 23(18), 24652468.Google Scholar
Porazinska, D.L., Fountain, A.G., Nylen, T.H., Tranter, M., Virginia, R.A. and Wall, D.H.. 2004. The biodiversity and biogeochemistry of cryoconite holes from McMurdo Dry Valley glaciers, Antarctica. Arct. Antarct. Alp. Res., 36(1), 8491.Google Scholar
Säwström, C., Mumford, P., Marshall, W., Hodson, A. and Laybourn-Parry, J.. 2002. The microbial communities and primary productivity of cryconite holes in an Arctic glacier (Svalbard 79°N). Polar Biol., 25(8), 591596.Google Scholar
Segawa, T., Takeuchi, N., Ishida, K., Kanda, H. and Kohshima, S.. 2010. Altitudinal changes in a bacterial community on Gulkana Glacier in Alaska. Microbes Environ., 25(3), 171182.Google Scholar
Stibal, M., Šabacká, M. and Kaštovská, K.. 2006. Microbial communities on glacier surfaces in Svalbard: impact of physical and chemical properties on abundance and structure of cyanobacteria and algae. Microbial Ecol., 52(4), 655–654.CrossRefGoogle ScholarPubMed
Stibal, M., Tranter, M., Benning, L.G. and Rehák, J.. 2008. Microbial primary production on an Arctic glacier is insignificant in comparison with allochthonous organic carbon input. Environ. Microbiol., 10(8), 21722178.Google Scholar
Takeuchi, N. 2009. Temporal and spatial variations in spectral reflectance and characteristics of surface dust on Gulkana Glacier, Alaska Range. J. Glaciol., 55(192), 701709.Google Scholar
Takeuchi, N., Kohshima, S., Goto-Azuma, K. and Koerner, R.M.. 2001a. Biological characteristics of dark colored material (cryoconite) on Canadian Arctic glaciers (Devon and Penny ice caps). Mem. Natl Inst. Polar Res., Special Issue 54, 495505.Google Scholar
Takeuchi, N., Kohshima, S. and Seko, K.. 2001b. Structure, formation, and darkening process of albedo-reducing material (cryoconite) on a Himalayan glacier: a granular algal mat growing on the glacier. Arct. Antarct. Alp. Res., 33(2), 115122.Google Scholar
Takeuchi, N., Nishiyama, H. and Li, Z.. 2010. Structure and formation process of cryoconite granules on Ürümqi glacier No. 1, Tien Shan, China. Ann. Glaciol., 51(56), 914.Google Scholar
Figure 0

Fig. 1. Hourly time series of ACA (black curve) including uncertainty (greyed area), incident radiation at the site (dotted curve) and air temperature at similar elevation (data from Gruvefjellet AWS: grey curve).

Figure 1

Fig. 2. Difference images highlighting spatial changes in cryoconite distribution at 48 hour intervals showing substantial reworking of the deposits. Note that downslope is parallel to the Y-axis, top to bottom.

Figure 2

Fig. 3. Plots of the change in (a) COM for cryoconite within the ∼0.04 m2 ROI, with each day highlighted with differing marker styles, and (b) the XY covariance indicating the degree of cryoconite dispersion over the ROI.

Figure 3

Fig. 4. Mean particle displacement as a function of time for sample of 50 random cryoconite granules; the coincident series of IR at the site is shown for comparison (grey dotted curve). Note that the velocity series is not continuous: data were examined at 24 hour intervals (n = 23 for 1 day), and images for which rescaling uncertainty occurred were omitted to avoid introduction of bias.

Figure 4

Fig. 5. (a) Compass plots illustrating net daily motion (magnitude and direction) for 50 cryoconite granules for DOY 203–211. Outer plot limits are 25 mm, with a 12.5 mm axis marked; downslope is oriented at 180°. (b) Daily sum of hourly-averaged IR data (bars) plotted alongside daily mean particle displacement (curve).

Figure 5

Fig. 6. Scatter plot of hourly particle motion (from a sample of 50 granules) vs shortwave radiation flux predicted using an adjusted Brock and Arnold (2000) energy-balance model.

Figure 6

Fig. 7. Scatter plot of ACA vs albedo (the ratio of reflected/incident radiation) corrected for the effect of surface slope.

Figure 7

Fig. 8. Cryoconite distribution at 0000 and 2300 h on DOY 207, showing initial positions (a) and motion tracks (b) of selected cryoconite granules. Of these, a subset (labelled 2:n and circled in black) show small, random ‘shuffling’ over the 23 hour period. Others (labelled 1:n) are associated with significant reworking of cryoconite hole morphology.