Monthly gravity field recovery using data from the GRACE and GRACE Follow-On missions includes errors limiting the spatial and temporal resolution of the estimated gravity fields. The major error contributions, besides the noise of the accelerometer instruments, arise from temporal aliasing errors due to imperfections in the non-tidal atmospheric and oceanic de-aliasing models and ocean tide models. We derive uncertainty information for the eight major tidal constituents from five different ocean tide models and introduce it into the gravity field recovery process in terms of a constrained normal equation system while expanding the parameter space by additional tidal parameters to be adjusted. We prove the effectiveness of the ocean tide variance-covariance information through realistic simulations and we assess its potential based on microwave and laser interferometer observations from the GRACE Follow-On mission. We show that errors are reduced by more than 20% ocean wRMS for a Gaussian filter radius of 300 km if uncertainty information for ocean tides is considered and stochastic modeling of instrument errors is applied, compared to the latest GFZ release 6.1. Our results also show the limited visibility of the effectiveness of the ocean tide variance-covariance information due to the dominating error contribution of non-tidal atmospheric and oceanic mass variations. Additionally, we investigate the option of estimating ocean tide parameters over a 1-year period while including ocean tide uncertainty information in order to improve ocean tide background modeling.
