ABSTRACT Seismic hazard expressed in terms of annual rate of exceedance and is used to calculate the earthquake insurance premium. Annual rate of exceedance is a complicated function of magnitudes, distances from site to earthquake sources and attenuation. Due to its complexity, determination of exact sampling distribution of earthquake insurance premium is not an easy task. This research proposes Monte Carlo simulation approach to determine the sampling distribution of earthquake insurance premium. Annual rate of exceedance was simulated first and then the insurance premium calculated based on simulation of annual rate of exceedance. The simulation involves quantifying synthetic catalogue similar to historical catalogue. Its simulation is conducted in order to construct annual rate of exceedance as an indicator of earthquake risk used in earthquake insurance. The simulation of 25 iteration and sample size of 100 shows that the sampling distribution of insurance premium is skewed to the right. The idea of using Monte Carlo simulation to study the sampling distribution of annual rate of exceedance is the originality of the study and the study contributes to the methodology of earthquake insurance.