Design of fair surfaces over irregular domains is a fundamental problem in computer-aided geometric design (CAGD), and has applications in engineering sciences (in aircraft, automobile, ship science etc.). In the design of fair surfaces over irregular domains defined over scattered data, it was widely accepted till recently that the classical Delaunay triangulation be used because of its global optimum property. However, in recent times it has been shown that for continuous piecewise linear surfaces, improvements in the quality of fit can be achieved if the triangulation pattern is made dependent upon some topological or geometric property of the data set or is simply data dependent. The fair surface is desired because it ensures smooth and continuous surface planar cuts, and these in turn ensure smooth and easy production of the surface in CAD/CAM, and favourable resistance properties. In this paper, we discuss a method for construction of C1 piecewise polynomial parametric fair surfaces which interpolate prescribed ℜ3 scattered data using spaces of parametric splines defined on H3 triangulation. We show that our method is more specific to the cases when the projection on a 2-D plane may consist of triangles of zero area, numerically stable and robust, and computationally inexpensive and fast. Numerical examples dealing with surfaces approximated on plates, and on ships have been presented.