Purpose: Model-assisted probability of detection (MAPOD) is an important approach used as part of assessing the reliability of nondestructive testing systems. The purpose of this paper is to apply the polynomial chaos-based Kriging (PCK) metamodeling method to MAPOD for the first time to enable efficient uncertainty propagation, which is currently a major bottleneck when using accurate physics-based models. Design/methodology/approach: In this paper, the state-of-the-art Kriging, polynomial chaos expansions (PCE) and PCK are applied to "a^ vs a"-based MAPOD of ultrasonic testing (UT) benchmark problems. In particular, Kriging interpolation matches the observations well, while PCE is capable of capturing the global trend accurately. The proposed UP approach for MAPOD using PCK adopts the PCE bases as the trend function of the universal Kriging model, aiming at combining advantages of both metamodels. Findings: To reach a pre-set accuracy threshold, the PCK method requires 50 per cent fewer training points than the PCE method, and around one order of magnitude fewer than Kriging for the test cases considered. The relative differences on the key MAPOD metrics compared with those from the physics-based models are controlled within 1 per cent. Originality/value: The contributions of this work are the first application of PCK metamodel for MAPOD analysis, the first comparison between PCK with the current state-of-the-art metamodels for MAPOD and new MAPOD results for the UT benchmark cases.