In the Jutulgryta area of Dronning Maud Land, Antarctica, subsurface melting ofthe ice sheet has been observed. The melting takes place during the summermonths in blue-ice areas under conditions of below-freezing air and surfacetemperatures. Adjacent snow-covered regions, having the same meteorological andclimatic conditions, experience little or no subsurface melting. To help explainand understand the observed melt-rate differences in the blue-ice andsnow-covered areas, a physically based numerical model of the coupledatmosphere, radiation, snow and blue-ice system has been developed. The modelcomprises a heat-transfer equation which includes a spectrally dependentsolar-radiation source term. The penetration of radiation into the snow and blueice depends on the solar-radiation spectrum, the surface albedo and the snow andblue-ice grain-sizes and densities. In addition, the model uses a completesurface energy balance to define the surface boundary conditions. It is run overthe full annual cycle, simulating temperature profiles and melting and freezingquantities throughout the summer and winter seasons. The model is driven andvalidated using field observations collected during the Norwegian AntarcticResearch Expedition (NARE) 1996–97. The simulations suggest that theobserved differences between subsurface snow and blue-ice melting can beexplained largely by radiative and heat-transfer interactions resulting fromdifferences in albedo, grain-size and density between the two mediums.