Extensive vertebrate tracksite from the Upper Red Formation (middle–late Miocene), west Zanjan, northwestern Iran

An extensive vertebrate tracksite from the middle– late Miocene Upper Red Formation in western Zanjan Province, northwestern Iran, provides new records of paleobiogeographical significance. These are records of common footprints of felids referred to as Felipeda lynxi Panin and Avram, canids referred to as Canipeda longigriffa Panin and Avram and less common bird footprints referred to as Iranipeda abeli Lambrecht. The Canipeda record establishes the late Miocene presence of canids on the Iranian Plateau as part of the Eurasia-wide “Eucyon event”. The Felipeda footprints are consistent with body fossil and footprint records elsewhere in Eurasia that indicate a widespread distribution of felids by late Miocene time. An unusual trace associated with the footprints is a large, shallow grazing or locomotion trace similar to Megaplanolites in some features but distinctive in various features and its occurrence in nonmarine facies, and it will be the subject of further study.


Introduction
The Upper Red Formation (URF) is known as the main and most extensive Miocene rock unit in Central Iran. The Lower Red Formation (LRF) and URF are the two main continental rock units, composed of marl, sandstone and gypsum red beds, that form the lower and upper boundaries, respectively, of the marine sediments of the Oligo-Miocene Qom Formation. All of these formations were deposited during the Oligocene to late Miocene after the late Eocene Pyrenean orogeny in the Central Iran basin. Both the LRF and URF are poor in index fossils, and their ages have been estimated by their lithostratigraphic position around the Qom Formation (Rupelian to Tortonian), so the LRF is considered early Oligocene in age, and the URF is assigned to the middle to late Miocene (Rahimzadeh, 1994;Gansser, 1955;Aghanabati, 2004). Continental facies of the URF are a good candidate for vertebrate tracksite studies (Amini, 2001;Rafiei et al., 2011;Amini, 1997), and numerous vertebrate footprints have been reported from this formation (Abbassi, 2010;Abbassi and Amini, 2008;Abbassi and Shakeri, 2005). Vertebrate footprints of the URF, however, are very diverse and abundant in its outcrops in northwestern Iran, and new tracksites have also been discovered by one of us (Nasrollah Abbassi) during the last decade in Zanjan Province. The main aim of the present research is to report and introduce one of these tracksites, which includes an extensive assemblage of abundant and well-preserved mammal and bird footprints on a large cliff in the Chehrabad area in western Zanjan Province (Fig. 1). The location of this large slab (12 × 10 m) is in the north Hamzalu village; GPS coordinates are 36 • 56.045 latitude and 47 • 50.262 longitude. Kavir area in north-central Iran show its maximum thickness (> 6000 m) with numerous evaporite diapiric structures. The thickness and rock facies of the URF decrease towards the northwest (1000-2000 m). Geologists of the NIOC divided the URF into two rock units in the Qom area as follows: M1 -lower evaporite subunit (300-400 m), composed of dark brown to red salt, anhydrite, pebbly terrigenous layers and thin shale beds, and M2 -upper subunit (3000-5000 m), which includes gypsiferous marls, brown to red sandstones, and alternations of light red marls with sandstones.
Additional data show that M1 is the uppermost member of the Qom Formation, and M2 is classifiable into three subunits in the Qom area as follows: m1-dark red gypsiferous sandstone, shale and siltstone with the green key bed in the lower part (2000 m thick), m2 -mainly eroded sandstone (about 1000 m thick), and m3 -pale yellow, strongly gypsiferous siltstone and marl with intercalations of soft calcareous sandstone (200-500 m thick) (Mostofi and Frei, 1959;Abaie et al., 1964).
The URF has very extensive outcrops in the southern and western parts of Zanjan Province. Stöcklin and Eftekhar-Nezhad (1969) divided the URF into two subunits in the southern parts of the province. The lower subunit (M1, 600-700 m thick) is characterized by its alternating beds that are characterized by colorfully banded white-green, white-pink, purple, violet and brown strata and by a relatively high content of evaporites. Colored marls with thin siltstone, sandstone, salt and gypsum beds and veins are the dominant rocks. There are gypsum beds, however, in a few diapiric structures in the lower part of the subunit. An andesitic crystal tuff (4 m thick) is also present in the lower subunit M1. The uniform, red-colored upper subunit (M2, 1500-2000 m thick) consists of clayey, marly and sandy shales with sandstone intercalations. Local to reginal conglomerate layers and lenses occur in this subunit. The gypsum layers are relatively less than in M1.
Outcrops of M1 in the western part of Zanjan Province around the Qezel Ozan River (1000-2000 m thick) include light red-to brown-colored marls with a distinctive gypsum layer at the lower boundary with the Qom. Vertebrate footprints have been found in two horizons in this subunit in the study area in the upper most layers of M1, which are green marls. M2 consists of dark brown-green marls with sandstone layers. Only one vertebrate-footprint-bearing horizon has been identified in M2 in western Zanjan (Fig. 2).
The tracksite studied here is north of Hamzalu village, in the Chehrabad area, western Zanjan, in green and brown marls with sandstone and siltstone intercalations, and it belongs to the second footprint-bearing layer of M1 (Fig. 2). Vertebrate footprints were discovered in detrital layers preserved in convex hyporelief with non-biogenic structures that include mud cracks, flutes, grooves and prod casts.
No age diagnostic fossils have been found in the URF in the Zanjan area. The URF is generally dated as middle to late Miocene in age because it overlies the Qom Formation and is overlain by Pliocene conglomerate. The uppermost layer of the Qom Formation is dated as Burdigalian, ca. 17 Ma (Schuster and Wielandt, 1999;Daneshian and Ramezani Dana, 2007). Magnetostratigraphic data show that the URF adjacent to the southern Alborz Mountains ranges in age between 17.5 and 7.5 Ma (Ballato et al., 2008). To refine the chronology of the URF in the Qezel Ozan area in the west of Zanjan, Ballato et al. (2017) determined an 40 Ar/ 39 Ar age of 15.0 ± 0.3 to 25.0 ± 1.3 Ma and a zircon U-Pb age of 10.7 ± 0.2 to 19.8 ± 0.9 Ma.
The uppermost reefal limestone of the Qom Formation (Member f) is overlain by evaporite and terrigenous layers, and it is early Miocene (Burdigalian) in age (Daneshian and Ramezani Dana, 2007;Mohammadi, 2020;Okhravi, 1998;Reuter et al., 2009). Thus, gypsum beds in the study area are assignable to the Qom Formation; the URF begins with red marls in the study area with two lithostratigraphic units, M1 and M2. M1 was deposited on the f member of the Qom Formation, immediately above its evaporites. Thus, the lower parts of the URF basin began in the Burdigalian, and the zircon U-Pb age of 10.7 ± 0.2 Ma is the youngest radiometric age date for the URF in western Zanjan (Ballato et al., 2017). Therefore, the continental basin of the URF on the southern flank of the Alborz Mountains in the Zanjan area was present during the late Burdigalian to early Tortonian. During this time interval (about 16-10 Ma), the prevailing climate of the URF basin has been estimated as warm and dry and affected by the surrounding topography such as the evolution of the Alborz Mountains (Ballato et al., 2010;Cronin, 2009).
The URF was deposited in continental environments by braided to meandering rivers, playa lakes and ephemeral ponds. Paleocurrent indicators show that sediments were transported from several points into the URF basin, particularly from the north and south (Amini, 1997). The fluvial depositional systems, including distal braided river, meandering river, and proximal braided river, were distinguished in the Avaj area, Hamadan Province, northwest Iran (Rafiei et al., 2011). The basin margins are delineated by a combination of active, strike-slip and normal fault systems (Alavi, 2004;Morley et al., 2009). These fault systems controlled the subsidence in small foredeep basins of the URF during the early-middle Miocene.

