The earthquake magnitude probability distribution is one of the underlying input data for certain earthquake analyses, such as probabilistic seismic hazard analysis. Nowadays, the method proposed by McGuire and Arabasz (1990) is commonly used for obtaining the (simulated) earthquake magnitude probability distributions. However, based on the observed earthquake data in 5 regions (Taiwan, Japan, California, Turkey, and Greece), the model did not fit the observation well. Instead, all of the case studies show that using the newly proposed gamma distribution can improve the simulation significantly compared to the conventional method.