Iceberg Calving Research Articles

A 3.3‐Million‐Year Record of Antarctic Iceberg Rafted Debris and Ice Sheet Evolution Quantified by Machine Learning

AbstractOver the last 3.3 million years, the Antarctic Ice Sheet (AIS) has undergone phases of ice sheet growth and decay, impacting sea level and climate globally. Presently, the largely marine‐terminating AIS loses mass primarily by iceberg calving and basal melt of ice shelves. Quantifying past rates and timing of AIS melt is vital to understanding future cryosphere and sea level changes. One proxy for past ice sheet instabilities is iceberg rafted debris (IRD) fluxes. However, traditional methods of IRD quantification are labor‐intensive. Here, we present a new method of identifying IRD grains in sediment core X‐ray images using a convolutional neural network machine learning algorithm. We present a 3.3‐million‐year record of AIS IRD melt events using sediment cores from International Ocean Discovery Program Sites U1536, U1537, and U1538 in the Southern Ocean's “Iceberg Alley.” We identify two increases in the IRD fluxes throughout this period, at ∼1.8 and 0.43 Ma. We propose that after 1.8 Ma, the AIS expanded and transitioned from a primarily terrestrial‐terminating to a primarily marine‐terminating ice sheet. Therefore, after 1.8 Ma, glacial terminations and AIS iceberg discharge are associated with variations in global ice volume, presumably through the mechanism of sea level and, therefore, grounding line change. The second AIS regime change occurs during the Mid‐Brunhes Event (∼0.43 Ma). After this time, there are heightened and continuous IRD fluxes at each glacial termination, indicating increased AIS size and instability after this time.

Calving front monitoring at a subseasonal resolution: a deep learning application for Greenland glaciers

Abstract. The mass balance of the Greenland Ice Sheet is strongly influenced by the dynamics of its outlet glaciers. Therefore, it is of paramount importance to accurately and continuously monitor these glaciers, especially the variation in their frontal positions. A temporally comprehensive parameterization of glacier calving is essential for understanding dynamic changes and constraining ice sheet modeling. However, many current calving front records are limited in terms of temporal resolution as they rely on manual delineation, which is laborious and not appropriate considering the increasing amount of satellite imagery available. In this contribution, we address this problem by applying an automated method to extract calving fronts from optical satellite imagery. The core of this workflow builds on recent advances in the field of deep learning while taking full advantage of multispectral input information. The performance of the method is evaluated using three independent test datasets. For the three datasets, we calculate mean delineation errors of 61.2, 73.7, and 73.5 m, respectively. Eventually, we apply the technique to Landsat-8 imagery. We generate 9243 calving front positions across 23 outlet glaciers in Greenland for the period 2013–2021. Resulting time series not only resolve long-term and seasonal signals but also resolve subseasonal patterns. We discuss the implications for glaciological studies and present a first application for analyzing the effect of bedrock topography on calving front variations. Our method and derived results represent an important step towards the development of intelligent processing strategies for glacier monitoring, opening up new possibilities for studying and modeling the dynamics of Greenland's outlet glaciers.

A Boulder Beach Formed by Waves From a Calving Glacier Revisited: Multidecadal Tsunami–Controlled Coastal Changes in Front of Eqip Sermia, West Greenland

ABSTRACTThe calving of glaciers regularly produces tsunami‐like waves that pose a serious threat to coastal environments. Those strong waves are not only able to move ice mélange and redistribute icebergs, growlers, or sea ice across a fjord but also flood and remodel neighbouring cliffs and beaches. Here, we analyze over 90 years (1929–2023) of coastal zone changes that occurred in front of Eqip Sermia. We show that calving waves play a dominant role in transforming the lateral moraine and forming a beach and spit system south of the glacier front. Part of the former moraine has transformed into a boulder‐dominated spit, which closed the lagoon over the years. By multidecadal analysis, we also detected a significant erosion of unconsolidated cliffs located on the opposite side of the bay (~0.53 m per year between 1985 and 2023). In addition, we demonstrate that even a single event (one calving wave) can remodel a beach surface by entrainment of up to 1.8‐m‐diameter boulders and the erosion of the beach surface by washing away sand and gravel from rocky outcrops. Our study constitutes important progress toward modes of paraglacial coastal evolution in regions characterized by rapidly retreating calving glaciers.