Study method
The footprint-bearing layer crops out as a large scarp cliff (12 × 10 m), and only its lower part is accessible. This layer is coarse-grained, dark brown sandstone and includes numerous carnivore and bird tracks with very large burrows or trails, all of them preserved in convex hyporelief. A quadcopter unmanned aerial vehicle (UAV), model DJI Phantom 3 GL300C, was used to prepare photographs of all parts of the layer (Fig. 3).
This research has used the UAV photogrammetry method to produce an accurate 3D model and orthomosaic of the bedding surface. The wall convergent imaging, at four different height rings, was undertaken using a 12 megapixel Sony Exmor 1/2.3 4K camera mounted on the DJI Phantom 3 UAV. Five scale bars consisting of coded targets have been used to develop a real-scale 3D model. The 3D point clouds were extracted from the captured images using the structure from motion (SFM) method. Then, mesh and texture were built on a point cloud. In the developed close-range photogrammetry network, the average accuracy of estimating the terrestrial coordinates of points per check scale bars is 1.4 mm. An orthomosaic with a resolution of 1 mm per pixel was produced by connecting these orthophotos, from which the effects of relief displacement and tilt have been removed. Using an orthomosaic similar to a map, the objects on the bedding surface are drawn with an accuracy of about 2 mm (Fig. 3a).
Footprints of the lower part of the cliff, however, were also traced on transparent plastic (Fig. 4). A Canon EOS M2 camera, the software package Agisoft Photo Scan Professional (Educational License), and Cloud Compare software were used (Falkingham, 2012;Falkingham et al., 2018).
No samples were collected, and the geometry of the footprints was measured on transparent paper, including footprint length (FL), footprint width (FW), pace angulation (PA), stride (S), trackway width (TW) and digit length (DL), following general conventions (Leonardi, 1987).

