Nonlinear variations in the coordinate time series of global navigation satellite system (GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects, including atmospheric, hydrological, and nontidal ocean loading. Continuous improvements in the accuracy of surface mass loading products, performance of Earth models, and precise data-processing technologies have significantly advanced research on the effects of environmental loading on nonlinear variations in GNSS coordinate time series. However, owing to theoretical limitations, the lack of high spatiotemporal resolution surface mass observations, and the coupling of GNSS technology-related systematic errors, environmental loading and nonlinear GNSS reference station displacements remain inconsistent. The applicability and capability of these loading products across different regions also require further evaluation. This paper outlines methods for modeling environmental loading, surface mass loading products, and service organizations. In addition, it summarizes recent advances in applying environmental loading to address nonlinear variations in global and regional GNSS coordinate time series. Moreover, the scientific questions of existing studies are summarized, and insights into future research directions are provided. The complex nonlinear motion of reference stations is a major factor limiting the accuracy of the current terrestrial reference frame. Further refining the environmental load modeling method, establishing a surface mass distribution model with high spatiotemporal resolution and reliability, exploring other environmental load factors such as ice sheet and artificial mass-change effects, and developing an optimal data-processing model and strategy for reprocessing global reference station data consistently could contribute to the development of a millimeter-level nonlinear motion model for GNSS reference stations with actual physical significance and provide theoretical support for establishing a terrestrial reference frame with 1 mm accuracy by 2050.