Heinrich event ice discharge and the fate of the Atlantic Meridional Overturning Circulation.

During Heinrich events, great armadas of icebergs episodically flooded the North Atlantic Ocean and weakened overturning circulation. The ice discharges of these episodes constrain the sensitivity of overturning circulation to iceberg melting. We reconstructed these ice discharges to be as high as 0.13 sverdrup (Sv) (1 Sv = 1 million cubic meters per second) during Heinrich event 4 and to average 0.029 Sv over all episodes. The present-day Greenland Ice Sheet calving of icebergs is comparable to that of a mid-range Heinrich event. As the future Greenland Ice Sheet recedes from marine-terminating outlets, its iceberg calving likely will not persist long enough for icebergs alone to cause catastrophic disruption to the Atlantic overturning circulation, although the accelerating Greenland runoff and continued global warming remain threats to the circulation stability.

A large-scale high-resolution numerical model for sea-ice fragmentation dynamics

Abstract. Forecasts of sea-ice motion and fragmentation are of vital importance for all human interactions with sea ice, ranging from those involving indigenous hunters to shipping in polar regions. Sea-ice models are also important for simulating long-term changes in a warming climate. Here, we apply the Helsinki Discrete Element Model (HiDEM), originally developed for glacier calving, to sea-ice breakup and dynamics. The code is highly optimized to utilize high-end supercomputers to achieve an extreme time and space resolution. Simulated fracture patterns and ice motion are compared with satellite images of the Kvarken region of the Baltic Sea from March 2018. A second application of HiDEM involves ice ridge formation in the Gulf of Riga. With a few tens of graphics processing units (GPUs), the code is capable of reproducing observed ice patterns that in nature may take a few days to form; this is done over an area of ∼100km×100km, with an 8 m resolution, in computations lasting ∼10 h. The simulations largely reproduce observed fracture patterns, ice motion, fast-ice regions, floe size distributions, and ridge patterns. The similarities and differences between observed and computed ice dynamics and their relation to initial conditions, boundary conditions, and applied driving forces are discussed in detail. The results reported here indicate that the HiDEM has the potential to be developed into a detailed high-resolution model for sea-ice dynamics at short timescales, which, when combined with large-scale and long-term continuum models, may form an efficient framework for forecasts of sea-ice dynamics.

A 3D glacier dynamics–line plume model to estimate the frontal ablation of Hansbreen, Svalbard

Abstract. Frontal ablation is responsible for a large fraction of the mass loss from tidewater glaciers. The main contributors to frontal ablation are iceberg calving and submarine melting, with calving often being the largest. However, submarine melting, in addition to its direct contribution to mass loss, also promotes calving through the changes induced in the stress field at the glacier terminus, so both processes should be jointly analysed. Among the factors influencing submarine melting, the formation of a buoyant plume due to the emergence of fresh subglacial water at the glacier grounding line plays a key role. In this study we used Elmer/Ice to develop a 3D glacier dynamics model including calving and subglacial hydrology coupled with a line plume model to calculate the calving front position at every time step. We applied this model to the Hansbreen–Hansbukta glacier–fjord system in southern Spitsbergen, Svalbard, where a large set of data are available for both the glacier and the fjord from September 2008 to March 2011. We found that our 3D model reproduced the expected seasonal cycle of advance–retreat. Besides, the modelled front positions were in good agreement with the observed front positions at the central part of the calving front, with longitudinal differences, on average, below 15 m for the period from December 2009 to March 2011. But there were regions of the front, especially the eastern margin, that presented major differences.

Environment‐Modulated Glacial Seismicity Near Dålk Glacier in East Antarctica Revealed by Deep Clustering

