 This paper proposes a new method to remove GPS systematic errors from K1 and K2 tidal constituents. The authors use quadratic fitting to determine smooth functions of tidal frequencies for tidal admittances, allowing them to interpolate K1 and K2 tides from GPS data. This method was tested with three GPS stations, and the results showed that the corrected K1 and K2 tides were much closer to the modeled K1 and K2 tides from FES 2014, one of the most accurate tide models. Additionally, the proposed method can be used to correct GPS systematic errors in other areas due to its validity for global tidal signals. This article was authored by Haidong Pan, Xiaoqing Su, Huaizhang, and others.