## 2 Introduction  
Paleoseismology studies with a geomorphologic viewpoint on tectonically active areas are essential for modern seismic hazard assessment and determining the seismic potential of seismogenic faults over much more extended periods than instrumental measurements [1]; moreover, these studies provide tightly constrained data for the location, size, and recurrence of large earthquakes associated with surface rupture in recent geological time, which are critical for the assessment of seismic potential and hazards along active faults [2; 3; 4; 5]. Strike-slip faults have played a critical role in the development of paleoseismology for several reasons. First, strike-slip faults are often the longest faults on continental landmasses and typically have conspicuous geographic expressions [6; 7]. Second, many of these faults have long records of seismicity because they pass through populated continental regions and have experienced surface ruptures during large and great historical earthquakes. Third, becauseoseismic deformation along strike-slip faults is horizontal, subsequent earthquakes do not deeply bury or expose erosion or traces of earlier events [5]. Using earthquake geology, paleoseismology, and archaeology in regions with long historical records, we can determine earthquake construction, potentially slippery conditions, and slip allocation over multiple earthquake cycles [8]. Paleoseismic research aims to ascertain past earthquake activity by locating fault movements, rupture sites, and slip/ground shaking strength distributions. One of the most fundamental principles of paleoseismology is reconstructing the sequence of past earthquakes in relatively young geological units [9]. One of the primary goals of earthquake geology and seismology is to determine the earthquake recurrence interval, which permits probabilistic predictions of the time and magnitude of future earthquakes along a fault [10; 11; 12; 13]. Iran, located in the Alpine-Himalayan belt, is one of the most seismically active areas on the planet. Its seismic activity is distributed along the same active fault lines (Figure 1), with individual faults rupturing every few hundred to thousands of years during strong earthquakes with a magnitude higher than 7.4 [14; 15; 16; 17; 18; 19; 20; 21; 22]. Iran experienced 145 earthquakes with a magnitude higher than 5.5 between 1900 and 2022, which resulted in 190,300 fatalities. More than 1000 people were killed in 17 of these earthquakes, 11 of which had a magnitude of 7 or higher [23]. The magnitude and time of the last seismic events are crucial to figuring out the earthquake's recurrence. Understanding a fault's past behavior can help predict how it might behave in the future, which can help with loss mitigation and seismic hazard estimation. Figure 1: Map of the main active faults of Iran [14; 15; 16; 17; 18; 19; 20; 24]: **(a)** The position of major Iran’s structural zones among the Arabian and Eurasian plates. AZ = Azerbaijan; KP = Kopeh Dagh; TP = Turan plate;  
S = Sistan suture zone; M = Makran; arrows are GPS vectors. (**b**) Distribution of active faults throughout Iran. Zagros Main Recent Fault and Zagros Reverse Fault separate the Zagros mountain range from central Iran. The blue frame shows the location of the study area. Doroud Fault, as one of the main segments of the Main Recent Fault, is located in the study area. KB = Karehbas Fault; S = Sarvestan Fault; SP = SabzPushan Fault. The Zagros mountain range (Figure 1a), one of the main structural zones of SW Iran, which includes a large number of Iran's devastating earthquakes, is a fold-and-thrust belt at the leading edge of the Arabia-Eurasia continental collision zone, which accounts for about one-third of the collision zone's overall 20 to 30 mm/year convergence rate [25; 26].