Abstract. While a lot of attention has been given to the numerical implementation of grounding lines and basal friction in the grounding zone, little has been done about the impact of the numerical treatment of ocean-induced basal melting in this region. Several strategies are currently being employed in the ice sheet modeling community, and the resulting grounding line dynamics may differ strongly, which ultimately add significant uncertainty to the projected contribution of marine ice sheets to sea level rise. We investigate here several implementations of basal melt parameterization on partially floating elements in a finite element framework, based on the Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP) setup: (1) melt applied only to entirely floating elements, (2) melt applied over the entire elements that are crossed by the grounding line, and (3) melt integrated partially over the floating portion of a finite element using two different sub-element integration methods. All methods converge towards the same state when the mesh resolution is fine enough. However, (2) and (3) will systematically overestimate the rate of grounding line retreat in coarser resolutions, while (1) converges faster to the solution in most cases. The differences between sub-element parameterizations are exacerbated for experiments with large melting rates in the vicinity of the grounding line and for a Weertman sliding law. As most real-world simulations use horizontal mesh resolutions of several hundreds of meters at best, and large melt rates are generally present close to the grounding lines, we recommend using (1) to avoid overestimating the rate of grounding line retreat.