Dry Snow Zone Research Articles

Validation of airborne and satellite altimetry data by Arctic Truck citizen science

The elevation of ice sheets changes due to climate change, and satellite altimetry is the preferred tool for measuring ice sheet-wide height changes. In situ validation is needed to ensure the quality of the observed elevation changes, but the cost often limits the amount of in situ data which can be collected. As more tourists are accessing the ice sheets, citizen science might provide in situ data in an environmentally friendly and cost-efficient way. Here, we investigate the opportunistic kinematic global positioning system (GPS) profiles across the Greenland ice sheet, collected by the American-Icelandic expedition on the Greenlandic icecap 2018. The collected GPS data are in good agreement with the widely used NASA’s Operation IceBridge Airborne LiDAR data measured within ± 10 days, with an average difference of 10.7 cm ± 11.7 cm. The main difference is attributed to changes in the compaction of the snow while driving and changes in the tires’ pressure. The kinematic GPS data are then compared with data from the European Space Agency’s CryoSat-2 mission. Here, an average bias of 92.3 cm ± 65.7 cm in the two records is observed between the spring CryoSat-2 and the truck GPS data obtained largely in the dry-snow zone. This suggests that the surface penetration of Ku-band radar on the Greenland ice sheet and the observed magnitude are consistent with the literature. Finally, we compared the 2018 GPS data to a profile obtained in 2005 near Kangerlussuaq, West Greenland. Here, the records show an average ice-elevation decrease of 9 m, with peaks at 26 m. These results show that the citizen science kinematic GPS data can provide high-resolution data necessary for the validation of satellite altimetry, with the added benefit of potential direct sampling properties of the surface and firn. Linking up with citizen-science expeditions is a beneficial way of providing cost-effective satellite validations and may also have a societal impact by involving more people in the climate monitoring of ice sheets.

Evaluation of Surface Melt on the Greenland Ice Sheet Using SMAP L-Band Microwave Radiometry

Monitoring melt extent and timing on the Greenland ice sheet is important for tracking the ice sheet's mass and energy balance as well as the global and Arctic climate variability and change. In this study, we use L -band (1.4 GHz) brightness temperature observations collected by NASA's soil moisture active passive (SMAP) mission to investigate the extent, duration, and intensity of melt events on the Greenland ice sheet from 2015 to 2021. SMAP provides nearly all-weather surface monitoring over all of Greenland twice daily with morning and evening overpasses at approximately 40-km spatial resolution. We applied empirical threshold and geophysical-model-based algorithms using horizontally and vertically polarized microwave brightness temperature differences to quantify both the intensity and extent of surface melting. Analysis of the melt seasons shows that Greenland experienced unusually strong melt events at the end of July 2019 and on August 14, 2021, which extended the melt area across much of the dry snow zone over a period of one and two days, respectively. In situ temperatures measured at Greenland's Summit station confirm the above freezing temperatures during these extreme events.

Facies detection of surge-type glaciers, Karakoram Himalaya using Sentinel-1 images

Glacier facies or radar glacier zones (RGZs) of the surge-type glaciers of the Karakoram Himalaya were classified based on the differences in the backscatter values utilizing multi-temporal SAR data from Sentinel-1 satellite of the years 2018 and 2019. Distinct zones of glaciers such as dry snow zones, percolation zones, wet snow zone, bare ice zone and debris-covered area can be discernible on the false-color composite image made from the early summer-late summer-winter (R-G-B) of SAR images obtained from Sentinel-1. The result suggests an interesting blended pattern of facies distribution on surge-type glaciers showing the impression of accumulation zones at the lower most parts of the glaciers contrary to typical non surge-type glaciers. We found that the facies of the surge-type glaciers of the Karakoram Himalaya did not follow the relationship with the elevation in terms of their distribution on the surface of glaciers. We also compared our results with the surface temperature of the glaciers derived from the thermal band of Landsat ETM+.

