The Global Navigational Satellite System (GNSS) technique is naturally sensitive to the geocenter motion, similar to all satellite techniques. However, the GNSS-based estimates of the geocenter used to contain more orbital artifacts than the geophysical signals, especially for the Z component of the geocenter coordinates. This contribution conveys a discussion on the impact of solar radiation pressure (SRP) modeling on the geocenter motion estimates. To that end, we process 3 years of GPS, GLONASS, and Galileo observations (2017–2019), collected by a globally distributed network of the ground stations. All possible individual system-specific solutions, as well as combinations of the available constellations, are tested in search of characteristic patterns in geocenter coordinates. We show that the addition of a priori information about the SRP-based forces acting on the satellites using a box-wing model mitigates a great majority of the spurious signals in the spectra of the geocenter coordinates. The amplitude of the 3 cpy (about 121 days) signal for GLONASS has been reduced by a factor of 8.5. Moreover, the amplitude of the spurious 7 cpy (about 52 days) signal has been reduced by a factor of 5.8 and 3.1 for Galileo and GPS, respectively. Conversely, the box-wing solutions indicate increased amplitudes of the annual variations in the geocenter signal. The latter reaches the level of 10–11 mm compared to 4.4 and 6.0 mm from the satellite laser ranging observations of LAGEOS satellites and the corresponding GNSS series applying extended empirical CODE orbit model (ECOM2), respectively. Despite the possible improvement in the GLONASS-based Z component of the geocenter coordinates, we show that some significant power can still be found at periods other than annual. The GPS- and Galileo-based estimates are less affected; thus, a combination of GPS and Galileo leads to the best geocenter estimates.