Two-dimensional interpolation – or surface fitting – is an approximation tool with applications in geodetic datum transformations, terrain modelling and geoid determination. It can also be applied to many other forms of geographic point data, including rainfall, chemical concentrations and noise levels. The problem of fitting of a smooth continuous interpolant to a bivariate function is particularly difficult if the dataset of control points is scattered irregularly. A typical approach is a weighted sum of data values where the sum of the weights is always unity. Weighting by inverse distance to a power is one approach, although a power greater than 1 is needed to ensure smooth results. One advantage over other methods is that data values can be incorporated into the interpolated surface. One disadvantage is the influence of distant points. A simple cut-off limit on distance would affect continuity. This study proposes a transition range of accelerated decline by means of an adjoining polynomial. This preserves smoothness and continuity in the interpolating surface. Case studies indicate accuracy advantages over standard versions of inverse-distance weighting.