To remove atmospheric pressure loading (ATML) effect from GNSS coordinate time series, surface pressure (SP) models are required to predict the displacements. In this paper, we modeled the 3D ATML surface displacements using the latest MERRA-2 SP grids, together with four other products (NCEP-R-1, NCEP-R-2, ERA-Interim and MERRA) for 596 globally distributed GNSS stations, and compared them with ITRF2014 residual time series. The five sets of ATML displacements are highly consistent with each other, particularly for those stations far away from coasts, of which the lowest correlations in the Up component for all the four models w.r.t MERRA-2 become larger than 0.91. ERA-Interim-derived ATML displacement performs best in reducing scatter of the GNSS height for 90.3% of the stations (89.3% for NCEP-R-1, 89.1% for NCEP-R-2, 86.4% for MERRA and 85.1% for MERRA-2). We think that this may be possibly due to the 4D variational data assimilation method applied. Considering inland stations only, more than 96% exhibit WRMS reduction in the Up direction for all five models, with an average improvement of 3–4% compared with the original ITRF2014 residual time series before ATML correction. Most stations (> 67%) also exhibit horizontal WRMS reductions based on the five models, but of small magnitudes, with most improvements (> 76%) less than 5%. In particular, most stations in South America, South Africa, Oceania and the Southern Oceans show larger WRMS reductions with MERRA-2, while all other four SP datasets lead to larger WRMS reduction for the Up component than MERRA-2 in Europe. Through comparison of the daily pressure variation from the five SP models, we conclude that the bigger model differences in the SP-induced surface displacements and their impacts on the ITRF2014 residuals for coastal/island stations are mainly due to the IB correction based on the different land–sea masks. A unique high spatial resolution land–sea mask should be applied in the future, so that model differences would come from only SP grids. Further research is also required to compare the ATML effect in ice-covered and high mountainous regions, for example the Qinghai–Tibet Plateau in China, the Andes in South America, etc., where larger pressure differences between models tend to occur.