GPS applications such as Precise Point Positioning (PPP) require the availability of precise ephemeris at high rate. To support these applications, several institutions such as the International GNSS Service (IGS) have developed precise orbital service. Unfortunately, however, the data rate of such precise orbits is usually limited to 15 minutes. To overcome this problem, a number of orbital interpolation methods are proposed. This paper examines the performance of four interpolation methods for IGS precise GPS orbits, namely Lagrange, Newton Divided Difference, Cubic Spline and Trigonometric interpolation. In addition, the paper discusses a new approach, which utilizes the residuals between the broadcast and precise ephemeris to generate a high density precise ephemeris. It is shown that the new approach produces better results than previously reported orbital interpolation accuracy.