Late Neolithic Vinča communities, spread over much of central and northern Balkans during the late sixth to mid-fifth millennium BC and characterised by unusually large and densely population centres, would have required highly organised food production systems. Zooarchaeological analysis indicates that domesticate livestock were herded, but little is known about the seasonal husbandry practices that helped ensure a steady supply of animal products to Vinča farming communities. Here, we present new stable carbon (δ13C) and oxygen (δ18O) isotopic measurements of incremental bioapatite samples from the teeth of domesticated livestock and wild herbivore teeth from two late Neolithic Vinča culture sites: Vinča-Belo brdo and Stubline (Serbia). Our results show a low variation overall within sheep and goats in terms of pasture type that may have been composed of seasonal halophyte plant communities, which have higher δ13C values due to the saline rich growing environments. Cattle feeding strategies were more variable and provided with supplementary forage, such as cut branches or leafy hay, during winter. The sharp distinction in the management of cattle and sheep/goat may be associated with the development of herding strategies that sought to balance livestock feeding behaviours with available forage or, more provocatively, the emergence of household-based control over cattle–an animal that held a central economic and symbolic role in Vinča societies.
Citation: Gillis RE, Bulatović J, Penezić K, Spasić M, Tasić NN, Makarewicz CA (2021) Of herds and societies—Seasonal aspects of Vinča culture herding and land use practices revealed using sequential stable isotope analysis of animal teeth. PLoS ONE 16(10): e0258230. https://doi.org/10.1371/journal.pone.0258230
Editor: Peter F. Biehl, University at Buffalo - The State University of New York, UNITED STATES
Received: June 2, 2021; Accepted: September 21, 2021; Published: October 7, 2021
Copyright: © 2021 Gillis et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: The analysis of the material was undertaken as a part of REG’s post-doc funded by the Deutsche Forschungsgemeinschaft (DFG, https://www.dfg.de/) at the Graduate School, Kiel (GSC 208). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
During the Late Neolithic (5300/5200 to 4600/4500 cal BC), much of the central Balkans and southeast parts of the Vojvodina and Transylvania regions were occupied by Vinča communities which produced dark burnished ceramics and figurines often bearing unique triangular faces [1,2]. Compared with earlier Starčevo-Körös-Çriş cultural complex, Vinča groups occupied a diverse range of site types, including tells, with long occupation sequences . Excavations of settlements have revealed extensive organisation of social space, with designated central places and well-defined lay-out of houses including regular spaced thoroughfares [4–7]. Population estimates based on the number of ‘house’ structures in Vinča culture settlements suggests these sites were densely populated with up to several hundred inhabitants (c.f. Porčić ). The increase in settlements distribution and size suggests a population increase that was supported by herding and agricultural activities. Vinča settlements were often situated on the alluvial plain and lower terraces of the Danube, Sava, Morava, and their tributaries, and subsistence systems shaped by periodic flooding of the landscape which would have provided rich seasonally replenished alluvial deposits ideal for cropping production . Crop growing may have been limited to small plots that were intensively cultivated as suggested by unpublished analysis of weed seed assemblages (Filipović pers. comm). Heavier soils of the wooded hill slopes were probably unlikely to have been cultivated.
Vinča communities exploited a wide spectrum of domesticated plants including einkorn (Triticum monococcum), emmer (Triticum dicoccum) and barley (Hordeum vulgare, hulled and naked), as well as pulses such as lentils (Lens culinaris) and pea (Pisum sativum) . Domesticated cattle (Bos taurus) and pigs (Sus scrofa) increased in socio-economic importance over time within Vinča communities, likely due to intensification in meat and dairy production associated with a diachronic increase in settlement size and populations . The high proportion of adult cattle identified at the late Vinča sites, Divostin and Gomolava, for example, may indicate dairying and use of animal traction [12,13], while at other sites, cattle appear to have mainly exploited for their meat indicated by the selective slaughter of sub-adult animals (18–36 months) [11,14–19]. Provisioning Vinča communities with an adequate supply of animal products would have required a high level of animal maintenance to ensure herd growth and productivity, in particular a sufficient source of good pasture and forage. Vinča herders may have employed at various times within the seasonal cycle complementary strategies that served to further provision their hungry herds by leading animals to crop stubble, or provisioning them with cut branches or leafy hay  collected from woodlands as has been previously hypothesised for Vinča domesticated pigs and cattle [21,22]. Waterlogged pastures caused by seasonal river flooding may have encouraged the development of mobile herding strategies where animals were moved to different, drier pastures within the seasonal calendar [11,21].
Carbon and nitrogen stable isotopes can be used to investigate husbandry practices, including forage and pasture management strategies [23–27]. The δ13C signatures of animal tissues reflect the carbon component of protein (bone) and carbohydrate (fats) intake, and total diet (teeth) and, for herbivores, tissue carbon isotope ratios directly reflect plant stable isotopic composition . Thus can be used to reconstruct pasture and fodder resource environments [23–25]. The nitrogen isotopic composition of bone and dentine collagen can indicate the level of anthropogenic activity associated with their food source [22,23,26,27]. Previous stable isotopic analyses of bone collagen from faunal material recovered from Stubline and Vinča-Belo brdo suggest that different types of forage were provisioned to different domesticated species, with cattle exhibiting a wide range of δ13C values indicative of open and closed forest pasturing . However, stable isotope values measured from bone collagen represent an integrated signal representing a long-term ‘dietary average’ and, due to constant bone remodelling, do not capture seasonal-scale shifts in dietary intake. Sequential sampling of tooth enamel captures at high-resolutions carbon and oxygen isotopic change especially in species bearing hypsodont teeth [29,30]. At mid- and high- latitudes, the stable oxygen isotope ratios of precipitation are strongly correlated to mean annual ambient temperature with low δ18O values expressed during the winter months and higher values during the summer [29,31]. Consequently, oxygen isotopes measured from incrementally sampled teeth provide a seasonal-scale temporal framework that assists in the interpretation of δ13C values that may be related to seasonal variation in pasture and fodder sources.
We present isotopic results measured from the bioapatite fraction of incrementally sampled teeth from domesticated and wild herbivores to further explore the Late Vinča fodder and pasture management strategies. These results are further contextualized within published stable isotopic results from bone collagen samples .
Cultural and climatic context
The epicentre of the Vinča culture was located in Serbia and dispersed throughout the river valleys of the northern and central Balkans, including southern–most Hungary and parts of the Transylvania. Vinča culture communities also occupied hilltop settlements including, for example, Gradac, Valač, Pljosna Stijena, Gradina-Stapari , while numerous settlement located in wetland and marshland landscapes have been identified in western Serbia. Cave sites are also known, including the sites of Petnica in western Serbia and Nandru Vale in Romania :45]. As part of this study and Gillis, Bulatović , we have focused on faunal material recovered from assemblages that date to the terminal phase of Vinča culture, i.e. Vinčа D/Vinča–Pločnik II phase (between 4850 and 4650/4600 cal BC [2,34]). On the broader regional scale, this phase is characterised by several major shifts in settlement organisation, material culture production, such as new vessel forms [35–39] and standardisation of ceramic production processes , and specialised animal economies. The distribution of similar material culture objects across the Vinča macro-region indicates the establishment of large-scale social networks within existing inter-regional trade/exchange routes [1,41]. The late Vinča phase also witnessed the widespread manufacture of copper items within communities located in the eastern and southern parts of Serbia [42–44], as well as the long-distance transmission of finished copper objects which were distributed over the almost entire area of the Vinča cultural sphere [2,45]. Overall, late Vinča settlements appear to be much larger in size and more complex in organisation than the earlier ones, consisting of up to several hundred elaborately built and equipped houses [5,46,47], complex of ditches surrounding the settlement [1,5], and supporting several thousand inhabitants .
We analysed material from two settlements, Vinča-Belo brdo and Stubline, both located on the alluvial plains and lower terraces of the Danube and Sava respectively, and their tributaries in northern Serbia (Fig 1). The prominent tell site of Vinča-Belo brdo stands several meters high on the southern bank of the Danube. The site has suffered erosion due to the Danube and covers an area of approx. 12 ha . The Neolithic Vinča-Belo Brdo sequence represents a long occupational history , established by Starčevo groups (c. 5700–5300 cal BC) and followed by Vinča culture groups until c. 4550 cal BC when occupation ceased [51,52]. Recent excavations have focused on the final Vinča occupation, represented by two successive horizons of burnt rectangular buildings, interspaced with a horizon of unburnt buildings, oriented along a north-west axis, composed of two to three rooms, and housing at least one kiln per structure, with some buildings containing several kilns . During the Vinča D period, the external space between houses was up to 2 m in places. Several communication lines (“streets”) were placed between houses in some areas, with no open spaces within settlements have been identified, suggesting that animal pens, gardens or cultivation areas were external to settlements during the late Vinča phases. Next to Vinča-Belo brdo is the Bolečica river, a tributary of the Danube, which annually floods part of the landscape surrounding the site . Damp alluvial soils in the river floodplain would have supported gallery forests and wetland vegetation, such as willow sp. (Salix sp.) and reeds (Phragmites communis). Thermophilous taxa, particularly oak, grew on rich chernozem soils and dominated forests that spread over hilltops and slopes .
1. Vinča-Belo brdo; 2. Stubline; 3. Opovo; 4. Gomolava; 5. Divostin; 6. Petnica; 7. Valač; 8. Gradac; 9. Gradina-Stapari; 10. Pljosna Stijena. The base map is available under the CC4.0 at the U.S. Geological Survey. (Credit: Penezić/Pendić).
The site of Crkvine at Stubline (herein Stubline) is located circa 50 km east of Vinča–Belo Brdo (Fig 1). In a broader geographical context, Stubline lies in the micro-region called Drenski Vis, located in the Posavina-Tamnava region. Approximately 12.5 ha in size, Stubline is situated on an elongated plateau of about 23.6 ha (100 m a.s.l.) [5,6,46,53,54]. The plateau was surrounded by two watercourses to the north and south, that meet at the break of slope in the south-western part of the settlement. Based on the geophysical survey, the earliest settlement consisted of about 120 houses surrounded by a system of double ditches irregular in plan. Over time, the Vinča community grew to occupy almost the entire eastern portion of the plateau and may have consisted of more than two hundred houses . At least two occupation phases have been so far identified at Stubline: Stubline Ia (Vinča D1 phase in the eastern part of the plateau) and Stubline Ib (Vinča D2 phase throughout the whole plateau) . A relative chronology based on material culture indicates that the Vinča culture occupation spanned a period of around 200 to 250 years (4850/4800–4650/4600 BC) .
The climate today is characterised as humid sub-tropical (Cfa) according to the Köppen climate classification . The warm season lasts on average for 3–4 months beginning in May with the average temperature reaching 24°C while the cold season last for on average for 3 months with temperatures falling to -2°C. Precipitation falls mainly in late spring and late autumn (c. 70mm), with most of the precipitation falling as snow (20-8mm) during the cold season . The climate during the occupation of the sites was wetter and warmer in the winter than today, with reduced variation between seasons .
Material and methods
All animal bone and teeth from Vinča-Belo brdo and Stubline were sampled by REG, assisted by JB and KP, with the permission of the local authorities and excavators (MS, NT). The collections are housed in City Museum, Belgrade. Details of the sampling protocol for wild and domesticated faunal bones for stable isotopic analysis of bone collagen are outlined in Gillis, Bulatović . The complete sample list is provided in S1 Table. The Vinča-Belo brdo herbivore tooth specimens were recovered from archaeological features (e.g., pit-holes, demolition layers) and cultural layers dated to the latest occupation phase (Vinča D2 cultural phase; 4700–4550/4500 cal BC . Specimens from Stubline were from archaeological features (e.g. waste pits, ditches) and cultural layers dated to the Vinča D period.
To examine the seasonality of husbandry practices, we selected teeth from cattle, sheep/goat, and domesticated pig, based on state of occlusal wear and overall specimen preservation. The third mandibular molars of domestic cattle (Bos taurus; Vinča-Belo brdo n = 5, Stubline n = 1) and sheep/goats (Ovis aries/Capra hircus; Vinča-Belo brdo n = 4, Stubline n = 5) were selected for analysis. The third molar in both species begins development around 10 months and erupts around 24 months [58,59]. We focused on teeth with dental wear up to stage G for both sheep/goats and cattle , in order to increase the chances of capturing a full seasonal cycle in the isotopic record. The isotopic composition of the third molar reflects adult dietary intake only and is free of ontogenic effects, including suckling and weaning, that influence the carbon and oxygen isotopes of tooth bioapatite.
For domesticated pigs (Sus scrofa), 3 individuals from both sites were sampled with both the mandibular M2 and M3 selected. The bunodont teeth of pigs display short crowns in comparison to ruminants. Enamel mineralisation of the lower M2 in modern pigs is initiated at around 3 months and finishes around 11 months, while mineralization of the mandibular M3 begins at approximately 11 months and extends through ca. 18 months [61,62]. Pigs are usually weaned by the second month of life , and therefore we do not expect the carbon and oxygen isotopic composition of tooth bioapatite to be affected by a suckling signal.
To provide a reference for a ruminant grazing/browsing within a forest environment we sampled two red deer (Cervus elaphus) individuals (7237, 7238) from Vinča-Belo brdo, with two teeth from a single individual (7237-M2/M3). Red deer mandibular M2 mineralisation begins around 4 months, around when weaning occurs, with the crown complete at 9 months. For the M3, mineralisation begins around 13 months with completion of the crown at 26 months based on a study on modern red deer from Richmond Park, London .
Stable isotope sampling and analysis methodology
Prior to sequential enamel sampling, the tooth surface was cleaned using an abrasive drill bit. A sequence of enamel samples was then drilled on the buccal surface from the enamel-root junction (ERJ) to the apex of the crown. Purification protocol for stable isotopic analysis followed . Purified bioapatite samples of domesticated samples were weighed between 550–650μg were analysed on a Kiel IV device interfaced to a Delta V Advantage IRMS at SSMIM platform, Paris. The analytical precision estimated from average of the internal standards (Marbre LM normalized to NBS 19: between four and eight per run) for δ13C analyses was 0.02‰ (range: 0.009–0.05‰) and for those of δ18O was 0.04‰ (range: 0.03–0.06‰). The red deer samples were analysed at the Leibniz labor at Kiel University (Germany). The analytical precision was for δ13C analyses was 0.02 (range: 0.01 to 0.03) and δ18O was 0.03‰ (range: 0.01 to 0.04).
Stable isotopic enrichment between diet and bioapatite
Studies concerning the enrichment of δ13C values between diet and bioapatite have shown that it can vary between species depending on dietary composition, physiology, and body size [30,66–69]. The magnitude of the enrichment depends on the plant composition (C3 or C4 plants) of an animals’ diet [67–69], and intra-tooth variation in the enrichment factor has been observed where the diet composition changes during the tooth’s development . We have chosen to use a single enrichment value has been used for all species of 14.1‰  as applied previously in Balasse, Bălăşescu  for the analysis of Romanian Late Neolithic and Chalcolithic settlements.
Correction of modern stable isotopic data
The burning of fossil fuels and to a lesser extent, the removal of large stretches of forests, has contributed to changes in the carbon isotope composition of atmospheric CO2. We applied a correction factor of 1.7‰ to carbon isotope data obtained from modern tissues. This calculation is based on the δ13C value of CO2 at present (−8‰) and the given studied period. The δ13C value of CO2 during the Late Vinča period was −6.3‰ based on interpolated data from Antarctic ice-cores .
The summary statistics used were the average midpoint value, range, and amplitude per species. The midpoint was calculated (δxxXMax+ δxxXMin)/2, and then an average and standard deviation was calculated per species. The amplitude was calculated by subtracting the δxxXMin from the δxxXMax, an average was then calculated per species. The maximum and minimum value per species is true range and not an average.
Stable oxygen isotope results
The δ18O values of Vinča-Belo brdo cattle average mid-point was −7.2‰ (range: −9.6 to −4.2‰) (Table 1). In comparison, the single Stubline cattle averaged −8.5‰ (range: −10.3 to −6.8‰). Vinča-Belo brdo sheep/goat specimens exhibited an average of −6.1‰ (range: −10.5 to −1.1‰) like those from Stubline (average −6.5‰, range: −10.6‰ to −2.4‰). The δ18O values exhibited by domesticated pigs from both sites fall within the range of −9.1‰ to −6.3‰, with the greater variation evident in Vinča-Belo brdo pigs, particularly the M2s (Table 1). Red deer δ18O values ranged from −9.0‰ to −4.6‰, which is comparable to variation exhibited in cattle teeth.
Stable carbon isotope results
The δ13C values from bioapatite samples of domesticated animals ranged from −14.1‰ and −7.9‰ with the samples from Vinča-Belo brdo cattle exhibiting the greatest range in values observed from all species (Table 2). Cattle from Vinča-Belo brdo had an average mid-point carbon isotope value enriched +1.5‰ in 13C relative to other domesticated species. Values from a single cattle tooth from Stubline fell within the range exhibited in cattle from Vinča-Belo brdo. The δ13C values for Vinča-Belo brdo sheep/goat teeth range between −14.1‰ and −10.1‰, with the average midpoint value of −12.0‰, while at Stubline, sheep/goat teeth exhibited a slightly larger range in values (−14.5‰ to −8.9‰). The M2s from domestic pigs at Vinča-Belo brdo exhibit values that were on average 0.3‰ higher midpoints values than M3s. At Stubline, there was little difference between the average midpoints from M2s and M3s of domesticated pigs (M2: −12.5‰; M3: −12.6‰).
Results for individual teeth are presented in S3 Table.
The two red deer individuals exhibited a large range in values (−15.6‰ to −10.9‰; Table 2). The second molar from red deer 7237 yielded δ13C values ranging from −15.6‰ to −10.9‰, with the lowest midpoint of the entire dataset of −13.3‰ while the M3s from the same individual and 7238, exhibited a smaller range in values (−13.7‰ to −11.9‰) and an average midpoint of −12.7‰.
Amplitude of intra-tooth carbon isotope change provides insight into isotopic variation in forage ingested by animals. The average amplitude for cattle from Vinča-Belo brdo was 1.7‰ (range: 0.7 to 2.4‰) with some individuals (Bos 7240 and 7241) exhibiting an amplitude of intra-tooth isotopic change of less than 1‰ to 2.4‰ (Bos 7230, Fig 2). At Vinča-Belo brdo, sheep/goat teeth exhibit the largest amplitude in intra-tooth isotopic change, averaging 2.5‰ in Δ13C (range: 1.4 to 3.3‰). In comparison, the single goat identified in the Vinča-Belo brdo assemblage (O/C 7245) exhibited a low amplitude of intra-tooth change (1.4‰). Two sheep/goat individuals (7267 and 7287) had the largest intro-tooth amplitude of 4.7‰ and 4.9‰ respectively. Overall, sheep/goat samples from Stubline exhibit the largest average amplitude of change for domesticated animals in the dataset (3.4‰). Pig intra-tooth amplitude is greatest at Vinča-Belo Brdo (0.7–1.3‰) compared to those from Stubline (0.3–0.4‰). Red deer M3s exhibit an average amplitude of 1.6‰, while the single M2 has the largest amplitude in the dataset of 4.6‰.
Combined incremental results
For cattle, there is considerable variation in the relationship between δ18O and δ13C values. Most individuals follow a pattern where δ13C values increase during the summer months (Fig 3B and 3D–3F). BOS 7241 and the single tooth from Stubline (BOS 7268) exhibits little variation over the seasonal cycle (Fig 3C and 3F). While BOS 7244 and 7230, the δ13C values remain high during winter when δ18O values decrease. The red deer teeth also show a similar decrease in δ13C values during the winter; this is particularly evident in 7237-M2. In contrast, carbon isotopic values increase in 7237-M3 when δ18O values decrease, i.e., during winter (Fig 3G–3I). For domesticated pigs analysed, little variation in δ13C values was observed reflecting in part the bunodont tooth geometry.
Insets: a) Bos 7230, b) Bos 7240, c) Bos 7241, d) Bos 7242, e) Bos 7244 (Vinča-Belo brdo), f) Bos 7268 (Stubline), g) Cervus 7237-M2, h) Cervus 7237-M3 and i) Cervus 7238-M3. ERJ is the enamel root junction.
Within our study, sheep/goat individuals from both sites exhibited a strong correspondence between δ18O and δ13C values, high δ13C values during summer and vice a versa (Fig 4) except for the single goat (O/C 7245, Fig 4A). In this individual, low δ13C values (> −11.9‰) coincide with high summer season δ18O values.
Previous stable carbon and nitrogen isotopic analysis of bone collagen samples from the Late Vinča phases (4800 to 4600 BC) at Vinča-Belo brdo and Stubline pasture and foddering practices for sheep, goat and cattle . The large range of δ13C values from bone collagen from domestic cattle at Vinča-Belo brdo suggests they may have been grazed across different environments (water meadows, open grassland, open/closed forest). While the relative low bone collagen δ13C values at Stubline may indicate that they were kept mainly within wooded or waterlogged environments. At both sites, caprines and domestic pigs consumed a diet enriched in 15N, which strongly suggests these animals were raised on either crop residue from cultivars raised on manured soils or pastures repeatedly visited by livestock . We previously proposed this supports the hypothesis that caprines herds and pigs were fed crop by-products while penned within the settlement , with access to the seasonal flooded landscape limited by herders to reduce foot rot in livestock . The δ13C and δ18O values from incremental bioapatite samples presented here provides us with the opportunity to test these hypotheses through integration of seasonal short-term stable isotopic records with those from bone collagen.
Oxygen isotopes of obligate and non-obligate drinkers
There were clear differences between ruminants and omnivores in their oxygen isotope observed at both sites. Sheep/goat samples from Vinča-Belo brdo and Stubline exhibit the largest average intra-tooth amplitude observed in domestic livestock, averaging 5.6‰ and 6.1‰, respectively for Vinča-Belo brdo and Stubline. The average amplitude per tooth for cattle from both sites were similar (Vinča-Belo brdo: 3.4; Stubline: 3.5‰), and close to that observed in deer (M2: 3.8‰; M3: 2.7‰), with domestic pigs exhibiting the lowest amplitude (<1 ‰; Fig 5). The oxygen stable isotopic composition of bioapatite is in near isotopic equilibrium with the oxygen isotopes of body water . These isotopic ratios are influenced by those of imbibed water, which is in turn affected by meteoric water δ18O values, which is influenced by altitude, local humidity, and temperature [29,73]. Drinking water source type, also influences body water oxygen isotope ratios. Groundwater typically represents a local meteoric water δ18O average, while rivers will vary seasonally depending on the input of melt water. Leaf water is enriched in 18O in comparison to the local average due to evapotranspiration . Here we discuss the differences in the midpoint value and intra-tooth amplitude between the maximum and minimum δ18O values between species. These will provide us an insight into water consumption habits of domesticated animals.
Cattle (B), pigs (SS), sheep/goat (OC) and red deer (CE) teeth from Vinča-Belo brdo (VBB) and Stubline (ST).
The large amplitude of intra-tooth oxygen isotopic change observed in sheep/goat may be due to ingestion of 18O enriched leaf water which comprises a large proportion of total body water in caprines . Leaf water oxygen isotope ratios are influenced by the local precipitation d18O values, aridity level, and evaporation rates so that leaf water is enriched in 18O relative to precipitation . Consequently, the δ18O record of semi-obligate species, which include caprines, is more sensitive to aridity . Serbia is characterised by a marked continental semi-arid climate, with high summer and low winter temperatures. Sheep and goat ingesting leaf water in arid and semi-arid environments may exhibit a large intra-tooth amplitude between high and low values within a single tooth due to animals obtaining water during the summer from leaf water . Our results chime with those previously seen in sheep/goat (4.8‰) from Popină Borduşani (5th millennium BC, Romania) located on the Danube .
Cattle, deer, and pigs are obligate drinkers, and it has been shown that leaf-water has limited effect on the δ18O values expressed in these species . The composition and seasonal variation in δ18O values of water sources is dependent on the frequency of a body of water is fed by sources, such as melt water or tributaries . The seasonal variation within δ18O values of standing water, such as lakes and wells, is dampened in comparison to local δ18O precipitation because of the effects of evaporation and influx water cancel each other out . In Laurentian great lakes, the seasonal shift in δ18O values was less than 2‰ . The amplitude between high and low δ18O values within single cattle and deer teeth was greater than 3‰, which is greater than observed by Jasechko, Gibson . It may suggest that cattle and deer ingested water from sources with a greater seasonal variation than observed in Laurentian lakes, like the Danube and the Sava. For pigs, the narrow range may reflect the tooth structure or formation or may reflect pigs ingesting water from standing water sources, such as wells.
The geometry of teeth between species will influence the stable isotopic record. Pig molars have different internal and external dental geometry in comparison to cattle and deer. The mineralisation process is faster in pigs and as with most species the timing of the process is not uniform across a crown , which may lead to a dampening of seasonal stable isotopic signal. This may explain the narrow range in oxygen values and low amplitude observed in pigs from this study. Frémondeau, Cucchi  demonstrated that incisors are a more suitable tooth for the purposes of capturing seasonality in δ18O values, perhaps due to a thinner enamel layers than is seen in molars. Future studies should incorporate incisors to directly investigate pig watering habits. This could potentially open-up new research avenue into analysing differences between species in terms of δ18O values to understand how communities may have managed water sources for livestock.
Sheep/Goat seasonal herding practices
Predicted dietary δ13C values from sheep/goat enamel bioapatite samples ranged from −28.6‰ to −23‰, considerably greater than the range previously observed in collagen samples (−26.4‰ to −23.4‰) . A strong synchroneity between the curves of δ13C and δ18O values, which indicates that that caprines consumed fodder sources that varied considerably in 13C over the seasonal cycle. Environmental factors, such as light levels and soil/atmospheric humidity, impact photosynthesis efficiency and respiration in C3 plants impact the isotopic fractionation of 12C/13C [82–84]. We propose three main hypothesis to explain the high variability in livestock δ13C values, 1) grazing in marsh land environments; 2) winter grazing in dense non-deciduous forests; 3) supplementary foddering during winter with leafy hay collected from either closed canopy forests during summer when the ‘canopy effect’ would be greatest, or collected from waterlogged galley forests, which also may exhibit low δ13C values .
The seasonally flooded landscape surrounding the sites may explain the large amplitude of intra-tooth change in δ13C values (Vinča-Belo brdo: 2.5; Stubline: 3.4)., which is greater than the values observed in hair samples from modern herbivores grazing on open pastures throughout the year (~1.5‰) . Plants growing in areas subject to seasonal flooding exhibit a strong seasonal variation in foliar δ13C values  with low carbon isotopes expressed in leaf tissue during flooding. This is caused by a decrease in CO2 intake by plants coupled with increased photosynthesis activity results in 12C being discriminated over 13C . In general, plant stomata during dry periods, particular in regions with hot summers, will close to minimizing water loss and reduce the discrimination against 13C, resulting in higher carbon isotopes values in plant tissues. For example, in semi-arid regions where the mean annual precipitation is less than 500mm, the mean value of C3 plants, adjusted for the Suess effect, were between −24.7‰ to −23.9‰ . We propose that animals pastured all-year round on a landscape with seasonal flooding under a strong continental climate will exhibit carbon values that reflect a strong seasonal pattern.
The Danube and Sava are active meandering rivers, with paleo-channels that are visible from satellite images. Focusing on Vinča-Belo Brdo, at present the Danube runs close to the site, however in the past it may have flowed further east. There is no evidence of flooding during the occupation levels, which suggests that the site was not as close to the river as today . Areas of standing water or vernal pools may form during flooding within extinct paleochannels . Plant communities in these pools vary in composition over a growing season due to changes in the groundwater levels. Coupled with changes in temperature from spring to summer, these environments may exhibit a large variation in foliar δ13C values. These environments contribute to osmotic stress in plants, thus reducing stomatal activity leading to an increase in the uptake of 13C increasing overall the foliar δ13C values . Therefore, animals grazing on these communities would exhibit δ13C values with strong seasonal variation with relatively high values in the summer.
Stubline and Vinča-Belo brdo lies on the edge of the Vojvodina region, which is part of the Pannonian Basin known for its salt steppe, pans, marshes, and lakes [89,90]. Most of these environments are of natural origin and support halophytes plants, such as succulents, xerophytic grasses, and other monocytes, that can tolerate high salt concentrations [89,90]. Drainage and intensive agricultural practices since the 18th and 19th century have led to loss of these environments . Halophyte plant communities exhibit high δ13C values due to the high salinity and low soil moisture, which can affect stomatal activity . Experiments have shown that C3 plants (Salicornia europaea ssp. and Puccinellia nuttalliana) native to highly saline marsh environments of Western Canada exhibit high δ13C values around −21.8‰ (adjusted for Suess effect), and vary between 6.5‰ and 4‰ across a 30m transect . Halophyte plant communities can be highly nutritious sources of forage for domesticated livestock by providing a source of salt and other minerals that may be absent in traditional pastures . We propose the strong synchronicity between δ13C and δ18O values may be a consequence of sheep/goat grazing on halophyte and/or vernal plant communities within the alluvial plain under a continental climatic regime.
The δ15N values of plants are also influenced by saline environments due to enrichment in soil 15N due to higher rates of denitrification that lead to 15N enriched soil nitrates . Moderately high δ15N values (8.1±1‰) were found in caprine collagen in comparison to cattle (7.1±0.9‰). Enriched δ15N may be an indication of sheep/goat individuals grazing within saline soils, as has been seen in other studies . On the other hand, it is possible that the 15N enrichment of bone collagen reflects animals visiting long-term pastures or grazing on crop by products under a manuring regime. Future studies are necessary to test this by determining the δ15N and δ13C values of remains of cereals and crop by-products.
During winter, the caprines from both sites exhibited low dietary δ13C values between −28.6‰ to −27.1‰, except for individual 7245 from Vinča-Belo Brdo. These values are comparable to ruminants dwelling in deciduous forested environments (−29.5‰ to −25.1‰)  and deer sampled from Vinča-Belo brdo (−29.7‰ to −25‰). Low plant foliar δ13C values can occur in forest with dense leaf cover caused by the canopy effect, which is caused by 13C depletion of atmospheric CO2 under the canopy due to uptake of recycled CO2 respired by 13C depleted organic matter , a decrease in light intensity near the forest floor decreasing efficiency of photosynthesis, discriminating against the transfer of 13C [83,84]. The net result is lower floral δ13C values in by 1 to 6‰ relative to floral growth in open environments, depending on canopy density which is dependent on tree species present in forests, the timing of the bud break, and leaf expansion [93,95,96]. In deciduous woodlands, canopy density is greatest during late spring/early summer, while non-deciduous forest stands retain its canopy throughout the year.
Charcoal analyses of Vinča-Belo brdo material indicate two main types of deciduous woodland environments were present near the settlement: riverine and hill slopes, with no large stands of non-deciduous trees . The loss of leaves during autumn reduces the canopy, thus caprines grazing in woodland during winter will not result in low carbon values. However, provisioning animals with winter leaf hay collected during the summer months when the trees are in full leaf would result low δ13C values registering in caprine tissues. Another possibility is that leafy hay was collected from arboreal species native to riverine environments, such as ash (Fraxinus excelsior). Waterlogging caused by seasonal flooding has a similar impact on leaf carbon isotopes values as the canopy effect, (86]. Provisioning animals with leafy hay during winter may explain the large intra-tooth amplitude in δ13C values observed in caprines. Compound specific stable isotope analysis (CSIA) of amino acids (AA) extracted from dentine can test for the consumption woody plants. This is due to the difference in δ15N values of glutamine (Glx) and phenylamine (Phe) AA, which has been found to be significantly different between herbaceous and woody plants . This method is independent of canopy effect principle concerning carbon isotope. When integrated with bioapatite analysis, it can verify winter leafy hay provisioning in bovids [98,99].
Variation in cattle pasture and foddering practices
Overall, the results from the sequential analysis indicate strong inter-individual variation in cattle, unlike that previously observed in sheep teeth. Previous analysis proposed that domestic cattle from Vinča-Belo brdo grazed across different types of pastures while animals at Stubline were kept mainly on wooded or waterlogged pastures . At Vinča-Belo brdo, dietary δ13C values from −28.2‰ to −22.1‰, based on bioapatite δ13C values exhibit a wider ranger greater than values previously observed for collagen (−26.4‰ to −23.0‰). This may reflect the attenuation of bone collagen stable isotope values due to bone turnover. The results from the sequential analysis indicate strong inter-individual variation in cattle from Vinča-Belo brdo, unlike that previously observed in sheep teeth.
Bos 7241 and Bos 7268 exhibit low δ13C values −26.1 to −24.5‰ with limited seasonal variation. However, three cattle individuals from Vinča-Belo brdo (Bos 7230, 7242 and 7244) exhibit δ13C values (>−24.1‰) as observed in caprine teeth. For Bos 7242 and 7244, these values occur in summer, as observed in caprine teeth. We propose this is the result of cattle grazing pastures supporting halophyte/vernal plant communities close to the sites, similar as proposed for caprines. For Bos 7230, high δ13C values reaching −21.1‰ during the period of winter season tooth formation may be evidence provisioning with hay collected from these vernal pastures in summer. There appears to be evidence for winter leaf hay provisioning in a single cattle, Bos 7240. This animal exhibits δ13C values as low as −28.2‰ during winter. Previous analyses of cattle from the LBK site of Bischoffsheim (Alsace ) and the 5th millennium BC site of Bercy (France ) found cattle exhibited winter season δ13C values lower than −27.0‰. Considering that cattle from Vinča-Belo brdo and Stubline would have experienced weather conditions of a strong continental character (hot summers, cold winters), this lack of evidence for supplementary fodder is surprising given that during winter the snow cover as well as spring floods may have restricted access to fresh pastures.
Our study benefits from paleoenvironment evidence from the teeth of deer, which are typically forest dwelling. Dietary δ13C values for red deer between −29.7‰ and −25.0‰ are similar to the range of values for modern forest dwelling ruminants  and mid-Holocene aurochs from Boreal environments . Carbon isotope values from individual 7240 fall within that of the deer dietary range. Leafy hay for cattle was collected from the same environment where the deer grazed would exhibit the same range of carbon isotopic values. It is also possible that Bos 7240 was also pastured within a woodland environment for the remaining part of the year. The dietary values from domesticated pigs (−27.2‰ to −25.6‰) were also within the deer range, thus potentially supporting the idea swine were managed under an open woodland pannage system for part of the year. The use of forest within the vicinity of Vinča-Belo brdo for pasture and pannage would support previous hypothesis proposed by Filipović, Marić  highlighting the diverse land use strategies practices by the Vinča communities.
Vinča communities’ land-use
The results from incremental bioapatite samples from domesticated and wild animals of Vinča-Belo brdo and Stubline allow for nuanced investigation of the land use by late Vinča communities. Previous paleoenvironmental and zooarchaeological studies have proposed that these communities exploited a wide range of niches for wild resources and to raise their herds, particularly cattle [9,11,21,22]. The analysis presented here and elsewhere  demonstrates the wide range in environments exploited by communities living at Stubline and Vinča-Belo brdo (Fig 6). The seasonally available floodplain meadows appear to have played an important role in caprine management systems, and suggests sheep were not kept within settlement all year round as previously suggested [11,21]. This hypothesis was based on the principal that sheep grazing on riverine landscapes leaves animals susceptible to foot rot [11,21] caused by bacteria (Fusobacterium necrophorum/Dichelobacter nodosus). It leads to lameness severely limiting animal access to food . Warm damp conditions found within permanent shelters promotes this condition, as do wet summers and warm winters, leading to a year-round endemic problem. Consequently, animals need to be grazed in large pastures and during winter, herders may have restricted grazing by penning animals within the settlement and provisioning animals leafy hay collected during late spring/summer to prevent foot rot and loss of sheep during floods.
The δ13C values indicate (from left to right): Deer and domesticated pigs sampled as part of this study, and Suess effected adjusted values for halophyte plants  to illustrate the gradient from dense to light forest canopies and salt marshes. The base images were made available on the ‘Nature scenes with fields’ https://www.vecteezy.com/vector-art/605564-nature-scenes-with-fields-and-mountians licensed to RG. The animal silhouettes are available under CC 4.0 from ArcheoZoo.org.
The use of crop by-products and stubble for animal pasture and fodder may have also occurred. Archaeobotanical analysis of carbonised cereals and weed assemblages from Opovo, a contemporaneous site north of Vinča-Belo brdo, suggests that the majority of crops were sowed in the winter and harvested in the summer . Allowing animals access after harvest would allow for manure to be directly distributed on cultivation plots. Crop plants may experience water stress in the final months prior to harvest and exhibit high δ13C values due to water stress , which may explain the δ13C values observed in sheep/goat individuals. The relatively high δ15N values found in sheep/goat bone collagen could reflect a foddering strategy where cereals were enriched with 15N due to manuring regimes [104,105].
Stable isotopic analysis of cattle teeth indicate several different husbandry strategies were applied. Some cattle exhibit δ13C and δ18O values similar to suggesting multi-species herding for all or part of the year (Fig 6), while other cattle exhibiting values more similar with cattle from contemporaneous LBK sites in Central Europe suggest both open grasslands or woodland grazing [106,107]. The large inter-individual variation seen in cattle may be due to different herding strategies applied to animals according to age and sex, use of more extensive pasturing strategies that evolved to cope with large herds, or grazing of small herds by individual households. Previous analysis of bone collagen values of sub-adult cattle (12 to 36 months) from at Popină Borduşani identified high variation in δ13C values between individuals. If these animals were raised for their meat, they may have been fed special diets to increase their weight or allowed to extensively range to enhance access to better quality pasture to increase their weight. The stable isotope values of cattle reflect diet during the formation of M3 (12–24 months). It is possible some of the individuals analysed here were male and this variation in δ13C values indicates animals being fattened for meat production.
Another source of carbon isotopic variation in cattle may be the emergence of extensive herding strategies where animals were able to roam between different habitats within the vicinity of the sites. This may represent a flexible and sustainable herding strategy where cattle were move naturally between pasture, reducing overgrazing. This type of management may only require one or two herdsmen that corrals individual household herds daily and moves them to pastures within the vicinity of the settlement. Collective herding strategies supports social cohesion and cooperation between households. Similar strategies are evident in traditional African societies, where village herds are graze together and returned to their owners in the evening for milking . The similarity in carbon isotopes evident within caprines suggests a collective management of small stock herds from different households while the diversity in cattle stable isotopic results may also indicate different households/groups managing their animals individually.
Continual herd growth likely took place alongside increased intensity of settlement and population size in during the Late Vinca phases and would have placed significant pressure on local pasture resources. Due to the proximity of Vinča with other settlements and coupled with the increase in specialisation of animal products , herders may have negotiated access to pastures in exchange for consumables (milk/butter/cheese/meat) and trade goods as well as beasts. These individual efforts on herding of cattle may be motivated by increasing symbolic status of cattle associated with their considerable productive capacity to generate large quantities of meat, milk and skins, as well their ability to transform landscapes via traction . This exchange may have resulted in the variation in stable isotope values found within the Vinča cattle. The negotiation and exchange of goods would of strengthened social relationships between Vinča communities  and the symbolic and economic importance for cattle. At the site of Opovo, Russell  proposes that over time, households may have built up their cattle herds, and this, coupled with the increasing representation of cattle in material culture, suggests the emergence of wealthy households . Furthermore, she proposes that it may not have been possible to function in the Vinca social sphere without owning cattle due to its importance as bride wealth, blood money, sacrifices, or for feast [Russell 199:159]. Only through the integrated stable isotopic and zooarchaeological analysis of other Vinča settlements across time and space can we begin to explore cattle-based wealth systems and its potential impact on Vinča social relationships.
Livestock herds were an important social and economic component of Vinča world and here we have uncovered the complex herding strategies employed for animal pasture and fodder. The combination of stable δ13C isotope values from different animal tissues has provided a means to reveal multi-faceted perspectives on herding strategies at different temporal scales (seasonal versus several years) thus highlighting the importance of multi-tissue analysis. Our study has uncovered the use by sheep flocks of seasonally available halophyte pastures and leafy hay, providing a new perspective on the sheep/goat management by Vinča communities. Future studies should consider including isotope analyses of botanical material to generate an integrated and detailed picture of herding and agricultural practices. There is also a need for detailed organic residues analysis of ceramic sherds to provide more information about processing of animal products.
The low variation in herding strategies for sheep/goat indicates a pooling of resources carried out a community level, while in contrast, the diversity in cattle pasturing strategies may be an indication of the emergence of household economies associated with the growing social and economic importance of cattle. This contrast in herding strategies between the two species highlights a level of complexity within the organisation of food procurement strategies for individual households and for the community in general. This study is limited to two sites, however, there is enormous potential for future stable isotopic studies considering the excellent preservation of faunal material from Vinča culture sites. More detailed analysis at other sites would enable a deeper understanding of Vinča culture herding practices as well as the landscape inhabited by these communities and their herds. Integrated with existing archaeological studies would provide a greater insight into the cultural evolution and the development of complex societies in Europe.
S1 Table. Sample list detailing context numbers, elements samples and age-at-death.
S2 Table. Detailed carbon isotope results for individual teeth.
Many thanks go to Dr Marie Balasse and Denis Fiorillo at the SMIMMS platform (MNHN, Paris) and Dr Nils Andersen at the Leibniz Labor (Kiel University) for assistance with mass spectrometry. Also, thanks go to Dr Delphine Frémondeau at UCL for her advice on the formation of and stable isotopic analysis of pig teeth.
- 1. Borić D, Hanks B, Šljivar D, Kočić M, Bulatović J, Griffiths S, et al. Enclosing the Neolithic World: A Vinča Culture Enclosed and Fortified Settlement in the Balkans. Current Anthropology. 2018;59(3):336–46.
- 2. Borić D. Absolute Dating of Metallurgical Innovations in the Vinča Culture of the Balkans. In: Kienlin TL, Roberts B, editors. Metals and Societies: Studies in honour of Barbara S Ottaway. Bonn: Verlag Dr Rudolf Habelt GMBH; 2009. p. 191–245.
- 3. Orton D, Gaastra J, Linden MV. Between the Danube and the Deep Blue Sea: Zooarchaeological Meta-Analysis Reveals Variability in the Spread and Development of Neolithic Farming across the Western Balkans. Open Quaternary. 2016;2.
- 4. Perić S, Bajčev O, Obradović Đ, Stojanović I. Neolitsko naselje Pavlovac-Gumnište: rezultati zaštitinih arheoloških iskopavanja 2011. In: Perić S, Bulatović A, editors. Arheološka istraživanja na autoputu E75 (2011–2014). Beograd: Arheološki institut; 2016. p. 221–74.
- 5. Spasić M. Metahousing: Neolithic and modern dwelling in Belgrade. Belgrade: Belgrade City Museum; 2012.
- 6. Crnobrnja AN. Group identities in the Central Balkan Late Neolithic. Documenta Praehistorica. 2012;XXXIX:155–65.
- 7. Tripković B. Domaćinstvo i zajednica: kućne i naseobinske istorije u kasnom neolitu centralnog Balkana. Beograd: Filozofski fakultet.; 2013.
- 8. Porčić M. An exercise in archaeological demography> estimating the population size of Late Neolithic settlements in the Central Balkans. Documenta Praehistorica. 2011;XXXVIII:323–32. pmid:29974304
- 9. Filipović D, Challinor D, Andrič M. Vinča tell in southeast Europe: Multi-proxy palaeobotanical evidence from Late Neolithic levels and the implications for the environment and economy. Quaternary International. 2017;429:13–23.
- 10. Filipović D, Obradović Ð. Archaeobotany At Neolithic Sites In Serbia: A Critical Overview Of The Methods And Results. In: Miladinović-Radmilović N, Vitezović S, editors. Bioarchaeology in the balkans balance and perspectives. Belgrade: Sremska Mitrovic; 2013. p. 25–56.
- 11. Orton D. Herding, Settlement, and Chronology in the Balkan Neolithic. European Journal of Archaeology. 2012;15(1):5–40.
- 12. Bökönyi S. The Neolithic fauna of Divostin. In: McPherron A, Srejović D, editors. Divostin and the Neolithic of Central Serbia. Pittsburg: University of Pittsburgh; 1988. p. 419–45.
- 13. Greenfield HJ. The origins of secondary product exploitation and the zooarchaeology of the Late Neolithic, Eneolithic, and Middle Bronze Age at Vinča-Belo Brdo, Serbia: the 1982 excavations. In: Greenfield HJ, editor. Animal Secondary Products Domestic Animal Exploitation in Prehistoric Europe, the Near East and the Far East. Oxford: Oxbow books; 2014. p. 274–334.
- 14. Bökönyi S. Tierknochenfunde der neuesten ausgrabungen in Vinča. In: Srejović D, Tasić NN, editors. Vinča and Its World: International Symposium The Danubium Region from 6000 to 3000 BC. Belgrade: Smederevska Palanka; 1990. p. 49–54. pmid:1966865
- 15. Clason A. The farmers of Gomolava in the Vinča and La Tène period. Rad Vojvođanskih Muzeja. 1979;25:60–114.
- 16. Bulatović J. Arheozoološki aspekti društvenih i kulturnih promena na centralnom Balkanu u petom milenijumu pre nove ere. Belgrade: University of Belgrade.; 2018.
- 17. Dimitrijević V. Vertebrate fauna of Vinča-Belo Brdo (excavation campaigns 1998–2003). Starinar. 2008;56:245–69.
- 18. Greenfield HJ. Fauna from the Late Neolithic of the Central Balkans: issues in subsistence and land use. Journal of Field Archaeology. 1991;18:161–86.
- 19. Russell N. Hunting, Herding and Feasting: human use of animals in Neolithic Southeast Europe. Berkeley: University of California; 1993.
- 20. Halstead P, Tierney J, Butler S, Mulder Y. Leafy Hay: an Ethnoarchaeological Study in NW Greece. Environmental Archaeology. 1996;1(1):71–80.
- 21. Filipović D, Marić M, Challinor D, Bulatović J, Tasić N. Natural environment and resources, and the long life of the Neolithic settlement at Vinča, southeast Europe. Archaeological and Anthropological Sciences. 2018;11(5):1939–60.
- 22. Gillis RE, Bulatović J, Penezić K, Spasić M, Tasić NN, Makarewicz CA. Herding and hunting at Vinča–Belo Brdo and Stubline during the Late Neolithic, a stable isotopic perspective. In: Bulatović J, Marković N, editors. Animal husbandry and hunting in the Southeast Europe and adjacent areas through time. Oxford: Archaeopress; 2020.
- 23. Makarewicz CA. Winter is coming: seasonality of ancient pastoral nomadic practices revealed in the carbon (δ13C) and nitrogen (δ15N) isotopic record of Xiongnu caprines. Archaeological and Anthropological Sciences. 2015;9(3):405–18.
- 24. Mainland I, Towers J, Ewens V, Davis G, Montgomery J, Batey C, et al. Toiling with teeth: An integrated dental analysis of sheep and cattle dentition in Iron Age and Viking–Late Norse Orkney. Journal of Archaeological Science: Reports. 2016;6:837–55.
- 25. Balasse M, Boury L, Ughetto-Monfrin J, Tresset A. Stable isotope insights (δ18O,δ13C) into cattle and sheep husbandry at Bercy (Paris, France, 4th millennium BC): birth seasonality and winter leaf foddering. Environmental Archaeology. 2013;17(1):29–44.
- 26. Fraser RA, Bogaard A, Schäfer M, Arbogast R, Heaton THE. Integrating botanical, faunal and human stable carbon and nitrogen isotope values to reconstruct land use and palaeodiet at LBK Vaihingen an der Enz, Baden-Württemberg. World Archaeology. 2013;45(3):492–517.
- 27. Gillis RE, Eckelmann R, Filipovic D, Muller-Scheessel N, Cheben I, Furholt M, et al. Stable isotopic insights into crop cultivation, animal husbandry, and land use at the Linearbandkeramik site of Vrable-Velke Lehemby (Slovakia). Archaeological and Anthropological Sciences. 2020;12(11). pmid:33123298
- 28. Tieszen LT, Fagre T. Effect of diet quality and composition on the isotopic composition of respiratory CO2, Bone collagen, Bioapatite, Soft tissues. In: Lambert JB, Grupe G, editors. Prehistoric human bone: Archaeology at the molecular level. Berlin: Springer-Verlag; 1993. https://doi.org/10.1111/j.1464-5491.1993.tb00121.x pmid:8365094
- 29. Gat JR. Oxygen and hydrogen isotopes in the hydrologic cycle. Annual Revue Earth Planet Science 1996;24:225–62.
- 30. Balasse M. Reconstructing dietary and environmental history from enamel isotopic analysis: time resolution of intra-tooth sequential sampling. International Journal of Osteoarchaeology. 2002;12(3):155–65.
- 31. Koch PL, Fisher AD, Dettman D. Oxygen isotopes variation in the tusks of extinct proboscideans: A measure of season of death and seasonality. Geology. 1989;17:515–9.
- 32. Ристић-Опачић Ј. Топографско-хронолошке карактеристике насеља винчанске културе на територији Србије. Гласник Српског археолошког друштва. 2005;21:75–112.
- 33. Chapman JC. Vinča Culture of South-East Europe: Studies in Chronology, Economy and Society. Oxford: Archaeopress; 1981.
- 34. Whittle A, Bayliss A, Barclay A, Gaydasrka B, Bánffy E, Borić D, et al. A Vinča potscape: formal chronological models for the use and development of Vinča ceramics in south-east Europe. Documenta Praehistorica. 2016;43:1–60.
- 35. Crnobrnja AN. Arrangement of Vinca culture figurines a study of social structure and organization. Documenta Praehistorica. 2011;XXXVIII:131–47.
- 36. Spasić M, Živanović S, Stojić D. Bowls With Everted Rim From Stubline—An Emblematic Vessel Type Of Late Vinča Culture. Studia Universitatis Hereditati. 2018;6:65–87.
- 37. Tasić NN. Anthropomorphic figurines from Vinča excavations 1998–2009. Documenta Praehistorica. 2011;XXXVIII:149–57.
- 38. Spasić M. A Group Find of Neolithic Figurines of the Vinča Culture from Stubline, Serbia. Les Carnets de l’ACoSt. 2014(12).
- 39. Hansen S. Bilder vom Menschen der Steinzeit. Untersuchungen zur anthropomorphen Plastik der Jungsteinzeit und Kupferzeit in Sudosteuropa. Mainz: Philipp von Zabern; 2007.
- 40. Vuković J. Late Neolithic pottery standardization: application of statistical analyses. Starinar. 2011;LXI: 81–100.
- 41. Tasić N. New evidence on salt use in the Neolithic of Southeast Europe. In: Nikolov V, Bacvarov K, editors. Salz und Gold: die Rolle des Salzes im prähistorischen Europa/Saltand Gold: The Role of Salt in Prehistoric Europe. Provadia: Veliko Tarnovo; 2012. p. 213–8.
- 42. Šljivar D, Kuzmanović-Cvetković J, Jacanović D. New contributions regarding the copper metallurgy in the Vinča culture. In: Tasić NN, Grozdanov C, editors. Homage to Milutin Garašanin. 251–266. Belgrade 2006.
- 43. Šljivar D, Kuzmanović-Cvetković J, Živković J. Pločnik: on copper metallurgy in the Vinča culture. Zbornik Narodnog Muzeja. 2012;20:27–46.
- 44. Radivojević M, Rehren T, Pernicka E, Šljivar D, Brauns M, Borić D. On the origins of extractive metallurgy: new evidence from Europe. Journal of Archaeological Science. 2010;37(11):2775–87.
- 45. Ottaway BS. Innovation, production, and specialization in early prehistoric copper metallurgy. European Journal of Archaeology. 2001;4:87–112.
- 46. Crnobrnja AN. The (E)neolithic Settlement Crkvine at Stubline, Serbia. In: Schier W, Draşovean F, editors. The Neolithic and Eneolithic in Southeast Europe; New approaches to dating and cultural Dynamics in the 6th to 4th Millennium BC. Prähistorische Archäologie in Südosteuropa. 28. Leidorf: Rahden/Westf.; 2014. p. 173–86.
- 47. Borić D. The End of the Vinča World: Modelling the Neolithic to Copper Age Transition and the Notion of Archaeological Culture. In: Hansen S., Raczky P., Anders A., Reingruber A, editors. Neolithic and Copper Age between the Carpathians and the Aegean Sea Chronologies and Technologies from 6th to 4th millennium BCE. Archäologie in Eurasien. Bonn: Verlag Marie Leidorf; 2015. p. 177–237.
- 48. Porčić M. An exercise in archaeological demography> estimating the population size of Late Neolithic settlements in the Central Balkans. Documenta Praehistorica. 2011;XXXVIII:323–32. pmid:29974304
- 49. Penezić K. Vinča-Belo Brdo settlement size. In: Urem-Kotsou D, Tasić NN, Burić M, editors. Making Spaces into Places The North Aegean, the Balkans and Western Anatolia in the Neolithic. Oxford: British Archaeological Reports; forthcoming.
- 50. Tasić NN, Marić M, Filipović D, Penezić K, Dunbar E, Reimer P, et al. Interwoven Strands for refining the Chronology of the Neolithic tell of Vinča-Belo Brdo, Serbia. Radiocarbon. 2016;58:795–831.
- 51. Tasić N, Marić M, Bronk Ramsey C, Kromer B, Barclay A, Bayliss A, et al. Vinča-Belo Brdo, Serbia: The times of a tell. Germania. 2015;93:1–76.
- 52. Tasić N, Marić M, Penezić K, Filipović D, Borojević K, Russell N, et al. The end of the affair: formal chronological modelling for the top of the Neolithic tell of Vinča-Belo Brdo. Antiquity. 2015;89(347):1064–82.
- 53. Crnobrnja AN. Investigations of Late Vinca House 1/2010 at Crkvine in Stubline. Starinar. 2012;LXII:45–64.
- 54. Spasić M. Neolitsko naselje u Stublinama. Godišnjak grada Beograda. 2013;LX:11–42.
- 55. Peel MC, Finlayson BL, McMahon TA. Updated world map of the Köppen-Geiger climate classification. Hydrology and Earth System Sciences Discussions, European Geosciences Union,. 2007;11(5):1633–44.
- 56. Cedar Lake Ventures I. Weather Spark [Available from: https://weatherspark.com/y/85769/Average-Weather-in-Belgrade-Serbia-Year-Round.
- 57. Mauri A, Davis BAS, Collins PM, Kaplan JO. The climate of Europe during the Holocene: a gridded pollen-based reconstruction and its multi-proxy evaluation. Quaternary Science Reviews. 2015;112:109–27.
- 58. Brown WAB, Christofferson PV, Massler M, Weiss Mb. Postnatal Tooth Development in Cattle. American Journal of Veterinary Researc. 1960;21(80):7–34. pmid:13805043
- 59. Ewbank JM, Phillipson DW, Whitehouse RD, Higgs ES. Notes: Sheep in the Iron Age: A Method of Study. Proceedings of the Prehisotric society. 1964;30:423–6.
- 60. Grant A. The use of tooth wear as a guide to the age of domestic ungulates. In: Wilson R, Grigson C, Payne S, editors. Ageing and sexing animal bones from archaeological sites. 109. Oxford: BAR International series; 1982. p. 91–108.
- 61. Carter RJ, Magnell O. Age estimation of wild boar based on molariform mandibular tooth development and its application to seasonality at the Mesolithic site of Ringkloster, Denmark. In: Albarella U, Dobney K, Ervynck A, Rowley-Conwy P, editors. Pigs and Humans: 10,000 Years of Interaction. Oxford: Oxford University press; 2007. p. 197–217.
- 62. Robinson C, Kirkman J, Weatherall JA, Richards A, Josephen K, Fejeskov O. Development stages in Permanent Porcine enamel. Acta anatomica. 1987;128:1–10. pmid:3825482
- 63. Frémondeau D, Cucchi T, Casabianca F, Ughetto-Monfrin J, Horard-Herbin MP, Balasse M. Seasonality of birth and diet of pigs from stable isotope analyses of tooth enamel (δ18O, δ13C): a modern reference data set from Corsica, France. Journal of Archaeological Science. 2012;39(7):2023–35.
- 64. Brown WAB, Chapman NG. Age assessment of red deer (Cervus elaphus): from a scoring scheme based on radiographs of developing permanent molariform teeth. Journal of Zoology Society of London. 1991;225:85–97.
- 65. Balasse M. Potential biases in sampling design and interpretation of intra-tooth isotope analysis. International Journal of Osteoarchaeology. 2003;13(1–2):3–10.
- 66. Cerling TE, Harris JM. Carbon isotope fractionation between diet and bioapatite in ungulate mammals and implications for ecological and paleoecological studies. Oecologia. 1999;120:347–63. pmid:28308012
- 67. Passey BH, Robinson TF, Ayliffe LK, Cerling TE, Sponheimer M, Dearing MD, et al. Carbon isotope fractionation between diet, breath CO2, and bioapatite in different mammals. Journal of Archaeological Science. 2005;32(10):1459–70.
- 68. Zazzo A, Balasse M, Passey B, Moloney AP, Monahan FJ, Schmidt O. The isotope record of short- and long-term dietary changes in sheep tooth enamel: Implications for quantitative reconstruction of paleodiets. Geochimica et Cosmochimica Acta. 2010;74:3571–86.
- 69. Zazzo A, Balasse M, Patterson WP. High-resolution δ13C intra-tooth profiles in bovine enamel: Implications for mineralization pattern and isotopic attenuation. Geochimica et Cosmochimica Acta. 2005;69(14):3631–42.
- 70. Balasse M, Bălăşescu A, Tornero C, Fremondeau D, Hovsepyan R, Gillis RE, et al. Investigating the scale of herding in Chalcolithic pastoral communities settled along the Danube River in the 5th millennium BC: A case study at Borduşani-Popină and Hârşova-tell (Romania). Quat Int. 2017;436:29–40.
- 71. Ferrio JP, Araus JL, Buxó R, Voltas J, Bort J. Water management practices and climate in ancient agriculture: inferences from the stable isotope composition of archaeobotanical remains. Vegetation History and Archaeobotany. 2005;14(4):510–7.
- 72. Luz B, Kolodny Y, Horowitz M. Fractionation of oxygen isotopes between mammalian bone phosphate and environmental drinking water. Geochimica et Cosmochimica Acta. 1984;48:1689–93.
- 73. Dansgaard W. Stable isotopes in precipitation. Tellus. 1964;16:436–68.
- 74. Pendazani S, Britton K. Oxygen isotopes in bioarchaeology: Principles and applications, challenges and opportunities. Earth-Science Reviews. 2019;188:77–107.
- 75. Levin NE, Cerling T, Passey B, Harris DA, Ehleringer J. A stable isotope aridity index for terrestrial environments. PNAS. 2006;103:11201–5. pmid:16840554
- 76. Cernusak LA, Barbour MM, Arndt SK, Cheesman AW, English NB, Feild TS, et al. Stable isotopes in leaf water of terrestrial plants. Plant Cell Environ. 2016;39:1087–102. pmid:26715126
- 77. Makarewicz CA, Pendazani S. Oxygen (δ18O) and carbon (δ13C) isotopic distinction in sequentially sampled tooth enamel of co-localized wild and domesticated caprines: Complications to establishing seasonality and mobility in herbivores. Palaeogeography, Palaeoclimatology, Palaeoecology. 2017.
- 78. Balasse M, Tresset A, Balasescu A, Blaise E, Tornero C, Gandois H, et al. Animal Board Invited Review: Sheep birth distribution in past herds: a review for prehistoric Europe (6th to 3rd millennia BC). Animal. 2017;11(12):2229–36. pmid:28532521
- 79. Gat JR. Stable isotopes of fresh and Saline lakes. Physics and Chemistry of Lakes. Berlin, Heidelberg: Springer 1995. p. 139–65.
- 80. Jasechko S, Gibson JJ, Edwards TWD. Stable isotope mass balance of the Laurentian Great Lakes. Journal of Great Lakes Research. 2014;40(2):336–46.
- 81. Jasechko S, Kirchner JW, Welker JM, McDonnell JJ. Substantial proportion of global streamflow less than three months old. Nature Geoscience. 2016;9(2):126–9.
- 82. Arens NC, Jahren AH, Amundson R. Can C3 plants faithfully record the carbon isotope composition of atmospheric carbon dioxide? Paleobiology. 2000;26(1):137–64.
- 83. Farquhar GD, Ehleringer J, Hubick KT. Carbon Isotope Discrimination And Photosynthesis. Annual Review of Plant Biology. 1989;40:503–37.
- 84. Van der Merwe NJ, Medina E. The Canopy Effect, Carbon Isotope Ratios and Foodwebs in Amazonia. Journal of Archaeological Science. 1991;18:249–59.
- 85. Fan R, Morozumi T, Maximov TC, Sugimoto A. Effect of floods on the δ13C values in plant leaves: a study of willows in Northeastern Siberia. PeerJ. 2018;6:e5374. pmid:30258705
- 86. Dungait JA, Docherty G, Straker V, Evershed RP. Seasonal variations in bulk tissue, fatty acid and monosaccharide delta(13)C values of leaves from mesotrophic grassland plant communities under different grazing managements. Phytochemistry. 2010;71(4):415–28. pmid:20079504
- 87. Reade H, Stevens RE, O’Connell TC, Barker G. Pleistocene and Holocene herbivore diets and palaeoenvironments in the Gebel Akhdar (Libya): Implications for past human populations. Palaeogeography, Palaeoclimatology, Palaeoecology. 2016;449:62–78.
- 88. Muller SD, Bruneton H, Soulie-Marsche I, Rey T, Thiery A, Waterkeyn A, et al. Long-Term Dynamics of a Mediterranean Alkaline Vernal Pool (Rhone Delta, Southern France). Wetlands. 2008;28(4):951–66.
- 89. Šefferová Stanová V, Janák M, Ripka J. Management of Natura 2000 habitats. 1530 *Pannonic salt steppes and salt marshes. European Commission; 2008.
- 90. Luković M, Stevanovic ZD, editors. Saline habitats as undiscovered ecotourism attractions. Tourism in function of development of the Republic of Serbia: Tourism and Rural Development (TISC 2020)–Thematic proceedings II; 2020.
- 91. Dajic-Stevanovic Z, Pecinar I, Kresovic M, Vrbnicanin S, Tomovic L. Biodiversity, utilization and management of grasslands of salt affected soils in Serbia. Community Ecology. 2008;9:107–14.
- 92. Guy RD, Reid DM, Krouse HR. Shifts in Carbon Isotope Ratios of Two C3 Halophytes Under Natural and Artificial Conditions. Oceologia Plantarum. 1980;44:241–7. pmid:28310563
- 93. Heaton TH. The 15N/14 N ratios of plants in South Africa and Namibia: relationship to climate and coastal/saline environments. Oecologia. 1987;74:236–46. pmid:28311996
- 94. Drucker DG, Bridault A, Hobson KA, Szuma E, Bocherens H. Can carbon-13 in large herbivores reflect the canopy effect in temperate and boreal ecosystems? Evidence from modern and ancient ungulates. Palaeogeography, Palaeoclimatology, Palaeoecology. 2008;266(1–2):69–82.
- 95. Tieszen LT. Natural variations in the carbon isotope values of plants: Implications for Archaeology, Ecology and Paleoecology. Journal of Archaeological Science. 1991;18:227–48.
- 96. Leuschner C, Meier IC. The ecology of Central European tree species: Trait spectra, functional trade-offs, and ecological classification of adult trees. Perspectives in Plant Ecology, Evolution and Systematics. 2018;33:89–103.
- 97. Kendall IP, Woodward P, Clark JP, Styring AK, Hanna JV, Evershed RP. Compound-specific δ15N values express differences in amino acid metabolism in plants of varying lignin content. Phytochemistry. 2019;161:130–8. pmid:30826700
- 98. Kendall IP. Development of an amino acid δ15N value-based proxy for the elucidation of the diets and habitats of Neolithic cattle: University of Bristol; 2017.
- 99. Gillis RE, Kendall IK, Balasse M, Evershed RP. Conduite de l’élevage dans un environnement forestier à Bischoffsheim: une approche biogéochimique. In: Arbogast R-M, Lefranc P, editors. Bischoffsheim: un village du Néolithique ancien en basse Alsace (France). Archéologiques. Luxembourg Musée national d’Histoire et d’Art du Grand Duché du Luxembourg; in press.
- 100. Noe-Nygaard N, Price TD, Hede SU. Diet of aurochs and early cattle in southern Scandinavia: evidence from 15N and 13C stable isotopes. Journal of Archaeological Science. 2005;32(6):855–71.
- 101. Bennett GN, Hickford JG. Ovine footrot: new approaches to an old disease. Vet Microbiol. 2011;148(1):1–7. pmid:20926208
- 102. Borojević K. Terra and Silva in the Pannonian Plain: Opovo agro-gathering in the Late Neolithic. Oxford: Archaeopress; 2006.
- 103. Winter-Schuh C, Schmütz K, Eberle N, Makarewicz CA. Restricted pasturing of domesticated cattle at a Late Neolithic settlement in Central Germany. Journal of Archaeological Science: Reports. 2018;22:285–97.
- 104. Fraser R, Bogaard A, Heaton TH, Charles M, Jones G, Christensen BT, et al. Manuring and stable nitrogen isotope ratios in cereals and pulses: towards a new archaeobotanical approach to the inference of land use and dietary practices. Journal of Archaeological Science. 2011;38(10):2790–804.
- 105. Styring AK, Fraser RA, Bogaard A, Evershed RP. The effect of manuring on cereal and pulse amino acid δ15N values. Phytochemistry. 2014;102:40–5. pmid:24631496
- 106. Kendall IP, Gillis RE, Balasse M, Evershed RP. Hodowla bydła w świetle analizy izotopów trwałych/Stable isotope perspectives of cattle husbandry practices. In: Pyzel J, editor. Ludwinowo, stanowisko 7 Osada neolityczna na Kujawach/Ludwinowo, site 7 Neolithic settlement in Kuyavia. Ocalone Dziedzictwo Archeologiczne. Pękowice–Gdańsk: PROFIL-ARCHEO; 2019. p. 277–88.
- 107. Berthon R, Kovačiková L, Tresset A, Balasse M. Integration of Linearbandkeramik cattle husbandry in the forested landscape of the mid-Holocene climate optimum: Seasonal-scale investigations in Bohemia. Journal of Anthropological Archaeology. 2018;51:16–27.
- 108. Dahl G, Hjort A. Having Herds. Stockholme: University of Stockholme; 1976.
- 109. Bogucki P. Animal traction and household economies in Neolithic Europe. Antiquity. 1993;67:492–503.
- 110. Russell N. Cattle as wealth in Neolithic Europe: where’s the beef? In: Bailey D, editor. The archaeology of value: essays on prestige and the process of valuation. 730. Oxford 1998.