Trivariate Joint Distribution Modelling of Compound Events Using the Nonparametric D-Vine Copula Developed Based on a Bernstein and Beta Kernel Copula Density Framework

Low-lying coastal communities are often threatened by compound flooding (CF), which can be determined through the joint occurrence of storm surges, rainfall and river discharge, either successively or in close succession. The trivariate distribution can demonstrate the risk of the compound phenomenon more realistically, rather than considering each contributing factor independently or in pairwise dependency relations. Recently, the vine copula has been recognized as a highly flexible approach to constructing a higher-dimensional joint density framework. In these, the parametric class copula with parametric univariate marginals is often involved. Its incorporation can lead to a lack of flexibility due to parametric functions that have prior distribution assumptions about their univariate marginal and/or copula joint density. This study introduces the vine copula approach in a nonparametric setting by introducing Bernstein and Beta kernel copula density in establishing trivariate flood dependence. The proposed model was applied to 46 years of flood characteristics collected on the west coast of Canada. The univariate flood marginal distribution was modelled using nonparametric kernel density estimation (KDE). The 2D Bernstein estimator and beta kernel copula estimator were tested independently in capturing pairwise dependencies to establish D-vine structure in a stage-wise nesting approach in three alternative ways, each by permutating the location of the conditioning variable. The best-fitted vine structure was selected using goodness-of-fit (GOF) test statistics. The performance of the nonparametric vine approach was also compared with those of vines constructed with a parametric and semiparametric fitting procedure. Investigation revealed that the D-vine copula constructed using a Bernstein copula with normal KDE marginals performed well nonparametrically in capturing the dependence of the compound events. Finally, the derived nonparametric model was used in the estimation of trivariate joint return periods, and further employed in estimating failure probability statistics.