The reliability of a vertical breakwater is calculated using direct integration methods based on joint density functions.The horizontal and uplifting wave forces on the vertical breakwater can be well fitted by the lognormal and the Gumbel distributions,respectively.The joint distribution of the horizontal and uplifting wave forces is analyzed using different probabilistic distributions,including the bivariate logistic Gumbel distribution,the bivariate lognormal distribution,and three bivariate Archimedean copulas functions constructed with different marginal distributions simultaneously.We use the fully nested copulas to construct multivariate distributions taking into account related variables.Different goodness fitting tests are carried out to determine the best bivariate copula model for wave forces on a vertical breakwater.We show that a bivariate model constructed by Frank copula gives the best reliability analysis,using marginal distributions of Gumbel and lognormal to account for uplifting pressure and horizontal wave force on a vertical breakwater,respectively.The results show that failure probability of the vertical breakwater calculated by multivariate density function is comparable to those by the Joint Committee on Structural Safety methods.As copulas are suitable for constructing a bivariate or multivariate joint distribution,they have great potential in reliability analysis for other coastal structures. The reliability of a vertical breakwater is calculated using direct integration methods based on joint density functions. Horizontal and uplifting wave forces on the vertical breakwater can be well fitted by the lognormal and the Gumbel distributions, respectively. The joint distribution of the horizontal and uplifting wave forces is analyzed using different probabilistic distributions, including the bivariate logistic Gumbel distribution, the bivariate lognormal distribution, and three bivariate archimedean copulas functions constructed with different marginal distributions simultaneously. We use the fully nested copulas to construct multivariate distributions taking into account related variables. Different goodness fitting tests are carried out to determine the best bivariate copula model for wave forces on a vertical breakwater. We show that a bivariate model constructed by Frank copula gives the best reliability analysis, using marginal distributions of Gumbel and lognormal to account f or uplifting pressure and horizontal wave force on a vertical breakwater, respectively. These results show that failure probability of the vertical breakwater calculated by multivariate density function is comparable to those by the Joint Committee on Structural Safety methods. As copulas are suitable for constructing a bivariate or multivariate joint distribution, they have great potential in reliability analysis for other coastal structures.