Most probabilistic seismic-hazard analysis procedures require that at least three seismic source parameters be known, namely the mean seismic activity rate λ, the Gutenberg-Richter b-value, and the area-characteristic (seismogenic source) maximum possible earthquake magnitude m max . In almost all currently used seismic-hazard assessment procedures that utilize these three parameters, it is explicitly assumed that all three remain constant over time and space. However, closer examination of most earthquake catalogs has indicated that significant spatial and temporal variations existed in the seismic activity rate λ, as well as in the Gutenberg-Richter b-value. In this study, the maximum likelihood estimation of these earthquake hazard parameters considers the incompleteness of the catalogs, the uncertainty in the earthquake magnitude determination, as well as the uncertainty associated with the applied earthquake-occurrence models. The uncertainty in the earthquake-occurrence models is introduced by assuming that both the mean seismic activity rate λ and the Gutenberg-Richter b-value are random variables, each described by the gamma distribution. This approach results in the extension of the classic frequency-magnitude Gutenberg-Richter relation and the