Delta-hedging longevity risk under the M7–M5 model: The impact of cohort effect uncertainty and population basis risk

In a recent project commissioned by the Institute and Faculty of Actuaries and the Life and Longevity Markets Association, a two-population mortality model called the M7–M5 model is developed and recommended as an industry standard for the assessment of population basis risk. In this paper, we contribute a delta hedging strategy for use with the M7–M5 model, taking into account of not only period effect uncertainty but also cohort effect uncertainty and population basis risk. To enhance practicality, the hedging strategy is formulated in both static and dynamic settings, and its effectiveness can be evaluated in terms of either variance or 1-year ahead Value-at-Risk (the latter is highly relevant to solvency capital requirements). Three real data illustrations are constructed to demonstrate (1) the impact of population basis risk and cohort effect uncertainty on hedge effectiveness, (2) the benefit of dynamically adjusting a delta longevity hedge, and (3) the relationship between risk premium and hedge effectiveness.     

Hedging mortality/longevity risks of insurance portfolios for life insurer/annuity provider and financial intermediary

In this paper, we propose two risk hedge schemes in which a life insurer (an annuity provider) can transfer mortality (longevity) risk of a portfolio of life (annuity) exposures to a financial intermediary by paying the hedging premium of a mortality-linked security. The optimal units of the mortality-linked security which maximize hedge effectiveness for a life insurer (an annuity provider) can be derived as closed-form formulas under the risk hedge schemes. Numerical illustrations show that the risk hedge schemes can significantly hedge the downside risk of loss due to mortality (longevity) risk for the life insurer (annuity provider) under some stochastic mortality models. Besides, finding an optimal weight of a portfolio of life and annuity business, the financial intermediary can reduce the sensitivity to mortality rates but the model risk; a security loading may be imposed on the hedge premium for a higher probability of gain to compensate the financial intermediary for the inevitable model risk.     

A common age effect model for the mortality of multiple populations

We introduce a model for the mortality rates of multiple populations. To build the proposed model we investigate to what extent a common age effect can be found among the mortality experiences of several countries and use a common principal component analysis to estimate a common age effect in an age–period model for multiple populations. The fit of the proposed model is then compared to age–period models fitted to each country individually, and to the fit of the model proposed by Li and Lee (2005).Although we do not consider stochastic mortality projections in this paper, we argue that the proposed common age effect model can be extended to a stochastic mortality model for multiple populations, which allows to generate mortality scenarios simultaneously for all considered populations. This is particularly relevant when mortality derivatives are used to hedge the longevity risk in an annuity portfolio as this often means that the underlying population for the derivatives is not the same as the population in the annuity portfolio.     

Multi-population mortality models: A factor copula approach

Modeling mortality co-movements for multiple populations have significant implications for mortality/longevity risk management. A few two-population mortality models have been proposed to date. They are typically based on the assumption that the forecasted mortality experiences of two or more related populations converge in the long run. This assumption might be justified by the long-term mortality co-integration and thus be applicable to longevity risk modeling. However, it seems too strong to model the short-term mortality dependence. In this paper, we propose a two-stage procedure based on the time series analysis and a factor copula approach to model mortality dependence for multiple populations. In the first stage, we filter the mortality dynamics of each population using an ARMA–GARCH process with heavy-tailed innovations. In the second stage, we model the residual risk using a one-factor copula model that is widely applicable to high dimension data and very flexible in terms of model specification. We then illustrate how to use our mortality model and the maximum entropy approach for mortality risk pricing and hedging. Our model generates par spreads that are very close to the actual spreads of the Vita III mortality bond. We also propose a longevity trend bond and demonstrate how to use this bond to hedge residual longevity risk of an insurer with both annuity and life books of business.     

Semi-parametric extensions of the Cairns–Blake–Dowd model: A one-dimensional kernel smoothing approach

Over the last few decades, there has been an enormous growth in mortality modeling as the field of mortality risk and longevity risk has attracted great attention from academic, government and private sectors. In this paper, we propose a time-varying coefficient (TVC) mortality model aiming to combine the good characteristics of existing models with efficient model calibration methods. Nonparametric kernel smoothing techniques have been applied in the literature of mortality modeling and based on the findings from Li et al.’s (2015) study, such techniques can significantly improve the forecasting performance of mortality models. In this study we take the same path and adopt a kernel smoothing approach along the time dimension. Since we follow the model structure of the Cairns–Blake–Dowd (CBD) model, the TVC model we propose can be seen as a semi-parametric extension of the CBD model and it gives specific model design according to different countries’ mortality experience. Our empirical study presented here includes Great Britain, the United States, and Australia amongst other developed countries. Fitting and forecasting results from the empirical study have shown superior performances of the model over a selection of well-known mortality models in the current literature.     

