First spectral measurements of light attenuation in Greenland Ice Sheet bare ice suggest shallower subsurface radiative heating and ICESat-2 penetration depth in the ablation zone

Abstract. Light transmission into bare glacial ice affects surface energy balance, bio-photochemical cycling, and light detection and ranging (LiDAR) laser elevation measurements but has not previously been reported for the Greenland Ice Sheet. We present in-ice solar irradiance measured over the spectral range 350–900 nm and 12–77 cm depth collected at a site in the western Greenland ablation zone. The acquired spectral irradiance measurements are used to calculate flux attenuation coefficients using an exponential decay Bouguer law model and are compared to values calculated from two-stream radiative transfer theory. Relative to asymptotic two-stream theory, our empirical attenuation coefficients are up to one order of magnitude larger in the range 350–530 nm, suggesting light absorbing particles embedded in ice enhance visible light absorption at our field site. The empirical coefficients accurately describe light attenuation in the ice interior but underestimate light attenuation near the ice surface. Consequently, Bouguer’s law overestimates transmitted flux by up to 50 % depending on wavelength. Refraction is unlikely to explain the discrepancy. Instead, vertical variation in the ice microstructure and the concentration of light absorbing particles appears to enhance near-surface attenuation at our field site. The magnitude of this near-surface attenuation implies that optical penetration depth is lower by up to 19 cm (28 %) at wavelengths relevant to visible-wavelength lidar altimetry of ice surface elevation (e.g. 532 nm for the Ice, Cloud, and Land Elevation Satellite-2) than is suggested by e-folding depths inferred from two stream theory for optically pure glacier ice. This enhanced near-surface attenuation implies shallower light transmission and therefore lower subsurface light availability for subsurface radiative heating and bio-photochemical cycling. We recommend radiative transfer models applied to bare ice in the Greenland Ice Sheet ablation zone account for vertical variation in light attenuation due to the vertical distribution of light absorbing particles and ice microstructure, and we provide new values of flux attenuation, absorption, and scattering coefficients to support model validation and parameterization.