Mt Etna volcanic activity is broadly divided into flank eruptions and summit paroxysms. Here, building on previously-available literature and data on the start time of these events, we collate two separate catalogs of the two activity types. Then we separately model their temporal occurrence. The catalog of flank eruptions, spanning the last 400 years, has been modelled by means of the most widely used renewal models, among which the best one (through Akaike Information Criterion) is the Brownian Passage Time. The catalog of summit paroxysms, covering the period 1986–2022, according to our cluster analysis is best characterized by 12 clusters of paroxysms. We separately analyze the inter-event times between onset times of successive clusters of paroxysms (inter-cluster inter-event times) and the inter-event times between successive paroxysms within clusters (intra-cluster inter-event times). Again, the Brownian Passage Time is the best-fitting model, obviously with very different parameters in the two cases. We test the best-fitting models by checking their ability to reproduce features of the real catalogs. Finally, we provide an example of how to use in practice such temporal models in the context of probabilistic hazard assessment, showing a possible use in the case of tephra fallout hazard from summit paroxysms.