Landfast ice constitutes a significant portion of sea ice mass in the central Arctic and Canadian Archipelago. It is instrumental in protecting coastlines from erosion, limiting coastal shipping, and modulating Arctic Ocean hydrography.  I present two different models for landfast ice simulation in E3SM;  both models differ from existing landfast ice schemes in order to circumvent significant uncertainties in the stress of keels against the ocean floor.  The InteRFACE Empirical Derichlet Model applies a static boundary condition surrounding grounded ice grid cells based on an analysis of 64,000 ridges off the Alaskan coast, using an empirical relationship between maximum keel depth and deformed rubble.  This approach uses results from a variational model to ground sea ice, typically within the 20m isobath, and then constrains ice movement in the vicinity of grounded keels.  The empirical model differs from previous landfast ice schemes by using advanced ridge analysis to enforce grounding.  The second landfast ice model takes this approach a step further to inversely derive basal stress acting on keels using a full variational model. The models are currently being coded into E3SM, applying global mesh constraints to the native unstructured sea ice grid.