A computationally efficient depression-filling algorithm for digital elevation models, applied to proglacial lake drainage
-
Published:2016-12-15
Issue:12
Volume:9
Page:4451-4460
-
ISSN:1991-9603
-
Container-title:Geoscientific Model Development
-
language:en
-
Short-container-title:Geosci. Model Dev.
Author:
Berends Constantijn J.,van de Wal Roderik S. W.
Abstract
Abstract. Many processes govern the deglaciation of ice sheets. One of the processes that is usually ignored is the calving of ice in lakes that temporarily surround the ice sheet. In order to capture this process a “flood-fill algorithm” is needed. Here we present and evaluate several optimizations to a standard flood-fill algorithm in terms of computational efficiency. As an example, we determine the land–ocean mask for a 1 km resolution digital elevation model (DEM) of North America and Greenland, a geographical area of roughly 7000 by 5000 km (roughly 35 million elements), about half of which is covered by ocean. Determining the land–ocean mask with our improved flood-fill algorithm reduces computation time by 90 % relative to using a standard stack-based flood-fill algorithm. This implies that it is now feasible to include the calving of ice in lakes as a dynamical process inside an ice-sheet model. We demonstrate this by using bedrock elevation, ice thickness and geoid perturbation fields from the output of a coupled ice-sheet–sea-level equation model at 30 000 years before present and determine the extent of Lake Agassiz, using both the standard and improved versions of the flood-fill algorithm. We show that several optimizations to the flood-fill algorithm used for filling a depression up to a water level, which is not defined beforehand, decrease the computation time by up to 99 %. The resulting reduction in computation time allows determination of the extent and volume of depressions in a DEM over large geographical grids or repeatedly over long periods of time, where computation time might otherwise be a limiting factor. The algorithm can be used for all glaciological and hydrological models, which need to trace the evolution over time of lakes or drainage basins in general.
Funder
Ministerie van Onderwijs, Cultuur en Wetenschap
Publisher
Copernicus GmbH
Reference23 articles.
1. Amante, C. and Eakins, B. W.: ETOPO1 Arc-Minute Global Relief Model: Procedures, Data Sources and Analysis, National Oceanic and Atmospheric Administration Technical Memorandum NESDIS NGDC-24, https://doi.org/10.7289/V5C8276M, 2009. 2. Arnold, N.: A new approach for dealing with depressions in digital elevation models when calculating flow accumulation values, Prog. Phys. Geog., 34, 781–809, https://doi.org/10.1177/0309133310384542, 2010. 3. Barber, D. C., Dyke, A., Hillaire-Marcel, C., Jennings, A. E., Andrews, J. T., Kerwin, M. W., Bilodeau, G., McNeely, R., Southon, J., Morehead, M. D., and Gagnon, J.-M.: Forcing of the cold event of 8,200 years ago by catastrophic drainage of Laurentide lakes, Nature, 34, 781–809, https://doi.org/10.1038/22504, 1999. 4. Broecker, W. S., Kenett, J. P., Flower, B. P., Teller, J., Trumbore, S., Bonani, G., and Wolfli, W.: Routing of meltwater from the Laurentide Ice Sheet during the Younger Dryas cold episode, Nature, 341, 318–321, https://doi.org/10.1038/341318a0, 1989. 5. Clark, P. U., Marshall, S. J., Clarke, G., Hostetler, S. W., Licciardi, J. M., and Teller, J.: Freshwater Forcing of Abrupt Climate Change During the Last Glaciation, Science, 293, 283–287, https://doi.org/10.1038/341318a0, 2001.
Cited by
9 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献
|
|