MATLAB code for firn thermodynamic and hydrological modelling

• Main Author: Marshall, Shawn
• Format: Dataset
• Language: unknown
• Published: 2021
• Subjects: firn; Earth and Environmental Sciences; code; thermodynamics; ice sheet; hydrology; glacier; energy balance; Ice Sheet
• Online Access: https://search.dataone.org/view/sha256:6ca269e83ff5d6cccc175c5dc6651d36a8e3829c8d466e4dab7617b82f101edf

This repository contains MATLAB source code for a model of glacier surface energy balance, coupled with a subsurface (snow and firn) thermodynamic and hydrological evolution. The default configuration of the subsurface model is for a 35-m firn column, with 0.1-m thick layers from 0-0.6 m, 0.2-m thick layers from 0.6-2 m, 0.4-m thick layers from 2-10 m, and 1-m thick layers below that. The model requires forcing from mean daily automatic weather station, climate model, or climate reanalysis data, with the default inputs being: minimum, mean, and maximum daily air temperature, mean daily relative or specific humidity, wind speed, air pressure, and longwave radiation, and mean and maximum daily incoming incoming shortwave. The model internally calculates snow albedo and snow surface temperature (hence, outgoing longwave radiation), as well as conductive heat flux to the snow surface, based on temperature gradients in the top three layers of the subsurface snow/firn model. See the referenced papers for further details on the model physics and parameterizations.