Radiative transfer model of snow for bare ice regions

Modeling a radiative transfer (RT) for coupled atmosphere-snow-bare ice systems is of fundamental importance for remote sensing applications to monitor snow and bare ice regions in the Greenland ice sheet and for accurate climate change predictions by regional and global climate models. Recently, the RT model for atmosphere-snow system was implemented for our regional and global climate models. However, the bare ice region where recently it has been expanded on the Greenland ice sheet due to the global warming, has not been implemented for these models, implying that this region leads miscalculations in these climate models. Thus, the RT model of snow for bare ice regions is needed for accurate climate change predictions. We developed the RT model for coupled atmosphere-snow-bare ice systems, and conducted a sensitivity analysis of the RT model to know the effect of snow, bare ice and geometry parameters on the spectral radiant quantities. The RT model considers snow and bare-ice inherent optical properties (IOPs), including snow grain size, air bubble size and its concentration and bare ice thickness. The conventional light scattering theory, Mie theory, was used for IOP calculations. Monte Carlo method was used for the multiple scattering. The sensitivity analyses showed that spectral albedo for the bare ice increased with increasing the concentration of the air bubble in the bare ice for visible wavelengths because the air bubble is scatterer with no absorption. For near infrared wavelengths, spectral albedo has no dependence on the air bubble due to the strong light absorption by ice. When increasing solar zenith angle, the spectral albedo were increased for all wavelengths. This is the similar trend with spectral snow albedo. Cloud cover influenced the bare ice spectral albedo by covering direct radiation into diffuse radiation. The purely diffuse radiation has an effective solar zenith angle near 50°. Converting direct into diffuse radiation reduces the effective solar zenith angle, resulting in reducing the spectral albedo. This is also the similar trend with spectral snow albedo. Further work should focus on the validation of the RT model using in situ measurement data through field and laboratory experiments.