The Greenland ice sheet (GrIS) is at present the largest single contributor to global-mass-induced sea-level rise, primarily because of Arctic amplification on an increasingly warmer Earth^1–5. However, the processes of englacial water accumulation, storage and ultimate release remain poorly constrained. Here we show that a noticeable amount of the summertime meltwater mass is temporally buffered along the entire GrIS periphery, peaking in July and gradually reducing thereafter. Our results arise from quantifying the spatiotemporal behaviour of the total mass of water leaving the GrIS by analysing bedrock elastic deformation measured by Global Navigation Satellite System (GNSS) stations. The buffered meltwater causes a subsidence of the bedrock close to GNSS stations of at most approximately 5 mm during the melt season. Regionally, the duration of meltwater storage ranges from 4.5 weeks in the southeast to 9 weeks elsewhere. We also show that the meltwater runoff modelled from regional climate models may contain systematic errors, requiring further scaling of up to about 20% for the warmest years. These results reveal a high potential for GNSS data to constrain poorly known hydrological processes in Greenland, forming the basis for improved projections of future GrIS melt behaviour and the associated sea-level rise⁶.

Mascon products derived from Gravity Recovery and Climate Experiment satellite gravimetry data are widely used to study the Greenland ice sheet mass balance. However, the products released by different research groups—JPL, CSR, and GSFC—show noticeable discrepancies. To understand them, we compare those mascon products with mascon solutions computed in-house using a varying regularization parameter. We show that the observed discrepancies are likely dominated by differences in the applied regularization. Furthermore, we present a numerical study aimed at an in-depth analysis of regularization-driven biases in the solutions. We demonstrate the ability of our simulations to reproduce 60%–80% of biases observed in real data, which proves that our simulations are sufficiently realistic. After that, we demonstrate that the quality of mascon-based estimates can be increased by a proper modification of the applied regularization: no correlation between mascons is assumed when they belong to different drainage systems. Using both simulations and real data analysis, we show that the improved regularization mitigates signal leakage between drainage systems by 11%–56%. Finally, we validate various mascon solutions over the SW drainage system, using trends from (i) the GOCO-06S model and (ii) the Input-Output Method as control data. In general, the in-house computed trend estimates are consistent with the trends from CSR and JPL solutions and the trends from the control data.