Systematic ichnology
Ichnoorder Carnivoripedida (Vialov, 1966) Ichnogenus Felipeda (Panin and Avram, 1962) Emended diagnosis Plantigrade to semi-plantigrade footprints exhibiting four digits (II to V), each with a spheroidal to ovoid or elongate digital pad. The pads form a semicircle in front of, or in about the front portion of, the metatarsal pads. Digital pad III is often the most anterior, and there never appears to be lateral symmetry on an anterior/posterior (AP) line between digital pads III and IV. Digital pads may be of equal or similar size or may show limited dimensional variation. Impressions of claw tips may be present but are usually absent (emended by Sarjeant et al., 2002).

Discussion
Felipeda is known as four-digit imprints of carnivore footprints, which are distinguished by their round digit imprints without claw impressions. The type specimen of Felipeda lynxi is incomplete, and it does not include complete sole imprints (Panin and Avram, 1962, fig. 15). Thus, the assessment of Felipeda as plantigrade or digitigrade is important. Vialov (1965) introduced Bestiopeda as imprints of four digits, the rear pad 53 mm wide and 62 mm long in B. bestia (type ichnospecies), each digit imprint an elongate oval, digits close together and almost contiguous, imprint of first digit has a noticeably narrower anterior edge, and digits II and III project anteriorly in front of digits I and IV (Sarjeant et al., 2002). Aramayo and Manera de Bianco (1987), on the other hand, described felid footprints as digitigrade and medium sized with imprints of a subtrapezoidal to subtriangular palmar pad or plant and four digital pads. They named these footprints Pumaeichnum and attributed them to the living puma (Felis concolor) based on size and morphology. The diagnosis of Pumaeichnum can be emended to include plantigrade and semi-plantigrade to digitigrade, obligate quadruped tetradactyl feloid footprint impressions exhibiting four digits (II-V), each with a spheroidal to ovoidal digital pad forming a single semicircular arc in front of the interdigital pad, sole and heel impressions; digital pads are of equal or similar size, and impressions of claw tips are usually absent (Remeika, 2001). There are no significant differences between this diagnosis and the diagnosis of Felipeda by Sarjeant et al. (2002). Thus, Pumaeichnum has been considered a junior synonym of the ichnogenus Felipeda (Agnolin et al., 2019). Sarjeant et al. (2002) considered Pumaeichnum and Bestiopeda as junior synonyms of Felipeda.
Pantheraichnus is a large plantigrade felid footprint from the Late Pleistocene in Germany. It has four round-oval digit II-V pads imprinted in front of the large pad. In the best preserved imprint, this large pad is oval-shaped in the manus and more kidney-shaped or rectangular in the pes (Diedrich, 2011). Nevertheless, there are no morphological differences between Pantheraichnus and Felipeda, so it is another junior synonym of Felipeda.
Quiritipes Sarjeant et al. (2002) from the Eocene in Wyoming was named for semi-digitigrade to semiplantigrade mammalian footprints of ovoidal to lanceolate outline and a small to moderate size. The pes is somewhat larger than the manus. Four digits are imprinted; all are acuminate (though not acute) without claws. The lack of claw imprints suggests a non-canid trackmaker and differs from Felipeda by size and the age of the footprints, which are older than the first occurrence of Felidae in the Oligocene.  The earliest felids appeared sometime between ca. 35 Ma (age of the sister group) and 28.5 Ma (minimum age of the earliest fossils) (Werdelin et al., 2010). The name Felipeda, however, has been used for Eocene carnivore footprints as a subichnogenus of Bestiopeda (Felipeda) by Scrivner and Bottjer (1986). Felipeda differs from Pycnodactylopus Sarjeant et al. (2002) in lacking broadly ovoid and broad digital impressions with two phalangeal pads. Abbassi and Shakeri (2005) reported felid footprints from the URF of the Mushampa area and for the first time from Iran and assigned them to Bestiopeda isp. In contrast, Abbassi and Amini (2008) reported poorly preserved felid footprints from the URF of the Eyvanekey section in eastern Tehran as Chelipus isp. Based on morphology and the lack of claw imprints, these reports are revised as Felipeda isp. Specimen Numerous footprints on the large slab, studied in the field and traced on transparent paper. There are 27 footprints arranged in six trackways measured on the lower part of the slab using transparent paper.