Prototype wireless sensors for monitoring subsurface processes in snow and firn

ABSTRACTThe detection and monitoring of meltwater within firn presents a significant monitoring challenge. We explore the potential of small wireless sensors (ETracer+, ET+) to measure temperature, pressure, electrical conductivity and thus the presence or absence of meltwater within firn, through tests in the dry snow zone at the East Greenland Ice Core Project site. The tested sensor platforms are small, robust and low cost, and communicate data via a VHF radio link to surface receivers. The sensors were deployed in low-temperature firn at the centre and shear margins of an ice stream for 4 weeks, and a ‘bucket experiment’ was used to test the detection of water within otherwise dry firn. The tests showed the ET+ could log subsurface temperatures and transmit the recorded data through up to 150 m dry firn. Two VHF receivers were tested: an autonomous phase-sensitive radio-echo sounder (ApRES) and a WinRadio. The ApRES can combine high-resolution imaging of the firn layers (by radio-echo sounding) with in situ measurements from the sensors, to build up a high spatial and temporal resolution picture of the subsurface. These results indicate that wireless sensors have great potential for long-term monitoring of firn processes.

Hyper-spectral data based investigations for snow wetness mapping

Spatial information on snow wetness content (SWC) is important for hydrology, climatology applications. Limited work is available on estimation of SWC using optical sensors. In present work, spectral signature characteristics of snow (~145 samples) acquired in winters of three years, using field spectral-radiometer (350–2500 nm) were correlated with synchronized SWC measurements. Correlation is found stronger in Near-Infra-Red (NIR) and Short-Wave-Infrared (SWIR) regions than Visible (VIS). Spectral peak width at 905 and 1240 nm is found negatively correlated with SWC, while positively correlated at 1025 nm. Asymmetry tends towards right as SWC increases and has stable positive correlations as compared to other characteristics. Sensitivity of widely used snow-related indices to SWC is also analyzed. Based on analysis, new ratio method at selected wavelengths is proposed to discriminate dry and wet snow zones using air/ground borne sensors. Proposed methodology is evaluated on air-borne hyper-spectral (AVIRIS-NG) data and 88% overall accuracy with kappa coefficient 77.6 observed after validation with reference observations.

500–2000-MHz Brightness Temperature Spectra of the Northwestern Greenland Ice Sheet

An ultra-wideband radiometer has been developed to measure subsurface properties of the cryosphere including ice sheets and sea ice. The radiometer measures brightness temperature spectra from 0.5 to 2 GHz using 12 channels, each of which measures scene brightness temperatures over an ~88-MHz bandwidth resolved into 0.24-MHz intervals. The instrument was flown over northwestern Greenland in September 2016 and acquired the first, wideband, low-frequency brightness temperature spectra over the ice sheet and coastal region. The results reveal strong spatial and spectral variations that correlate well with the physical properties of the surface encountered along the flight path, which started over ocean, then passed the rock near the coast, and then up onto the ablation, wet, percolation, and dry snow zones of the interior ice sheet. In particular, strong spectral responses in percolation and dry snow zones are observed and plausibly explained by varying the distribution of horizontal density layers and isolated icy bodies in the upper portion of the firn. The success of the airborne deployment of the instrument and subsequent implementation of algorithms to limit radio frequency interference in unprotected bands is motivating continued airborne investigations as well as stimulating research into the feasibility of a spaceborne instrument.

The darkening of the Greenland ice sheet: trends, drivers, and projections (1981–2100)

