Mass loss from the Greenland Ice Sheet (GrIS) has accelerated over the past two decades, coincident with rapid Arctic warming and increasing moisture transport over Greenland by atmospheric rivers (ARs). Summer ARs affecting western Greenland trigger GrIS melt events, but the physical mechanisms through which ARs induce melt are not well understood. This study elucidates the coupled surface–atmosphere processes by which ARs force GrIS melt through analysis of the surface energy balance (SEB), cloud properties, and local- to synoptic-scale atmospheric conditions during strong summer AR events affecting western Greenland. ARs are identified in MERRA-2 reanalysis (1980–2017) and classified by integrated water vapor transport (IVT) intensity. SEB, cloud, and atmospheric data from regional climate model, observational, reanalysis, and satellite-based datasets are used to analyze melt-inducing physical processes during strong, >90th percentile “AR₉₀₊” events. Near AR “landfall,” AR₉₀₊days feature increased cloud cover that reduces net shortwave radiation and increases net longwave radiation. As these oppositely signed radiative anomalies partly cancel during AR₉₀₊events, increased melt energy in the ablation zone is primarily provided by turbulent heat fluxes, particularly sensible heat flux. These turbulent heat fluxes are driven by enhanced barrier winds generated by a stronger synoptic pressure gradient combined with an enhanced local temperature contrast between cool over-ice air and the anomalously warm surrounding atmosphere. During AR₉₀₊events in northwest Greenland, anomalous melt is forced remotely through a clear-sky foehn regime produced by downslope flow in eastern Greenland.

We use the CloudSat 2006–2016 data record to estimate snowfall over the Greenland Ice Sheet (GrIS). We first evaluate CloudSat snowfall retrievals with respect to remaining ground-clutter issues. Comparing CloudSat observations to the GrIS topography (obtained from airborne altimetry measurements during IceBridge) we find that at the edges of the GrIS spurious high-snowfall retrievals caused by ground clutter occasionally affect the operational snowfall product. After correcting for this effect, the height of the lowest valid CloudSat observation is about 1200 m above the local topography as defined by IceBridge. We then use ground-based millimeter wavelength cloud radar (MMCR) observations obtained from the Integrated Characterization of Energy, Clouds, Atmospheric state, and Precipitation at Summit, Greenland (ICECAPS) experiment to devise a simple, empirical correction to account for precipitation processes occurring between the height of the observed CloudSat reflectivities and the snowfall near the surface. Using the height-corrected, clutter-cleared CloudSat reflectivities we next evaluate various Z–S relationships in terms of snowfall accumulation at Summit through comparison with weekly stake field observations of snow accumulation available since 2007. Using a set of three Z–S relationships that best agree with the observed accumulation at Summit, we then calculate the annual cycle snowfall over the entire GrIS as well as over different drainage areas and compare the derived mean values and annual cycles of snowfall to ERA-Interim reanalysis. We find the annual mean snowfall over the GrIS inferred from CloudSat to be 34±7.5 cm yr⁻¹ liquid equivalent (where the uncertainty is determined by the range in values between the three different Z–S relationships used). In comparison, the ERA-Interim reanalysis product only yields 30 cm yr⁻¹ liquid equivalent snowfall, where the majority of the underestimation in the reanalysis appears to occur in the summer months over the higher GrIS and appears to be related to shallow precipitation events. Comparing all available estimates of snowfall accumulation at Summit Station, we find the annually averaged liquid equivalent snowfall from the stake field to be between 20 and 24 cm yr⁻¹, depending on the assumed snowpack density and from CloudSat 23±4.5 cm yr⁻¹. The annual cycle at Summit is generally similar between all data sources, with the exception of ERA-Interim reanalysis, which shows the aforementioned underestimation during summer months.

Remote-sensing observations are needed to estimate the regional and global impacts of snow. However, to retrieve accurate estimates of snow mass and rate, these observations require augmentation through additional information and assumptions about hydrometeor properties. The Precipitation Imaging Package (PIP) provides information about precipitation characteristics and can be utilized to improve estimates of snowfall rate and accumulation. Here, the goal is to demonstrate the quality and utility of two higher-order PIP-derived products: liquid water equivalent snow rate and an approximation of volume-weighted density called equivalent density. Accuracy of the PIP snow rate and equivalent density is obtained through intercomparison with established retrieval methods and through evaluation with colocated ground-based observations. The results confirm the ability of the PIP-derived products to quantify properties of snow rate and equivalent density, and demonstrate that the PIP produces physically realistic snow characteristics. When compared to the National Weather Service (NWS) snow field measurements of six-hourly accumulation, the PIP-derived accumulations were biased only +2.48% higher. Additionally, this work illustrates fundamentally different microphysical and bulk features of low and high snow-to-liquid ratio events, through assessment of observed particle size distributions, retrieved mass coefficients, and bulk properties. Importantly, this research establishes the role that PIP observations and higher-order products can serve for constraining microphysical assumptions in ground-based and spaceborne remotely sensed snowfall retrievals.