Abstract. As surface melt is increasing on the Greenland ice sheet (GrIS), quantifying the retention capacity of the firn layer is critical to link meltwater production to meltwater runoff. Firn-densification models have so far relied on empirical approaches to account for the percolation-refreezing process, and more physically based representations of liquid water flow might therefore bring improvements to model performance. Here we implement three types of water percolation schemes into the Community Firn Model: the tipping bucket approach, the Richards Equation in a single-domain and the Richards Equation in a dual-domain, which accounts for partitioning between matrix and fast preferential flow. We investigate their impact on firn densification at four locations on the GrIS and compare model results with observations. We find that for all of the flow schemes, significant discrepancies remain with respect to observed firn density, particularly the density variability in depth, and that inter-model differences are large. The simple bucket scheme is as efficient in replicating observed density profiles as the single-domain Richards Equation. The most physically detailed dual-domain scheme does not necessarily reach best agreement with observed data. However, we find that the implementation of preferential flow does allow for more frequent ice layer formation and for deeper percolation. We also find that the firn model is more sensitive to the choice of densification scheme than to the choice of water percolation scheme. The disagreements with observations and the spread in model results demonstrate that progress towards an accurate description of water flow in firn is necessary. The numerous uncertainties surrounding firn micro- and macro-structure, its hydraulic properties, and the one dimensionality of firn models render the implementation of physically based percolation schemes difficult. An improved understanding of the parameters affecting evolution of polar firn, of the effects of the climatic forcing on the densification process and more accurate treatment of liquid water would benefit further developments.;
