The Swarm mission consists of three satellites orbiting the Earth at low orbital altitudes. The onboard GPS receivers, star cameras, and laser retro-reflectors make the Swarm mission an interesting candidate to explore the contribution of Swarm GPS data to the recovery of both the static and time-variable gravity fields. We use 1.5 years of Swarm GPS and attitude data to generate kinematic positions of high quality to perform gravity field determination using the Celestial Mechanics Approach. The generated gravity fields reveal severe systematic errors along the geomagnetic equator. Their size is correlated with the ionospheric density and thus strongly varying over the analyzed time period. Similar to the findings of the GOCE mission, the systematic errors are related to the Swarm GPS carrier phase data and may be reduced by rejecting GPS data affected by large ionospheric changes. Such a measure yields a strong reduction of the systematic errors along the geomagnetic equator in the gravity field recovery. Long wavelength signatures of the gravity field may then be recovered with a similar quality as achieved with GRACE GPS data, which makes the Swarm mission well suited to bridge a potential gap between the current GRACE and the future GRACE Follow-On mission.

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only ≈10 %; in the absence of the polar data gap, SLR residuals agree by ≈96 %. From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.