A step-by-step guide to building two-population stochastic mortality models

Two-population stochastic mortality models play a crucial role in the securitization of longevity risk. In particular, they allow us to quantify the population basis risk when longevity hedges are built from broad-based mortality indexes. In this paper, we propose and illustrate a systematic process for constructing a two-population mortality model for a pair of populations. The process encompasses four steps, namely (1) determining the conditions for biological reasonableness, (2) identifying an appropriate base model specification, (3) choosing a suitable time-series process and correlation structure for projecting period and/or cohort effects into the future, and (4) model evaluation.For each of the seven single-population models from Cairns et al. (2009), we propose two-population generalizations. We derive criteria required to avoid long-term divergence problems and the likelihood functions for estimating the models. We also explain how the parameter estimates are found, and how the models are systematically simplified to optimize the fit based on the Bayes Information Criterion. Throughout the paper, the results and methodology are illustrated using real data from two pairs of populations.     

Modeling multi-country mortality dependence and its application in pricing survivor index swaps—A dynamic copula approach

This paper introduces mortality dependence in multi-country mortality modeling using a dynamic copula approach. Specifically, we use time-varying copula models to capture the mortality dependence structure across countries, examining both symmetric and asymmetric dependence structures. In addition, to capture the phenomenon of a heavy tail for the multi-country mortality index, we consider not only the setting of Gaussian innovations but also non-Gaussian innovations under the Lee–Carter framework model. As tests of the goodness of fit of different dynamic copula models, the pattern of mortality dependence, and the distribution of the innovations, we used empirical mortality data from Finland, France, the Netherlands, and Sweden. To understand the effect of mortality dependence on longevity derivatives, we also built a valuation framework for pricing a survivor index swap, then investigated the fair swap rates of a survivor swap numerically. We demonstrate that failing to consider the dynamic copula mortality model and non-Gaussian innovations would lead to serious underestimations of the swap rates and loss reserves.     

Risk-minimization for life insurance liabilities with basis risk

In this paper we study the hedging of typical life insurance payment processes in a general setting by means of the well-known risk-minimization approach. We find the optimal risk-minimizing strategy in a financial market where we allow for investments in a hedging instrument based on a longevity index, representing the systematic mortality risk. Thereby we take into account and model the basis risk that arises due to the fact that the insurance company cannot perfectly hedge its exposure by investing in a hedging instrument that is based on the longevity index, not on the insurance portfolio itself. We also provide a detailed example within the context of unit-linked life insurance products where the dependency between the index and the insurance portfolio is described by means of an affine mean-reverting diffusion process with stochastic drift.     

Abundant bursting patterns of a fractional-order Morris–Lecar neuron model

In this paper, a fractional-order Morris–Lecar (M–L) neuron model with fast-slow variables is firstly proposed. The fractional-order M–L model is a generalization of the integer-order M–L model with fast-slow variables, where the fractional-order derivative is used to characterize the memory effect and power law of membranes. Then the bursting patterns of the new model are investigated by using the bifurcation theory of fast-slow dynamical systems. Numerical simulation shows that the new model exhibits some bursting patterns that appear in some common neuron models with properly chosen parameters but do not exist in the corresponding integer-order M–L model. Further, on the basis of a comparison of the nonlinear dynamics between the fractional-order M–L model and the integer-order M–L model, we show that the fractional-order derivative can activate the slow potassium ion channel faster and play an important role to modulate the firing activity of the new model.     

Delta–Gamma hedging of mortality and interest rate risk

One of the major concerns of life insurers and pension funds is the increasing longevity of their beneficiaries. This paper studies the hedging problem of annuity cash flows when mortality and interest rates are stochastic. We first propose a Delta–Gamma hedging technique for mortality risk. The risk factor against which to hedge is the difference between the actual mortality intensity in the future and its “forecast” today, the forward intensity. We specialize the hedging technique first to the case in which mortality intensities are affine, then to Ornstein–Uhlenbeck and Feller processes, providing actuarial justifications for this selection. We show that, without imposing no arbitrage, we can get equivalent probability measures under which the HJM condition for no arbitrage is satisfied. Last, we extend our results to the presence of both interest rate and mortality risk. We provide a UK calibrated example of Delta–Gamma hedging of both mortality and interest rate risk.     

