Four new bivariate distributions with gamma type marginals are introduced. Various representations are derived for their joint densities, product, moments, conditional densities and conditional moments. Some of these representations involve special functions such as the complementary incomplete gamma and Whittaker functions. Construction of multivariate generalizations is discussed. Finally, all application to rainfall data from Florida is provided.