Abstract. The surface energy balance and meltwater production of the Greenland ice sheet (GrIS) are modulated by snow and ice albedo through the amount of absorbed solar radiation. Here we show, using space-borne multispectral data collected during the 3 decades from 1981 to 2012, that summertime surface albedo over the GrIS decreased at a statistically significant (99 %) rate of 0.02 decade−1 between 1996 and 2012. Over the same period, albedo modelled by the Modèle Atmosphérique Régionale (MAR) also shows a decrease, though at a lower rate ( ∼ −0.01 decade−1) than that obtained from space-borne data. We suggest that the discrepancy between modelled and measured albedo trends can be explained by the absence in the model of processes associated with the presence of light-absorbing impurities. The negative trend in observed albedo is confined to the regions of the GrIS that undergo melting in summer, with the dry-snow zone showing no trend. The period 1981–1996 also showed no statistically significant trend over the whole GrIS. Analysis of MAR outputs indicates that the observed albedo decrease is attributable to the combined effects of increased near-surface air temperatures, which enhanced melt and promoted growth in snow grain size and the expansion of bare ice areas, and to trends in light-absorbing impurities (LAI) on the snow and ice surfaces. Neither aerosol models nor in situ and remote sensing observations indicate increasing trends in LAI in the atmosphere over Greenland. Similarly, an analysis of the number of fires and BC emissions from fires points to the absence of trends for such quantities. This suggests that the apparent increase of LAI in snow and ice might be related to the exposure of a "dark band" of dirty ice and to increased consolidation of LAI at the surface with melt, not to increased aerosol deposition. Albedo projections through to the end of the century under different warming scenarios consistently point to continued darkening, with albedo anomalies averaged over the whole ice sheet lower by 0.08 in 2100 than in 2000, driven solely by a warming climate. Future darkening is likely underestimated because of known underestimates in modelled melting (as seen in hindcasts) and because the model albedo scheme does not currently include the effects of LAI, which have a positive feedback on albedo decline through increased melting, grain growth, and darkening.

Drivers of ASCAT C band backscatter variability in the dry snow zone of Antarctica

ABSTRACTC band backscatter parameters contain information about the upper snowpack/firn in the dry snow zone. The wide incidence angle diversity of the Advanced Scatterometer (ASCAT) gives unprecedented characterisation of backscatter anisotropy, revealing the backscatter response to climatic forcing. TheA(isotropic component) andM2(bi-sinusoidal azimuth anisotropy) parameters are investigated here, in conjunction with data from atmospheric and snowpack models, to identify the backscatter response to surface forcing parameters (wind speed and persistence, precipitation, surface temperature, density and grain size). The long-term meanAparameter is successfully recreated with a regression using these drivers, indicating strong links between theAparameter and precipitation on long timescales. While the ASCAT time series is too short to determine which factors drive observed trends, factors influencing the seasonal and short timescale variability are revealed. On these timescales,Astrongly responds to the propagation of surface temperature cycles/anomalies downward through the firn, via direct modulation of the dielectric constant. The influence of precipitation onAis small at shorter timescales. TheM2parameter is controlled by wind speed and persistence, through modification of monodirectionally-aligned surface roughness. This variability indicates that throughout much of coastal Antarctica, a microwave ‘snapshot’ is generally not representative of longer-term conditions.

Neither dust nor black carbon causing apparent albedo decline in Greenland's dry snow zone: Implications for MODIS C5 surface reflectance

AbstractRemote sensing observations suggest Greenland ice sheet (GrIS) albedo has declined since 2001, even in the dry snow zone. We seek to explain the apparent dry snow albedo decline. We analyze samples representing 2012–2014 snowfall across NW Greenland for black carbon and dust light‐absorbing impurities (LAI) and model their impacts on snow albedo. Albedo reductions due to LAI are small, averaging 0.003, with episodic enhancements resulting in reductions of 0.01–0.02. No significant increase in black carbon or dust concentrations relative to recent decades is found. Enhanced deposition of LAI is not, therefore, causing significant dry snow albedo reduction or driving melt events. Analysis of Collection 5 Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data indicates that the decline and spectral shift in dry snow albedo contains important contributions from uncorrected Terra sensor degradation. Though discrepancies are mostly below the stated accuracy of MODIS products, they will require revisiting some prior conclusions with C6 data.