Diagnosis
Felidae footprint with four oval digit imprints separate from sole imprints. The footprint is 52 mm wide, and footprint length cannot be measured because the trace is not complete. The digit II and III imprints are closer to each other and have a frontal position in contrast to the lateral digits I and IV, so digit III is 3 mm forward of digit IV, and digit II is 2-3 mm forward of digit I. Digit imprints are 17 × 12 mm in dimensions. Claw imprints are absent, suggesting that they were retractable. Footprints of Felis chaus and Lynx caracal are the closest analogies to the studied footprints (from Panin and Avram, 1962; translated from Romanian).

Description
These plantigrade footprints are Felidae footprints; fourdigital impressions are spheroidal to ovoid in shape, arranged in front of the metacarpal and metatarsal pad imprints. Metacarpal and metatarsal pads show conspicuous posterior indentations. A gap separates digits II-V from the metacarpal and metatarsal pads. The imprints of digits III-IV are larger and closer to each other. No claw imprints are visible. Table 1 shows dimensions of these footprints.

Discussion
Shape of the digit imprints, morphology and size of metatarsal and/or metacarpal imprints, span between digits and between metatarsal and/or metacarpal imprints, size comparison with living felids, and symmetry or asymmetry of the footprints are the main characters for ichnospecies identification of the ichnogenus Felipeda. Based on the new combination of Bestiopeda and Pumaeichnum as Felipeda, and the newly reassigned ichnospecies, there are eight valid ichnospecies of Felipeda: F. lynxi Panin 1965, F. biancoi Aramayo and Manera de Bianco, 1987, F. felis Panin 1965, F. milleri Remeika 1999, F. scrivneri Sarjeant et al. 2002, F. bottjeri Sarjeant et al. 2002, F. stouti (= Pumaeichnum stouti) (Remeika, 1999) and F. paryula (Anton et al., 2016). Kordos (1985) introduced Bestiopeda maxima for mammalian carnivore footprints of great size with the sole impression proximodistally flattened with five separate toe prints that are oval in plain view. Sarjeant et al. (2002) considered this as Felipeda maxima, a new combination. Indeed, this footprint comprises five digit imprints and it should be assigned to Felipeda.
Bestiopeda gracilis was considered as the type ichnospecies of Chelipus by Sarjeant and Langston (1994) and attributed to the Canidae. Vialov (1966) described B. gracilis as a small tetradactyl track, 35 mm long, 29 mm wide, with digits not in contact but fanned out, and with wide digits that are sharply pointed (clawed), and he identified it as the tracks of a small canid (Lucas, 2007). Pehuencoichnum gracilis Aramayo and Manera de Bianco (1987) has been combined with Chelipus (Remeika, 2001) because it has a footprint structure strikingly like that of a canid rather than of a felid, as originally assigned, and it is certainly referable to Chelipus as a member of this ichnogenus. A reexamination of photos of Pehuencoichnum shows that it has claw imprints and should be assigned to Canipeda (Melchor et al., 2019). Sarjeant et al. (2002) considered Bestiopeda gracilis as Felipeda gracilis in a new combination. The photo of the holotype of B. gracilis shows exactly tapered digit imprints with a small metatarsal and/or metapodium and a gap between the digit and metatarsal and/or metapodium pad imprints. Felipeda lacks claw imprints or tapered digit imprints, so B. gracilis should not be combined with Felipeda.
In the last report of vertebrate footprints from the lower part of the slab studied here, 25 of the footprints were identified as Bestiopeda isp. (Khoshyar et al., 2016a). These footprints have spheroidal digital and metatarsal and/or metapodium pad impressions without claw imprints, so we assign them to Felipeda lynxi.
Ichnogenus Canipeda (Panin and Avram, 1962) Emended diagnosis Digitigrade to semi-digitigrade, tetradactyl, paraxonic, longer than wide footprints; arranged in quadrupedal and homopodial trackways. Elliptical, similar-sized and clawed digital pads form an arc in front of, and are well separated from, a large metapodial pad. Metapodial pad is rounded to triangular or trapezoidal (emended by Melchor et al., 2019).

