Surface mass balance model evaluation from satellite and airborne lidar mapping

We present estimates of Greenland Ice Sheet (GrIS) surface elevation change from a novel combination of satellite and airborne laser altimetry measurements. Our method combines measurements from the Airborne Topographic Mapper (ATM), the Land, Vegetation and Ice Sensor (LVIS) and ICESat-1 to generate elevation change rates at high spatial resolution. This method allows to extend the records of each instrument, increases the overall spatial coverage compared to a single instrument, and produces high-quality, coherent maps of surface elevation change. In addition by combining the lidar datasets, we are able to investigate seasonal and interannual surface elevation change for years where Spring and Fall Operation IceBridge campaigns are available. We validate our method by comparing with the standard NSIDC elevation change product calculated using overlapping Level-1B ATM data. We use the altimetry-derived mass changes to evaluate the uncertainty in surface mass balance, particularly in the runoff component, from two Regional Climate Models (RCM's), the Regional Atmospheric Climate Model (RACMO) and the Modéle Atmosphérique Régional (MAR), and one Global Climate Model (GCM), MERRA2/GEOS-5. We investigate locations with low ice sheet surface velocities that are within the estimated ablation zones of each regional climate model. We find that the surface mass balance outputs from RACMO and MAR show good correspondence with mass changes derived from surface elevation changes over long periods. At two sites in Northeast Greenland (NEGIS), the MAR model has better correspondence with the altimetry estimate. We find that the differences at these locations are primarily due to the characterization of meltwater refreeze within the ice sheet.