Ice Height and Backscattering Coefficient Variability over Greenland Ice Sheets Using SARAL Radar Altimeter

Ice sheets investigation is important with regard to climate change and contribution to the sea level rise or fall. Radar altimetry in complement with laser altimetry can serve as a suitable candidate for precise monitoring of ice sheet evaluations. SARAL due to higher observation into the polar region (up to 82.5°N) can cover nearly 100% of the Greenland ice sheet. Continuous ice tracking mode retracker can provide useful information about ice surfaces, that is, determining the snow coverage, ice sheet transaction margin, and the evolution of snow depth during winter more accurately. This study present the results obtained with SARAL satellite Altika radar altimeter over the Greenland ice sheet region. The altimeter high rate waveforms products are used for utilizing the full capability of the instrument. High resolution DEM (1 km) generated using ICESAT/GLAS altimeter has been used for selecting the good quality data over the study region. Four different retrackers—Ocean, ICE-1, ICE-2, and Sea-Ice—were tested on the SARAL altimeter data set and compared with the DEM extracted ice sheet elevations. Three different data analysis—region of interest (ROI), track analysis, and cross-over analysis—were performed for in-depth analysis of the ice height changes and back scattering coefficient variability. ROI's (1° × 0.5°) were selected based on accumulation dry snow zone, percolation zone, wet snow zone, and ablation zone. Finally to observe the effect of Ka band, SARAL results has been compared with the Envisat altimeter in terms of back scatter and error in the height retrieval due to penetration problem within the ice sheet layer. The new SARAL data set confirms the potential of ice altimetry and provides a new opportunity to monitor the ice sheet surface topography evolution.

Assessing the role of sublimation in the dry snow zone of the Greenland ice sheet in a warming world

AbstractMeteorological and glaciological data obtained over an intensive 2 year measurement period (2000–2002) are used to run a physically based climatic mass balance model to characterize a seasonal variability in mass and energy exchanges at Summit, Greenland. The model resolves the full surface energy balance and the subsurface temperature profile by inclusion of energy release from penetrating shortwave radiation. A Monte Carlo approach using 1000 different parameter combinations is adopted to assess model uncertainty, with output compared to measured surface and subsurface temperatures, changes in surface height, and eddy correlation data. The heat exchanges associated with the change in phase of water are very small in all seasons, with the average turbulent latent heat flux equal to 0.4 (±0.2) W m−2. This suggests that the mean annual water vapor gradient is toward the surface, resulting in a mass gain of 4.1 mm WE yr−1. The mass gain represents only a small fraction of the total accumulation (<2%), in part because of the change in sign of the water vapor flux from winter (deposition) to summer (sublimation), but if assumed to be typical of the entire dry snow zone (40% of the total ice sheet area) is equivalent to approximately 5.5 Gt yr−1. A simple experiment based on 2012 atmospheric conditions suggests that mass turnover from water vapor exchanges will likely be enhanced in a warming climate, with sublimation increasing more than deposition. Should the sign of the mean turbulent latent heat flux change due to warming, the present mass gain in the dry snow zone could easily become a mass loss of equal proportion, which would further enhance the negative mass balance of the Greenland ice sheet.

Considerations for Ku-Band Scatterometer Calibration Using the Dry-Snow Zone of the Greenland Ice Sheet

Postlaunch calibration of satellite-borne scatterometers using backscatter data from natural land targets helps to maintain scatterometer accuracy. Due to its temporal stability, the dry-snow zone of the Greenland ice sheet has been proposed in previous studies as a calibration target. Using QuikSCAT data, this letter examines the backscatter properties of the dry-snow zone that are relevant to Ku-band scatterometer calibration including temporal and spatial variabilities, and azimuth-angle and polarization and incidence-angle dependences. The backscatter is found to seasonally vary throughout the dry-snow zone by 0.4 dB on average. Small interannual variations (less than 1.5 dB over a nine year period) in the backscatter are also present in some regions. Azimuth modulation is generally less than 0.4 dB in magnitude and is not significant in some regions of the dry-snow zone. Melting and refreezing appear to cause the quasi-polarization ratio to temporarily decrease. Spatially consistent and relatively temporally stable regions that are well within the interior of the dry-snow zone are best suited as calibration sites.