AbstractOver the past decades, seismic monitoring has been increasingly used to track glacial activities associated with ice loss. Many seismological studies focus on West Antarctica, whereas glacial seismicity in East Antarctica is much less studied. Here, we apply unsupervised deep learning to a dense nodal seismic array near Dålk Glacier, East Antarctica, operating from 6 December 2019 to 2 January 2020. An autoencoder is used to automatically extract event features, which are then input into a Gaussian mixture model for clustering. We divide the data into 50 clusters and merge them according to their temporal and spectral characteristics. The results reveal five main types of seismic signals: two groups with monochromatic and high frequencies, two groups with broadband frequency and short duration, and a group with mainly low frequency and long duration. By comparing the environmental conditions (wind, temperature and tides), we infer that the two monochromatic groups are wind‐induced vibrations of the near‐station flag markers and topography; the two broadband groups are likely thermal contractions on the blue ice surface and stick‐slip events at the ice base; and the low‐frequency events are water‐filled basal crevassing and iceberg calving. In particular, we observe one type of low‐frequency event preceded by high‐frequency onset, which is likely basal crevassing near the grounding line of Dålk Glacier and predominantly occurred during rising tides. Our findings show that deep clustering is effective in identifying a wide range of glacial seismic events and can contribute to the rapid growth of passive glacier seismic monitoring.

Uncrewed aerial vehicle with onboard winch system for rapid, cost-effective, and safe oceanographic profiling in hazardous and inaccessible areas

Interactions between coastal waters and marine-terminating glaciers in the Polar Regions play a significant role in global sea level rise fueled by a rapidly warming Arctic.The risk of glacier calving, and the abundance of ice, can make it impossible for surface vessels to access the waters near glacier termini. Alternative methods using manned aircraft are expensive. As a result, oceanographic measurements are limited near glacier termini.We present an uncrewed aerial vehicle (UAV) with an on-board winch system that allows oceanographic profiling in remote, hazardous areas using a commercial conductivity, temperature, and depth (CTD) sensor payload. The UAV is optimized for easy handling and deployment and is capable of high-speed and efficient cruise flight. An autopilot system provides pilot assistance and autonomous flight capabilities. The total weight of the UAV including payload is 6.5 kg with an endurance of 24 min.Testing of the system was conducted in South Greenland during winter conditions in March 2023 with successful profiles collected near a glacier terminus (<5 m) and in small openings in ice mélange (2.2 m). The system proved capable, reliable, and efficient. Further development of the system will allow other sensors for an even more flexible measurement suite.

Discrimination between icequakes and earthquakes in southern Alaska: an exploration of waveform features using Random Forest algorithm

SUMMARY This study examines the feature space of seismic waveforms often used in machine learning applications for seismic event detection and classification problems. Our investigation centres on the southern Alaska region, where the seismic record captures diverse seismic activity, notably from the calving of marine-terminating glaciers and tectonic earthquakes along active plate boundaries. While the automated discrimination of earthquakes and glacier quakes is our nominal goal, this data set provides an outstanding opportunity to explore the general feature space of regional seismic phases. That objective has applicability beyond ice quakes and our geographic region of study. We make a noteworthy discovery that features rooted in the spectral content of seismic waveforms consistently outperform statistical and temporal features. Spectral features demonstrate robust performance, exhibiting resilience to class imbalance while being minimally impacted by factors such as epicentral distance and signal-to-noise ratio. We also conduct experiments on the transferability of the model and find that transferability primarily depends on the appearance of the waveforms. Finally, we analyse misclassified events and find examples that are identified incorrectly in the original regional catalogue.

Stability of Ice Shelves and Ice Cliffs in a Changing Climate

The largest uncertainty in future sea-level rise is loss of ice from the Greenland and Antarctic Ice Sheets. Ice shelves, freely floating platforms of ice that fringe the ice sheets, play a crucial role in restraining discharge of grounded ice into the ocean through buttressing. However, since the 1990s, several ice shelves have thinned, retreated, and collapsed. If this pattern continues, it could expose thick cliffs that become structurally unstable and collapse in a process called marine ice cliff instability (MICI). However, the feedbacks between calving, retreat, and other forcings are not well understood. Here we review observed modes of calving from ice shelves and marine-terminating glaciers, and their relation to environmental forces. We show that the primary driver of calving is long-term internal glaciological stress, but as ice shelves thin they may become more vulnerable to environmental forcing. This vulnerability—and the potential for MICI—comes from a combination of the distribution of preexisting flaws within the ice and regions where the stress is large enough to initiate fracture. Although significant progress has been made modeling these processes, theories must now be tested against a wide range of environmental and glaciological conditions in both modern and paleo conditions. ▪Ice shelves, floating platforms of ice fed by ice sheets, shed mass in a near-instantaneous fashion through iceberg calving.▪Most ice shelves exhibit a stable cycle of calving front advance and retreat that is insensitive to small changes in environmental conditions.▪Some ice shelves have retreated or collapsed completely, and in the future this could expose thick cliffs that could become structurally unstable called ice cliff instability.▪The potential for ice shelf and ice cliff instability is controlled by the presence and evolution of flaws or fractures within the ice.

