In this paper, we introduce a family of bivariate generalized Gompertz-power series distributions. This new class of bivariate distributions contains several models such as: bivariate generalized Gompertz -geometric, -Poisson, - binomial, -logarithmic, -negative binomial and bivariate generalized exponental-power series distributions as special cases. We express the method of construction and derive different properties of the proposed class of distributions. The method of maximum likelihood and EM algorithm are used for estimating the model parameters. Finally, we illustrate the usefulness of the new distributions by means of application to real data sets.