Future projections of the Greenland ice sheet energy balance driving the surface melt

Abstract. In this study, simulations at 25 km resolution are performed over the Greenland ice sheet (GrIS) throughout the 20th and 21st centuries, using the regional climate model MAR forced by four RCP scenarios from three CMIP5 global circulation models (GCMs), in order to investigate the projected changes of the surface energy balance (SEB) components driving the surface melt. Analysis of 2000–2100 melt anomalies compared to melt results over 1980–1999 reveals an exponential relationship of the GrIS surface melt rate simulated by MAR to the near-surface air temperature (TAS) anomalies, mainly due to the surface albedo positive feedback associated with the extension of bare ice areas in summer. On the GrIS margins, the future melt anomalies are preferentially driven by stronger sensible heat fluxes, induced by enhanced warm air advection over the ice sheet. Over the central dry snow zone, the surface albedo positive feedback induced by the increase in summer melt exceeds the negative feedback of heavier snowfall for TAS anomalies higher than 4 °C. In addition to the incoming longwave flux increase associated with the atmosphere warming, GCM-forced MAR simulations project an increase of the cloud cover decreasing the ratio of the incoming shortwave versus longwave radiation and dampening the albedo feedback. However, it should be noted that this trend in the cloud cover is contrary to that simulated by ERA-Interim–forced MAR for recent climate conditions, where the observed melt increase since the 1990s seems mainly to be a consequence of more anticyclonic atmospheric conditions. Finally, no significant change is projected in the length of the melt season, which highlights the importance of solar radiation absorbed by the ice sheet surface in the melt SEB.

Present status and variations in the Arctic energy balance

The total solar irradiance (TSI, or solar constant) acquired a new value: 1361 W m−2 instead of 1365 W m−2. However a long-term variation of TSI was not detected. The solar irradiance at the earth's surface is considerably smaller (170 W m−2) than previously believed (e.g. 198 W m−2 of IPCC AR4). The previous overestimation is due to the underestimation of the absorption of solar radiation in the atmosphere. The absorption of solar radiation in the atmosphere at about 90 W m−2, or 25–28% of the primary solar radiation from space. The global mean atmospheric downward terrestrial radiation is much larger (345 W m−2) than previously assumed (325 W m−2 of IPCC AR4). The Arctic has regions of negative annual net radiation, a very rare phenomenon on the globe. These regions are the Central Arctic Ocean with its multi-year ice coverage and the accumulation area of the Greenland ice sheet. The energy balance of these regions is presented. Long-wave incoming radiation has been increasing in the Arctic at a rate of 4–5 W m−2/Decade. The Greenland ice sheet exhibits a large vertical difference in net radiation from the ablation area to the dry snow zone in summer. It ranges from 80 W m−2 in the ablation area to 20 W m−2 at the equilibrium line and to 10 W m−2 in the dry snow zone. This gradient determines the melt gradient on the ice sheet, and is mainly caused by the altitude variation in atmospheric long-wave radiation, seconded by the albedo variation. The effect of albedo in summer for various surfaces is discussed. Simulation capabilities of radiation for many GCMs are investigated.

Melt regimes, stratigraphy, flow dynamics and glaciochemistry of three glaciers in the Alaska Range

