A spatially resolved estimate of High Mountain Asia glacier mass balances from 2000 to 2016

High Mountain Asia hosts the largest glacier concentration outside the polar regions. These glaciers are important contributors to streamflow in one of the most populated areas of the world. Past studies have used methods that can provide only regionally averaged glacier mass balances to assess the glacier contribution to rivers and sea level rise. Here we compute the mass balance for about 92% of the glacierized area of High Mountain Asia using time series of digital elevation models derived from satellite stereo-imagery. We calculate a total mass change of -16.3 +/- 3.5 Gt yr(-1) (-0.18 +/- 0.04 m w.e. yr(-1)) between 2000 and 2016, which is less negative than most previous estimates. Region-wide mass balances vary from 4.0 +/- 1.5 Gt yr(-1) (-0.62 +/- 0.23 m w.e. yr(-1)) in Nyainqentanglha to +1.4 +/- 0.8 Gt yr(-1) (+0.14 +/- 0.08 m w.e. yr(-1)) in Kunlun, with large intra-regional variability of individual glacier mass balances (standard deviation within a region similar to 0.20m w.e. yr(-1)). Specifically, our results shed light on the Nyainqentanglha and Pamir glacier mass changes, for which contradictory estimates exist in the literature. They provide crucial information for the calibration of the models used for projecting glacier response to climatic change, as these models do not capture the pattern, magnitude and intra-regional variability of glacier changes at present.