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...