AbstractWe used ground-penetrating radar (GPR), GPS and glaciochemistry to evaluate melt regimes and ice depths, important variables for mass-balance and ice-volume studies, of Upper Yentna Glacier, Upper Kahiltna Glacier and the Mount Hunter ice divide, Alaska. We show the wet, percolation and dry snow zones located below ~2700ma.s.l., at ~2700 to 3900ma.s.l. and above 3900ma.s.l., respectively. We successfully imaged glacier ice depths upwards of 480 m using 40-100 MHz GPR frequencies. This depth is nearly double previous depth measurements reached using mid-frequency GPR systems on temperate glaciers. Few Holocene-length climate records are available in Alaska, hence we also assess stratigraphy and flow dynamics at each study site as a potential ice-core location. Ice layers in shallow firn cores and attenuated glaciochemical signals or lacking strata in GPR profiles collected on Upper Yentna Glacier suggest that regions below 2800ma.s.l. are inappropriate for paleoclimate studies because of chemical diffusion, through melt. Flow complexities on Kahiltna Glacier preclude ice-core climate studies. Minimal signs of melt or deformation, and depth-age model estimates suggesting ~4815 years of ice on the Mount Hunter ice divide (3912ma.s.l.) make it a suitable Holocene-age ice-core location.

Spatially extensive estimates of annual accumulation in the dry snow zone of the Greenland Ice Sheet determined from radar altimetry

Abstract. We present estimates of accumulation rate along a 200 km transect ranging in elevation from 2750 to 3150 m in the dry snow zone on the western slope of the Greenland Ice Sheet. An airborne radar altimeter is used to estimate the thickness of annual internal layers and, in conjunction with ground based snow/firn density profiles, annual accumulation rates between 1998 and 2003 are derived. A clear gradient in the thickness of each layer observed by the radar altimeter and in the associated estimates of annual accumulation is seen along the transect, with a 33.6% ± 16% mean decrease in accumulation from west to east. The observed inter-annual variability is high, with the annual mean accumulation rate estimated at 0.359 m.w.e. yr−1 (s.d. ± 0.049 m.w.e. yr−1). Mean accumulation rates modelled using meteorological models overestimate our results by 16% on average, but by 32% and 42% in the years 2001 and 2002. The methodology presented here demonstrates the potential to obtain accurate and spatially extensive accumulation rates from radar altimeters in regions of ice sheets where field observations are sparse, and accumulation rates greater than several tens of cm.

ESA'S POLarimetric Airborne Radar Ice Sounder (POLARIS): design and first results

The Technical University of Denmark has developed and tested a P-band ice sounding radar for European Space Agency (ESA). With the recent by the International Telecommunication Union (ITU) allocation of a radar band at 435 MHz, increased interest in space-based sounding of the Earth's ice caps has been encountered. ESA's POLarimetric Airborne Radar Ice Sounder (POLARIS) is intended to provide a better understanding of P-band scattering and propagation through ice sheets and to verify novel surface clutter suppression techniques in preparation for a potential space-based ice sounding mission. POLARIS is a nadir-looking, fully polarimetric radar featuring aperture synthesis, a multi-aperture antenna for surface clutter suppression and a large dynamic range based on a shallow/deep sounding approach. The system is installed in a De Haviland DHC-6 Twin Otter aircraft, and in May 2008, a proof-of-concept campaign was conducted in Greenland. This study outlines the design and implementation of the system, and based on first results it is concluded that in the central dry snow zone of Greenland, POLARIS can resolve shallow and deep internal ice layers, penetrate the thickest ice encountered and detect the bedrock.

Geographic variability of nitrate deposition and preservation over the Greenland Ice Sheet

