This study introduces an improved hypocentral version of the space‐time Epidemic‐Type Aftershock Sequence (ETAS) model, entitled as the 3D‐finite source (3D‐FS) ETAS model, and applies it to the analysis of the Southern California earthquake catalog. By stochastic reconstruction, we are able to reconstruct the patterns of aftershock productivity density along the mainshock ruptures. Detailed analysis of the productivity patterns reveals that: (1) Directly triggered aftershocks make up 21% to 41% of all earthquakes within the mainshock rupture areas, and show significant spatial heterogeneity and temporal migrations; (2) Major aftershocks tend to locate in low productivity areas, at the edges of clusters formed by small aftershocks; (3) Large slip areas are depleted of aftershocks, over 60% of all productivity distributes in areas with slip less than 0.3 times of the maximum slip, and the trajectory of the productivity pattern on the fault plane demonstrates apparent compensation to coseismic slip. We relate the difference in triggering abilities of four mainshocks to the heat flow in corresponding regions. Simulation results suggest that the 3D‐FS ETAS model has apparent advantages of improving the performance of short‐term aftershock forecast. Moreover, the later aftershocks are more correlated with the locations of subsequent events than earlier aftershocks, suggesting that the migration of aftershocks is important for mitigating aftershock hazard.