The migration characteristics of nanoplastics (NPs) in the natural freezing process are complex and have attracted increasing attention in simulating natural freezing in recent years. However, simulated freezing conditions often fall short of replicating natural freezing processes, and studies on the vertical distribution of NPs remain inadequate. This study established a more realistic simulation of the natural freezing process in surface water by controlling both the air temperature (T1) and the water temperature (T2). Additionally, we introduced a new parameter, the local distribution coefficient (Kiw1), to compare with the effective distribution coefficient (Kiw2). The values of Kiw1 and Kiw2 for PS-500 nm were 0.18 and 0.21, respectively, at T1 = -20 °C and T2 = 1 °C. The results revealed the NPs concentration differed in ice, near-ice liquid, and far-ice liquid. Both properties of NPs and environmental factors could regulate the vertical distribution of NPs. The findings underscored the importance of freezing temperature regulated by T1 and T2, elucidating the roles of various influencing factors on the vertical distribution characteristics of NPs and unraveling the mechanisms of NPs distribution in the ice-water system. This study can provide valuable insights for understanding the migration of NPs in surface water in cold regions.