Discussion
Canidae footprints are known as elongate footprints with claw marks that suggest non-retractable claws, a feature of canids. The digit imprints are close to each other with a symmetrical pattern along the AP line, and metatarsal and/or metacarpal imprints may be asymmetrical in the posterior lobe (Iliopoulos et al., 2012;Sarjeant et al., 2002).  The tetradactyl, quadrupedal, digitigrade, paraxonic and homopodial tracks include Canipeda, Creodontipus, Quiritipes, Tetrastoibopus, and Bestiopeda (Melchor et al., 2019).
Tetrastoibopus was introduced for asymmetrical, plantigrade to semi-plantigrade footprints exhibiting four single small pads on the digits (II to V) and short claw imprints (Sarjeant and Langston, 1994). Carnotipus has large, asymmetrical tetra-to pentadactyl footprints with elongate imprints with tapered distal ends and without claw imprints (López et al., 1992). Sarjeant and Langston (1994) established Chelipus as a canid footprint and replaced it with Bestiopeda gracilis (Vialov, 1965). Canipeda differs from Bestiopeda either by larger size or by closely set, ovoid digit prints. Metatarsal pad imprints of Canipeda are larger, and the ratio of metatarsal pad length to footprint length (MPL / FL) is about 0.50. Bestiopeda, however, has a MPL / FL value of 0.36. Abbassi and Shakeri (2005) reported Creodontipus from URF in western Zanjan as a tetradactyl pes overprinted on the manus footprints without claw imprints: first Canidae footprints from the Aghajari Formation (late Miocene-Pliocene) in the Zagros Mountains, southern Iran, had been reported, and because of the presence of distinctive claw imprints, these were attributed to Canipeda (Abbassi, 2020b). Specimens Numerous footprints in the studied slab. The geometry of six footprints in the lower part have been studied and measured.

Diagnosis
Canipeda footprints with four oval digital pads, closely spaced with distinctive claw footprints. Footprint length/width ratio is 1.6 to 1.9, with low divarication of outer digits (II-V) less than 30 • .

Description
Tetradactyl footprints with ellipsoidal or sub-spheroidal digit imprints, spread from large metapodial imprints. Middle digit imprints (III-IV) are longer than the lateral digit (II and V) imprints. Metatarsal/metapodium imprints are of the same size and heart shaped with two to three posterior lobes. Distinctive claw imprints are visible in front of the digit imprints (Fig. 7). The contact line of the claw imprint with the digit imprints is well preserved in some footprints. The positions of the claw imprints on digits III and IV are straight and inclined on lateral digits II and V. The claw imprint of digit III is aligned in front of the digit, and the claw imprint of digit IV, however, is inclined to the midline of the footprint. Positions of pes and manus imprints are not constant, and many pes imprints are overlapped on the manus imprints or separated from each other, or they are isolated footprints. Pes and manus imprints are of the same size. Table 1 shows the dimensions of these footprints.

Discussion
At first, Canipeda was established as a monospecific C. longigriffa, which includes elongate tetradactyl footprints with long claw imprints. Its claw imprints converge in the middle digit imprints and are straight in lateral digit imprints. Based on the presence or absence of claw imprints as a distinctive diagnostic feature, some of the Bestiopeda ichnospecies were revised in the new combination ichnospecies of Canipeda by Sarjeant et al. (2002). In addition to the type ichnospecies (C. longigriffa), there are three valid ichnospecies of the ichnogenus Canipeda (Melchor et al., 2019): C. gracilis (Vialov, 1965), C. sanguinolenta (Vialov, 1965) and C. therates (Remeika, 1999).
Canipeda gracilis is known as having larger manus than pes imprints and a footprint length/width ratio between 1.0 and 1.3, less than in C. longigriffa. The ichnospecies C. sanguinolenta (= Bestiopeda sanguinolenta Vialov, 1966) has relatively large, tetradactyl carnivore footprints with tapered digit imprints. C. therates is a small-sized footprint with small metatarsal pad imprints.
Hyaena footprints differ from canid footprints by their larger and closer digit imprints (Abi-Said and Abi-Said, 2007). Iliopoulos et al. (2012) reported tetradactyl footprints with large claw imprints from the late Miocene of Crete (Greece) that might be referred to a large-sized hyaena. Anton et al. (2004), however, believe that assigning C. longigriffa to Canidae is incorrect because of the large size of the interdigital pad relative to the toe pads, as well as the detailed morphology of the interdigital pad. They attributed C. longigriffa to a herpestid (mongoose). To assess this idea requires the reexamination of the holotype of C. longigriffa, which is lost, so only its photo is available. Ichnoorder Avipedia Vialov 1961 Morphofamily Gruipedidae Sarjeant and Langston, 1994 Ichnogenus Iranipeda (Lambrecht, 1938) Emended diagnosis Large tetradactyl avian track with prominent broad, consistently diverging and acutely tapering digit impressions. Prominent heel pad and hallux impressions present. No traces of webbing present (Doyle et al., 2000).

