The Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA's Terra and Aqua satellites conducts continuous monitoring of the Earth's land surface and oceans. Recently, a sharp discontinuity (averaging 1.9°C) has been noticed at 60° N in both MODIS daytime and night-time land surface temperature (LST) products. This linear artefact arises because the CO₂ high cloud test in the operational code for the generation of the MODIS cloud mask product is used only in the non-polar region (between 60° N and 60° S). The resulting discontinuity clearly has negative implications for any statistical applications of these temperature data. In this technical note we present a new algorithm, which minimizes this discontinuity. The method uses edge detection and elimination based on a mixture of Sobel and non-linear Laplacian filters (edge detection and quantification), cubic splines (edge modelling), and a controllable power function for image restoration. The implementation of this algorithm is demonstrated on an image of average minimum night-time LST between 2001 and 2008.