Intensity-duration-frequency (IDF) curves are widely used for peak discharge estimation in designing hydraulic structures. The traditional Gumbel probability method entails selecting annual maximum rainfall depths (intensities) conditioned on a fixed time window width (which in general will not coincide with the rainfall event duration) from a continuous record to perform a frequency analysis in terms of the marginal distribution. The digitized database contains annual maximum intensities for selected discrete durations. This method presents problems when intensities are required for arbitrary durations which are not part of the selected durations. Accurate interpolated and especially extrapolated intensity values are hard to obtain. The present study offers two methods both involving a joint probability approach to overcome the deficiencies inherent in the traditional method of IDF analysis. The first joint probability approach employs Box-Cox and modulus transformations to transform original data to near bivariate normality. The second method does not require such a transformation. Instead, it uses the closed-form bivariate Burr III cumulative distribution to fit the data. Another advantage of the joint probability approach is that it allows one to gauge the rarity of certain extreme events, such as probable maximum precipitation, in terms of the joint occurrence of its extremely high intensity and a sufficiently long duration (e.g. 24 hours). The joint probability approach is applied to three data sets. The resulting conditional probability intensity estimates are quite close to those obtained by traditional Gumbel IDF analysis. In addition, reliable interpolated and extrapolated intensities are available because the approach essentially fits a flexible surface to the discrete data with the capability of providing a complete probabilistic structure. / Master of Science