Geocenter motion estimated from GRACE orbits: The impact of F10.7 solar flux

We assess the impact of orbit modeling on the origin offsets between GRACE kinematic and reduced-dynamic orbits. The origin of the kinematic orbit is the center of IGS network (CN), whereas the origin of the reduced-dynamic orbit is assumed to be the center of mass of the Earth (CM). Theoretically, the origin offset between these two orbits is associated with the geocenter motion. However, the dynamic property of the reduced-dynamic orbit is highly related to orbit parameterizations. The assessment of the F10.7 impact on the geocenter motion is implemented by using different orbit parameterization setups in the reduced-dynamic method. We generate two types of reduced-dynamic orbits using 15 and 240 empirical parameters per day from 2005 to 2012. The empirical parameter used in Bernese GNSS Software is called piece-wise constant empirical acceleration (PCA) and is mainly to absorb the non-gravitational forces mostly related to the atmospheric drag and solar radiation pressure. The differences between kinematic and dynamic orbits can serve as a measurement for geocenter. The RMS value of the geocenter measurement in the 15-PCA case is approximately 3.5 cm and approximately 2 cm in the 240-PCA case. The correlation between the orbit difference and F10.7 is about 0.90 in the 15-PCA case and −0.10 to 0 in the 240-PCA case. This implies that the reduced-dynamic orbit modeled with 240 PCAs absorbs the F10.7 variation, which aliases to the 15-PCA orbit solution. The annual amplitudes of the geocenter motion are 3.1, 3.1 and 2.5 mm in the 15-PCA case, compared to 0.9, 2.0 and 1.3 mm in the 240-PCA case in the X, Y and Z components, respectively. The 15-PCA solution is thus closer to the geocenter motions derived from other space-geodetic techniques. The proposed method is limited to the parameterizations in the reduced-dynamic approach.