Discussion
Large tetradactyl bird footprints have long been known from Miocene sediments. They include the Ardeipeda gigantea of Panin and Avram (1962), which is 23 cm long and 16.3 cm wide, and Iranipeda abeli, which is about 24 cm long and 22 cm wide (Lambrecht, 1938). Distinctively, Ardeipeda differs from Iranipeda by jointed digit imprints to the metapodium, and the hallux imprint is aligned to the digit III imprints. Also, the ratio of footprint length to hallux length of Iranipeda abeli (FL / DL hallux = 5.7) is greater than that of Ardeipeda gigantea (FL / DL hallux = 2.8). Avipeda filiportatis is large, four-digit tracks up to 19 cm long (Vialov, 1965), and Abbassi et al. (2020) reported large bird footprints from the late Miocene-Pliocene in northern Iraq with 22 cm length and 17 cm width. Lockley and Harris (2010), however, transferred Avipeda filiportatis to Ardeipeda filiportatis. Among these footprints, Iranipeda is large, tetradactyl avian footprints with prominent, broad, consistently diverging and acutely tapering digit impressions with circular heel imprints (Abbassi et al., 2020). Ardeipeda gigantea is characterized by connected and forward-directed digit imprints and backward-directed, long and isolated digit I imprints.
Ardeipeda filiportatis shows distinctive circular heel imprints with connected, forward-directed digit imprints and has a shorter digit III than the lateral digit, and the hind digit I is much shorter (Vialov, 1989). Sarjeant and Langston (1994) considered Iranipeda as a synonym of the ichnogenus Gruipeda Avram, 1962. However, Doyle et al. (2000) and Abbassi et al. (2015) concluded that Iranipeda is distinctive from Gruipeda and is a valid ichnogenus. Iranipeda are attributed to large birds with unwebbed feet, such as Ardeidae.
Ichnospecies Iranipeda abeli Lambrecht, 1938 Fig. 8 Specimen Numerous bird footprints on the large slab, about 35 footprints mapped in the lower part of the slab.

Diagnosis
Large tetradactyl bird footprints with circular heel imprint and isolated and wide digit imprints. The digit III imprint is larger than the lateral digit imprints, and the digit I imprint is distinct. No interdigital webbing visible.

Description
Tetradactyl footprints with three straight, large, forwarddirected digit imprints and a smaller hind digit I imprint.
Most of the footprints show isolated digit imprints with a distinctive circular heel imprint. Imbricated scale imprints are visible on some heel imprints. Three-digit pad imprints are visible in some of the footprints. Tips of digits are sharp without separate claw imprints.

Discussion
Three ichnospecies have been assigned to Iranipeda: I. abeli, I. millumi and I. intermedia (= Gruipeda intermedia Panin 1965, new combination by Doyle et al., 2000). I. abeli has wide digit imprints but a relatively shorter digit I than in Ardeipeda. The imprint of digit I of I. millumi is not clear. Doyle et al. (2000) reported tridactyl digit imprints like Iranipeda isp. with broad digits, identifiable heel impression and a relatively small size. The hallux (digit I) print is absent in this single tridactyl Iranipeda print. Based on the diagnosis of Iranipeda, it is difficult to assign it to Iranipeda. Panin (1965) introduced Gruipeda intermedia as large tetradactyl bird footprints (FL = 15.5 to 16 cm) with slender digit imprints. Based on the size and the position of the digit imprints, it has been considered as an ichnospecies of Iranipeda intermedia in a new combination. Iranipeda intermedia definitely has isolated digit imprints with a circular heel imprint with a short digit I imprint. Iranipeda intermedia was mostly made by crane-like birds and are distinct from heron footprints, which show all the digits joined at a single spot and digit I that is short and rarely registers.