Updating glacier inventories on the periphery of Antarctica and Greenland using multi-source data

Abstract Melting and calving of glaciers and ice caps in Antarctica and Greenland could potentially contribute significantly to global sea level rise. Updates to existing outlines that provide critical glacier baseline information in both regions could help in the analysis of particular changes in glacier parameters such as area and volume from time-series inventories. Here we synthesize previously established techniques and apply new multi-source datasets to update glacier outlines in selected test areas of Antarctica and Greenland, as well as to reduce uncertainties and errors during the mapping process. The workflow includes mapping glacier boundaries, subdividing glaciers by watersheds and assigning glacier attributes. Complicated glacier scenarios and updating challenges in polar regions are discussed and demonstrated by representative case studies. For the first time in Antarctica, we analyze the effect of terminus types on mapped glacier areas, and in Greenland we compare the differences with glacier mapping results using Landsat OLI and ETM+. With new data sources, the methods described in this study might help to create glacier outlines on a larger scale in Antarctica and Greenland. Although data sources can be substituted, the enormous amount of manual labor required to update glacier inventories remains a significant challenge.

Sedimentary stratigraphy and provenance off Dronning Maud Land (East Antarctica) during the mid-Pleistocene transition: Implications for paleoclimate and ice dynamics

Glacial-marine sediments from the Antarctic continental margin provide a record of depositional environment, oceanographic variability and ice dynamics that is tapped with scientific ocean drilling. This study focuses on Ocean Drilling Program Core 693A-2R, a 9.7 m sediment core retrieved from near the continental margin of the Archean Grunehogna Craton in Dronning Maud Land (DML), East Antarctica. The results contribute to a better understanding of ice-shelf behavior in DML during the mid-Pleistocene transition (MPT), a well-known transition from 40-kyr to 100-kyr cycle periods. The age model, constructed based on Sr isotope stratigraphy and geomagnetic reversals, indicates that the core spans 1.20 to 0.65 Ma. The dynamic behavior of DML ice shelves with periodic iceberg calving is revealed by the glacial–interglacial variation in sedimentation patterns, with interglacials characterized by higher concentrations of ice-rafted debris (IRD) associated with enhanced paleoproductivity than glacial intervals. The responses of DML ice shelves to warm climates are represented by a prolonged interglacial period at 1.0–1.1 Ma (MIS 31–27) and significant interglacial expressions during MIS 19 and 17. The 40Ar/39Ar ages of individual ice-rafted hornblende grains are compared with the on-land geology of DML and neighboring regions to determine the provenances of IRD. Specifically, 40Ar/39Ar results record primarily late Neoproterozoic to Cambrian ages (600–400 Ma) with a predominant peak of 520–480 Ma. This Pan-African/Ross orogeny signature is very common in East Antarctica but is not found in the most proximal margin of the Grunehogna Craton, and is instead associated with the region of DML several hundred kilometers east of the deposition site. This indicates that significant discharges of icebergs occurred in the remote DML, which were then transported by the westward-flowing Antarctic Coastal Current to deposit IRD at the studied site during the MPT. This study establishes a confirmed MPT sedimentary sequence off DML, against which future MPT proxy records from the Weddell Sea embayment and other sectors in Antarctica can be compared and correlated, and provides a basis for more detailed analyses of the response of DML ice sheet to Pleistocene climate variations.

Array processing in cryoseismology: a comparison to network-based approaches at an Antarctic ice stream

