We construct 3D models to simulate South American subduction since 100 Ma.
The effect of seafloor age results in slab dip angle distribution opposite to observed.
Dynamic suction from the overriding plate has a long-wavelength and prolonged effect on reducing slab dip angle.
Subducting oceanic plateau and aseismic ridges effectively flatten the slab at local scales.
Formation of South American flat slabs represents a multi-scale and multi-physics process.