An analysis of 96 snow pit and ice cores distributed over the Greenland ice sheet is used to determine the main drivers of variability in the preserved records of nitrate concentration. The data set provides samples from spatially distributed locations, allowing us to investigate the effect of snow accumulation rate, temperature, and sublimation on nitrate concentration. The mean ice sheet concentration in the dry snow zone (2000 ≥ mean annual sea level (masl)) is 132 ng g−1, ranging between 47 and 265 ng g−1 with a standard deviation of ±37 ng g−1. Nitrate flux varies between 1.1 and 14.7 μg cm−2 a−1 with a mean of 4 ± 2 μg cm−2 a−1. Large‐scale spatial variability exists as a result of accumulation gradients, with concentration 5% greater in the northern plateau, yet flux over the northern plateau is 30% lower than the dry snow zone as a whole. While spatially, flux appears to be more dependent on accumulation, preservation of flux shows increasing dependence on concentration with increasing accumulation. The relationship between concentration and accumulation is nonlinear, showing less dependence in the low‐accumulation regions versus high‐accumulation regions. Accumulation alone is insufficient to account for the observed variability in nitrate flux in the low‐accumulation regions, and evidence supports the need for additional components to a transfer function model for nitrate that includes photochemistry, temperature, and sublimation. Spatial variability across the ice sheet is nonuniform, and changes in nitrate concentration have occurred in some regions at a greater rate than others. While the data support that overall the ice sheet acts as an archive of paleoatmospheric concentration despite the effects of postdepositional processing, one needs to consider spatial variables to properly account for trends and variability in the records. This is tested by evaluating past spatial relationships and yields the result that the significant geographic shifts with respect to reactive N concentrations have occurred over the ice sheet in the past century.

WindSat Passive Microwave Polarimetric Signatures of the Greenland Ice Sheet

WindSat has systematically collected the first global fully polarimetric passive microwave data over both land and ocean. As the first spaceborne polarimetric microwave radiometer, it was designed to measure ocean surface wind speed and direction by including the third and fourth Stokes parameters, which are mostly related to the asymmetric structures of the ocean surface roughness. Although designed for wind vector retrieval, WindSat data are also collected over land and ice, and this new data has revealed, for the first time, significant land signals in the third and fourth Stokes parameter channels, particularly over Greenland and the Antarctic ice sheets. The third and fourth Stokes parameters show well-defined large azimuth modulations that appear to be correlated with geophysical variations, particularly snow structure, melting, and metamorphism, and have distinct seasonal variation. The polarimetric signatures are relatively weak in the summer and are strongest around spring. This corresponds well with the formation and erosion of the sastrugi in the dry snow zone and snowmelt in the soaked zone. In this paper, we present the full polarimetric signatures obtained from WindSat over Greenland, and use a simple empirical observation model to quantify the azimuthal variations of the signatures in space and time.

Empirical Retrieval of Surface Melt Magnitude from Coupled MODIS Optical and Thermal Measurements over the Greenland Ice Sheet during the 2001 Ablation Season.

Accelerated ice flow near the equilibrium line of west-central Greenland Ice Sheet (GIS) has been attributed to an increase in infiltrated surface melt water as a response to climate warming. The assessment of surface melting events must be more than the detection of melt onset or extent. Retrieval of surface melt magnitude is necessary to improve understanding of ice sheet flow and surface melt coupling. In this paper, we report on a new technique to quantify the magnitude of surface melt. Cloud-free dates of June 10, July 5, 7, 9, and 11, 2001 Moderate Resolution Imaging Spectroradiometer (MODIS) daily reflectance Band 5 (1.230-1.250μm) and surface temperature images rescaled to 1km over western Greenland were used in the retrieval algorithm. An optical-thermal feature space partitioned as a function of melt magnitude was derived using a one-dimensional thermal snowmelt model (SNTHERM89). SNTHERM89 was forced by hourly meteorological data from the Greenland Climate Network (GC-Net) at reference sites spanning dry snow, percolation, and wet snow zones in the Jakobshavn drainage basin in western GIS. Melt magnitude or effective melt (E-melt) was derived for satellite composite periods covering May, June, and July displaying low fractions (0-1%) at elevations greater than 2500m and fractions at or greater than 15% at elevations lower than 1000m assessed for only the upper 5 cm of the snow surface. Validation of E-melt involved comparison of intensity to dry and wet zones determined from QSCAT backscatter. Higher intensities (> 8%) were distributed in wet snow zones, while lower intensities were grouped in dry zones at a first order accuracy of ∼ ±2%.