In seismology, an important goal is to attain a better understanding of the earthquake process. In this study of the physics of aftershock generation, I couple statistical analysis with modelling of physical processes in the postseismic period. I present a theoretical formulation for the distribution of interevent times for aftershock sequences obeying the empirically well established Omori law. As opposed to claims by other authors, this work demonstrates that the duration of the time interval between two successive earthquakes cannot be used to identify whether or not they belong to the same aftershock sequence or occur as a result of the same underlying process. This implies that a proper understanding of earthquake interevent time distributions is necessary before conclusions regarding the physics of the earthquake process are drawn. In a discussion of self-organised criticality (SOC) in relation to empirical laws in seismology, I find that Omori's law for aftershocks cannot be used as evidence for the theory of SOC. Instead, I consider that the occurrence of aftershocks in accordance with Omori's law is a result of a physical process that can be modelled and understood. I analyse characteristic features in the spatiotemporal distribution of aftershocks in the south Iceland seismic zone, following the two M6.5 June 2000 earthquakes and a M4.5 earthquake in September, 1999. These features include an initially constant aftershock rate, whose duration is larger following a larger main shock, and a subsequent power law decay that is interrupted by distinct and temporary deviations in terms of rate increases and decreases. Based on pore pressure diffusion modelling, I interpret these features in terms of main shock initiated diffusion processes. I conclude that thorough data analysis and physics-based modelling are essential components in attempts to improve our understanding of processes governing the occurrence of earthquakes.