Abstract. Seismicity at glaciers, ice sheets, and ice shelves provides observational constraint on a number of glaciological processes. Detecting and locating this seismicity, specifically icequakes, is a necessary first step in studying processes such as basal slip, crevassing, imaging ice fabric, and iceberg calving, for example. Most glacier deployments to date use conventional seismic networks, comprised of seismometers distributed over the entire area of interest. However, smaller-aperture seismic arrays can also be used, which are typically sensitive to seismicity distal from the array footprint and require a smaller number of instruments. Here, we investigate the potential of arrays and array-processing methods to detect and locate subsurface microseismicity at glaciers, benchmarking performance against conventional seismic-network-based methods for an example at an Antarctic ice stream. We also provide an array-processing recipe for body-wave cryoseismology applications. Results from an array and a network deployed at Rutford Ice Stream, Antarctica, show that arrays and networks both have strengths and weaknesses. Arrays can detect icequakes from further distances, whereas networks outperform arrays in more comprehensive studies of a particular process due to greater hypocentral constraint within the network extent. We also gain new insights into seismic behaviour at the Rutford Ice Stream. The array detects basal icequakes in what was previously interpreted to be an aseismic region of the bed, as well as new icequake observations downstream and at the ice stream shear margins, where it would be challenging to deploy instruments. Finally, we make some practical recommendations for future array deployments at glaciers.

Monitoring glacier calving using underwater sound

Abstract. Climate shifts are particularly conspicuous in glaciated areas. Satellite and terrestrial observations show significant increases in the melting and breakup of tidewater glaciers and their influence on sea level rise and ocean mixing. Increasing melt rates are creating an urgency to better understand the link between atmospheric and oceanic conditions and glacier frontal ablation through iceberg calving and melting. Elucidating this link requires a combination of short- and long-timescale measurements of terminus activity. Recent work has demonstrated the potential of using underwater sound to quantify the time and scale of calving events to yield integrated estimates of ice mass loss (Glowacki and Deane, 2020). Here, we present estimates of subaerial calving flux using underwater sound recorded at Hansbreen, Svalbard, in September 2013 combined with an algorithm for the automatic detection of calving events. The method is compared with ice calving volumes estimated from geodetic measurements of the movement of the glacier terminus and an analysis of satellite images. The total volume of above-water calving during the 26 d of acoustical observation is estimated to be 1.7±0.7×107 m3, whereas the subaerial calving flux estimated by traditional methods is 7±2×106 m3. The results suggest that passive cryoacoustics is a viable technique for long-term monitoring of mass loss from marine-terminating glaciers.

Annual mass budget of Antarctic ice shelves from 1997 to 2021.

Antarctic ice shelves moderate the contribution of the Antarctic Ice Sheet to global sea level rise; however, ice shelf health remains poorly constrained. Here, we present the annual mass budget of all Antarctic ice shelves from 1997 to 2021. Out of 162 ice shelves, 71 lost mass, 29 gained mass, and 62 did not change mass significantly. Of the shelves that lost mass, 68 had statistically significant negative mass trends, 48 lost more than 30% of their initial mass, and basal melting was the dominant contributor to that mass loss at a majority (68%). At many ice shelves, mass losses due to basal melting or iceberg calving were significantly positively correlated with grounding line discharge anomalies; however, the strength and form of this relationship varied substantially between ice shelves. Our results illustrate the utility of partitioning high-resolution ice shelf mass balance observations into its components to quantify the contributors to ice shelf mass change and the response of grounded ice.

Soundscape during summer in Kongsfjorden, Svalbard, Arctic

This paper summarises research on ambient noise data collected during the summers of 2016 and 2017 by an autonomous Ambient Noise Measurement System (ANMS) deployed in the Kongsfjorden, Arctic. The primary goals of this study are to identify the various soundscape components and the variation of the soundscape in the Arctic environment during the summers of 2016 and 2017. This work mainly quantifies the ambient noise sources during each month of the summer. After analysing the data, five distinct sources of soundscape were identified, out of which three geophysical sources, including iceberg bubbling, glacier calving, and rain noise, as well as two anthropogenic sources, such as shipping noise and ice ramming noise, were identified. No biological sources were identified either in the summer of 2016 or in 2017. A comparison of average Sound Pressure Levels (SPL) in the summers of 2016 and 2017 shows that the difference in average SPL below 8 kHz is 3 dB re 1 μPa2/Hz.

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

