Abstract Following the method of Kijko and Sellevoll ∗ ∗ (1987). the maximum likelihood estimation of seismic hazard parameters (maximum regional intensity, mean activity rate and b-value) is obtained for several seismic source regions in Spain, combining data from large historical events together with complete instrumental catalogues. Uncertainties in seismic hazard parameters estimates are also given. The mean activity rate ( I ≥ 4) ranges from about two events/yr in the Western Pyrenees, the Gulf of Cadiz and the Atlantic Ocean region southwest of Cap San Vicente to values around 0.5 in the Catalan Coastal Ranges and the Guadalquivir basin. The b-value of the frequency intensity distribution ranges between 0.4 and 0.6 in all the seismic source zones considered. The maximum regional intensity estimate has its largest value in the Pyrenees, the Betic Cordillera and the At antic Ocean region southwest of Cap San Vicente. The lowest value corresponds to the western part of the Alboran Sea.