We study an iterated temporal and contemporaneous aggregation of N independent copies of a strongly stationary subcritical Galton-Watson branching process with regularly varying immigration having index α ∈ (0, 2). Limits of finite dimensional distributions of appropriately centered and scaled aggregated partial sum processes are shown to exist when first taking the limit as N → ∞ and then the time scale n → ∞. The limit process is an α-stable process if α ∈ (0, 1) ∪ (1, 2), and a deterministic line with slope 1 if α = 1.Recently, Puplinskaitė and Surgailis [21,22] studied iterated aggregation of random coefficient autoregressive processes of order 1 with common innovations and with so-called idiosyncratic innovations, respectively, belonging to the domain of attraction of an α-stable law. Limits 2010 Mathematics Subject Classifications 60J80, 60F05, 60G10, 60G52, 60G70.