Abstract
Footprints represent a relevant vestige providing direct information on the biology, locomotion, and behaviour of the individuals who left them. However, the spatiotemporal distribution of hominin footprints is heterogeneous, particularly in North Africa, where no footprint sites were known before the Holocene. This region is important in the evolution of hominins. It notably includes the earliest currently known Homo sapiens (Jebel Irhoud) and the oldest and richest African Middle Stone Age hominin sites. In this fragmented ichnological record, we report the discovery of 85 human footprints on a Late Pleistocene now indurated beach surface of about 2800 m2 at Larache (Northwest coast of Morocco). The wide range of sizes of the footprints suggests that several individuals from different age groups made the tracks while moving landward and seaward across a semi-dissipative bar-trough sandy beach foreshore. A geological investigation and an optically stimulated luminescence dating of a rock sample extracted from the tracksite places this hominin footprint surface at 90.3â±â7.6 ka (MIS 5, Late Pleistocene). The Larache footprints are, therefore, the oldest attributed to Homo sapiens in Northern Africa and the Southern Mediterranean.
Similar content being viewed by others
Introduction
Since the discovery of the 3.66Â Ma hominin trackways at Laetoli in Tanzania in 19761,2,3 paleoanthropologists have emphasized the importance of footprint studies. Footprints represent snapshots of life, and provide direct information on the biological (stature, age, body mass), biomechanical (speed, gait) and behavioural features of hominin groups3,4,5,6,7,8,9.
Over the past 2Â decades, a large number of hominin footprint sites have been discovered on different continents7,8,9,10,11,12,13,14,15,16,17,18,19. Despite an increase in the number of findings, their spatiotemporal distribution is still highly heterogeneous. For instance, in Africa they are mainly concentrated in the east1,2,3,7,12,20 and south of the continent15,21,22,23. In contrast, evidence from North Africa found in Holocene deposits is sparse4,24. This fact is relevant considering that this region is central to an understanding of the evolution of hominins. Occupied since 2.4Â Ma25, North Africa also includes the Jebel Irhoud site where the earliest Homo sapiens remains were discovered26,27. For more recent periods, the rich Mousterian and Aterian industries28,29 testify to various human settlements during the Upper Pleistocene, particularly in coastal environments. In 2022, in the course of a field mission as part of a research project investigating the origin and dynamics of the coastal boulders that litter this rocky coastline in NW Morocco30, we identified new hominin footprints on a rocky beach in Larache (NW of Morocco). The site covers an area of about 2800Â m2 over which have been preserved a large number of tracks and trackways that were documented by both aerial and terrestrial photogrammetry. This discovery presents an opportunity to enhance our knowledge on the patchy ichnological context of the North African coast through recourse to a multiproxy research approach. A prior analysis of the geology, geomorphology, and stratigraphy of the site sets the context for a description of these footprints and inferences on their biological characteristics from morphometric data. The geochronology of the site and the importance of these footprints from a geographical and temporal perspective are then discussed.
Results
The rocky beach of Larache (NW of Morocco)
Larache is located on the northwestern Atlantic coast of Morocco, south of the Loukkos River estuary and the city of Larache (Fig. 1). The sandy beach to the north of the estuary is made up of fine and medium quartz and bioclastic sand, sometimes covered by dunes up to 10 m high or limited inland by a Quaternary sandstone cliff that constrains the mouth of the Loukkos. The rocky shore platform beach of Larache is situated south of this estuary. The platform is part of the Rmel Plateau which is bounded seaward by actively retreating cliffs with an elevation ranging from 20 to 30 m above sea level. Older inactive cliffs are present 1â3 kms inland. The active cliffs form an almost straight shoreline with a few promontories, including Punta Nador and Laghdira (the northern limit of the study area). At the Larache site, the intertidal width of the platform varies from 20 m at the promontories to 100 m in the widest sectors, with an average width of 60 m. The surface of the platform dips gently seaward or landwards30 (Fig. 1). The cliffs consist of a succession of poorly consolidated Quaternary deposits with large foreset beds of aeolianites, thin red sandy edaphogenic levels with continental gastropods, and metric cross-bedded sets of beach deposits30 (Fig. 2).
Geological sketch map of the Larache area and panoramic photograph of the rocky shore platform showing the location of the indurated beach. The area delimited by the dotted red line corresponds to the footprint discovery zone. The vertical white line corresponds to the location of the sediment log analyzed in this study. The geological sketch map was generated by QGIS software (v.3.28.4) software (https://www.qgis.org/fr/site/forusers/download.html) based on the Geological Map of Morocco (Scale: 1/1,000,000 N.M.S.Gâ260â1985).
As these quartz-arenite cliffs retreat, the sandstones on the shore platform surface are exposed (Fig. 1). Overall, the stratigraphy of the Larache site is characterized by the alternation of different types of sedimentary deposits and the ongoing erosion of the cliffs, which has led to the exposure of the underlying sandstones on the shore platform surface30. The rocky shore platform sector where we discovered the footprints is geologically interpreted as a fossil beach foreshore and appears to be better preserved from erosion and weathering than the adjoining area, since it is partially protected from the lateral impact of swell on the south coast by a large boulder and frontally by the higher elevation of the platform (Fig. 1). On the north coast, the site is highly eroded, with partial collapse of the platform, likely leading to future weakening and demise of the footprint site (Fig. 1). The footprints (Figs. 3 and 4) are stratigraphically located in the lower part of the Pleistocene succession of the Larache cliffs, being part of its basal unit (Unit 1). The tracked surfaces appear at the top of two sheets of the lowermost visible stratum. The tracked layers are seaward-inclined plane beds, exhibiting isolated or aligned footprints. The sediment sample used for OSL dating was taken from the lower tracked surface to the northeast of the tracks site (Figs. 2 and 3, Fig. S1).
Images of some of the hominin tracks in Larache. (a) Two tracks side by side at the bottom of the photo, which also depicts a ground control point (chequered cardboard) for differential GPS surveying. (b) Two cross trackways and photography for 3D footprint modelling. (c) to (g) Detailed view of some footprints. White Scale barsâ=â20 cm.
The Larache footprint site. (a) Position of the footprints (white crosses) plotted on an aerial photo from a drone flight (orthophotography). (b) A Digital Elevation Model (DEM) of the site with the orientation of the footprints. The green square represents the point at which the sediment sample was taken for OSL dating. The coordinates of the sampling point are: (X: 428,456,313âY: 506,198,258âZ: 1,7). The wind rose shows the preferential direction of the footprints (orientation of the tracks) relative to the coastline. The maps in (a) and (b) were created using QGIS software (v.3.28.4)(https://www.qgis.org/fr/site/forusers/download.html). (a) was based on the photogrammetric orthomosaic produced using Agisoft Metashape software (version.2.0.1, https://www.agisoft.com/forum/index.php?topic=14904.0). (b) is a Digital Elevation Model (DSM) with a resolution of 1.09Â cm per pixel which was computed using Agisoft Metashape software (version.2.0.1) and imported to a GIS environment using QGIS software (v.3.28.4). (a) was obtained by a low-altitude programmed flight using an Unmanned Aerial Vehicle (UAV) on July the 5th, 2022. A multirotor Mavic Pro Platinum quadcopter (DJI, Nanshan District, Shenzhen, China) equipped with a 12.35 MP camera was collected at nadirs at an elevation of 12 m.
The local geology characterized by a steep coastal cliff is divided into 4 units using sedimentological criteria (Fig. 2). Unit 1 consists of alternating sets of fine sandstones displaying landward-inclined cross-bedded and seaward-inclined parallel laminations. The sandstones are composed of quartz grains of medium-fine sand in a 90% (Fig. S2), but also include some bioclastic grains of bivalve shell fragments, echinoderms and particular benthic foraminifera (Quinqueloculina sp.) (Fig. S2). The rock is consolidated by two different types of cements: a primary calcareous cement and a secondary one of siliceous nature. The discovered footprints are part of tracked surfaces inclined seaward but separated by a set of landward-dipping crossbeds. This cross-bedded set wedges towards the sea so that it disappears horizontally by about 6 m (lower picture in Fig. 2). Thus, the tracked surfaces appear superimposed in the westernmost area and would be considered as two sheets of the same stratum. The lowermost part of this unit comprises of a 2 m-thick layer of more compacted horizontal laminated shelly sandstones, with low porosity. The levels of Unit 1 under these shelly sandstones are not visible. This unit includes the surfaces tracked by hominins. The OSL-dated tracked layer yielded an age of 90.3â±â7.6 ka (see Table 1 and âMethodsâ section). Unit 2 starts with a metric edaphic level of red silty sands. However, most of the unit consists of a 7 m-thick set of complex cross-bedded deposits that culminates in a strongly root-bioturbated level. This complete unit can be interpreted as an aeolianite. Unit 3 is similar to unit 2. It starts with a meter-thick pedogenic level of red sandy silts and continues with 6 m of several sets of curved-base cross beds. Many of the cross-bedding sets are deformed, possibly by seismic activity. The last of these levels is also strongly bioturbated by roots. This unit has also been interpreted as an aeolianite. Unit 4 consists of a metric stratum of red silty sands of pedogenic origin and is overlain by Holocene dunes.
Ichnological assemblage
The laterality of 91% of the 85 footprints identified at Larache (Fig. 4), yields 30 (35%) as left and 39 (46%) as right footprints (Table S1). This assemblage contains isolated footprints as well as trackways, including 2â6 footprints.
The footprints are preferentially oriented in a WNW-ESE direction (Table S1), relatively perpendicular to the present shoreline, and secondarily in a SWâNE direction (Fig. 4).
The morphology of the Larache footprints is variable, a phenomenon regularly observed particularly in coastal or dune contexts7,9,13,16,18. Indeed, footprint morphology is not only impacted by the biological and biomechanical features of the track-makers but also by the nature of the substrate and the taphonomic agents4,31,32. Despite this extramorphological variation, the Larache footprints reflect human foot characteristics such as a rounded heel and an adducted hallux. The toe impressions can be difficult to discern from the rest of the footprint, which has also been observed at other sites with similar depositional conditions9,13,18. However, the hallux is distinguishable on some footprints (Fig. 5âe.g., C1-01). The depth distribution is also variable, although trends can be observed: the area associated with the heel and the forefoot are the deepest. The deepest area of the heel is generally located in its middle. The deepest area of the forefoot is usually located in the areas associated with the most medial metatarsal heads or toes (hallux and second toe). This depth distribution is consistent with the distribution of plantar pressures during human walking. During the stance phase, the heel is the first part to contact the ground (âheel strikeâ). The weight is then transferred to the forefoot, and in particular, the metatarsal heads and then the toes (especially the most medial) during the propulsion33,34. Finally, a difference in depth can be observed depending on the footprint dimensions. The smallest footprints are shallower than the longest ones (Fig. S3).
Hominin footprints from Larache. (a) Natural views and shaded 3D elevation of (a) C01-01, (b) C03-01, (c) C02-01 and (d) C03-21). Scale bar, 2 cm. 3D models were generated with Agisoft Metashape software (version.2.0.1, https://www.agisoft.com/forum/index.php?topic=14904.0) using between 18 and 20 images taken with a Nikon D7500 (20.9 MP, Nikon AF-S DX35mm f/1.8 G) from a height of approximately 0.5â1.5Â m.
The length was measured for 81 of the 85 footprints (Table S1). It ranges from 12.7 to 30.0 cm (mean: 22.7 cm) following an asymmetric normal distribution (ShapiroâWilk: p valueâ>â0.05). Footprints longer than 20 cm are more numerous (80%) than those shorter than 20 cm (20%). These lengths correspond to estimated statures between 120.8 and 189.0 cm (average: 160.0 cm). Most of these footprints (89%) have an estimated stature higher than 140 cm (Fig. 6).
The lengths correspond to the three age classes: children, adolescents and adults. The distribution of the footprints according to these three age classes is also relatively balanced (Fig. 7). There are slightly more footprints attributed to children (31) than to adolescents (26) or adults (24). However, the model used to estimate age classes has some uncertainties based on the average foot length for each age (cf. âMethodsâ section). Thus, footprints associated with older children may have been made by adolescents, and vice versa. The same uncertainty exists for the footprints close to the boundary between adolescents and adults (cf. Discussion).
These palaeobiological estimates of stature and age are made for each footprint and not for each individual. Although part of the assemblage is composed of isolated footprints, it is possible to estimate a minimum number of individuals from the experimental knowledge of intraindividual variation (cf. âMethodsâ section). The 81 measured footprints correspond to at least 5 individuals. Among this minimum of 5 individuals, there is a young child (between 1 and 4Â years), an older child (between 4 and 8Â years), an adolescent or small adult (148.8â161.8Â cm), a medium-sized adult (162.9â175.9Â cm) and a tall adult (176.2â189.0Â cm) who is likely male based on the sexual dimorphism of height known in fossil and extant Homo sapiens35,36.
Discussion
The sedimentological features of the tracked layers are evidence informing about the environment in which footprints were made. Unit 1 is composed of a rhythmic alternation of two different primary sandy facies. Facies of land-dipped cross-beds and facies of sea-inclined parallel lamination. The sediment that constitutes both facies is a medium-to-fine sand, mainly of quartz, but also including bioclastic grains. The grains are partially compacted by calcareous and silica cements. The grain size, gentle slope, and conjunction of these facies in the form of a sequence is the same as that identified on the foreshore zone of a semi-dissipative beach with transverse landward migrating bars and troughs37,38. Thus, the landward-inclined crossbedding represents the migration of the bars. In contrast, once a bar becomes attached to the beachface, a seaward-inclined slope is generated, wherein the waves build up parallel sheets of sand inclined in the same sense until the arrival of a new bar.
In the Larache case, the footprints were probably left: (1) in the course of fair-weather wave conditions (no storm surge), at the limits of swash flow (2) and the landward limits of the spring high tidal zone over the beach foreshore bar (3), three conditions that would have corresponded to limited swash flow duration and very thin laminar swash flow, all potentially the most favorable conditions in the course of the spring-neap tidal cycle to salt-crusting and inception of preservation. These conditions correspond to a situation where the footprints left on the upper beach were not washed out by vigorous swash on the beach face and have time for the salt-crusting that we deem as initiating the preservation process. It is interesting to note that the tracked surfaces are seaward-inclined whereas these surfaces are covered and preserved by sets of landward cross-beds. The presence of bioclastic shell fragments indicates a marine contribution to the essentially terrigenous beach quartz sand. The presence of a primary calcareous cement suggests the possibility of the occurrence of an early cementation that ensued in the same sedimentary environment and forming a hardened beach. This process could contribute to the final preservation of the footprints. The dependence of footprints preservation essentially on taphonomic processes that potentially preserve bedforms, such as the afore-mentioned early cementation (of the sediment or within the microbial mat), rapid covering by sediment, and overgrowth by microbial mats has been well-described39. According to Thulborn40, footprints are most commonly preserved in environments that experience periodic or cyclic sediment accumulation as observed in our case. In the Larache case, the mobility of the bars covered up footprint-bearing surfaces located within water-saturated troughs separating the bars. This is observed in other coastal environments with similar dynamics41,42,43,44,45.
Considering the OSL dating, the estimated age for this track site is 90.3â±â7.6 thousand years (ka) (MIS 5c, formerly known as the "Ouljian" throughout Morocco). Studies on sea-level fluctuation, geomorphological and sedimentary context of the Larache coastal zone are sparse. Pierre et al.46 referenced an Ouljian outcrop situated approximately 1 m above the present sea level on a cliff north of the Loukkous estuary. In contrast, in a description the currently active cliffs south of Larache, Adil47 mentioned the presence of the Ouljian at the base of these cliff but did not specify a date or chronological context. To the north of Larache (Tangier Region), El Kadiri et al.48, and El Abdellaoui et al.49, have predominantly focused on dating MIS 5 marine terraces (between MIS 5a and MIS 5e), while underscoring the influence of neotectonics in uplifting sedimentary facies of the Gibraltar arc by approximately 0.06 mm/year to 0.1 mm over the last 100,000 years. This corresponds to an elevation of 7â10 m. South of Larache, between Rabat and Casablanca, MIS 5 sea-level oscillations have been extensively documented. Weisrock50 and Chahid et al.51 placed sea level during MIS 5c at around ââ20 m above the present level. This MIS 5c level is observed at various elevations between present sea level and +â6 m. It was marked in its final phase by the formation of a high continuous aeolian sandy ridge dated between 103.8â±â8.2 and 93.3â±â7.1 ka51. The chronostratigraphic evidence both north and south of Larache is consistent with the dating of the sedimentary level on which the tracks were established on a beach foreshore environment, subsequently preserved in the medium and long terms by the formation of a dune. However, the impact of neotectonics and the rate of sedimentary facies uplift at our site remain to be determined. Using this dating for estimating the coastal position when the tracks were recorded is feasible. Precise dating of the indurated platform would throw more light on the paleoenvironmental context and deposition of coastal and marine sediments in this region.
The discovery of the Larache footprints represents further evidence of the importance of North Africa, and the Moroccan region in particular, during hominin evolution. Occupied by hominins since at least 2.4Â Ma25, North Africa has notably yielded the earliest occurrences of Homo sapiens: the 315,000-year-old fossils of Jebel Irhoud26,27. Occupations in Morocco have remained important during the most recent periods, including the Middle Stone Age, during which the Larache footprints were left. Fossil remains, and notably, the rich Mousterian/Aterian industries testify to the presence of hominins in the Atlantic and Mediterranean coasts during this period52,53. The contribution of recent studies has considerably clarified our knowledge of the chronology of these records. The western Atlantic coast between Rabat and Tangier, where the Larache footprints are located, was occupied in several places. A large part of these occupations is in the Rabat-Temara region, where paleoanthropological remains and archaeological artefacts attest to occupations by Homo sapiens between 150 and 50Â ka ago and notably during MIS 554,55,56, when the Larache footprints were made.
From an ichnological point of view, the Larache footprints represent an important discovery. Indeed, no other site in North Africa has yielded footprints dating from the Pleistocene or Pliocene. They are, therefore, the oldest human footprints in this region and among the oldest footprints attributed to Homo sapiens worldwide. Only two other regions have yielded older Homo sapiens footprints: South Africa, with several coastal sites, the oldest of which is dated to 150 ka ago15,23, and the Arabian Peninsula, where footprints dated to 120â110 ka ago were discovered at Alathar16. The site of Larache is also important for the number of footprints that have been discovered. Indeed, the ichnological assemblage, including 85 footprints, is important. Most sites, including those mentioned above, have yielded, in the best of cases, less than a few dozen footprints57. The absence of archaeological material (animal bones, tools, occupation structures) or paleoanthropological material (hominin skeletal remains) is another feature of the Larache ichnological assemblage. However, such an association is extremely rare within the fossil hominin footprint record57. Footprints from Melka Kunture7, Schöningen14 and Le Rozel9 are among the only ones to be associated with archaeological material. In the vast majority of cases, ichnological assemblages are not associated with any other type of material.
One of the most original aspects of footprints is that they provide access to very brief periods of life, like snapshots20,58. They can provide information not only on the composition of groups but also on their behaviour6,8,9,18,59. In Larache, the morphometric study shows that a group of young children, adolescents and adults left these prints with a minimum number of 5 individuals. The 81 prints analysed represent a balanced distribution between children, adolescents and adults. Although there is no doubt that the track-making group was composed of different age classes (children, adolescents and adults), it is necessary to exercise caution when making inferences about the group's composition from estimated statures and age classes. First of all, these estimates are based on the morphology of the footprints, which varies even if a single individual makes them. The intra-individual morphological variation can be significant, particularly in dune environments where footprints are left in soft substrates32,60,61. It is, therefore, necessary to be cautious when estimating biological characteristics from footprints, especially when they are isolated and intra-individual variation cannot be quantified. Indeed, in the case of isolated footprints, the difference between estimated and actual stature can exceed 10%61. In addition to this uncertainty linked to intra-individual morphological variation, there is also uncertainty linked to taphonomic modifications. Although footprints discovered in the fossil record are rapidly covered, they can be eroded by taphonomic agents, thus modifying their morphology and consequently potentially biasing palaeobiological inferences62. Secondly, age class estimates are made from an average variation between foot length and age obtained for different modern populations. This average variation is not perfect since the attribution of age classes for some individuals may be erroneous, particularly for individuals whose footprint lengths are close to the limits between two age classes. In the case of the Larache footprints, many are close to the limit between children and adolescents. Thus, some footprints attributed to adolescents could have been left by children and vice versa. An identical uncertainty is present for the limit between adolescents and adults. Furthermore, the relative distribution obtained may also overestimate the number of children for biometric reasons. Indeed, children need to make more strides to cover the same distance as an adult, leaving more footprints. Lastly, the composition reflected by the footprints is not necessarily that of the entire social group that occupied the area. Only part of the social group may be present, as the performance of certain tasks may have led to the selection of individuals5,6,59.
A key question when studying footprints is what the individuals were doing at the site. Since no occupation structures were found, this site may correspond to a passage and/or foraging site. While Pleistocene Homo sapiens were hunter-gatherers, individuals likely left the Larache footprints while probably searching for resources. Numerous archaeological discoveries, particularly in Morocco and notably in the Rabat-Temara region, have shown the importance of coastal areas for access to resources, whether raw materials, prey or even plants52,63,64,65,66. In this context, the preferential orientation of the Larache footprints towards the offshore could maybe indicate the search for marine resources. The presence of young children, in Larache, possibly contributing to the search for these resources, could provide unique information on the social behaviour of Pleistocene Homo sapiens populations. However, further studies will be needed to validate this hypothesis.
Conclusion
The Larache coastal hominin tracksite reported here is dated to 90.3â±â7.6 ka and represents one of the worldâs largest and best-preserved Late Pleistocene tracks sites, and the sole documented site in North Africa and Southern Mediterranean. Morphological characteristics undoubtedly relate these tracks to hominins, and due to their geological age and geographical location, they belong to Homo sapiens. According to the morphometric comparisons made with the complete footprints, the statures of the track-makers vary from 120.8 and 189.0 cm (mean: 160 cm), and the tracks to a minimum of 5 individuals. The distribution of the footprints with respect to the shoreline provides a snapshot of the movements of multigenerational Homo sapiens individuals (children, adolescents, and adults) in this site. The distribution supports the ecological relationship between Homo sapiens populations and coastal areas where no other type of anthropological or archaeological evidence has yet been found. The protection of the Larache site from intense marine erosion until now has enabled the intact preservation of the traces. However, the ongoing collapse of the rocky shore platform in the northern zone of the site could lead to its eventual demise, and that of the tracks it has preserved thus far. The loss of the site in the medium and long term is likely to be caused by rising sea level and storm events. In the short term, other new footprints will be discovered as sediments are eroded. It would be interesting to monitor this erosion to expose new, complementary tracks that could be used to clarify hominin group size. The various caves along Larache's southern coastline should also be explored for any fossils or lithic traces that might be present.
Methods
Geological context and sedimentary analysis
The geological context was determined by standard field mapping of exposures and outcrops. The detailed sedimentary column of the sequence in which the trackâs surface is preserved has been recorded, including dip direction and facies description through analysis of the grain characteristics (Fig. 2).
Footprint analysis
Footprint identification
The identification of the footprints was based on the anatomical characteristics of the human foot: a rounded heel, a plantar arch, relatively short toes and an adducted hallux4,67.
Footprint recording
Descriptions, measurements, and photos of footprints were made directly in the field (Fig. 2). We carried out a low-altitude UAV flight of the footprint surface on the indurated shore platform and a series of DGPS measurements of the most visible tracks and the limits of footprint surfaces (clusters). UAV images were collected at nadirs along calibrated flight paths parallel to the coast, at an elevation of 12 m using a Mavic Pro Platinum quadcopter (DJI, Nanshan District, Shenzhen, China) equipped with a 12.35 MP camera. DJI GS Pro Flight Planner software⢠was used for programming the automatically and autonomously executed flight plan. High-overlap image sets ensure sufficient target points in adjacent images as the standard procedure; we planned the flight paths to provide 80% overlap between successive images and 70% sidelap between adjacent paths. A total of 461 photographs were collected over the study area. We used Agisoft Metashape Professional software (v.1.6.2) to process the images. During the flight survey, 12 ground control points (GCPs) consisting of black-and-white chequered cardboard, essential to produce georeferenced aerial images, were deployed at different locations on the shore platform to provide the widest possible coverage (Fig. 3a). Their positions were recorded using a Trimble differential GPS in real-time kinematics mode (Fig. 3a). Digital elevation models (DSM) with a resolution of 1.09 cm per pixel and orthophotographs with a resolution of 5.43 mm per pixel were computed using Agisoft Metashape software and imported to a GIS environment using QGIS software (v.3.28.4) to characterize the spatial distribution and orientation of the footprints. The data collected during this flight were used to characterize the spatial distribution and orientation of the footprints using QGIS software (v.3.28.4). In addition, 28 footprints were scanned by close-up photography and vertical mosaic of high-resolution digital images taken with a Nikon D7500 (20.9 MP, Nikon AF-S DX35mm f/1.8 G) from a height of approximately 0.5â1.5 m (Fig. 3b). Photographs were photogrammetrically modelled using Agisoft Metashape software (v.2.0.1) to obtain 3D models of the best-preserved footprints.
Morphometry
The lengths were measured in 3D along the longitudinal axis of the footprints. For a complete footprint, this length corresponds to the distance between the base of the heel and the tip of the second toe. The measurements taken from the 3D models were compared to those obtained directly in the field and those measured on the orthomosaics obtained with the drone (Table S1). No significant difference was observed between the measurement techniques (KruskalâWallis: p valueâ>â0.05; mean differenceâ<â0.1 cm).
Stature estimates
Stature (S) is a variable commonly estimated from footprint length (L). This estimate is based on the strong relationship between foot length and stature (between 13 and 16%). One stature was estimated for each Larache footprint whose length was measured using two published experimental regressions. Both experimental studies were carried out under substrate conditions close to those found at Larache. The first regression equation (Sâ=â3.8âÃâLâ+â73.3, râ=â0.78)68.The second regression Eq. (4.1âÃâLâ+â67.6, râ=â0.69)68. The stature presented in the results section corresponds to the average of the estimates obtained by these two equations. The differences in stature between these two regressions for the Larache footprints sample are minor (mean: 1.1 cm) (Table S1).
Age class estimates
An age was also estimated from the length of the Larache footprints. To do this, a curve representing the growth in foot length (Fig. 7, Fig. S4) in modern populations was realised from published data69,70,71,72,73,74 including almost 12,000 individuals from different geographical areas (America, Europe, Africa, Asia). The Larache footprints were positioned on the curve according to their length. An age class was then estimated for each footprint based on its position on the curve.
Minimum number of individual estimates (MNI)
Although the footprint assemblage is partly composed of isolated footprints, it was possible to determine a Minimum Number of Individuals. For this, we used the intra-individual variation of footprints. Actually, footprints of a single individual in a soft substrate have different sizes and shapes. However, this intra-individual morphological variation is limited; it is possible to use this limit for a sample of isolated footprints to determine a minimum number of individuals. Several studies on experimental or Holocene footprints have shown that this intra-individual variation can be important; the deviation from the individual mean can exceed 10%4,32,61,68. Therefore, we used the same assumption as Ashton et al.13 in their study of Happisburgh footprints by considering that a single individual could have made footprints falling within 屉10% of each other.
Luminescence dating
Sampling method
A rock sample was extracted from a zone of sedimentary continuity on the surface of the footprints for OSL dating. This sampling was carried out with all possible precautions and following the sampling protocols for OSL dating. A rectangular sample measuring approximately 15 cmâÃâ15 cm was extracted by hand using a hammer and chisel and packed directly into an opaque PVC box. The sample was collected at the lower tracked surface (See Fig. 2, Fig. S1). The coordinates of the sampling point are: (X: 428,456,313âY: 506,198,258âZ: 1,7). A fragmentation of this sample was used for analysis by transmitted light microscopy (Fig. S2).
Sample preparation procedure
The outer layer of the sample was removed to a depth of 25 mm on each side to ensure that the material exposed to sunlight was effectively eliminated. Subsequently, the core of the sample was carefully crushed and wet sieved to obtain grain fractions of 90â180 μm and 180â250 μm. The 180â250 μm fraction was chemically treated to isolate the quartz grains. Initially, it was exposed to a 10% HCl solution until the carbonates were completely dissolved. Following this step, H2O2 was employed to remove any organic matter present. Two density separation processes were then carried out using sodium polytungstate to separate quartz from feldspar and other heavier minerals effectively. Finally, a 40% HF solution was employed to remove any remaining feldspar residuals and for etching the outer layer of quartz, which might have been influenced by alpha radiation. Lastly, the resulting quartz grains were carefully dried and sieved once more, with the specific fractions falling within the 180â250 μm range selected for subsequent measurements.
SAR procedure
Samples were measured using the SAR (Single Aliquot Regenerative dose) procedure69. A portion of 180â250 μm quartz grains from the sample was divided into aliquots that were measured individually. Table 2 shows the SAR protocol used. The OSL signal was measured for 40 s at 0.1 s per data point, giving 400 data points (Fig. S5). OSL signal was measured using a blue LED light source. OSL signal was collected with a reading temperature of 125 °C. The signal was dominated by the fast component, considering the first 5 data points (0.5 s) for measurement. The background was calculated considering the end part of the spectrum (50 data points) and subtracted from the measured signal. No feldspar contamination was detected using IR stimulation. Figs. S6 and S7 have been chosen as representative of the OSL measurement process. They show the OSL decay and doseâresponse curves obtained from the sample.
The resulting populations are normally distributed with overdispersion valuesâ<â20% (Fig. S8). Central Age Model (CAM)70 has been applied to calculate the equivalent dose (i.e., the accumulated dose due to the ionising radiation received by the quartz grains over the period they have been buried) of each sample.
Pre-heat dose recovery test
A pre-heat temperature of 220 °C was considered for the OSL measurements. A pre-heat dose recovery test was conducted for the sample to determine this temperature. Fifteen aliquots were bleached using a daylight simulator for 5 h and subsequently irradiated with 51.6 Gy, previously determined by a dose range test. Five sets of 3 aliquots each were preheated at 5 different temperatures (180 °C, 200 °C, 220 °C, 240 °C and 260 °C). Then, the given/recovered dose ratio was determined for each preheat temperature. The best averaged given/recovered dose ratio was obtained for 220 °C.
OSL measurement and analysis
OSL measurements were conducted using the TL/OSL reader Risø TL-DA 20. This instrument is equipped with a calibrated 90Sr/90Y beta source, which delivers an approximate dose rate of 0.09 Gy/s at the location of the sample disc.
The OSL signal was collected for each measurement from small quartz multi-grain aliquots (30â60 grains) ranging from 1 to 2Â mm in diameter. A total of 48 aliquots were measured to ensure accurate and reliable results (Fig. S8).
Dose rate calculations
The dose rate for the sample was calculated from the mean activity concentrations of 40K, 232Th and 238U, which were measured by high-resolution gamma spectrometry. Conversion factors were used according to Guérin et al.71. The water content considered is indicated in the results (Table 1). It was calculated as the average value between the sample's water content measured at the laboratory and a reference value for the saturation of the sand layers (25%), as the sample was submerged by the effect of the tides for part of the year. The dose attenuation contribution was determined considering this value. The contribution of cosmic radiation to the total dose rate has been calculated as a function of latitude, longitude, altitude, burial depth and mean density of the coating based on data from Prescott and Hutton, 199472. Alpha grain-size attenuation was determined according to Brennan et al.73. According to Guérin et al.71 beta grain-size attenuation was determined. The minimum and maximum etch depths were 8 and 10 microns, respectively. Beta etch depth attenuation factor was used according to Bell74.
Dating results
The age of the sample was determined by applying the following expression:
where the equivalent dose is expressed in Gray (Gy) and the dose rate in Gray per kilo-year (Gy/ka). The age is therefore expressed in kiloyears (ka) (Table 1). Luminescence dating basics can be found in Aitken75.
Dating measurement uncertainty
To estimate the measurement uncertainty, the contributions of the standards, the measurement method, the environmental conditions and those derived from the testing process were considered. The values obtained correspond to the time at which the measurements were made. The expanded measurement uncertainty was expressed with a coverage factor kâ=â2, which, in a normal distribution, corresponds to a coverage probability of approximately 95%. The measurement uncertainty was determined according to JCGM 100:200876.
Data availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].
References
Leakey, M. D. & Hay, R. L. Pliocene footprints in the Laetolil Beds at Laetoli, northern Tanzania. Nature 278, 317â323 (1979).
Masao, F. T. et al. New footprints from Laetoli (Tanzania) provide evidence for marked body size variation in early hominins. eLife 5, e19568 (2016).
McNutt, E. J. et al. Footprint evidence of early hominin locomotor diversity at Laetoli, Tanzania. Nature 600, 468â471 (2021).
Bennett, M. R. & Morse, S. A. Human Footprints: Fossilised Locomotion? (Springer, 2014).
Hatala, K. G., Demes, B. & Richmond, B. G. Laetoli footprints reveal bipedal gait biomechanics different from those of modern humans and chimpanzees. Proc. R. Soc. B 283, 20160235 (2016).
Hatala, K. G. et al. Footprints reveal direct evidence of group behavior and locomotion in Homo erectus. Sci. Rep. 6, 28766 (2016).
Altamura, F. et al. Archaeology and ichnology at Gombore II-2, Melka Kunture, Ethiopia: Everyday life of a mixed-age hominin group 700,000 years ago. Sci. Rep. 8, 2815 (2018).
Bustos, D. et al. Footprints preserve terminal Pleistocene hunt? Humanâsloth interactions in North America. Sci. Adv. 4, eaar7621 (2018).
Duveau, J., Berillon, G., Verna, C., Laisné, G. & Cliquet, D. The composition of a Neandertal social group revealed by the hominin footprints at Le Rozel (Normandy, France). PNAS 116, 19409â19414 (2019).
Mietto, P., Avanzini, M. & Rolandi, G. Human footprints in Pleistocene volcanic ash. Nature 422, 133â133 (2003).
Webb, S., Cupper, M. L. & Robins, R. Pleistocene human footprints from the Willandra Lakes, southeastern Australia. J. Hum. Evol. 50, 405â413 (2006).
Bennett, M. R. et al. Early hominin foot morphology based on 1.5-million-year-old footprints from Ileret, Kenya. Science 323, 1197â1201 (2009).
Ashton, N. et al. Hominin footprints from early Pleistocene deposits at Happisburgh, UK. PLoS ONE 9, e88329 (2014).
Altamura, F. et al. Fossil footprints at the late Lower Paleolithic site of Schöningen (Germany): A new line of research to reconstruct animal and hominin paleoecology. Quat. Sci. Rev. 310, 108094 (2023).
Helm, C. W., Lockley, M. G., Cole, K., Noakes, T. D. & McCrea, R. T. Hominin tracks in southern Africa: A review and an approach to identification. Palaeontol. Afr. 53, 8196 (2019).
Stewart, M. et al. Human footprints provide snapshot of last interglacial ecology in the Arabian interior. Sci. Adv. 6, eaba8940 (2020).
Bennett, M. R. et al. Evidence of humans in North America during the Last Glacial Maximum. Science 373, 1528â1531 (2021).
Mayoral, E. et al. Tracking late Pleistocene Neandertals on the Iberian coast. Sci. Rep. 11, 1â12 (2021).
Zhang, D. D. et al. Earliest parietal art: Hominin hand and foot traces from the middle Pleistocene of Tibet. Sci. Bull. 66, 2506â2515 (2021).
Hatala, K. G. et al. Snapshots of human anatomy, locomotion, and behavior from Late Pleistocene footprints at Engare Sero, Tanzania. Sci. Rep. 10, 7740 (2020).
Roberts, D. & Berger, R. Last interglacial (c. 117 kyr) human footprints from South Africa. S. Afr. J. Sci. 93, 349â350 (1997).
Roberts, D. L. Last interglacial hominid and associated vertebrate fossil trackways in coastal Eolianites, South Africa. Ichnos 15, 190â207 (2008).
Helm, C. W. et al. Dating the Pleistocene hominin ichnosites on South Africaâs Cape south coast. Ichnos 30, 49â68 (2023).
Mafart, B. Ãtude anthropologique dâempreintes de pas humains Holocènes de la sebkra El Azrag (République islamique de Mauritanie). LâAnthropologie 110, 766â787 (2006).
Sahnouni, M. et al. 1.9-million-and 2.4-million-year-old artifacts and stone toolâcutmarked bones from Ain Boucherit, Algeria. Science 362, 1297â1301 (2018).
Hublin, J.-J. et al. New fossils from Jebel Irhoud, Morocco and the pan-African origin of Homo sapiens. Nature 546, 289â292 (2017).
Richter, D. et al. The age of the hominin fossils from Jebel Irhoud, Morocco, and the origins of the middle stone age. Nature 546, 293â296 (2017).
Bouzouggar, A. et al. 90,000 year-old specialised bone technology in the Aterian middle stone age of North Africa. PLoS ONE 13, e0202021 (2018).
Sehasseh, M. et al. Early middle stone age personal ornaments from Bizmoune Cave, Essaouira, Morocco. Sci. Adv. 7, eabi8620 (2021).
Sedrati, M. et al. Using UAV and structure-from-motion photogrammetry for the detection of boulder movement by storms on a rocky shore platform in Laghdira, Northwest Morocco. Remote Sens. 14(16), 4102 (2022).
DâAoût, K., Meert, L., Gheluwe, B. V., Clercq, D. D. & Aerts, P. Experimentally generated footprints in sand: Analysis and consequences for the interpretation of fossil and forensic footprints. Am. J. Phys. Anthropol. 141, 515â525 (2010).
Morse, S. A. et al. Holocene footprints in Namibia: The influence of substrate on footprint variability. Am. J. Phys. Anthropol. 151, 265â279 (2013).
Elftman, H. & Manter, J. Chimpanzee and human feet in bipedal walking. Am. J. Phys. Anthropol. 20, 69â79 (1935).
Alexander, R. M. Principles of Animal Locomotion (Princeton University Press, 2003).
Gray, J. P. & Wolfe, L. D. Height and sexual dimorphism of stature among human societies. Am. J. Phys. Anthropol. 53, 441â456 (1980).
Carretero, J. M. et al. Stature estimation from complete long bones in the Middle Pleistocene humans from the Sima de los Huesos, Sierra de Atapuerca (Spain). J. Hum. Evol. 62, 242â255 (2012).
Holman, R. A. & Bowen, A. J. Bars, bumps and holes: Models for the generation of complex beach topography. J. Geophys. Res. 87, 457â468 (1982).
Walker, R. G. & Plint, A. G. Wave- and storm-dominated shallow marine systems. In Facies Models: Response to Sea Level Change (eds Walker, R. G. & James, N. P.) 219â238 (Geological Association of Canada, 1992).
Marty, D., Strasser, A. & Meyer, C. A. Formation and taphonomy of human footprints in microbial mats of present-day tidal-flat environments: implications for the study of fossil footprints. Ichnos 16, 127â142 (2009).
Thulborn, T. Dinosaur Tracks (Chapman and Hall, 1990).
Roehl, P. Stony Mountain (Ordovician) and Interlake (Silurian) facies analogs of recent low energy marine and subaerial carbonates, Bahamas. AAPG Bull. 51, 1979â2032 (1967).
Monty, C. L. V. Recent algal stromatolitic deposits, Andros Island, Bahamas, preliminary report. Geol. Rundsch. (Int. J. Earth Sci.) 62, 742â783 (1972).
Hardie, L. A. Sedimentation on the modern carbonate tidal flats of northwest Andros Island, Bahamas. Johns Hopkins Univ. Stud. Geol. 22, 95â96 (1977).
Boss, S. K. & Neumann, A. C. Impacts of Hurricane Andrew on carbonate platform environments, northern Great Bahama Bank. Geology 21, 897â900 (1993).
Bourrouilh-Le Jan, F. G. Very high sedimentation (supratidal hurricane deposits) and mid-Holocene highstand on carbonate platforms, Andros, Bahamas: An alternative view. Sediment. Geol. 199, 29â49 (2007).
Pierre, G., Gigout, M. & Rampnoux, J. P. Sur le Quaternaire du Littoral Marocain de Larache à Tétouan et ses Déformations. 2ème Réunion Annuelle des SCIENCES de la Terre. (Société Géologique de France Ed., 1975).
Adil, S. Les Formations Quaternaires de Larache: Sédimentologie, Paléoenvironnement, Datation Radiochimique et Problème Dâaménagement. 3rd. Cycle Thesis, Université Mohammed V Faculté des Sciences, Rabat, Morocco (1996).
El Kadiri, K. et al. Eustatic and tectonic controls on Quaternary Ras Leona marine terraces (Strait of Gibraltar, northern Morocco). Quat. Res. 74, 277â288 (2010).
El Abdellaoui, J. E., Petit, F., Ghaleb, B. & Ozer, A. Sea level fluctuation during MIS 5e and geomorphological context on the southern coast of the Strait of Gibraltar (Morocco). Geomorphol. Relief Process. Environ. 22(3), 287â301 (2016).
Weisrock, A. Niveaux marins du Maroc atlantique durant le dernier Interglaciaire (SIM 5.5, SIM 5.3 et SIM 5.1). Geomorphol. Relief Process. Environ. 22(3), 245â251 (2016).
Chahid, D. et al. Nouvelles données morpho-stratigraphiques et géochronologiques sur le cordon littoral externe (SIM 5-c) de Rabat-Témara, Maroc. Geomorphol. Relief Process. Environ. 22(3), 253â264 (2016).
McBrearty, S. & Brooks, A. S. The revolution that wasnât: A new interpretation of the origin of modern human behavior. J. Hum. Evol. 39(5), 453â563 (2000).
Aumassip, G. Préhistoire du Sahara et de ses abords (LâHarmattan, 2019).
Barton, R. N. E., Bouzouggar, A., Collcutt, S. N., Schwenninger, J. L. & Clark-Balzan, L. OSL dating of the Aterian levels at Dar es-Soltan I (Rabat, Morocco) and implications for the dispersal of modern Homo sapiens. Quat. Sci. Rev. 28(19â20), 1914â1931 (2009).
Schwenninger, J. L. et al. A new Luminescence Chronology for Aterian Cave Sites on the Atlantic Coast of Morocco. South-Eastern Mediterranean Peoples Between 130,000 and 10,000 Years Ago 1â19 (Garcea, 2010).
Ben Arous, E. et al. An improved chronology for the middle stone age at el mnasra cave, Morocco. PLoS ONE 17(2), e0261282 (2022).
Hatala, K. G., Roach, N. T. & Behrensmeyer, A. K. Fossil footprints and what they mean for hominin paleobiology. Evol. Anthropol. 32, 39â53 (2023).
Duveau, J. The Le Rozel footprints: Snapshots of Neandertal groups in the Late Pleistocene. A combined morphometric and experimental approach. Bull. Mem. Soc. Anthropol. Paris 33(2), 53â64 (2021).
Roach, N. T. et al. Pleistocene footprints show intensive use of lake margin habitats by Homo erectus groups. Sci. Rep. 6(1), 1â9 (2016).
Bennett, M. R. et al. Walking in mud: Remarkable Pleistocene human trackways from white sands national park (New Mexico). Quat. Sci. Rev. 249, 106610 (2020).
Duveau, J. From footprint morphometrics to the stature of fossil hominins: A common but uncertain estimate. LâAnthropologie 126, 103067 (2022).
Wiseman, A. L. & De Groote, I. A three-dimensional geometric morphometric study of the effects of erosion on the morphologies of modern and prehistoric footprints. J. Archaeol. Sci. Rep. 17, 93â102 (2018).
Nespoulet, R. et al. Palaeolithic and Neolithic occupations in the Témara region (Rabat, Morocco): Recent data on hominin contexts and behavior. Afr. Archaeol. Rev. 25, 21â39 (2008).
El Hajraoui, M. A., Nespoulet, R., Debénath, A. & Dibble, H. L. La Préhistoire de la région de Rabat-Témara Vol. 3 (Institut National des Sciences de lâArchéologie et du Patrimoine, 2012).
Marean, C. W. The origins and significance of coastal resource use in Africa and Western Eurasia. J. Hum. Evol 77, 17â40 (2014).
Marquer, L. et al. The first use of olives in Africa around 100,000 years ago. Nat. Plants 8(3), 204â208 (2022).
Kim, J. Y., Kim, K. S., Lockley, M. G. & Matthews, N. Hominid ichnotaxonomy: An exploration of a neglected discipline. Ichnos 15, 126â139 (2008).
Dingwall, H. L., Hatala, K. G., Wunderlich, R. E. & Richmond, B. G. Hominin stature, body mass, and walking speed estimates based on 1.5 million-year-old fossil footprints at Ileret, Kenya. J. Hum. Evol. 64, 556â568 (2013).
Murray, A. S. & Wintle, A. G. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiat. Meas. 32, 57â73 (2000).
Galbraith, R. F., Roberts, R. G., Laslett, G. M., Yoshida, H. & Olley, J. M. Optical dating of single and multiple grains of quartz from Jinmium rock shelter, northern Australia: Part I, experimental design and statistical models. Archaeometry 41, 339â364 (1999).
Guérin, G., Mercier, N. & Adamiec, G. Dose-rate conversion factors: Update. Anc. TL 29, 5â8 (2011).
Prescott, J. R. & Hutton, J. T. Cosmic ray contributions to dose rates for luminescence and ESR dating: Large depths and long-term time variations. Radiat. Meas. 23, 497â500 (1994).
Brennan, B. J., Lyons, R. G. & Phillips, S. W. Attenuation of alpha particle track dose for spherical grains. Int. J. Rad. Appl. Instrum. D 18, 249â253 (1991).
Bell, W. T. Attenuation factors for the absorbed radiation dose in quartz inclusions for thermoluminescence dating. Anc. TL 8, 1â12 (1979).
Aitken, M. J. An Introduction to Optical Dating (Oxford University Press, 1998).
BIPM et al. JCGM 100. Evaluation of Measurement DataâGuide to the Expression of Uncertainty in Measurement (2008).
Acknowledgements
The authors are grateful to Iliass Sedrati, Abdelkarim Tadibaght, Mrani-Alaoui Beni for their untiring help and support during the field surveys, and Reda Hawat for operating the drone. We also thank the Provincial Directorate in Larache of the Ministry of Equipment, Transport, Logistics, and Water, as well as the Moroccan Royal Armed Forces, for authorizations and support in the field. This work is part of the project Gof-Boulders âGeomorphological control on offshore boulders transport along a wave-dominated coastâ funded by the University South Brittany and ISblue-the Interdisciplinary graduate school for the blue planet (ANR-17-EURE-0015)âand co-funded by a grant from the French government under the program âInvestissements dâAvenirâ. J.D. research is funded by the FYSSEN foundation. I.D.M. is supported by a Ramón y Cajal fellowship (RYC-2022) and by the Ministry of Science and Innovation of Spain. This article has also received support from Project PID2019-104625RB-100, funded by MCIN/AEI/10.13039/501100011033, by the Andalusian Government to the Research Group RNM276 and by the Centro CientÃfico-Tecnológico de Huelva (CCTH).
Author information
Authors and Affiliations
Contributions
M.S., J.A.M. and J.D. studied the conception and design. M.S., J.A.M., J.D., A.E.M., E.M., I.D.M., E.J.A., G.B., A.S., R.L.G., A.S. and J.R.S. Compiled the data acquisition. Methodology: M.S., J.A.M, J.D., J.R.S.; Analysed data: M.S., J.A.M, J.D., A.E.M., E.M., I.D.M., E.J.A., G.B., A.S., R.L.G., A.S. and J.R.S. Writing-original draft: M.S., J.A.M., J.D., A.E.M., E.M., I.D.M., E.J.A., G.B., A.S., R.L.G., A.S. and J.R.S. Writing-review and editing: M.S., J.A.M. and J.D.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sedrati, M., Morales, J.A., Duveau, J. et al. A Late Pleistocene hominin footprint site on the North African coast of Morocco. Sci Rep 14, 1962 (2024). https://doi.org/10.1038/s41598-024-52344-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-52344-5