Ice sheets are covered by a firn layer, which is the transition stage between fresh snow and ice. Accurate modelling of firn density properties is important in many glaciological aspects. Current models show disagreements, are mostly calibrated to match specific observations of firn density and lack thorough uncertainty analysis. We use a novel calibration method for firn models based on a Bayesian statistical framework, which results in improved model accuracy and in uncertainty evaluation.