A stacked sequence of Pliocene Gilbert-type fan deltas in the Loreto basin was shed from the footwall of the dextral-normal Loreto fault and deposited at the margin of a marine basin during rapid fault-controlled subsidence. Fan-delta parasequences coarsen upward from marine siltstone and sandstone at the base, through sandy bottomsets and gravelly foresets, to gravelly nonmarine topsets. Each topset unit is capped by a thin shell bed that records marine flooding of the delta plain. Several mechanisms may have produced repetitive vertical stacking of Gilbert deltas: (1) autocyclic delta-lobe switching; (2) eustatic sea-level fluctuations; (3) climatically controlled fluctuations in sediment input; and (4) episodic subsidence produced by temporal clustering of earthquakes. We favor hypothesis 4 for several reasons, but hypotheses 2 and 3 cannot be rejected at this time. Earthquake clustering can readily produce episodic subsidence at spatial and temporal scales consistent with stratigraphic trends observed in the Loreto basin. This model is supported by comparison with paleoseismological studies that document clustering on active faults over a wide range of time scales. Earthquake clustering is a new concept in basin analysis that may be helpful for understanding repetitive stratigraphy in tectonically active basins.