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 a pairwise dependency. Recently vine copula has been recognized as the highly flexible approach to constructing a higher dimensional joint density framework. In such construction, parametric class copula with parametric univariate marginal distributions is often involved. Such incorporation can lack flexibility due to parametric functions with 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 is applied to 46 years of flood characteristics collected on the west coast of Canada. The univariate flood marginal distribution is modelled using nonparametric kernel density estimation (KDE). The 2-D Bernstein estimator and Beta kernel copulas estimator are 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 is selected using goodness-of-fit (GOF) test statistics. The performance of the nonparametric vine approach is also compared with the vine constructed in the parametric and semiparametric fitting procedure. Investigation reveals that the D-vine constructed using Bernstein copula with normal KDE marginals nonparametrically performed well in capturing dependence of the compound events. Finally, the derived nonparametric model is used in the estimation of trivariate OR- and AND-joint return periods, further employed in estimating failure probability (FP) statistics. The trivariate return periods for the AND-joint case are higher than for the OR-joint case for the same flood combination. Also, the trivariate flood hazard results in a high-value FP than bivariate and univariate events. Ignoring the trivariate dependence could result in the underestimation of FP