The effect of the non-tidal ocean load (NTOL) and the continental water load (CWL) on crustal deformation are calculated for global satellite laser ranging (SLR) stations and on 4!×4! grids (only over the land). For the regions most severely affected, the peak-to-peak displacements due to the NTOL are found to be as large as 3 mm for the horizontal components and 10 mm for the vertical component. The peak-to-peak displacements due to the CWL reach 3 mm for the horizontal components and 15 mm for the vertical component. We apply the time series of NTOL and CWL to precise SLR analysis. The LAGEOS orbit analysis reveals that the Estimating the Circulation and Climate of the Ocean (ECCO) model makes the root mean square (RMS) of the range residual 0.2% smaller, and that the CWL makes it 0.8% smaller, compared with the case where loading displacement is neglected. On the other hand, with the NTOL derived from Topex/Poseidon altimetry data, the SLR orbit fit is not improved.