Precise point positioning (PPP) technique is practical in estimating the zenith tropospheric delay (ZTD). To comprehensively study the factors affecting its accuracy in the polar region, different elevation cutoff angles for each single system and multi-system combination, different tropospheric mapping functions, tropospheric gradients and also seasonal variations are assessed. Meanwhile, European Centre for Medium-Range Weather Forecasts (ECMWF) products are adopted to validate the accuracy of PPP-derived ZTD. The results show that the root mean square (RMS) value of the estimated ZTD and ECMWF products is around 4 mm when the elevation cutoff angle is set to 10 degrees or less. From the view of system combination, it improves the ZTD accuracy in comparison with each individual system. The dual-system combination has 21.8%, 19.0% and 19.1% improvements in its ZTD accuracy compared with the Global navigation satellite system (GLONASS), Beidou navigation satellite system (BDS) and Galileo satellite navigation (GAL) systems. The four-system combination can reach 7.7% further improvement in ZTD estimation than that of the dual-system. For mapping functions, Vienna mapping function 3 (VMF3) function has the best accuracy among the four mapping functions of Niell mapping function (NMF), global mapping function, VMF1 and VMF3, while NMF does not perform as well as others. Considering the tropospheric gradients could get higher accuracy in ZTD estimation than ignoring them. The tropospheric delay estimated in winter is more accurate than that in summer and specifically the RMS value in winter is about 0.9 mm lower than that in summer. Overall, from investigation the estimation of ZTD in the polar region can be optimally handled with a multi-system and VMF3 mapping function at an elevation cutoff angle of 7°, considering the tropospheric gradients.