Eurasian griffon vulture (Gyps fulvus) as a bone modifying agent and its implications for archaeology

Neo-taphonomic studies have allowed us to detect bone damage patterns linked to carnivore preferences and behavioral traits as well as to improve our understanding of the origin of different alterations on vertebrate fossil faunas. However, taphonomically speaking vultures are among the least studied of all common, obligate scavengers. The research reported here contributes to characterise Eurasian griffon vulture (Gyps fulvus) behavior from a taphonomic perspective describing bone damage on 12 small-sized ungulate carcasses. The combination of observational data from photo/video-trap together with taphonomic analyses allowed us to manage factors like feeding behavior or time of consumption, as well as to accurately record bone modified items. Some bone-modifying effects are described here for the first time as vulture-made bone-damage distinctiveness. Still, some others may pose equifinality problems especially regarding small carnivores. This taphonomic conundrum leaves an interpretation problem particularly in archaeological sites in which those agents are present and consequently, an individualization dilemma about the taphonomic actors involved in bone modified assemblages.

Currently, the menaces facing bird of prey populations around the world are countless and the majority of species are under pressure to a lesser or greater extent 32,33 .In this sense, vultures as species with larger body sizes, low fecundity and high specialization are more likely to be endangered since their populations are more threatened with extinction if external forces affect their demographic parameters increasing their mortality rates 34,35 .This is the case of Old World vultures with 68% of species threatened 32,36 .8][39][40][41][42][43][44][45] ).For instance, forensic taphonomy has also enlarged this inquiry for the benefit of forensic anthropology 46,47 .The effects of vulture modification on human remains and by analogy on pigs (their internal anatomy, fat distribution, general lack of fur, and omnivorous diet are similarly analogous to those of humans) have been documented so as to examine timing and sequence OPEN 1 Institut català de Paleoecologia Humana i Evolució Social (IPHES-CERCA), Zona Educacional 4, Campus Sescelades URV (Edifici W3), 43007, Tarragona, Spain. 2 Departament d'Història i Història de l'Art, Universitat Rovira i Virgili (URV), Avinguda de Catalunya 35, 43002 Tarragona, Spain. 3 Instituto Pirenaico de Ecología (CSIC), Jaca, Huesca, Spain. 4 Instituto de Investigación en Recursos Cinegéticos (IREC), CSIC-UCLM-JCCM, Ciudad Real, Spain. 5Departament d'Economia Aplicada, Universitat de Les Illes Balears, Palma, Mallorca, Spain. 6Models for Information Processing and Fuzzy Information (MOTIBO) Research Group, Health Research Institute of the Balearic Islands, Idisba, 07120 Palma, Mallorca, Spain.* email: rblascolopez@gmail.com;rblasco@iphes.cat of processes such as the rate of skeletonization, disarticulation, and dispersal of bone remains as well as crucial factors being the postmortem interval estimations (PMI) or time since death 26,[48][49][50] .In this manner, vultures emerge in the role of taphonomic indicators that can help to better understand postmortem events and interval estimations at those scavengers' modified scenes.Unfortunately, these studies are primarily focused on soft tissue decomposition and spatial dispersion without regard to bone-modifying effects.
Some preliminary examples of bone alteration by vultures were showed by Cáceres 25 without delving too deep into the subject.The first researcher in describing the taphonomic effects of New World vulture scavenging was Reeves 26 , who documented relatively shallow scratches mostly on the crania and mandibles and linear surface scratches with no depth.Those types of marking were so shallow that intentional cleaning or weather exposure would distort or eliminate their presence.In this study, American black vultures (Coragyps atratus) and turkey vultures (Cathartes aura) waited approximately 24 h before starting to scavenge and between 3 and 27 h of feeding to completely skeletonize the pig carcasses.The main focus of the research was to determining accurate postmortem intervals (PMI) so as to better interpreting taphonomic events and animal scavenging on human remains.The results showed that the raptors would modify the scene, altering and dispersing bones as well as changing the rate of decomposition.In this case, domesticated pigs (Sus scrofa) were used since their internal anatomy, general lack of fur and omnivorous diet, are generally accepted as a human alternative.Taphonomically speaking, the results reported by Domínguez-Solera and Domínguez-Rodrigo 27 were more accurate.The main goal of this study was to show the patterns of bone surface modification according to their frequency and anatomical distribution on a carcass scavenged by griffon vultures (Gyps fulvus).A complete female adult deer carcass was feed on to vultures in order to record their feeding behavior and the resulting bone damage.The results showed considerable bone dispersion while the rib cage, the skull and the limbs remained attached by the skin.As previously reported by Reeves 26 , the disarticulated process was a consequence of the competition to gain access to edible tissues.According to the authors, a diagnostic feature of these raptors is the location of most of the bone marks on the diaphysis of long bones.Almost 22% of the anatomical elements bore marks on their surface, mostly in the form of shallow irregular scores and striae.Several types of marks were described as punctures and circular or oval pits containing crushed bone inside which raised an obvious case of equifinality with mammalian carnivore tooth marks.Trying to enhance the body of research, Fetner and Sołtysiak 28 exposed a bull head to a griffon vulture with the aim of recording the marks on the cranium and mandible, specifically their distribution, length and morphology.The marks were described as shallow single lines although sometimes clustered in groups forming V or L-shape which could eventually disappear due to different diagenetic agents.The mark breadth was less than 1 mm in all cases while 0.3 mm was the maximum depth.The results of all those works show particular marks, mostly in the form of irregular scores and sometimes forming specific shapes.However, all modifications are mostly described as shallow scratches without depth.The mere fact of cleaning the bones to proceed with their analysis or even their exposure to different diagenetic agents could make these types of marks disappear.Some of the features described in these works make such marks prone to being confused with trampling-related effects.Only in Domínguez-Solera and Domínguez-Rodrigo 27 the problem of equifinality with carnivore tooth marks arises although focused on specific modifications such as pits bearing some crushing points.Unfortunately, those works remain insufficient to build the griffon vulture-made pattern either because of the scarce/narrow sample together with the lack of variables considered such as age categories or prey species variety.Here, we present the characterization of scavenger behavior from a neo-taphonomic perspective describing European griffon vulture taphonomic bone-signature on ungulate carcasses in the Catalan Pyrenees (Val d' Aran, Pallars Sobirà and Pallars Jussà).

Study area and subject population
The study was carried out between summer 2022 and winter 2023 in the Pallars Jussà (Pre-Pyrenees, Catalonia).The region displays a diverse range of altitudes, spanning from approximately 500 to nearly 3000 m a.s.l.This variation in altitude contributes to significant climatic differences, resulting in contrasting landscapes and a rich array of vegetation.This region is sparsely populated by humans, who subsist there largely on a mixed economy, based mostly on animal husbandry and tourism.Mediterranean shrubland and mixed woodlands consisting of oak species (Quercus ilex, Q. faginea) and Scots pine (Pinus sylvestris) are the dominant vegetation types in the cultivated lands of the lower valleys, which are situated below 700 m a.s.l.Montane forest zones occur at the middle altitudes (700-1600 m a.s.l.), mainly dominated by red Scots pine, European beech (Fagus sylvatica) and a wide variety of scrublands, open grasslands and pasturelands.In the upper mountain areas (1600-1800 m a.s.l.upwards), mountain pines (Pinus uncinata) are abundant, and above the tree line (2300 m a.s.l.), a mosaic of different types of alpine pastures dominate 51 .The area is inhabited by the four obligate vulture species: Griffon vulture, Bearded vulture, Egyptian vulture (Neophron percnopterus) and Cinereous vulture (Aegypius monachus).Extensive and semi-extensively reared livestock (cattle, sheep and horses) provide most of the biomass for the scavenger guild 52 .Extensive livestock are kept outdoors from April to October, depending on the weather conditions, and some transhumant herds migrate to high summer pastures.Significant populations of herbivorous ungulates, such as red deer (Cervus elaphus), fallow deer (Dama dama), roe deer (Capreolus capreolus), and Pyrenean chamois (Rupicapra pyrenaica), inhabit the area.Of the obligate scavengers within a 30 km radius of the central study area, almost 800 pairs of griffon, 16 pairs of cinereous (Aegypius monachus), 21 pairs of bearded and 50 pairs of Egyptian vultures breed in the area 53 .Regarding facultative scavengers, the region exhibits a noteworthy population of avian (e.g., corvids and raptors such red kite Milvus milvus and black kite Milvus migrans) and mammal species.Additionally, shepherd dogs and hunting dogs that roam freely are frequently observed 37,53 .
The Griffon vulture is a cliff-nesting, socially monogamous colonial species that may nest in large colonies of over 150 pairs 54 .Both sexes incubate the eggs and feed the chicks with an equal division during the breeding season 55 .In the study area, the griffon vulture starts their breeding period in January-February coinciding with egg-laying.The incubation period is about 54-58 days and the chick (only one chick per pair is raised) remains at nest during 1 months (fledging take place in June-August).Both sexes have a similar role in parental roles incubating, and sharing chick rearing tasks proportionally as occurs on other vulture species 56 .Foraging areas are very extensive (hundreds of km) and individuals monitored with GPS transmitters in the study area can visit areas of Extremadura or Alps, as have been documented recently 45,57 .

Vulture-derived sample
The remains of 12 small-sized ungulate carcasses (two Equus caballus, two Bos taurus, six Sus scrofa, one Capreolus capreolus and one Ovis aries) scavenged by vultures (henceforth abbreviated as OB[s], as each of them comes from a different observation) in the region of Pallars Jussà during summer and autumn months of 2022 and winter 2022/23 were analyzed.All carcasses used in our observations weight less than 100 kg and were deposited whole in the selected places.The carcasses were made available to the vultures by: (1) vulture scavenging of natural deaths; (2) road kills/required culling provided by the rangers and transported by ourselves to official bird-watching points.In the first case we were able to document vulture scavenging of carcasses on domestic animals because livestock owners notified rangers' teams of those feeding bouts.Only in one observation (OB1) the carcass was unsupervised for a few hours at dawn so the possibility of some secondary scavengers exists.In the second case we were able to record exclusively vulture consumption since photo-and video trap systems were installed in the vicinity of carcasses.The ungulates were only exposed during daily time, from morning to sunset in order to prevent carnivore intervention at night.In a few cases, vultures did not access the remains the same day they were exposed, therefore, we covered the carcass during the night and uncover it the next morning.Thus, non-vulture scavenging was null and the whole process allowed us to record only vulture bone-damaging.Table 1 summarizes the main features of the neo-taphonomic study organized by individual observations.

Taphonomic analyses
All the skeletal remains were recovered and cleaned to expose bone surfaces using non-aggressive procedure, such as boiling in a mixture of water and non-enzymatic detergent.More resistant tissues were removed by using hair brushes and wooden/soft pointed tools that did not compromise the bone surface.Only ungulate carcasses weighting less than 100 kg were included in this case study; which includes different age groups, from neonate to senile.
Bone damage on the skeletal remains was analyzed using an Olympus SZ 11 stereoscopic (magnification up to 110), together with an ESEM (FEI QUANTA 600 Environmental Scanning Electron Microscope) when necessary.The bone surface modifications were documented combining the criteria defined by Maguire et al. 58 , Haynes 59 and Blumenschine and Selvaggio 60 for mammalian carnivores with those used by Armstrong and Avery 20 and Lloveras et al. 23 for birds of prey-on previous studies of Andrews 61 .In the case at hand, the vulture-induced bone damage was classified as follows: (1) Pits and punctures are marks caused when the tooth/beak (or teeth) of the taphonomic agent penetrates into (pitting) or through (puncturing) the cortical bone surface.(2) Scores are essentially elongated marks usually U-shaped and generated by dragging a tooth/beak along the surface of the bone.www.nature.com/scientificreports/ (3) Notches are angled marks usually without depth cut in the edge or top of the bone, usually recorded on flat bones.(4) Grooves described as a transversal narrow cut or depression in the surface of the bone that are divided here into: (1) general grooves as described right in the above line.(2) V-grooves which are channelshaped notches usually located vertically/transversally on the ends of flat bones such as scapulae, ilium and ischium, mandibular angle (gonium) and, distal shafts of ribs, and especially in non-adult individuals belonging to large size (e.g.E. caballus).Their angular shape begins in the bone cortex and continues through the cancellous tissue with considerably depth and width angles.(3) U-grooves are described like the previous ones and located in the same skeletal elements, portions and individuals, and have been considered as notches in thinner cortical tissue.The difference between the two alterations (V-U) remains in their U-shape.This circular shape begins in the bone cortex and continues vertically/transversally through the trabeculae with considerably depth and relatively wide initial openings.( 5) Crushing is produced by the breakdown of concentrated areas of cortical tissue in the form of short cracks and splits.We add to this study longitudinal crushing as extended longitudinal fractures located on the borders/edges of the thin cortical bone (e.g.spinous/transverse processes of vertebrae and rib shafts) with considerable length but narrow fracture width and conspicuous splits along the modification.( 6) Crenulated edges often occur on flat bones in the form of notched and ragged fractures when the force applied overcomes the density and strength of the bone tissue.( 7) Classic and general peeling usually occurs in bone elements susceptible to fracture by bending, such as ribs or vertebrae, although it can be documented on other skeletal elements, such as pelvis or scapulae.It can be identified by roughened surfaces where cortical layers have been stripped away and where cracks are present at the ends of the fractures with irregular shapes 62 .Incipient peeling includes fraying and bent ends on the cortical layers 63 .(8) Longitudinal cracks are fractures following the lines of bone collagen usually produced by the expansion of a perforation.( 9) Fractures are complete breakages of cortical tissue.They appear as abrupt and clean interruptions of the bone structure, in such a way that the complementary ones can be recovered.(10) Hollow out is described here as an emptying of cancellous tissue similar to scooping out on epiphyseal sections but located in flat bones such as the scapulae (medial border), pelvis (ilium and ischium) and distal shafts of the ribs and mostly recorded in large-sized non-adult individuals.It is caused with the beak by dragging away soft bone tissue in those specific skeletal portions and leaving orifices in the form of conspicuous hollow punctures.

Statistical analyses
We utilized a multinomial model to examine the bone taphonomic modifications, with a focus on predicting outcomes and emphasizing a comparative analysis among three distinct scavenger species namely Gyps fulvus, Ursus arctos and Vulpes vulpes.Although other classification techniques are available, such as regression trees and their derivatives or neural networks, to mention a few, we decided to use multinomial models since they are not black box models, as they allow an easy and clear interpretation of the results.Those models enable, in addition, quantification of changes in the probability of classification in the face of variations in the conditioning variables.That is, we are interested not only in classifying correctly, but also in learning what the determinants of classification are (see, for instance 64 ).Moreover, the accuracy (i.e.out of sample percentage of correct classification) of the estimated multinomial logistic model outperforms that of the regression tree in our empirical case-the sample accuracy for the multinomial logistic model is 81.2% whereas that of the pruned regression tree is 79.6%.Using the tenfold cross validation method, where the sample is randomly divided into 10 parts and in 10 iterations 90% of the sample is used to estimate the model and the remaining 10% is predicted, gives an accuracy of 77.8% for the multinomial model versus 74.4% for the regression tree.Using the Leave-One-Out cross-validation method, the multinomial model gives an accuracy of 77.7% whereas that of the pruned regression trees is 74.1%.Suppose we have a response variable Y that can take on one of K possible outcomes, with K > 2. In our case, we observed the three categories of scavengers, so that K = 3.Also, we observed a p-dimensional vector, X, of characteristics of a damaged bone (number of pits, scores, fractures, type of bone, damaged area, etc.).A multinomial model seeks to determine the probability of Y taking a specific value, given X.This is mathematically denoted as P(Y = k|X).In the case of multinomial logistic regression, one of the most common approaches, one models the log-odds as a linear function of X.In this model, the conditional probability of Y given X is: for k ∈ {1,…,K}, where β k is a vector of parameters conformable with the dimension of X i , i ∈ {1,…,n}, and n is the number of observations.Given that if we make β * k = β k + γ for all k we get the exactly the same probabilities, it is usual to normalize the model taking β 1 = 0 , the Gyps fulvus in our setting.Thus, we may write Then, the log of the odds is linear for k ∈ {2,3}, x ji is the value of variable/characteristic j (j ∈ {1,…,p}) in observation i (i ∈ {1,…,n}), β k,j is its associated coefficient of this variable in category k, and β k,0 is a constant term.For quantitative variables, the value of β k,j should be interpreted as the effect of an increase of one unit in the log-odds, whereas for qualitative variables, which are introduced in the model using dummy-variables, they should be interpreted as the effect on the log-odds from moving from the base category of the qualitative variable to the category indicated by the dummy variable.The base category is no-damage for all of the qualitative variables used in our model except Age and Skeletal group.If we want to compare the log-odds of categories 2 and 3, and we have to look at the difference among the coefficients.The most common way to estimate the parameters of a multinomial model is through the method of maximum likelihood.All the computations and data analysis were performed in R language version 4.3.0(R Core Team 2023; www.R-proje ct.org/) and using package nnet 65 and ggplot2 66 .
Here we present the overview of bone surface modification incidences and anatomical locations.Note the co-occurrence of more than one alteration on the same specimen: 1. Pits/punctures and scores are fairly common in the vulture-created sample with 63 affected-skeletal specimens (9.7%) and the most significant values recorded on the skull (66.7% related to recovered specimens of the same anatomical regions).These are the only alterations documented on long limb bone specimens with the femur being the most affected skeletal element (16.7%) followed by humerus and tibia (4.2% www.nature.com/scientificreports/both).With regard to pits/puncture beak marks we emphasize on the distinction between both marks as a way to individualize/differentiate different taphonomic agents in the case of equifinality problems.
Punctures are recorded on 39 bone elements (6%) whereas pits are documented on seventeen specimens (2.6%).Among the observed beak marks, a percentage of 58.3% cranial specimens are punctured while pits are recorded with values of 16.7% on the same anatomical region.The maxilla (n = 3), the hard palate (n = 3) and the sphenoid (n = 3) exhibit the highest number of punctures among these anatomical regions (Fig. 1d).Hemi-mandible (n = 2), sternum (n = 1) and femur (n = 3) bear only puncture marks on its ramus, body and diaphysis/epiphyses respectively.Six hemi-pelvises are punctured (n = 3) and pitted (n = 3) as well as nine (4.5%) and four (2%) ribs are punctured and pitted respectively.In this case, two pits and three punctures are documented on the costochondral joint whereas the rest (n = 5 punctures and n = 2 pits) occurs predominantly on the distal shafts of the ribs.Twelve vertebrae are punctured (seven thoracic, three lumbar and two cervical) and seven are pitted (four thoracic, two lumbar and one cervical) with most puncturing and pitting damage on the spinous processes and much more rarely on transverse processes.Only one sacral vertebra shows this type of modification (1 pit) on its spinous process and a sternum (1 puncture) on its body.Regarding exclusively pit marks we recorded a total of 0.8% of angular shaped pits (n = 5) located on three ribs, one lumbar vertebra and one hemi-pelvis (Fig. 2c).The size of those beak marks ranges from 1 × 0.7 mm to 5.3 × 5.26 mm with an average of 2.1 × 1.5 mm in the case of pits and 2 × 1.1 mm to 13.1 × 9.1 mm with the average of 5.2 × 3.7 mm regarding puncture marks (Fig. 3).All pits/ puncture beak marks were located on thin cortical tissue with the exception of one puncture recorded on the body of a sternum and three punctures documented on the proximal diaphysis (n = 1), distal diaphysis (n = 1) and distal epiphyses (n = 1) of one femur.2. Scores are detected on 2.3% (n = 15) of the sample with the cranial region showing the highest incidence at 8.3%.These scores are located on both the maxillary body (skull) and ramus (mandible).The appendicular region is also affected with a 4.2% of scoring specifically on one femur, one tibia and one humerus.Two cervical vertebrae and one sacral vertebra bear scores mainly on its transverse process, and the superior border of one scapula is also scored.Eventually, four ribs bear this alteration on the mid-shaft (n = 2), distal shaft (n = 1) and costochondral articulation (n = 1).Most of the skeletal elements exhibit shallow scores with random directions especially on appendicular specimens.Those marks are manly superficial with 2.7 mm and 0.2 mm in the maximum and minimum breadth size respectively and an average of 0.8 mm. 3. Notches are hardly represented with 1.4% (n = 9) of the sample especially on the maxilla of one skull (8.3%).
In relation to the axial skeleton, four ribs display notches on the distal shaft (n = 3) and costochondrale articulation (n = 1); two thoracic vertebrae (4.1%) bear notches on the transverse process and the superior articular facet and two lumbar vertebrae (1.6%) show this type-mark on the border of the spinous process.4. Grooves are moderately represented in the assemblage with a total of 37 (5.7%) modified items.V-groove beak marks are recorded on 21 (3.2%) skeletal elements and are included in the previous general percentage.This subtype is mostly documented on the angle of hemi-mandibles (n = 6) and the ilium of hemipelvises (n = 3) (Figs.1a, f, 2c and 4a).Five ribs contain this alteration basically on the distal shaft, and the transverse (n = 2) and the spinous (n = 3) processes of the vertebrae are also v-grooved.One scapula shows v-grooves on its medial border.The hooked tip's sharp drop of the vulture's beak due to its angular shape and the mobility of the neck would be the cause of these modifications (Fig. 4c).The u-grooves are observed on 14 items (2.2%) that are also included in the general groove percentage.Innominate bones bear the highest percentage (13.6%) of this inflicted mark mainly on its ilium and the angle, and the condyle of two hemi-mandibles are also u-grooved (8.3%).Two scapulae (8.7%) exhibit u-grooves on the medial border and one single rib shows this mark on its distal shaft (Fig. 1c).Regarding the spine, the damage is located on the spinous process, laminae and superior articular process of three thoracic vertebrae and one transverse process, one spinous process and one body of three lumbar vertebrae.This particular shape has to be given by the slight upwards beak's curvature.General grooves are registered on 6 specimens (0.9%).
Concerning the axial skeleton, three lumbar vertebrae (4.3%) are grooved on the transverse (n = 2) and spinous (n = 1) processes and the spinous process of one thoracic vertebra (0.8%) has also this beak mark.One hemi-mandible bears 4 grooves on the angle and one rib is grooved on its costochondrale joint.The dimension of v-groove marks ranges from 1.7 mm × 1.2 mm to 8.3 mm × 6.2 mm with an average range of 5.1 mm × 3.5 mm regarding to v-grooves and from 3. mainly longitudinally crushed (Fig. 2d).Two hemi-mandibles are commonly crushed on its angle and six hemi-pelvises show this modification particularly on the ilium.6. Fractures induced by griffon vultures in the form of crenulated edges are recorded on 13 (2%) specimens.
These modifications are observed on the girdles, with the scapulae being the only altered skeletal element (n = 5; 21.7%) which is affected on its medial border (Fig. 2a), spine and supraspinous and infraspinous fossa.One hemi-mandible (4.2%) is crenulated on its posterior border of ramus regarding to the cranial skeleton.Vertebrae are the only bone specimens crenulated in relation to the axial skeleton, including the sacral and lumbar vertebrae (4.8% and 4.3% respectively), as well as the cervical and thoracic (2% and 1.6% respectively) vertebrae.7.Although being registered mainly in primate produced faunas (e.g. 63) and bear-modified samples 5,10 , we were able to document a relatively modest yet unexpected percentage (3.7%) of peeling damage in the griffon vulture-modified assemblage.Classic peeling is the most common type (n = 10) primarily affecting lumbar vertebrae (n = 5) with substantial damage on the transverse processes (n = 5) (Fig. 1b) and, to a lesser extent, on the spinous processes (n = 3).Additionally, four ribs exhibit peeling damage on their distal shafts (Fig. 1e).General peeling is documented on four lumbar vertebrae specifically on transverse (n = 4) and spinous (n = 2) processes and two ribs bear this type on damage on its head (n = 2) and distal shaft (n = 1).Seven lumbar vertebrae are incipiently peeled (bent end) on the transverse process and two ribs show fraying on the distal and mid-shaft (Fig. 2e).8. Fifteen specimens are longitudinally cracked (2.3%) clustering on axial skeleton and girdles.Regarding the former, thoracic (n = 4) and lumbar (n = 1) vertebrae display this alteration exclusively on the spinous process and seven ribs (3.5%) are cracked on its distal (n = 4) and mid-shaft (n = 3).Concerning the latter, two scapulae (8.7%) show longitudinal cracks on the supraspinous fossa and, the ilium of one (4.5%)hemi-pelvis is also cracked.9. Transverse fractures are recorded mainly on the hyoid bone (n = 2; 50%) and spinous and transverse processes of lumbar vertebrae (n = 29; 41.4%) (Fig. 2b).Thoracic and cervical vertebrae are also fractured on its transverse and spinous processes (16.3% and 14.3% respectively) together with sacral vertebrae (19%).Ten rib specimens are fractured mostly on the distal shaft (n = 12), mid-shaft (n = 3), proximal shaft (n = 3), head (n = 4), neck (n = 1) and tubercle (n = 1).Fractured portions of vertebrae and ribs, especially transverse/spinous processes and distal/mid-shafts are also related to the intra-specific vulture competi- tion in foraging soft tissues quickly.The rest of the fractured specimens include one skull with a split on its styloid process, one breakage on the ilium of a hemi-pelvis and one fractured coronoid process of a hemi-mandible.10.Hollow out is registered on 3.1% of the sample (n = 20).This alteration is mainly observed on the axial skeleton, particularly on the transverse and spinous processes of thoracic (n = 6) and lumbar (n = 2) vertebrae, as well as on the distal shafts of the ribs (n = 8).Additionally, one head and one costochondrale articulation are also hollowed out (Fig. 4b).Regarding girdle elements, scapulae (n = 2) show the highest percentages of this modification recorded on the medial border, and one hemi-pelvis is hollowed out on its ilium (Fig. 1g).One skull (8.3%) bears also this vulture induced-damage on both mastoid processes.
When split by age categories, the highest percentages of bone induced-damage belong to infantile individuals (n = 195; 27.8%) followed by juvenile (n = 34; 14.3%), senile (n = 15; 10.9%) and neonate ungulates (n = 26; 9.4%).Although neonates have the minimum general percentage, its cranial region shows the highest induceddamage with 83.3% of alterations specifically grooves and scores on the mandibular angle and punctures on the eye socket.We recorded 73.7% and 73.1% of vulture inflicted-damage on the axial and girdle areas (n = 157 and n = 19 respectively) of infantile individuals.In relation to the axial region, 105 vertebrae, 51 ribs and one sternum are damaged, as well as ten scapulae and nine hemi-pelvises that are also modified.The cranial area shows a significant degree of bone damage (66.7%) regarding this same age category, and limbs bear 5.1% of scoring mainly.The alterations with the highest percentages of beak marks inflicted on skeletal specimens of this same age category are crushing (10.1%), fractures (9.8%) and pits/punctures/scores (6.6%).Bone damage reaches 42.9% (n = 3) on the girdles and 33.3% (n = 31) on the axial portions of juvenile ungulates whereas no cranial or limb bones are modified.All anatomical regions of senile individuals present moderate vulture activity with the cranial region displaying 33.3% (n = 1) of bone modifications followed by girdles (n = 1; 25%), axial skeleton (n = 12; 18.5%) and limbs (n = 1; 6.3%).Respecting ungulate carcass species, a notable more intensive percentage of modifications is observed on E. caballus (35.1%) and C. capreolus (31.4%) followed by S. scrofa (19.2%), whereas O. aries (10.9%) and B. taurus (9.4%) seem to have slightly less intensive damage.There seems to be a deeper scavenging on summer (35.1%) and winter (24%) than on autumn (15.7%), although when split by skeletal parts we found the axial region (n = 27; 96.4%) to be the most exploited area in winter followed by girdles (n = 6; 85.7%) in the same year period.We found a general decrease in scavenging activity during autumn for all skeletal regions, including cranial (45.8%), girdle (38.7%), axial (32.8%) and limb bones (0.8%).
The multinomial logistic regression model was applied to the complete dataset consisting of 1782 observations which include data from neotaphonomic studies on brown bears 5 and red foxes 6 .Table 4 presents the maximum likelihood estimates derived from this analysis, taking the griffon vulture as the reference category (i.e.Y i = 1 ), and 9 quantitative conditioning variables (number of pits, scores, fracture, etc.) as well as 43 qualitative dichotomous ones (categories of age, skeletal group, and damaged areas).It also shows the standard errors and the associated p-value of each estimate.
By emphasizing the significance of p-values in identifying pertinent variables for classification purposes, we observe that within the quantitative variables, the number of pits (N.Pits) plays a crucial role in distinguishing between Gyps fulvus and Vulpes vulpes exerting notably a positive effect that favors the latter scavenger.However, the number of scores (N.Scores) has a highly favorable effect for griffon vulture dismissing both, Ursus arctos and Vulpes vulpes.There is a slight negative effect observed (p=0.068) when considering the impact of fractures on the classification of Vulpes vulpes but the same bone alteration does have a relatively significant effect on the classification of Ursus arctos.Conversely, the number of crushing exhibits a favorable effect on the classification of that large carnivore.The occurrence of longitudinal crushing or bent ends provides additional evidence supporting the classification of Gyps fulvus while an increase in the number of general peeling reduces the probability of Vulpes vulpes being classified as the scavenger.Finally, when considering the quantitative variables, the amount of furrowing provides supporting evidence for both Ursus arctos and Vulpes vulpes.
In terms of the qualitative variables, the presence of classic peeling provides some evidence favoring the classification of Ursus arctos.On the other hand, the occurrence of fraying decreases the probability of a bone being damaged by Gyps fulvus.When considering the skeletal group, the inclusion of cervical and lumbar vertebrae along with ribs contributes to an increased probability of Vulpes vulpes being identified as the scavenger.Moreover, thoracic vertebrae give strong evidence in favor of this last small carnivore.The same skeletal element provides supporting evidence for the classification of Ursus arctos, however, the cranial skeleton and girdles exhibit a relatively low probability, suggesting they are less indicative of this large carnivore as the main scavenger.Lastly, the presence of small bones strengthens the evidence against Gyps fulvus as the main taphonomic agent.In relation to the dummy variables that are associated with the damaged areas of the bones, these simply collect changes in the starting point of the probability curve associated with the number of damages, that is, the case in which the bone is not damaged.For this reason, its interpretation is clearer when its effect is included together with the quantitative variable.Figures 5 and 6 show how the probability of classifying the three species of scavengers varies as the number of pits or scores increases in relation to the different bone damaged areas (dist shaft and other).We considered the bone being a rib and the age category being inf (infantile).Based on observations, in the absence of pits, Gyps fulvus is the most likely scavenger, followed by Ursus arctos, while the probability of Vulpes vulpes is very low.However, as the number of pits increases, the probability of Vulpes vulpes steadily rises and eventually approaches 1 for a large number of pits.In contrast, when examining the occurrence of scores, a small number of such types of modification are associated with the highest probability being attributed to Ursus arctos, followed by Vulpes vulpes.Nevertheless, as the number of scores increases, the probability of Gyps fulvus rapidly rises and approaches 1 when the number of scores exceeds 8.

Discussion
Studies on obligate scavengers depict their taphonomic signature grounded primarily on shallow bone damage, mainly scores and striae, with the majority of bone surface modifications located on the diaphysis of long limb bones.Slight evidence such as irregularly shaped shallow scratches or linear surface striae without depth has also been documented principally on the crania and mandibles 26,27 together with superficial single lines sometimes grouped and forming V or L-shapes 28 .However, the controversy emerges when describing a species-taphonomic hallmark based on a limited sample and, by extension, a limited data set.Our analysis results on 12 small-sized ungulate carcasses scavenged by griffon vultures point to a dissimilar scenario.
On the topic of bone induced-damage, ten different types of marks have been described and recorded in the present study without counting on subtypes and the well-founded division between pits and punctures as distinctiveness of vulture-made bone modification due to the help the latter split can provide when deciphering equifinality processes.As opposed to the low number of conspicuous pits and punctures that have been previously registered, our research shows an elevated frequency of these beak-marks ranging the third percentage-position only preceded by crushing and transversal fractures with a pit: puncture ratio of 17:39 (see Table 2).The location, shape and dimensions of isolated pits could mislead taphonomists on interpreting some specimens given that those marks located on the same anatomical areas are conventional in carnivore bone-altered assemblages.Therefore, the analysis of solely punctures and their combination with other type marks would help to individualize the responsible agency.
Some bone-unique made effects are listed here as v-grooves, u-grooves, longitudinal crushing and hollow out.Correspondingly, common grooves should be associated to v-grooves as a way to identify alleged scavenger especially on axial and girdles regions.Infantile individuals bear the highest percentage of those inflicted marks particularly on those regions together with the cranial skeleton of neonate specimens.The greater elasticity and less hardness of the compact bone tissue of these age groups compared to the pertinent in adult ones seem to facilitate this fact 67 .Additional singularities that would segregate small carnivores are, on the one hand, puncture, notches and groove-marks together with hollow out located on the costocondrale joint and recorded as individual damage.The precise forepaws and the elastic neck of the vulture family together with most of the vertebral column of the carcass articulated would be needed to access these particular points.Therefore, obvious physical reasons, for instance cranial shape and dentition, would prevent carnivores to inflict these peculiar marks neither on articulated thoracic cage nor on single individual ribs.On the other hand we registered pits in the medial area of the spinous process especially on thoracic vertebrae also retrieved from anatomical connected vertebral columns.The movements of the head together with the shape of the beak allow vultures to reach remote hidden small pieces of edible tissue and, unintentionally stamp their hallmark.Features such as this one and the fact that griffon vultures feed essentially on the softer viscera causes them to be called as gulpers 68 .Likewise, longitudinal crushing appears in remarkable percentages.A ratio of 4:7 (crushing: longitudinal crushing) is documented on the border of the ribs and a significant percentage is also documented on thoracic and lumbar vertebrae.We associate most longitudinal crushing with the eagerness of pecking soft tissue embedded between distal and mid-shafts of the ribs and vertebral transverse processes.From this perspective, the closeness towards the carcass seems to have an age-hierarchical character in the feeding bouts 69 which leaves sub-adult, immature and juvenile griffon vultures competing for the last scraps and generating most of the previous inflicted bone-damage.The morphology of the complex head-neck system of these obligate scavengers performs in a variety of adaptations regarding its feeding strategy.Related to this, we were able to recover a stillborn calf (OB12) hollowed out from a skin hole in the lower abdomen with the solely result of a fractured hyoid-bone which illustrates the easiness of their long/elastic neck to reach remote and concealed scraps of soft tissue without hindrance.This particular morphologic feature would be associated in some cases with beak marks on both dorsal and ventral sides of skeletal elements such as ribs and transverse processes of vertebrae.
Nevertheless, the inference regarding the identification of the taphonomic agent turns out to be unambiguous when analyzing isolated incidences of cause-and-effect in some instances.As pointed out by Domínguez-Solera and Domínguez-Rodrigo 27 the risk of potential equifinality problems increases as we widen the referential frameworks.In our case we identified some related-points such as pits/punctures together with scores on specimens of the axial skeleton and transversal fractures especially on the same skeletal region.Most of those alterations both, isolated or in association with different bone-type marks, have been recorded in small-sized ungulates scavenged by carnivores 6,7 .
In relation to the first point, as mentioned in the results section, we have observed that the distinguishing of pits and punctures located in certain skeletal regions can be beneficial when determining particular taphonomic agents.Moreover, and regarding pits/scores based on the quantitative variables, the multinomial model reveals that the number of pits plays a significant role in differentiating between griffon vulture and red fox, exerting a noticeable positive effect that favors the red fox as the main scavenger.On the other hand, the number of scores has a highly favorable effect for Gyps fulvus, while dismissing both Ursus arctos and Vulpes vulpes.Concerning the second aspect, we can observe a detrimental impact on the fox classification, with a higher probability of it with the brown bear and the vultures being the taphonomic agents as the number of fractures increases.There is currently insufficient evidence to differentiate between the last two scavengers, but expanding the dataset will facilitate obtaining more accurate outcomes.

Conclusions
Previous studies show that griffon vultures reveal certain patterns primarily characterized by shallow damage to the bones, consisting mainly of scores and grooves located on the shafts of limb bones.Additionally, there is some evidence of irregularly shaped superficial scratches or linear grooves on the skulls and jaws, often forming V or L-shapes, along with superficial single lines that can occur individually or grouped together.The results we Figure 1.(a) V-groove on the mandibular angle of a calf (OB2); (b) Classic peeling on the transverse process of a lumbar vertebra (OB3); (c) U-groove on the distal shaft of a rib (OB10); (d) Punctures on the maxilla of a wild boar (OB14); (e) Co-ocurrence of a pit and classic peeling (incipient) on the distal shaft of a rib (OB10) (f) V-groove on the ilium of a pelvis (OB6); (g) Hollow out on the medial border of a scapula (OB10).Scale: 1 cm.

Figure 2 .
Figure 2. (a) Crenulated edge on the medial border of a scapula (OB14); (b) Fractures on the spinous and transverse processes and puncture on the superior articular facet of a vertebra (OB14); (c) V-groove associated to an angular pit on the ilium region of a pelvis (OB13); (d) Longitudinal crushing on a rib of a roe deer (OB6); (e) Fraying on a rib of a calf (OB2).Scale: 1 cm.

Figure 3 .
Figure 3. Box plots show the measurements of pits/punctures (width and length), and the dimensions of Uand V-grooves (depth and breadth) occurring on thin cortical bone.

Figure 5 .
Figure 5.A graphical representation showing the likelihood distribution of pit quantities on ribs among the three distinct scavenger species.

Figure 6 .
Figure 6.Plot illustrating the likelihood distribution of the number of scores on ribs among the three distinct scavenger species.

Table 1 .
General characteristics of the experimental series organized by individual observations (OB).*Totaltime from initial scavenging to final removal.OB

Table 2 .
Recovered bones vs damaged bones by skeletal elements split by different bone damage type.Only the skeletal categories with specimens bearing inflicted-damage are registered in the table.Note that some specimens show co-occurrence of damage and therefore, the total number (summing up all types) can be higher in some categories.Rcv Bn, Recovered bones; Dmg Bn, Damaged bones; Pt, Pct, Sc, Pits, Punctures, Scores; Ntc, Notches; Grv, Grooves; TFrac, Transversal Fractures; Crs, Crushing; LCracks, Longitudinal cracks; Hollow out, Hllw; Pl, Peeling; CrEdg, Crenulated edges.