Contributions of terrestrial and GRACE data to the study of the secular geoid changes in North America

This paper tests and discusses different statistical methods for modelling secular rates of change of the geoid in North America. In particular, we use the method of principal component/empirical orthogonal functions (PC/EOF) analysis to model the geoid rates from Gravity Recovery and Climate Experiment (GRACE) satellite data. As demonstrated, the PC/EOF analysis is useful for studying the contributions from different signals (mainly residual hydrology signals and leakage effects) to the GRACE-derived geoid rates. The PC/EOF analysis leads to smaller geoid rates compared to the conventional least-squares fitting of a trend and annual and semi-annual cycles to the time series of the spherical harmonic coefficients. This is because we filter out particular spatiotemporal modes of the regional geoid changes.We apply the method of least-squares collocation with parameters to combine terrestrial data (GPS vertical velocities from the Canadian Base Network and terrestrial gravity rates from the Canadian Gravity Standardization Net) with the GRACE-derived vertical motion to obtain again the geoid rates. The combined model has a peak geoid rate of 1.4 mm/year in the southeastern area of Hudson Bay contrary to the GRACE-derived geoid rates that show a large peak of 1.6–1.7 mm/year west of Hudson Bay. We demonstrate that the terrestrial data, which have a longer time span than the GRACE data, are important for constraining the GRACE-derived secular signal in the areas that are well sampled by the data.