To offer more reliable and precise ultra-rapid orbits is very necessary for improving the performance of multi-constellation Global Navigation Satellite Systems GNSS (multi-GNSS) real-time (RT) precise positioning service. After the significance of high order ionospheric delays remaining in the ionosphere-free observations was confirmed, such corrections are suggested in precise GNSS data processing. However, in the precise orbit determination (POD) for generating ultra-rapid orbits, the higher order corrections are usually not considered most likely because a RT ionospheric model needed for calculating the higher order corrections is hardly available or the impact is believed rather small compared to the accuracy of the predicted orbit. In this contribution, we firstly analyze the characteristics in temporal-spatial variations of higher order ionospheric effects using data collected from International GNSS Service (IGS) global ground stations and fluctuations up to centimeter-level are detected during a high ionospheric activity period. Afterwards, we evaluate the higher order ionospheric effects on multi-GNSS ultra-rapid orbit determination. The improvement in GNSS ultra-rapid orbits by applying higher order ionospheric corrections is shown by the reduced orbit difference over the overlapping period. Furthermore, the improvement due to the higher order corrections is also validated by employing satellite laser ranging (SLR) data.