We present a global model of ocean tide loading (OTL)-induced surface deformation for the M2 harmonic constituent. The model quantifies three-dimensional displacement amplitudes and phases across the horizontal (west and south) and vertical (up) components. Computations were performed using the advanced VILMA-E software (Tanaka et al., 2021; Huang et al., 2021), which integrates the TPXO9.5 global tide model (Egbert & Erofeeva, 2002) with a three-dimensional anelastic Earth structure optimized for OTL observations (Huang et al., 2025). The final output is provided as a high-resolution global grid (0.1° × 0.1°) in netCDF format, compatible with standard geospatial processing tools such as NCO and CDO for efficient data access and analysis.
As a supplement to Huang et al. (2022) “The influence of sediments, lithosphere and upper mantle (anelastic) with lateral heterogeneity on ocean tide loading and ocean tide dynamics”, we provide for the advanced earth model LH-Lyon-3Dae [consisting of 3D elastic sediments, lithosphere and 3D anelastic upper mantle structures, see Huang et al.(2022) for details] the solutions of vertical ocean tide loading (OTL) displacement, self-attraction and loading (SAL) elevation, and ocean tides. Solutions for three tidal constituents, i.e., M2, K1 and Mf, are given. As a comparison, solutions based on the 1D elastic model PREM and the 1D anelastic LH-Lyon-1Dae are also presented. With these solutions, the primary results in Huang et al. (2022) such as the model amplitude differences, RMS differences and the predictions in GNSS stations can be reconstructed.