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.