Reverse mortgage pricing and risk analysis allowing for idiosyncratic house price risk and longevity risk

Reverse mortgages provide an alternative source of funding for retirement income and health care costs. The two main risks that reverse mortgage providers face are house price risk and longevity risk. Recent real estate literature has shown that the idiosyncratic component of house price risk is large. We analyse the combined impact of house price risk and longevity risk on the pricing and risk profile of reverse mortgage loans in a stochastic multi-period model. The model incorporates a new hybrid hedonic–repeat-sales pricing model for houses with specific characteristics, as well as a stochastic mortality model for mortality improvements along the cohort direction (the Wills–Sherris model). Our results show that pricing based on an aggregate house price index does not accurately assess the risks underwritten by reverse mortgage lenders, and that failing to take into account cohort trends in mortality improvements substantially underestimates the longevity risk involved in reverse mortgage loans.     

Estimation of time-varying mortality rates using continuous models for Daphnia magna

Structured population models that make the assumption of constant demographic rates do not accurately describe the complex life histories seen in many species. We investigated the accuracy of using constant versus time-varying mortality rates within discrete and continuously structured models for Daphnia magna. We tested the accuracy of the models we considered using density-independent survival data for 90 daphnids. We found that a continuous differential equation model with a time-varying mortality rate was the most accurate model for describing our experimental D. magna survival data. Our results suggest that differential equation models with variable parameters are an accurate tool for estimating mortality rates in biological scenarios in which mortality might vary significantly with age.     

Identifiability issues of age–period and age–period–cohort models of the Lee–Carter type

The predominant way of modelling mortality rates is the Lee–Carter model and its many extensions. The Lee–Carter model and its many extensions use a latent process to forecast. These models are estimated using a two-step procedure that causes an inconsistent view on the latent variable. This paper considers identifiability issues of these models from a perspective that acknowledges the latent variable as a stochastic process from the beginning. We call this perspective the plug-in age–period or plug-in age–period–cohort model. Defining a parameter vector that includes the underlying parameters of this process rather than its realizations, we investigate whether the expected values and covariances of the plug-in Lee–Carter models are identifiable. It will be seen, for example, that even if in both steps of the estimation procedure we have identifiability in a certain sense it does not necessarily carry over to the plug-in models.     

Discrete May–Leonard Competition Models III

A discretization method attributed to Kahan is used to approximate the May–Leonard (M–L) competition model for three species. The local dynamics of this discrete-time M–L model are analyzed. This model differs from the discrete M–L models being studied previously. This discrete model shows dynamical consistency with the continuous M–L model. Numerically, we showed that the discrete M–L model has a degenerate Hopf-bifurcation, which is consistent with the continuous M–L model.     

A more meaningful parameterization of the Lee–Carter model

A new Lee–Carter model parameterization is introduced with two advantages. First, the Lee–Carter parameters are normalized such that they have a direct and intuitive interpretation, comparable across populations. Second, the model is stated in terms of the “needed-exposure” (NE). The NE is the number required in order to get one expected death and is closely related to the “needed-to-treat” measure used to communicate risks and benefits of medical treatments. In the new parameterization, time parameters are directly interpretable as an overall across-age NE. Age parameters are interpretable as age-specific elasticities: percentage changes in the NE at a particular age in response to a percent change in the overall NE. A similar approach can be used to confer interpretability on parameters of other mortality models.     

Heath–Jarrow–Morton modelling of longevity bonds and the risk minimization of life insurance portfolios

