This study aims to model the joint drought duration and severity distribution using two-dimensional copulas. The method of inference function for margins (IFM method) is employed to construct copulas. Two separate maximum likelihood estimations of univariate marginal distributions are performed first, then followed by a maximization of the bivariate likelihood as a function of the dependence parameters. The drought duration and severity are assumed to be exponential and gamma distributions, respectively. Several copulas are tested to determine the best data fitted copula. Droughts, defined using the Standardized Precipitation Index (SPI), of Wushantou (Taiwan) are employed as an example to illustrate the proposed methodology. The copula fitting results for drought duration and severity are quite satisfactory. The bivariate drought analyses, including the joint probabilities and bivariate return periods, based on the derived copula-based joint distribution are also investigated to demonstrate the advantages of bivariate modeling of droughts.