Spatial and temporal patterns of large fire (>100 ha) incidence in Portugal over the period 1984–2004 were modeled using extreme value statistics, namely the Peaks Over Threshold approach, which uses the Generalized Pareto Distribution (GPD) as a model. The original dataset includes all fires larger than 5 ha (30 616 fires) that were observed in Portugal during the study period, mapped from Landsat satellite imagery. The country was divided into eight regions, considered internally homogeneous from the perspective of their fire regimes and respective environmental correlates. The temporal analysis showed that there does not appear to be any trend in the incidence of very large fires, but revealed a cyclical behavior in the values of the GPD shape parameter, with a period in the range of 3 to 5 years. Spatial analysis highlighted strong regional differences in the incidence of large fires, and allowed the calculation of return levels for a range of fire sizes. This analysis was affected by the presence of a few outlying observations, which may correspond to clusters of contiguous fire scars, resulting in artificially large burned areas. We discuss some of the implications of our findings in terms of consequences for fire management aimed at preventing the occurrence of extremely large fires, and present ideas for extending the present study.