This study focuses on the kriging based metamodeling for the prediction of parameter-dependent mode coupling instabilities. The high cost of the currently used parameter-dependent Complex Eigenvalue Analysis (CEA) has induced a growing need for alternative methods. Hence, this study investigates capabilities of kriging metamodels to be a suitable alternative. For this aim, kriging metamodels are proposed to predict the stability behavior of a four-degree-of-freedom mechanical system submitted to friction-induced vibrations. This system is considered under two configurations defining two stability behaviors with coalescence patterns of different complexities. Efficiency of kriging is then assessed on both configurations. In this framework, the proposed kriging surrogate approach includes a mode tracking method based on the Modal Assurance Criterion (MAC) in order to follow the physical modes of the mechanical system. Based on the numerical simulations, it is demonstrated by a comparison with the reference parameter-dependent CEA that the proposed kriging surrogate model can provide efficient and reliable predictions of mode coupling instabilities with different complex patterns.