Data Set ID:

High Mountain Asia Near-Global Multi-Decadal Glacial Lake Inventory, Version 1

This data set contains polygons of glacial lake extent on a near-global scale, averaged over five multi-year periods between 1990 and 2018.

This is the most recent version of these data.

Version Summary:

Initial release

STANDARD Level of Service

Data: Data integrity and usability verified

Documentation: Key metadata and user guide available

User Support: Assistance with data access and usage; guidance on use of data in tools

See All Level of Service Details

Data Format(s):
  • ESRI Shapefile
Spatial Coverage:
N: 82.6, 
S: -55.4, 
E: 180, 
W: -180
Spatial Resolution:
  • 30 m x 30 m
Sensor(s):ETM+, MSS, OLI, TIRS, TM
Temporal Coverage:
  • 1 January 1990 to 31 December 2018
Temporal ResolutionVariesMetadata XML:View Metadata Record
Data Contributor(s):Dan Shugar

Geographic Coverage

Other Access Options

Other Access Options


As a condition of using these data, you must cite the use of this data set using the following citation. For more information, see our Use and Copyright Web page.

Shugar, D., A. Burr, U. K. Haritashya, J. Kargel, C. S. Watson, M. C. Kennedy, A. R. Bevington, R. A. Betts, S. Harrison, and K. Strattman. 2020. High Mountain Asia Near-Global Multi-Decadal Glacial Lake Inventory, Version 1. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: [Date Accessed].

Literature Citation

As a condition of using these data, we request that you acknowledge the author(s) of this data set by referencing the following peer-reviewed publication.

  • Shugar, D., A. Burr, U. K. Haritashya, J. Kargel, C. S. Watson, M. C. Kennedy, A. R. Bevington, R. A. Betts, S. Harrison, and K. Strattman. 2020. Rapid worldwide growth of glacial lakes since 1990, Nature Climate Change. . .

16 January 2019
Last modified: 
3 September 2020

Data Description


The main parameter in this data set is glacial lake extent in m2 (Area_m2). The data files also contain a set of additional attributes, such as Lat, Lon, Year_Start, and Country. Figure 1 is a screenshot of the QGIS user interface displaying the contents of the file HMA_GLI_glacial_lakes_LatLonWGS84_1990_99.shp.

Figure 1. Left: locations of various glacial lake polygons. Right: parameters for the red-highlighted polygon.

File Information


The data files are provided as ESRI Shapefiles (.shp, .shx, .dbf, .prj, .sbx, .sbn, .cpg).

Naming Convention

Example file name:


The files are named according to the following convention, which is described in Table 1:


Table 1. File Naming Convention
Variable Description
HMA_GLI Data set ID
glacial_lakes Indicates the focus on glacial lakes
LatLonWGS84 Indicates the use of the WGS 84 coordinate reference system
YYYY_yy Years included in the data file. E.g.: 1990_99 means that the data include the years 1990 to 1999.
.ext Indicates ESRI shapefiles: .shp, .shx, .dbf, .prj, .sbx, .sbn, .cpg

Spatial Information


The spatial coverage is nearly global, but does not include Antarctica and various smaller, high-latitude islands.

Northernmost latitude: 82.6° N
Southernmost latitude 55.4° S
Easternmost longitude: 180.0° E
Westernmost longitude: 180.0° E


Individual polygons with varying spatial extents and distances from each other.


Table 2 provides geolocation information for this data set.

Table 2. Geolocation Details
Geographic coordinate system WGS 84
EPSG code 4326
PROJ4 string +proj=longlat +datum=WGS84 +no_defs

Temporal Information


01 January 1990 to 31 December 2018


The polygons represent glacial lake extents averaged over five multi-year periods: 1990-1999, 2000-2004, 2005-2009, 2010-2014, and 2015-2018.

Data Acquisition and Processing


The size of glacial lakes has been rapidly changing due to climate change and the associated retreat of glaciers. The implications of continued glacial melt for glacial lake outburst floods and local and regional water resources are important both from societal and ecological standpoints. This data set can be used to explore the changes in glacial lake extents on a near-global scale by contrasting five different time periods between 1990 and 2018. For more background information, see Shugar et al. (2020).

Acquisition and Processing

Using optical imagery from Landsat missions 4, 5, 7, and 8 available in Google Earth Engine (254,795 individual scenes), an initial mapping of water pixels on a near-global scale was performed. To this end, the calibrated surface reflectance data from the 254,795 scenes were stacked to produce a multi-sensor data cube. The cube was restricted to only include relatively (typically <20%) cloud-free scenes from the ablation season, in order to minimize the likelihood of capturing snow and frozen lake surfaces. Using a Normalized Difference Water Index-based model implemented in Google Earth Engine, surface water was identified and outlined and subsequently filtered by a set of variables to retain only glacial lakes (i.e., supraglacial, proglacial, and ice-marginal lakes). The images were then aggregated by epoch and verified for complete coverage of glacier proximal areas in order to avoid biases related to differing spatiotemporal image densities.

For more information on the data acquisition and processing steps, see the Methods section in Shugar et al. (2020).

Quality, Errors, and Limitations

The largest error sources arise from varying water and adjacent terrain properties, such as the color of the glacial lakes, frozen lake conditions, and shadows cast by terrain or clouds. Errors due to omissions or inaccuracies in the source data also affect the quality of this data set. For more information on the error assessments performed on this data set, see Shugar et al. (2020).

Software and Tools

The data files can be opened using software that recognizes the ESRI Shapefile file format, such as QGIS, command-line tools from the GDAL project, and ArcMap.

Related Data Sets

High Mountain Asia at NSIDC | Data Sets

Related Websites

High Mountain Asia at NSIDC | Overview
NASA High Mountain Asia Project
NASA Research Announcement: Understanding Changes in High Mountain Asia
GLIMS: Global Land Ice Measurements from Space


Dr. Dan H. Shugar
University of Calgary


Support for this work was provided by NASA (NNX16AQ62G and 80NSSC19K0653) to J. S. Karge, U. K. Haritashya, and D. H. Shugar, and from NSERC (Discovery Grant 2020-04207 and Discovery Accelerator Supplement 2020-00066) to D. H. Shugar. Without free access of the Landsat data archive, this and many other scientific efforts would not have been possible.


Shugar, D. H., Burr, A., Haritashya, U. K., Kargel, J. S., Watson, C. S., Kennedy, M. C., Bevington, A. R., Betts, R. A., Harrison, S., & Strattman, K. (2020). Rapid worldwide growth of glacial lakes since 1990. Nature Climate Change.
No technical references available for this data set.

How To

Programmatic Data Access Guide
Data from the NASA National Snow and Ice Data Center Distributed Active Archive Center (NSIDC DAAC) can be accessed directly from our HTTPS file system or through our Application Programming Interface (API). Our API offers you the ability to order data using specific temporal and spatial filters... read more
Filter and order from a data set web page
Many NSIDC data set web pages provide the ability to search and filter data with spatial and temporal contstraints using a map-based interface. This article outlines how to order NSIDC DAAC data using advanced searching and filtering.  Step 1: Go to a data set web page This article will use the... read more