Large grazing or locomotion traces Description
Numerous large grazing or locomotion traces are visible on the studied slab (Fig. 3). These trace fossils are preserved in convex hyporelief and were created first and then strongly affected by mud-cracks or trampled by vertebrate tracks, as well as, of course, by weathering after exposure. Thus, tracing the paths of these trails is difficult. These traces are straight and unbranched, curved, and with looping and selfcross-cutting and are subcylindrical in shape with large striations running along some of them. Sediment fill of the subcylinder is the same as the host rock sandstone. The width of the traces varies from 10 to 40 cm (mean 35 cm), and the longest trail has a course that is 18 m long. The depth (relief) of the traces reaches 2 to 5 cm. The rims of the traces are mostly ambiguous. One or two non-continues tubular bulges are preserved in the middle part or on both sides of the traces.

Discussion
The identification and interpretation of large burrows or trails is difficult, in either attributing them to distinct trace makers or distinguishing them from inorganic structures (Rodríguez-Tovar et al., 2008). The large burrows or trails in the Chehrabad tracksite are likely attributable to a relatively large animal grazing in or crawling on the sediment because of self-crossing or looping structures. Megapermichnus aliabadensis is an elongate, bifurcated and inclined structure from Permian sediments of the Alborz Mountains, which has been interpreted as being produced by a fish (Jenny and Jenny-Deshueese, 1978). Giant, tubular, unilobate trace fossils similar to Planolites in shape and design, but markedly larger in size, have been named Megaplanolites ibericus (Calvo et al., 1987).
Megagyrolithes asdescensis is a large spiral burrow attributed to crustaceans (Gaillard, 1980). Myers et al. (2019), in an aquaria-based study, reported sedimentary structures resembling planar lamina and ripple cross lamina with grain sizes ranging from fine sand to gravel made by a sand-sifting goby, and they noted that these structures may have a long history in the geological record reflecting feeding styles and bio-resuspension behaviors of macrofauna from different time periods (Myers et al., 2019). Numerous gutter-like furrows, up to 60 cm wide and up to 9 m long, have been reported from the shallow Middle Jurassic sea floor from Liesberg, Switzerland, and interpreted as feeding traces made by large marine vertebrates, most likely plesiosaurs and ichthyosaurs searching for food in the lime mud (Geister, 1998).
The large traces described here resemble Megaplanolites, to which they are tentatively assigned, and may have been created by a large snake-like animal. Megaplanolites has been reported from shallow marine carbonate shelves to littoral conditions; however, the Chehrabad Megaplanolites was created in a continental environment. Distinctive features of the somewhat similar but small-scale cylindrical burrows, Planolites and Palaeophycus, were well discussed by Pemberton and Frey (1982) and Keighley and Pickerill (1995): Planolites is diagnosed primarily on the basis that it is an unlined simple burrow that was actively backfilled with sediment having textural and fabric-like characters unlike those of the host rock. In contrast, Palaeophycus is a lined burrow and represents passive sedimentation, filled with sediments typically identical to those of the surrounding matrix.
Some of the large traces described here show both characteristics of Planolites and Palaeophycus, including no distinctive wall or lining, sediment fill the same as the surrounding matrix and striations on the surface of some of the traces (Fig. 9). These large traces may be the crawling traces of a snake-like animal on the surface during the first phases after sedimentation. A snake crawling on sediment produces closed to open trails with lateral sediment deformations and striations of variable width (Abbassi, 2020a) (Fig. 10). The width variations of the traces described here, presence of striations, shape and self-crossing suggest that these trails may have been made by a snake or snake-like animal. This compression is not outlined by lateral sediment deformation, so crawling traces show sharp and unique rims with surrounding sediments. Usually, forces generated by snakes are using variable numbers of points for lateral undulation movement across a flat surface (arrows in the sketch of Fig. 10). Thus, the relief of a crawling impression may differ along the trace, which may not be present in the studied traces. Further study of these traces is planned.

