The optimal combination of global positioning system (GPS) geometric heights with gravimetrically derived geoid undulations for the determination of orthometric heights above mean sea level (GPS/geoid leveling), or more precisely with respect to a vertical geodetic datum, requires the incorporation of a parametric corrector surface model. Such a parametric model is needed to absorb the datum inconsistencies and systematic distortions inherent among the different types of height data. An analysis for the achievable accuracy of relative GPS/geoid leveling is performed in this paper, using the covariance (CV) matrix of the estimated parameters in the corrector surface model, the standard measuring accuracy of GPS height differences, and the relative internal accuracy of a gravimetric geoid model. A test network of spirit leveled GPS control points located in the southwestern part of Canada has been used to perform a simulative integrated least-squares adjustment of all three types of height data, and to obta...