Previous studies on red blood cell (RBC) aggregation have elucidated the inverse relationship between shear rate and RBC aggregation under Poiseuille flow. However, the local parabolic rouleaux pattern in the arterial flow observed in ultrasonic imaging cannot be explained by shear rate alone. A quantitative approach is required to analyze the spatiotemporal variation in arterial pulsatile flow and the resulting RBC aggregation. In this work, a 2D RBC model was used to simulate RBC motion driven by interactional and hydrodynamic forces based on the depletion theory of the RBC mechanism. We focused on the interaction between the spatial distribution of shear rate and the dynamic motion of RBC aggregation under sinusoidal pulsatile flow. We introduced two components of shear rate, namely, the radial and axial shear rates, to understand the effect of sinusoidal pulsatile flow on RBC aggregation. The simulation results demonstrated that specific ranges of the axial shear rate and its ratio with radial shear rate strongly affected local RBC aggregation and parabolic rouleaux formation. These findings are important, as they indicate that the spatiotemporal variation in shear rate has a crucial role in the aggregate formation and local parabolic rouleaux under pulsatile flow.