Discussion
Despite the abundant plant fossils in the studied area, unfortunately there are no paleobotanical studies in the URF. A few and sparse herb, reed and tree trunk fossils have been identified in the study area (Fatemeh Vaez Javadi, personal communication, 2020). An early Miocene gymnosperm tree trunk fossil has been reported in the Ardabil area, northern Zanjan Province, northwest Iran (Mustoe et al., 2020). Thus, it seems that around the playa and the pond sedimentary environments of the URF, there was grassland meadow vegetation. Abundant middle Miocene carnivore and herbivores lived in the savanna plain of the URF environment.
Based on previous reports, Miocene trackmakers of the URF include carnivores (felids, Creodonta, mustelids, canids), birds (Charadriiformes, Anseriformes, Gruiformes, Ciconiiformes) and herbivores (artiodactyls, proboscideans) (Abbassi, 2010;Abbassi and Amini, 2008;Abbassi and Shakeri, 2005;Alavi et al., 2016;Khoshyar et al., 2016a, b). The tracksite of the URF in Chehrabad studied here comprises abundant imprints of canids and felids with sparse bird footprints. The mean length of the canid footprints in the lower part of this tracksite is 57 mm with a FL / FW ratio (57 / 58 mm) of 0.98, which is less than the known ichnospecies of Canipeda (Melchor et al., 2019). Today, canids are cosmopolitan mammals, but they were restricted to North America for all of their early history, not reaching the Old World until the Miocene at the latest. The oldest and most primitive canids are the Hesperocyoninae of the late middle Eocene (Duchesnean) to the middle Miocene (Barstovian) in North America (Rose, 2006;Wang, 1994). Canidae was not present in Europe until the Miocene and not in Asia and Africa until the late Miocene (Kemp and Kemp, 2005). Late Miocene Vulpes riffautae is the oldest Canidae fossil from Africa (De Bonis et al., 2007). The migration of Canidae from North America to Eurasia and their first appearance in Eurasia is known as the "Eucyon event", which occurred in the late Miocene (Sotnikova and Rook, 2010;Rook, 2009;Sotnikova, 2010). The late Miocene Maragheh fauna in northwest Iran is a well-known late Miocene faunal site (Bernor, 1985;Kostopoulos, 2009). Canidae fossils, however, have not been reported from this pre-eminent vertebrate fossil site. The numerous Canidae footprints from the Chehrabad tracksite are the first reported Canidae fossil remains from the Iranian Plateau, and they occur during the Tortonian.
The earliest records of felids are from the Oligocene; the genus Proailurus appears just after the end of the Eocene in Eurasia and expands to other continents in the late Miocene Johnson et al., 2006;Gaubert and Veron, 2003). The early-middle Miocene radiation of Felidae is complex, so Pseudaelurus is reported from Europe, Arabia, China and North America (Werdelin et al., 2010). Felis attica, Metailurus orientalis and Amphimachairodus aphanistus are Felidae from the middle and up- per parts of the Maragheh Formation (8.16-7.4 Ma), and the latter two felids are of large size (Bernor, 1986;Ataabadi et al., 2013). In contrast, the felid footprints documented here are relatively small in size: FL = 61 mm, FW = 59 mm, and FL / FW equals 1.03. These felid footprints are smaller in size then the Mushampa sample (= Bestiopeda isp.; Abbassi and Shakeri, 2005): FL = 125 mm, FW = 118 mm, and FL / FW equals 1.05. Overall, the felid footprints are fewer in number than the canid footprints in the studied slab. Nevertheless, they represent an important regional record of felids.

Conclusion
The tracksite of Chehrabad is one of the numerous vertebrate tracksites of the URF in the west of Zanjan Province. The URF was deposited in a continental playa-like environment, Figure 10. Recent snake trails on creek sediments. Arrows show lateral pressures during the crawling of snakes (Abbassi, 2020). Length of hammer equals 30 cm. with braided to meandering rivers and ephemeral ponds under a warm climate during the middle to late Miocene. The studied large slab contains abundant carnivore footprints mostly arranged in a track corridor. In addition to sparse large bird footprints, Iranipeda abeli, numerous canid footprints are attributable to Canipeda longigriffa and felid footprints to Felipeda lynxi. The presence of canid footprints confirms that the first appearance of the Canidae in Eurasia was no later than the late Miocene. In this track corridor, a grazing or crawling snake-like animal produced a trace tentatively attributed to Megaplanolites, which formed in a continental setting.
Author contributions. NA and MK collected the investigated material in the field. NA identified ichnotaxa and prepared the manuscript. SGL improved the ichnotaxonomy and edited the manuscript in the English language. FE performed the UAV photogrammetry and accurate 3D modeling and orthomosaic of the studied slab.