A model of the weathering crust and microbial activity on an ice-sheet surface

Shortwave radiation penetrating beneath an ice-sheet surface can cause internal melting and the formation of a near-surface porous layer known as the weathering crust, a dynamic hydrological system that provides home to impurities and microbial life. We develop a mathematical model, incorporating thermodynamics and population dynamics, for the evolution of such layers. The model accounts for conservation of mass and energy, for internal and surface-absorbed radiation, and for logistic growth of a microbial species mediated by nutrients that are sourced from the melting ice. It also accounts for potential melt-albedo and microbe-albedo feedbacks, through the dependence of the absorption coefficient on the porosity or microbial concentration. We investigate one-dimensional steadily melting solutions of the model, which give rise to predictions for the weathering crust depth, water content, melt rate, and microbial abundance, depending on a number of parameters. In particular, we examine how these quantities depend on the forcing energy fluxes, finding that the relative amounts of shortwave (surface-penetrating) radiation and other heat fluxes are particularly important in determining the structure of the weathering crust. The results explain why weathering crusts form and disappear under different forcing conditions and suggest a range of possible changes in behaviour in response to climate change.