Internal solitary waves in polar regions have attracted much interest recently. It is important to understand how sea ice affects them as this may have a profound influence on human activities and the environment. In this study, experiments on internal solitary waves with and without two types of sea ice (ice sheet and ice keel) are presented, as well as corresponding simulations using the Korteweg-de Vries (KdV) equation, the Benjamin-Ono (BO) equation, and the variable-coefficient Korteweg-de Vries (vKdV) equation, which is a derivation of the KdV equation. Comparison between experiments without sea ice and simulations using the KdV and BO equations proves the suitability of the former over the latter for this study. The experiments with sea ice and theoretical simulations using the vKdV equation provide evidence for wave deformation, oscillation occurring in the rear of the wave, and a decrease in amplitude. The latter suggests possibilities of energy dissipation or the emission of small amplitude linear waves. The sharp vertices of the ice result in occasional inconsistency with the vKdV predictions. Nonetheless, the vKdV equation is still suitable for modeling internal solitary waves under sea ice, giving generally accurate results that can assist further studies. This is the first time the vKdV equation has been applied to investigate the impacts of sea ice on internal solitary waves.