Ionospheric-free(IF)linear combination(LC)is usually used to eliminate the first-order ionospheric delay error of low Earth orbit(LEO)satellite reduced-dynamic precision orbit determi-nation(POD),ignoring the high-order ionospheric(HOI)delay error.With the development of LEO satellite POD technology,calculating the HOI delay at different orbital altitude and exploring its changes has become an important means to further improve POD accuracy.Firstly,using the international reference ionosphere-2016(IRI-2016)and international geomagnetic reference field:the 13th generation(IGRF-13)models,the ionospheric pierce point(IPP)position and geomag-netic field intensity are calculated.Secondly,the ionospheric slant total electron content(STEC)is calculated using smooth satellite-borne GNSS data.Then,the second-order and third-order ion-ospheric delay of GOCE,GRACE-A,and SWARM-A/B satellites are calculated,respectively.Fi-nally,the effects of HOI delay on LEO satellite overlapping orbit analysis,satellite laser ranging(SLR)verification,and precision science orbit(PSO)comparison results are evaluated.The ex-perimental results show that the effects of HOI delay on the reduced-dynamic POD of LEO satel-lites is approximately on the order of millimeters to centimeter level;the effects of HOI delay on the accuracy of reduced-dynamic POD for LEO satellites is 0.92,0.22,0.21,and 0.18 mm,re-spectively;as the orbital altitude of LEO satellites increases,the effects of HOI delay on the re-duced-dynamic POD of LEO satellites decreases.