Strong crustal earthquakes in Greece are typically followed by aftershocks, the properties of which are important factors in seismic hazard assessment. In order to examine the properties of earthquake sequences, we prepared an earthquake catalog comprising aftershock sequences with mainshocks of Mw ≥ 5.5 from 1995 to 2021. Regional aftershock parameters were estimated to highlight variations in aftershock decay and productivity among regions with similar seismotectonic characteristics. A statistically based method of estimating aftershock duration and a metric of relative aftershock productivity to examine the variations among the different cases were employed. From the detailed analysis of the selected seismic sequences, we attempt to unravel the physical mechanisms behind deviations in aftershock duration and productivity and resolve the relative contribution of background seismicity, the Omori–Utsu law parameters and the mainshock faulting properties. From our analysis, the duration of aftershock sequences depends upon the rupture process of the mainshock, independently of its magnitude. The same applies to aftershock productivity, however, other tectonic setting (e.g., seismic coupling) or source-related (e.g., focal depth, stress drop) parameters also contribute. The estimated regional parameters of the aftershock rate models could be utilized as initial ones to forecast the aftershock occurrence rates at the early stage following a mainshock.