This paper has two parts. In the first, we apply the Heath–Jarrow–Morton (HJM) methodology to the modelling of longevity bond prices. The idea of using the HJM methodology is not new. We can cite Cairns et al. [Cairns A.J., Blake D., Dowd K, 2006. Pricing death: framework for the valuation and the securitization of mortality risk. Astin Bull., 36 (1), 79–120], Miltersen and Persson [Miltersen K.R., Persson S.A., 2005. Is mortality dead? Stochastic force of mortality determined by arbitrage? Working Paper, University of Bergen] and Bauer [Bauer D., 2006. An arbitrage-free family of longevity bonds. Working Paper, Ulm University]. Unfortunately, none of these papers properly defines the prices of the longevity bonds they are supposed to be studying. Accordingly, the main contribution of this section is to describe a coherent theoretical setting in which we can properly define these longevity bond prices. A second objective of this section is to describe a more realistic longevity bonds market model than in previous papers. In particular, we introduce an additional effect of the actual mortality on the longevity bond prices, that does not appear in the literature. We also study multiple term structures of longevity bonds instead of the usual single term structure. In this framework, we derive a no-arbitrage condition for the longevity bond financial market. We also discuss the links between such HJM based models and the intensity models for longevity bonds such as those of Dahl [Dahl M., 2004. Stochastic mortality in life insurance: Market reserves and mortality-linked insurance contracts, Insurance: Math. Econom. 35 (1) 113–136], Biffis [Biffis E., 2005. Affine processes for dynamic mortality and actuarial valuations. Insurance: Math. Econom. 37, 443–468], Luciano and Vigna [Luciano E. and Vigna E., 2005. Non mean reverting affine processes for stochastic mortality. ICER working paper], Schrager [Schrager D.F., 2006. Affine stochastic mortality. Insurance: Math. Econom. 38, 81–97] and Hainaut and Devolder [Hainaut D., Devolder P., 2007. Mortality modelling with Lévy processes. Insurance: Math. Econom. (in press)], and suggest the standard pricing formula of these intensity models could be extended to more general settings.In the second part of this paper, we study the asset allocation problem of pure endowment and annuity portfolios. In order to solve this problem, we study the “risk-minimizing” strategies of such portfolios, when some but not all longevity bonds are available for trading. In this way, we introduce different basis risks.     

Statistical emulators for pricing and hedging longevity risk products

We propose the use of statistical emulators for the purpose of analyzing mortality-linked contracts in stochastic mortality models. Such models typically require (nested) evaluation of expected values of nonlinear functionals of multi-dimensional stochastic processes. Except in the simplest cases, no closed-form expressions are available, necessitating numerical approximation. To complement various analytic approximations, we advocate the use of modern statistical tools from machine learning to generate a flexible, non-parametric surrogate for the true mappings. This method allows performance guarantees regarding approximation accuracy and removes the need for nested simulation. We illustrate our approach with case studies involving (i) a Lee–Carter model with mortality shocks; (ii) index-based static hedging with longevity basis risk; (iii) a Cairns–Blake–Dowd stochastic survival probability model; (iv) variable annuities under stochastic interest rate and mortality.     

Stochastic portfolio specific mortality and the quantification of mortality basis risk

In the last decade a vast literature on stochastic mortality models has been developed. However, these models are often not directly applicable to insurance portfolios because:

(a) For insurers and pension funds it is more relevant to model mortality rates measured in insured amounts instead of measured in the number of policies.

(b) Often there is not enough insurance portfolio specific mortality data available to fit such stochastic mortality models reliably.

Therefore, in this paper a stochastic model is proposed for portfolio specific mortality experience. Combining this stochastic process with a stochastic country population mortality process leads to stochastic portfolio specific mortality rates, measured in insured amounts. The proposed stochastic process is applied to two insurance portfolios, and the impact on the Value at Risk for longevity risk is quantified. Furthermore, the model can be used to quantify the basis risk that remains when hedging portfolio specific mortality risk with instruments of which the payoff depends on population mortality rates.     

Longevity risk management for life and variable annuities: The effectiveness of static hedging using longevity bonds and derivatives

For many years, the longevity risk of individuals has been underestimated, as survival probabilities have improved across the developed world. The uncertainty and volatility of future longevity has posed significant risk issues for both individuals and product providers of annuities and pensions. This paper investigates the effectiveness of static hedging strategies for longevity risk management using longevity bonds and derivatives (q-forwards) for the retail products: life annuity, deferred life annuity, indexed life annuity, and variable annuity with guaranteed lifetime benefits. Improved market and mortality models are developed for the underlying risks in annuities. The market model is a regime-switching vector error correction model for GDP, inflation, interest rates, and share prices. The mortality model is a discrete-time logit model for mortality rates with age dependence. Models were estimated using Australian data. The basis risk between annuitant portfolios and population mortality was based on UK experience. Results show that static hedging using q-forwards or longevity bonds reduces the longevity risk substantially for life annuities, but significantly less for deferred annuities. For inflation-indexed annuities, static hedging of longevity is less effective because of the inflation risk. Variable annuities provide limited longevity protection compared to life annuities and indexed annuities, and as a result longevity risk hedging adds little value for these products.