The covariance of real-time global positioning system (GPS) orbits has been drawing attention in various fields such as user integrity, navigation performance improvement, and fault detection. The international global navigation satellite system (GNSS) service (IGS) provides real-time orbit standard deviations without correlations between the axes. However, without correlation information, the provided covariance cannot assure the performance of the orbit product, which would, in turn, causes significant problems in fault detection and user integrity. Therefore, we studied real-time GPS orbit covariance characteristics along various coordinates to effectively provide conservative covariance. To this end, the covariance and precise orbits are estimated by means of an extended Kalman filter using double-differenced carrier phase observations of 61 IGS reference stations. Furthermore, we propose a new method for providing covariance to minimize loss of correlation. The method adopted by the IGS, which neglects correlation, requires 4.5 times the size of the covariance to bind orbit errors. By comparison, our proposed method reduces this size from 4.5 to 1.3 using only one additional parameter. In conclusion, the proposed method effectively provides covariance to users.