Greenland Ice Sheet (GrIS) albedo has decreased over recent decades, contributing to enhanced surface melt and mass loss. However, it remains unclear whether GrIS darkening is due to snow grain size increases, higher concentrations of light‐absorbing impurities (LAIs), or a combination. Here, we assess albedo controls in the western GrIS percolation zone using in situ albedo, LAI, and grain size measurements. We find a significant correlation between albedo and snow grain size ( p < 0.01), but not with LAIs. Modeling corroborates that LAI concentrations are too low to significantly reduce albedo, but larger grain sizes could reduce albedo by at least ∼3%. Strong atmospheric blocking increases grain sizes and reduces albedo through increased surface temperature, fewer storms, and higher incoming shortwave radiation. These findings clarify the mechanisms by which anomalously strong blocking contributed to recent GrIS albedo decline and mass loss, highlighting the importance of improving projections of future blocking.

Abstract We present continuous estimates of snow and firn density, layer depth and accumulation from a multi-channel, multi-offset, ground-penetrating radar traverse. Our method uses the electromagnetic velocity, estimated from waveform travel-times measured at common-midpoints between sources and receivers. Previously, common-midpoint radar experiments on ice sheets have been limited to point observations. We completed radar velocity analysis in the upper ~2 m to estimate the surface and average snow density of the Greenland Ice Sheet. We parameterized the Herron and Langway (1980) firn density and age model using the radar-derived snow density, radar-derived surface mass balance (2015–2017) and reanalysis-derived temperature data. We applied structure-oriented filtering to the radar image along constant age horizons and increased the depth at which horizons could be reliably interpreted. We reconstructed the historical instantaneous surface mass balance, which we averaged into annual and multidecadal products along a 78 km traverse for the period 1984–2017. We found good agreement between our physically constrained parameterization and a firn core collected from the dry snow accumulation zone, and gained insights into the spatial correlation of surface snow density.

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow, and Ku-band radar and optical imagery acquired in 2014, 2016, and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88∘ S, an elevation bias validation site for ICESat-2 and potential validation site for CryoSat-2. We find significant small-scale variability (<10 km) in snow accumulation based on the snow radar subsurface stratigraphy, indicating areas of strong wind redistribution are prevalent at 88∘ S. In general, highs in snow accumulation rate correspond with topographic lows, resulting in a negative correlation coefficient of r2=-0.32 between accumulation rate and MSWD (mean slope in the mean wind direction). This relationship is strongest in areas where the dominant wind direction is parallel to the survey profile, which is expected as the geophysical surveys only capture a two-dimensional cross section of snow redistribution. Variability in snow accumulation appears to correlate with variability in MSWD. The correlation coefficient between the standard deviations of accumulation rate and MSWD is r2=0.48, indicating a stronger link between the standard deviations than the actual parameters. Our analysis shows that there is no simple relationship between surface slope, wind direction, and snow accumulation rates for the overall survey area. We find high variability in surface roughness derived from laser altimetry measurements on length scales smaller than 10 km, sometimes with very distinct and sharp transitions. Some areas also show significant temporal variability over the course of the 3 survey years. Ultimately, there is no statistically significant slope-independent relationship between surface roughness and accumulation rates within our survey area. The observed correspondence between the small-scale temporal and spatial variability in surface roughness and backscatter, as evidenced by Ku-band radar signal strength retrievals, will make it difficult to develop elevation bias corrections for radar altimeter retrieval algorithms.

The mass balance of the Greenland Ice Sheet (GrIS) in a warming climate is of critical interest in the context of future sea level rise. Increased melting in the GrIS percolation zone due to atmospheric warming over the past several decades has led to increased mass loss at lower elevations. Previous studies have hypothesized that this warming is accompanied by a precipitation increase, as would be expected from the Clausius–Clapeyron relationship, compensating for some of the melt-induced mass loss throughout the western GrIS. This study tests that hypothesis by calculating snow accumulation rates and trends across the western GrIS percolation zone, providing new accumulation rate estimates in regions with sparse in situ data or data that do not span the recent accelerating surface melt. We present accumulation records from sixteen 22–32 m long firn cores and 4436 km of ground-penetrating radar, covering the past 20–60 years of accumulation, collected across the western GrIS percolation zone as part of the Greenland Traverse for Accumulation and Climate Studies (GreenTrACS) project. Trends from both radar and firn cores, as well as commonly used regional climate models, show decreasing accumulation rates of 2.4±1.5 % a−1 over the 1996–2016 period, which we attribute to shifting storm tracks related to stronger atmospheric summer blocking over Greenland. Changes in atmospheric circulation over the past 20 years, specifically anomalously strong summertime blocking, have reduced GrIS surface mass balance through both an increase in surface melting and a decrease in accumulation rates.

We report annual snow accumulation rates from 1959 to 2004 along a 250 km segment of the Expéditions Glaciologiques Internationales au Groenland (EGIG) line across central Greenland using Airborne SAR/Interferometric Radar Altimeter System (ASIRAS) radar layers and high resolution neutron-probe (NP) density profiles. ASIRAS-NP-derived accumulation rates are not statistically different (95 % confidence interval) from in situ EGIG accumulation measurements from 1985 to 2004. ASIRAS-NP-derived accumulation increases by 20 % below 3000 m elevation, and increases by 13 % above 3000 m elevation for the period 1995 to 2004 compared to 1985 to 1994. Three Regional Climate Models (PolarMM5, RACMO2.3, MAR) underestimate snow accumulation below 3000 m by 16–20 % compared to ASIRAS-NP from 1985 to 2004. We test radar-derived accumulation rates sensitivity to density using modeled density profiles in place of NP densities. ASIRAS radar layers combined with Herron and Langway (1980) model density profiles (ASIRAS-HL) produce accumulation rates within 3.5 % of ASIRAS-NP estimates in the dry snow region. We suggest using Herron and Langway (1980) density profiles to calibrate radar layers detected in dry snow regions of ice sheets lacking detailed in situ density measurements, such as those observed by the Operation IceBridge campaign.

Abstract Accumulation is a key parameter governing the mass balance of the Greenland ice sheet. Several studies have documented the spatial variability of accumulation over wide spatial scales, primarily using point data, remote sensing or modeling. Direct measurements of spatially extensive, detailed profiles of accumulation in Greenland, however, are rare. We used 400 MHz ground-penetrating radar along the 1009 km route of the Greenland Inland Traverse from Thule to Summit during April and May of 2011, to image continuous internal reflecting horizons. We dated these horizons using ice-core chemistry at each end of the traverse. Using density profiles measured along the traverse, we determined the depth to the horizons and the corresponding water-equivalent accumulation rates. The measured accumulation rates vary from ~0.1 m w.e. a–1 in the interior to ~0.7 m w.e. a–1 near the coast, and correspond broadly with existing published model results, though there are some excursions. Comparison of our recent accumulation rates with those collected along a similar route in the 1950s shows a ~10% increase in accumulation rates over the past 52 years along most of the traverse route. This implies that the increased water vapor capacity of warmer air is increasing accumulation in the interior of Greenland.