<jats:title>Abstract</jats:title><jats:p>Firn describes the interstage product between snow and ice in cold regions of the earth, where annual snow fall exceeds the amount of snow melting. The continuing accumulation of snow leads to its densificiation due to overburden stress until it becomes ice. In the field of glaciology various attempts on simulating firn densification have been made and new models are still developed, as the knowledge of the firn column's density structure allows important derivations.</jats:p><jats:p>The presented study reassesses a model description for low density firn based on the process of grain boundary sliding presented by Alley in 1987 [1] using an optimisation approach. By comparing simulation results to 159 measured firn density profiles from Greenland and Antarctica it finds a possible additional dependency of the constitutive relation on the mean surface mass balance. This result is interpreted as an insufficient description of the stress regime.</jats:p>