PurposeThis paper aims to propose a method of the safety boundary protection for unmanned aerial vehicles (UAVs) in the icing conditions.Design/methodology/approachForty icing conditions were sampled in the continuous maximum icing conditions in the Appendix C of the Federal Aviation Regulation Part 25. Icing numerical simulations were carried out for the 40 samples and the anti-icing thermal load distribution in full evaporation mode were obtained. Based on the obtained anti-icing thermal load distribution, the surrogated model of the anti-icing thermal load distribution was established with proper orthogonal decomposition and Kriging interpolation. The weather research and forecasting (WRF) model was used for meteorological simulations to obtain the icing meteorological conditions in the target area. With the obtained icing conditions and surrogated model, the anti-icing thermal load distribution in the target area and the variation with time can be determined. According to the energy supply of the UAVs, the graded safety boundaries can be obtained.FindingsThe surrogated model can predict the effects of five factors, such as temperature, velocity, pressure, median volume diameter (MVD) and liquid water content (LWC), on the anti-icing thermal load quickly and accurately. The simulated results of the WRF mode agree well with the observed results. The method can obtain the graded safety boundaries.Originality/valueThe method has a reference significant for the safety of the UAVs with the limited energy supply in the icing conditions.