In this article, we explore the use of multivariate Archimedean copulas in modelling the mortality dependence between different countries and pricing mortality bonds. We study the fitting performance of multi-dimensional, fully nested, and partially nested Archimedean copulas and test 11 types of generators and two skewed distributions. To evaluate their practical usefulness, we adopt the fitted models to compute the market prices for some typical mortality bond structures. The results show that the copula assumption has a significant impact on the calculation of the prices of mortality-linked securities and the management of extreme mortality risks.

We propose a two-regime Markov switching copula to depict the evolution of mortality dependence. One regime represents periods of high dependence and the other regime represents periods of low dependence. Each regime features a regular vine (R-vine) copula that, built on bivariate copulas, provides great flexibility for modelling complex high-dimensional dependence. Our estimated model indicates that the years of recovery from extreme mortality deterioration and the years of health care reform more likely fall into the low regime, while the years in which extreme mortality deteriorating events break out and the peaceful years without major mortality-impacting events more likely fall into the high regime. We use a case study to illustrate how the regime-switching copula can be applied to assess the effectiveness of longevity risk hedge with different beliefs about future mortality dependence evolution incorporated.