Bivariate copula trees for gross loss aggregation with positively dependent risks

We propose several numerical algorithms to compute the distribution of gross loss in a positively dependent catastrophe insurance portfolio. Hierarchical risk aggregation is performed using bivariate copula trees. Six common parametric copula families are studied. At every branching node, the distribution of a sum of risks is obtained by discrete copula convolution. This approach is compared to approximation by a weighted average of independent and comonotonic distributions. The weight is a measure of positive dependence through variance of the aggregate risk. During gross loss accumulation, the marginals are distorted by application of insurance financial terms, and the value of the mixing weight is impacted. To accelerate computations, we capture this effect using the ratio of standard deviations of pre-term and post-term risks, followed by covariance scaling. We test the performance of our algorithms using three examples of complex insurance portfolios subject to hurricane and earthquake catastrophes.