Simulating the processes controlling ice-shelf rift paths using damage mechanics

Abstract Rifts are full-thickness fractures that propagate laterally across an ice shelf. They cause ice-shelf weakening and calving of tabular icebergs, and control the initial size of calved icebergs. Here, we present a joint inverse and forward computational modeling framework to capture rifting by combining the vertically integrated momentum balance and anisotropic continuum damage mechanics formulations. We incorporate rift–flank boundary processes to investigate how the rift path is influenced by the pressure on rift–flank walls from seawater, contact between flanks, and ice mélange that may also transmit stress between flanks. To illustrate the viability of the framework, we simulate the final 2 years of rift propagation associated with the calving of tabular iceberg A68 in 2017. We find that the rift path can change with varying ice mélange conditions and the extent of contact between rift flanks. Combinations of parameters associated with slower rift widening rates yield simulated rift paths that best match observations. Our modeling framework lays the foundation for robust simulation of rifting and tabular calving processes, which can enable future studies on ice-sheet–climate interactions, and the effects of ice-shelf buttressing on land ice flow.

Impacts of tidewater glacier advance on iceberg habitat

Abstract Icebergs in proglacial fjords serve as pupping, resting and molting habitat for some of the largest seasonal aggregations of harbor seals (Phoca vitulina richardii) in Alaska. One of the largest aggregations in Southeast Alaska occurs in Johns Hopkins Inlet, Glacier Bay National Park, where up to 2000 seals use icebergs produced by Johns Hopkins Glacier. Like other advancing tidewater glaciers, the advance of Johns Hopkins Glacier over the past century has been facilitated by the growth and continual redistribution of a submarine end moraine, which has limited mass losses from iceberg calving and submarine melting and enabled glacier thickening by providing flow resistance. A 15-year record of aerial surveys reveals (i) a decline in iceberg concentrations concurrent with moraine growth and (ii) that the iceberg size distributions can be approximated as power law distributions, with relatively little variability and no clear trends in the power law exponent despite large changes in ice fluxes over seasonal and interannual timescales. Together, these observations suggest that sustained tidewater glacier advance should typically be associated with reductions in the number of large, habitable icebergs, which may have implications for harbor seals relying on iceberg habitat for critical life-history events.

Proglacial river stage derived from georectified time-lapse camera images, Inglefield Land, Northwest Greenland

The Greenland Ice Sheet is a leading source of global sea level rise, due to surface meltwater runoff and glacier calving. However, given a scarcity of proglacial river gauge measurements, ice sheet runoff remains poorly quantified. This lack of in situ observations is particularly acute in Northwest Greenland, a remote area releasing significant runoff and where traditional river gauging is exceptionally challenging. Here, we demonstrate that georectified time-lapse camera images accurately retrieve stage fluctuations of the proglacial Minturn River, Inglefield Land, over a 3 year study period. Camera images discern the river’s wetted shoreline position, and a terrestrial LiDAR scanner (TLS) scan of riverbank microtopography enables georectification of these positions to vertical estimates of river stage. This non-contact approach captures seasonal, diurnal, and episodic runoff draining a large (∼2,800 km2) lobe of grounded ice at Inglefield Land with good accuracy relative to traditional in situ bubble-gauge measurements (r2 = 0.81, Root Mean Square Error (RMSE) ±0.185 m for image collection at 3-h frequency; r2 = 0.92, RMSE ±0.109 m for resampled average daily frequency). Furthermore, camera images effectively supplement other instrument data gaps during icy and/or low flow conditions, which challenge bubble-gauges and other contact-based instruments. This benefit alone extends the effective seasonal hydrological monitoring period by ∼2–4 weeks each year for the Minturn River. We conclude that low-cost, non-contact time-lapse camera methods offer good promise for monitoring proglacial meltwater runoff from the Greenland Ice Sheet and other harsh polar environments.