Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Using cementochronology to assess the seasonality of catastrophic events in medieval mass graves (Kutná Hora-Sedlec, Czechia, 14th century): Preliminary results

  • Eliška Zazvonilová ,

    Roles Conceptualization, Investigation, Writing – original draft, Writing – review & editing

    zazvonie@natur.cuni.cz

    Affiliations Institute of Archaeology of the Czech Academy of Sciences, Prague 1, Czech Republic, Faculty of Science, Department of Anthropology and Human Genetics, Charles University, Prague 2, Czech Republic

  • Hana Brzobohatá,

    Roles Conceptualization, Data curation, Project administration, Writing – original draft, Writing – review & editing

    Affiliation Institute of Archaeology of the Czech Academy of Sciences, Prague 1, Czech Republic

  • Jan Frolík,

    Roles Data curation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Institute of Archaeology of the Czech Academy of Sciences, Prague 1, Czech Republic

  • Petr Velemínský,

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation Department of Anthropology, National Museum, Prague 1, Czech Republic

  • Jaroslav Brůžek

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation Faculty of Science, Department of Anthropology and Human Genetics, Charles University, Prague 2, Czech Republic

Abstract

While season-of-death estimation using cementochronology is routine in archaeozoology, its use is much less frequent in bioarchaeology. Based on the character of the outermost increment (bright or dark), two seasons (spring/summer, autumn/winter) can be distinguished. Although many studies mention its potential and possible use in forensic anthropology or bioarchaeology, few exist with estimation results. This study aimed to apply cementochronology–a histological method based on counting and assessing regular circa-annual acellular cementum increments–to 42 individuals from medieval mass graves from Kutná Hora-Sedlec (Czechia, 14th century) to estimate the season-of-death. The mass graves belong to two stratigraphically distinct groups; written and archaeological sources relate them to two catastrophic events (the famine of 1318 and the plague epidemic of 1348–1350). Using cementochronology, we distinguished two distinct seasons corresponding to the two groups of graves, with individuals from the first group dying predominantly in spring/summer, while those from the second group died in autumn/winter. Taking into account the typical seasonal dynamics of epidemics, the results would be more in line with written sources. However, during the evaluation, we faced difficulties identifying the outermost increment and detecting the dark (thinner) increment; we recommend including only young and middle-aged adults in future studies, due to the difficulty of evaluation, and to consider the readability of the tissue (often affected by diagenesis). In conclusion, cementochronology has potential in the context of estimating the season-of-death, but the technical possibilities for enhancing the outermost increment need to be addressed, and the amount of data analysed expanded.

Introduction

Reconstructing seasonality in bioarchaeology

Seasonality in a human-environmental relationship is an important and complex issue in bioarchaeology [1]. The reconstruction of seasonal patterns finds application in topics such as the seasonality of particular activities, occupation of sites, or function of some particular sites [25]. Methods for investigating past seasonality, such as isotope analyses or the study of the microstructure of tooth cementum, can also provide a deeper understanding particularly of shorter-term activities in comparison to radiocarbon dating [1]. The study of the microstructure of tooth cementum using the method of cementochronology has a long tradition, and not only in archaeology [6]. This histological method is based on the seductive principle of the circa-annual cementum increments in the acellular tooth cementum [7]. Under light microscopy, each pair of light and dark bands represents one year [8,9]. The season-of-death can be assessed based on the outermost increment –a bright thicker increment represents spring/summer, and a dark thinner band the autumn/winter months [8,10]. Compared to archaeozoology, where cementochronology is routinely used for age and season-of-death estimation, in anthropology, the method is not often used, and if so, for age estimation, with a few studies on the season-of-death assessment [e. g., 11–13]. Preserved teeth allow season-of-death estimation not only in cases of the topics mentioned above, but also in another important topic related to seasonality–episodes of famine or infectious diseases.

Catastrophic events in medieval Kutná Hora-Sedlec

Catastrophic events such as famine or infectious diseases were not exceptional in Europe during the Middle Ages. Kutná Hora-Sedlec (Czechia) represents the most extensive set of medieval mass graves relating to catastrophic events on a European scale [14]. During the reconstruction of a church with an ossuary, 32 complete and incomplete mass burials (with around nine hundred skeletons) and 878 individual graves were excavated [14,15]. Kutná Hora-Sedlec has been intensively studied from the archaeological [15,16] and the bioarchaeological and palaeodemographic [14,17,18] perspectives. The chronology of the mass graves at the ossuary is based on an investigation of the stratigraphic configuration, with two levels of pits [15]. Younger mass graves have been preliminarily assigned to the fatalities of the 1348–1350 plague epidemic (based on coin recoveries from two of them), while the older level is provisionally linked to famine and the burial of victims before the gates of Sedlec in Kutná Hora, based on the 1318 report of the Zbraslav Chronicle (p. 248) [19]. The plague epidemic in the Czech lands is also mentioned in several other historical sources [20], but this area was less affected area than, for example, France or Italy [21]. This leads to fewer written records of the Bohemian plague outbreak. While written materials are invaluable sources of information that cannot be obtained from skeletal remains, the link to the plague is a preliminary result, and archaeogenetic analysis can provide a definitive confirmation. However, it can be argued that the catastrophic events which occurred at Kutná Hora-Sedlec in the 14th century caused increased mortality, manifested by a change in the funerary ritual (mass graves).

The seasonal dynamics of famine and plague

Famine and plague periods, mentioned in the historical documents, can be considered as probable causes of the mass dying at Kutná Hora-Sedlec. Estimating seasonality alone cannot solve the problem of the apparent causes of the mass graves, but it can help to contribute to understanding the differences between their two stratigraphically distinct levels. We can also get closer to interpreting the cause of mass graves by linking the results of seasonality estimation to the fact that several diseases have specific seasonal dynamics. Plague, in particular, is greatly influenced by seasonal factors [22,23]. Famine is usually defined as a discrete mass starvation event triggered by a food shortage [24,25]. The causes can be both exogenous and endogenous (e. g. extreme weather, bovine diseases, or factors internal to society, such as regulations and institutions) [26]. The complex scenario of the outbreak of famine can be demonstrated in the so-called Great Famine that afflicted Northern and Western/Central Europe between 1315 and 1322, affecting the Czech lands as well. This was recorded by chroniclers of the time, and narrative sources of Czech provenance also exist (e. g. the Chronicon Regiae Aulae, written at the Cistercian monastery of Zbraslav near Prague by Peter of Zittau (1275–1339) [19]. The Great Famine was caused by an agrarian crisis that happened in a period of climate transition: between the Medieval Climatic Optimum (lasting roughly until 1200) and the onset of a cooler period, the Little Ice Age, beginning in the early 14th century, when average temperatures declined [27]. As for the season of the peak of famine mortality, we rely mainly on demographic studies of post-medieval populations and the death registers of parishes [2830]. These sources show that demographic crises caused by widespread and severe food scarcity culminate in the Northern Hemisphere in the winter or early spring months. Plague is the only medieval epidemic for which the temporal and geographical propagation characteristics have been studied and described [23,31]. The temperature plays a significant role in the timing and growth of epidemics [32].

The data composed during the time of the Medieval Black Death suggest that mortality peaked mostly during the warmer weather months between April and October [33]. This peak mortality coincides with a seasonal peak in travel and trade during this period [23] investigated seasonal patterns and the epidemic peak timing of historical plague outbreaks (mainly from the 16-19th centuries) together with its association with latitude as a rough proxy measure for climatic conditions. They found that the typical plague season differed by latitude, and that the more northerly the location, the later the epidemic mortality peak occurred. Based on the findings of [23] and the city’s latitude (49°56′54″N), Kutná Hora would fall into the category of having August–September plague mortality peaks.

Our study aimed to apply cementochronology for the season-of-death estimation in 42 individuals from five mass graves coming from two stratigraphically different levels. Two distinct seasons were expected in the examined burials, hypothetically corresponding to two catastrophic events recorded in historical documents: plague and famine periods, with different seasonal peaks in mortality. We also sought to consider technical aspects of using cementochronology in bioarchaeology, including avoidable and unavoidable issues.

Material and methods

Material

Thirty-two complete or incomplete mass graves have been excavated at Kutná Hora-Sedlec. For our purposes, individuals from both levels of mass graves were chosen. Since the detection of the outermost line in older individuals could be problematic (due to the high number of cementum increments), only teeth from juveniles, young and middle-aged adults were included. In terms of tooth type selection, canines were preferentially chosen because of their lower susceptibility to caries [34]. Where canines were not present, first or second incisors were used. Only healthy teeth without visible caries were included. The ages of all individuals were previously estimated based on the degree of dental attrition, and the metamorphosis of sternal rib ends and articular facets of the ilium [3539]. Sampling took place at the place where the skeletal remains are deposited, in the Department of Anthropology of the National Museum in Prague (grave numbers are listed in the Material section and also in S1 Table). No permits were required for the described study, which complied with all relevant regulations.

Of the total number of 42 examined individuals, 22 belonged to two "plague" mass graves, and 20 individuals belonged to "famine" mass graves. The 5 mass graves were marked with the letters A-E (plague mass graves nos. 516 and 764 = A and B, famine mass graves nos. 765, 853 and 1578 = C, D, E) (Fig 1). Only one tooth was selected from each individual. The distribution of selected individuals relative to the total number of individuals within the mass grave is as follows: A 8/43, B 14/54, C 16/33, D 2/17, E 2/23.

thumbnail
Fig 1. The location of Kutná Hora-Sedlec (indicated by a red dot), based on Open Maps For Europe (https://www.mapsforeurope.org), © EuroGeographics [2023].

https://doi.org/10.1371/journal.pone.0295757.g001

Methods

For the preparation of thin sections, the standardized protocol [40] was used. According to protocol, teeth were at first cleaned with distilled water. Epoxy resin (Araldite 2020®) was used for embedding the teeth (without the crown). Subsequently, the tooth embedded into the resin was fixed in a holder to a precision saw with a diamond blade (Isomet Buehler 1000®), and five sequential cuts (from the middle part of the root) with a thickness of 100 μm were prepared at low speeds. Sections were mounted on slides and covered with coverslips (using Canada balsam). Each section was observed using a microscope (Olympus Bx-51), using polarized light (at 200x and 400x magnification). Segments with readable cementum were captured and saved as TIFF images. These were then edited and processed in Fiji (ImageJ) software. Prior to the age or season-of-death, the quality of the tissue affecting the results should be considered [9,41,42]. Samples were therefore included or excluded based on the readability index (measuring the distinctness of the tissue, scale 0–4; [7]). Samples (thin sections) with zero values in all five sections/tooth were excluded from the sample. Diagenetic apatite crystals can create false increments which should be distinguished from true cementum [5]. This can be done by using lambda plate, which is a frequently used approach (e. g., [4,5,43]).

Recording the season-of-death was performed visually, based on the saved images. After detecting the readable area of acellular cementum and distinguishing between acellular cementum and surrounding tissues, the character of the outermost increment was recorded (translucent or opaque). We assume that one increment (annual layer) consists of one pair of opaque (dark) and translucent (bright) lines [44] (Fig 2), and that the opaque (dark) line forms from late autumn to early spring, while the translucent (bright) line forms from spring to autumn [12,13,44]. When conflicting or unambiguous results were recorded, an “indeterminable” character was assigned. We then subjected our results to measures of intra- and inter-observer error, in both cases selecting ten random images that were evaluated by two observers and one observer several months apart. The results are expressed both by the verbal description in the table and by Cohen’s kappa (range of values: 0–1; 0 indicates no agreement, 1 perfect agreement).

thumbnail
Fig 2. Cross-section of a tooth illustrating cementum increments (C = cementum, D = dentine) under polarized transmitted light with lambda plate.

Arrows indicate opaque (dark) cementum increments.

https://doi.org/10.1371/journal.pone.0295757.g002

All data were analysed in Microsoft Excel, Fiji (ImageJ) and Zoner Photo Studio X.

Results

Sample inclusion & exclusion

Good readability of the tissue is essential for evaluating incremental lines for age or season-of-death estimation. Of the 42 teeth, 14 of them (33%) were completely excluded from sample, having zero values in all five thin sections. Mostly, prominent diagenetic changes destroyed the tissue and observing the outermost increment (or other increments) was impossible (Fig 3). In the initial phase, 14 teeth were excluded, and season-of-death estimations could be made for the remaining 28. Only in 4 teeth did all five thin sections have readable cementum enabling the estimation; of the remainder, 3/29 had four readable sections, 8/29 three readable sections, 6/29 two readable sections and 7/29 only one readable section (all results are listed in S1 and S2 Tables).

thumbnail
Fig 3. Prominent diagenetic changes affecting the readability of cementum increments.

Arrows indicate destroyed areas of the acellular cementum (D = dentine, C = cementum).

https://doi.org/10.1371/journal.pone.0295757.g003

Season-of-death estimation

The initial hypothesis was that two different seasons can be observed in the five examined mass graves. In the first examined group of mass graves (A and B, “plague” pits), the season-of-death estimation could be performed for 13 individuals. In 2 cases, the season was rather indeterminable, because of an unclear combination of bright and dark bands. In 11 individuals, the outermost line was bright, which corresponds with the spring/summer season (Table 1 and Fig 4). In the second groups of mass graves (C, D, E), the estimation could be performed for 15 individuals. In 3 of these, the season was indeterminable; in the remaining 12 individuals, the season-of-death estimation was autumn/winter, based on the dark outermost increment (Table 2).

thumbnail
Fig 4. Thin section with readable cementum increments; the arrow indicates the translucent outermost increment.

https://doi.org/10.1371/journal.pone.0295757.g004

thumbnail
Table 1. Season-of-death estimation based on cementochronology, mass graves A and B (“plague” mass graves/pits).

https://doi.org/10.1371/journal.pone.0295757.t001

thumbnail
Table 2. Season-of-death estimation based on cementochronology, mass graves C, D, E “famine” mass graves/pits).

https://doi.org/10.1371/journal.pone.0295757.t002

Inter- and intra-observer error

Ten sections were randomly selected to assess inter-observer agreement in the identification of the outermost increment. One more experienced evaluator and one with no previous experience with cementochronology scored the same photos. Surprisingly, the agreement rate is high, in this case Cohen’s kappa equalling 0.9 (80%) (Table 3).

Ten sections were also randomly selected to assess intra-observer agreement in the identification of the outermost increment. In this case, the estimation is different for the three cases, and the agreement is therefore moderate (0.4). As in [45], the two originally dark lines have now been identified as light (Table 4). We also include a comparison of the level of difficulty in detecting the last increment according to the age of the individual (Fig 5).

thumbnail
Fig 5. Comparison of the difficulty of assessing cementum increments in the context of age estimation and seasonality in individuals of different ages (left: 40–60 years, right: 17–19 years).

https://doi.org/10.1371/journal.pone.0295757.g005

Discussion

Reconstructing seasonality using cementochronology: Technical notes

In our study, the season-of-death in 42 individuals from two stratigraphically distinct levels of mass graves was assessed using cementochronology. Our aim was to better understand the origin of the graves with special regard to the historical documents mentioning catastrophic events and related archaeological findings. We faced challenges in the evaluation that make the method a useful tool, but progress is needed in the protocol for evaluation or sampling. Although many studies indicate the potential of cementochronology for season-of-death estimation in anthropology (e. g., [4649]), very few papers actually provide any results (e. g., [11,13,50]). The results of Wedel [12], with a 99% success rate for the method (distinguishing between two seasons: spring/summer and autumn/winter) sound very promising. Even so, since this study, there have been no other more numerous results that critically discuss evaluation procedures. Yet, there is a plethora of topics in bioarchaeology related to seasonality, and its assessment could enrich other indicators besides, for example, oxygen isotope analysis.

In archaeozoology, cementochronology is widely used for estimating the age and season of death; in bioarchaeology, the method of cementochronology generally meets with less involvement. Typically, the method serves as a useful tool for age-at-death estimation, although some critical studies have pointed out some limitations [7,5154]. However, tooth cementum can be an advantageous indicator, especially for lifelong continuous tissue growth or generally good preservation of teeth in the archaeological record. Incremental lines in acellular cementum are considered to correlate with seasonal growth in most species [44]. Histological differences in the orientation of collagen fibres and their degree of mineralization are thought to explain the alternating light and dark increments in dental cementum [44,55]. Seasonal changes in diet, nutritional intake, hormonal cycles, sunlight or mastication are considered to be related factors [56]. What remains unclear is why, contrary to wild mammals, most modern humans and domestic mammals lacking a seasonal diet [45,57,58] still have regularly deposited cementum increments. The unclear etiology of the increments is often cited as one of the reasons why the method is not more widely used.

The application of cementochronology to estimate seasonality on human teeth is associated with several problematic factors, which we faced in this study. The terminology describing the outermost increment differs among researchers–the terms translucent and opaque are used (but with respect to the type of microscope in which a light increment may appear dark and vice versa), but also growth zone/annulus and slow/rapid growth [55]. In our study, the terms bright and dark increments are used. Human teeth are sometimes many times smaller than animal teeth, which makes the analysis more difficult. Observation can also be complicated by the large number of increments in individuals of advanced calendar age as well (Fig 5) [7]. This is a problem in anthropology [7], but also in archaeozoology [59]. Also, Wedel & Wescott [13] included only relatively young individuals in their study. There is also a relationship between the number of increments and tissue readability [7,60], with higher numbers of increments increasing the difference between chronological and estimated age. Regarding the archaeological material, the application of the method is moreover limited by diagenetic changes, possibly leading to a high percentage of discarded samples [61]. Another issue is related to the difference in detecting the two types of increments (translucent and opaque)—even in animals; the opaque line is more difficult to detect, because it is very thin [62]. In our case, it was also more difficult to identify the dark increment, which very often rather merged into the edge of the tooth, or remnants of fibres on it. To date, there are only a very small number of studies, few of which discuss inter- and intra-observer error, too. If this method is to be adopted for routine practice, it is necessary to expand the results. Compared to Wedel’s [12] study, Ralston [58] showed 61.5 to 71.2% correct season identification, and also 28.8% inter-observer error. Similarly, Thomas & Corby [63] demonstrated 58.3 to 73.9% correct season identification. Moreover, both studies faced the problem of identifying the true outermost increment, due to optical edge effects and inconsistencies in the outermost increments’ appearance. We too faced such issues, and agree with Thomas & Corby [63] that a modified protocol enhancing the outermost increment is strongly needed.

Good tissue visibility is critical, but may already be affected by age. This can be partly balanced by the preparation protocol, which should ensure the best possible result–with no cutting and grinding marks, and especially no glue marks on the tooth margin. Technical defects can be considered as reason for excluding samples [62]. However, it may still be the case that remnants of other tissues persisting on the tooth margin and partially cover the cementum. Finally, although some studies have shown no influence of periodontal disease on the results of age-at-death estimation [64,65], a no longer functional tooth can have disrupted apposition of the tissue, which can affect the results. Some pathological conditions (e. g. a fractured tooth, coronal destruction, caries) do not affect the results of the age estimation using cementochronology, but this does not apply to damaged desmodontium and periodontium, where they may render the tooth unusable for analysis or cause greater variation in age estimation [66]. Problems related to the preservation of tissues are unavoidable in an archaeological context; in bioarchaeology and forensic anthropology, however, the potential of season-of-death estimation is valuable and deserves further research. This may include a validation study in documented samples, larger samples balancing a high number of excluded teeth with diagenetic changes, or advanced technical steps such as measuring the outermost increment using measurement tools or considering the intra- and inter-observed error. Also, developing protocols for recording and presenting increment data would prevent researchers from subjectivity and counting false lines (caused by using different microscopy).

Seasonality of famine and plague epidemics

The historical and archaeological evidence from Kutná Hora-Sedlec documents that some part of the population was impacted by serious catastrophes–episodes of famine (1318) and plague (1348–1350). Both events have specific seasonal dynamics, but their study in archaeology encompasses several limitations. Kutná Hora-Sedlec site is exceptional in terms of the number of individual and mass graves excavated during various excavation campaigns between 2010 and 2018, as part of restoration work on its famous ossuary chapel [15,16]. Individual graves (single inhumations) dated by radiocarbon dating belong to the 13-14th and 14-16th centuries, while the mass graves are dated to the 14th century [68]. This period was a time of dramatic climate change for Bohemia (and Europe in general). While the 11th-13th centuries were relatively warm and stable, around 1300 temperatures fell, and crop failures and famines occurred in the temperate zone. This gave rise to various negative phenomena, such as the mass migrations of people to cities and the emergence of diseases. The causes of the considerable die-off during the 14th century at Kutná Hora-Sedlec are written in historical documents, which mention the victims of the famine of 1318. The so-called Great Famine is the best-documented famine, impacting northern and western Europe between 1315 and 1322 [67]. As for seasonality, famine episodes mostly occurred during the winter or early spring months [2830]. Since famines were common among certain social groups in the Middle Ages, the famine must have been intense for it to appear in the chronicles, which are an invaluable source of information today [68]. Otherwise, it is true that hunger or the diet of the poor appeared rather rarely in medieval written sources [26]. Indeed, it is also very difficult to document famines on the basis of archaeological and osteological evidence alone [69]. Basically, two types of famine can be distinguished–acute and chronic. However, at the level of the skeletal remains, it is difficult to decide if an acute or chronic famine, or malnutrition, can be identified. In terms of osteological evidence, some signs associated with malnutrition (e. g. cribra orbitalia, Harris lines) [70] can be observed. These indicators are however non-specific, so they cannot be directly associated with famine (or e. g. infectious disease or other stress). The situation is also relativized by the osteological paradox, which draws attention to the important fact that an individual who dies of an acute disease, given that it is on the skeleton for a short duration, does not develop stress manifestations, and might be considered as healthy [71].

As regards the plague epidemic (1348–1350), we are faced with a lack of historical sources [21]. The main reason is that the plague epidemic was not as severe in Czechia as it was in, for example, in Italy or France. According to František Pražský, the plague reached Austria in 1349 and also began in Bohemia. Also, the first attack of the plague was pushed out by “fresh and cold air” [20]. This would be in line with information about known plague seasonal dynamics. On the basis of being related to latitude, temperature and precipitation, the epidemic growth of plague outbreaks was positive between 11.7°C and 21.5°C, and second pandemic outbreaks peaked in the summer [23,32]. The problem is that there is no completely accurate meteorological data, so we have to work with model data. Mild winters or short-term fluctuations are difficult to reflect. While there are possibilities to distinguish between famine and plague graves on the basis of osteological evidence alone (without archaeogenetic analyses), options are limited. The situation is further complicated by the fact that there are situations in which a population faced the famine and plague epidemics at the same time [21,72]. Catastrophes such as wars or disease may lead to a change in funerary rite, and also result in specific features in the skeletal population [73]. At the level of a group of individuals (or population), a specific (attritional) mortality profile that differs from the normal one can support the identification of a catastrophic event [7375]. Compared to attritional mortality, where the highest risk of death is for infants and older adults, catastrophic mortality should mimic the age structure of the living population, having a high risk of death in all age categories [7577]. In this case too, great attention should be given to the choice of age estimation methods, e.g. an excess number of middle-aged adults may affect the mortality curve and its interpretation [75]. Interestingly, in the study of Lépinau et al. [17], the examined mass graves (numbers 765, 853 and 1578, in our study: C, D, E) were assigned more to plague mortality, on the basis of the age structure and pattern of skeletal lesions. Trying to find a specific age pattern for catastrophic events (mainly plague epidemics) is the subject of research [7578], the mortality pattern of the 14th century Black Death was not very different from normal medieval mortality, and plague preferentially killed older adults and individuals with poor health. Together with our findings it is clear that additional research including archaeogenetic analyses [79] is needed to clear up the nature of mass graves in Kutná Hora-Sedlec. In conclusion, the results of season-of-death estimation using cementochronology identified two different seasons in the investigated mass graves. This makes the method a useful tool, for example, in cases where a question related to seasonality is being addressed on archaeological material.

Conclusions

Cementochronology has so far rarely been applied in bioarchaeology to estimate the season-of-death. The reasons may vary, but in the context of explanatory value, the method can potentially address issues related to seasonality. In our study, we managed to identify two different burial periods which can hypothetically be associated with catastrophic events known from historical and archaeological records. It is impossible to say whether a specific catastrophic event was directly involved based on season-of-death estimation alone. Still, our findings contribute to understanding the differences between the two groups of graves, and to discussing the results with a knowledge of the seasonality of famine figand plague. Thus, cementochronology could also be applied in bioarchaeology, but with inevitable (diagenetic changes) and avoidable limitations (selection of younger and middle-aged adults, emphasis on standardization of procedures, maximum removal of observer subjectivity).

Supporting information

S1 Table. Readability index values for the assessed thin sections.

https://doi.org/10.1371/journal.pone.0295757.s001

(DOCX)

S2 Table. Season-of-death estimations for all mass graves.

https://doi.org/10.1371/journal.pone.0295757.s002

(DOCX)

Acknowledgments

The authors would like to thank Alastair Millar for language advice, Dan Bursák for his help with creating maps, and to two anonymous reviewers for their helpful and inspiring comments.

References

  1. 1. Kwiecien O, Braun T, Brunello CF, Faulkner P, Hausmann N, Helle G, et al. What we talk about when we talk about seasonality–A transdisciplinary review. Earth-Sci Rev. 2022;225: 103843.
  2. 2. Knudson KJ, Aufderheide AE, Buikstra JE. Seasonality and paleodiet in the Chiribaya polity of southern Peru. J Archaeol Sci. 2007;34: 451–462.
  3. 3. Hallmann N, Burchell M, Schöne BR, Irvine GV, Maxwell D. High-resolution sclerochronological analysis of the bivalve mollusk Saxidomus gigantea from Alaska and British Columbia: techniques for revealing environmental archives and archaeological seasonality. J Archaeol Sci. 2009;36: 2353–2364.
  4. 4. Rendu W. Hunting behavior and Neanderthal adaptability in the Late Pleistocene site of Pech-de-l’Azé I. J Archaeol Sci. 2010;37: 1798–1810.
  5. 5. Schmaus TM, Chang C, Tourtellotte PA. A model for pastoral mobility in Iron Age Kazakhstan. J Archaeol Sci Rep. 2018;17: 137–143.
  6. 6. Weinand DC. Seasonality examination of White-tailed deer increment structures from coastal Georgia, USA. Archaeozoologia. 2000;11: 65–78.
  7. 7. Bertrand B, Cunha E, Bécart A, Gosset D, Hédouin V. Age at death estimation by cementochronology: Too precise to be true or too precise to be accurate? Am J Phys Anthropol. 2019;169: 464–481. pmid:31049939
  8. 8. Klevezal GA. Recording Structures of Mammals: Determination of Age and Reconstruction of Life History. 1st ed. Routledge; 2017. https://doi.org/10.1201/9780203741146
  9. 9. Bertrand B, Oliveira-Santos I, Cunha E. Cementochronology: a validated but disregarded method for age at death estimation. In: Adserias-Garriga J, editor. Age Estimation: A Multidisciplinary Approach. Elsevier; 2019. pp. 169–186. Available: https://linkinghub.elsevier.com/retrieve/pii/B9780128144916000121.
  10. 10. Pike-Tay A. Red deer hunting in the Upper Paleolithic of south-west France: a study in seasonality. BAR Publishing; 1991.
  11. 11. Klevezal GA, Shishlina NI. Assessment of the Season of Death of Ancient Human from Cementum Annual Layers. J Archaeol Sci. 2001;28: 481–486.
  12. 12. Wedel VL. Determination of Season at Death Using Dental Cementum Increment Analysis. J Forensic Sci. 2007;52: 1334–1337. pmid:17868272
  13. 13. Wedel VL, Wescott DJ. Using dental cementum increment analysis to estimate age and season of death in African Americans from an historical cemetery in Missouri. Int J Paleopathol. 2016;15: 134–139. pmid:29539548
  14. 14. Brzobohatá H, Frolík J, Zazvonilová E. Bioarchaeology of Past Epidemic- and Famine-Related Mass Burials with Respect to Recent Findings from the Czech Republic. Interdisci Archaeol. 2019;10: 79–87. https://doi.org/10.24916/IANSA.2019.1.6.
  15. 15. Frolík J. Pohřbívání ve vrcholném středověku a v novověku na Chrudimsku, Pardubicku a Kolínsku. Archaeol Hist. 2017;42: 187–205.
  16. 16. Horák J, Weiss-Krejci E, Frolík J, Velímský F, Šmejda L. The Cemetery and Ossuary at Sedlec near Kutná Hora: Reflections on the Agency of the Dead. Interdisciplinary Explorations of Postmortem Interaction. Cham: Springer Nature; 2022. pp. 269–295.
  17. 17. Lépinau A de, Castex D, Brzobohatá H, Frolík J, Velímský F, Brůžek J, et al. Entre peste et famine: caractérisation d’une crise de mortalité par l’étude de trois sépultures multiples du site de Kutná Hora–Sedlec (République tchèque, XIVe siècle). Bull Mém Société D’Anthropologie Paris. 2021;33.
  18. 18. Brzobohatá H, Frolík J, Velímský F. Wealth or Just Job Seekers: Medieval Skeletal Series from Kutná Hora-Sedlec (Czech Republic) with a Notable Surplus of Men. Interdisci Archaeol. 2023; 12: 79–92.
  19. 19. FRB IV: Petra Žitavského Kronika zbraslavská (Chronicon Aulae Regiae). Pragae: Fontes rerum Bohemicarum IV; 1884.
  20. 20. Maur E. Morová epidemie roku 1380 v Čechách. Hist Demogr. 1986;10: 37–71.
  21. 21. Nodl M. Impacts of the plague epidemic on the Kingdom of Bohemia in the second half of the XIVth century and at the beginning of the XVth century. Stud Hist Gedanensia. 2021;Tom 12 (2021)/2:s: 44–62.
  22. 22. Eckert EA. Seasonality of Plague in Early Modern Europe: Swiss Epidemic of 1628–1630. Clin Infect Dis. 1980;2: 952–959. pmid:7012992
  23. 23. Krauer F, Viljugrein H, Dean KR. The influence of temperature on the seasonality of historical plague outbreaks. Proc R Soc B Biol Sci. 2021;288: 20202725. pmid:34255997
  24. 24. Devereux S. Famine in the twentieth century. Institute of Development Studies, University of Sussex; 2000. Available: https://www.ids.ac.uk/publication/famine-in-the-twentieth-century.
  25. 25. Horocholyn K, Brickley M, Department of Anthropology, McMaster University, Hamilton, Ontario, Canada. Pursuit of Famine: Investigating Famine in Bioarchaeological Literature. Bioarchaeology Int. 2017;1: 101–115.
  26. 26. McCleery I. Getting Enough to Eat: Famine as a Neglected Medieval Health Issue. The Sacred and the Secular in Medieval Healing: Sites, Objects and Texts. London and New York: Routledge; 2016.
  27. 27. Fagan B. Malá doba ledová. Jak klima formovalo dějiny v letech 1300–1850. Praha: Academie; 2007.
  28. 28. Healey J. Famine and the female mortality advantage: sex, gender and mortality in northwest England, c. 1590–1630. Contin Change. 2015;30: 153–192.
  29. 29. Steinbachová L. Demografický vývoj za hladomoru v letech 1771–1772. Hist Demogr. 2001;25: 101–129.
  30. 30. Grulich J. Populační vývoj a životní cyklus venkovského obyvatelstva na jihu Čech v 16. až 18. století. České Budějovice: Jihočeská univerzita České Budějovice; 2008.
  31. 31. Christakos G, Olea RA, Yu H-L. Recent results on the spatiotemporal modelling and comparative analysis of Black Death and bubonic plague epidemics. Public Health. 2007;121: 700–720. pmid:17544041
  32. 32. Krauer F. Temperature, toponyms and thresholds: A modelling approach to understanding the spread of plague during the second pandemic. Dissertation presented for a degree of Philosophiae Doctor (PhD), Faculty of Mathematics and Natural Sciences. 2020.
  33. 33. Welford MR, Bossak BH. Validation of Inverse Seasonal Peak Mortality in Medieval Plagues, Including the Black Death, in Comparison to Modern Yersinia pestis-Variant Diseases. Carter DA, editor. PLoS ONE. 2009;4: e8401. pmid:20027294
  34. 34. Batchelor PA, Sheiham A. Grouping of tooth surfaces by susceptibility to caries: a study in 5–16 year-old children. BMC Oral Health. 2004;4: 2. pmid:15511295
  35. 35. İşcan MY, Loth SR, Wright RK. Metamorphosis at the sternal rib end: A new method to estimate age at death in white males. Am J Phys Anthropol. 1984;65: 147–156. pmid:6507605
  36. 36. Lovejoy CO. Dental wear in the Libben population: Its functional pattern and role in the determination of adult skeletal age at death. Am J Phys Anthropol. 1985;68: 47–56. pmid:4061601
  37. 37. Brooks S, Suchey JM. Skeletal age determination based on the os pubis: A comparison of the Acsádi-Nemeskéri and Suchey-Brooks methods. Hum Evol. 1990;5: 227–238.
  38. 38. Schmitt A. Une nouvelle méthode pour estimer l’âge au décès des adultes à partir de la surface sacro-pelvienne iliaque. Bull Mém Société Anthropol Paris. 2005;17: 89–101.
  39. 39. Calce SE. A new method to estimate adult age-at-death using the acetabulum. Am J Phys Anthropol. 2012;148: 11–23. pmid:22331613
  40. 40. Colard T, Bertrand B, Naji S, Delannoy Y, Bécart A. Toward the adoption of cementochronology in forensic context. Int J Legal Med. 2018;132: 1117–1124. pmid:25773917
  41. 41. Wittwer-Backofen U, Gampe J, Vaupel JW. Tooth cementum annulation for age estimation: Results from a large known-age validation study. Am J Phys Anthropol Off Publ Am Assoc Phys Anthropol. 2004;123: 119–129. pmid:14730646
  42. 42. Klevezal GA, Kirillova IV, Shishlina NI, Sokolov AA, Trunova YE. Growth layers in tooth dentin and cementum: problems and perspectives of their use in the study of fossil and subfossil mammal remains including humans. Doc Archaeobiol. 2006;4: 113–124.
  43. 43. Pike-Tay A, Cosgrove R, Garvey J. Systematic seasonal land use by late Pleistocene Tasmanian Aborigines. J Archaeol Sci. 2008;35: 2532–2544.
  44. 44. Lieberman DE. The Biological Basis for Seasonal Increments in Dental Cementum and Their Application to Archaeological Research. J Archaeol Sci. 1994;21: 525–539. https://doi.org/10.1006/jasc.1994.1052.
  45. 45. Meckel LA. The utility of dental cementum increment analysis for estimating season-of-death in naturally decomposed skeletons. M. A. Thesis, Texas State University. 2016.
  46. 46. Kaur P, Astekar M, Singh J, Arora K, Bhalla G. Estimation of age based on tooth cementum annulations: A comparative study using light, polarized, and phase contrast microscopy. J Forensic Dent Sci. 2015;7: 215. pmid:26816462
  47. 47. Le Cabec A, Tang NK, Ruano Rubio V, Hillson S. Nondestructive adult age at death estimation: Visualizing cementum annulations in a known age historical human assemblage using synchrotron X-ray microtomography. Am J Phys Anthropol. 2019;168: 25–44. pmid:30431648
  48. 48. Cerrito P, Bailey SE, Hu B, Bromage TG. Parturitions, menopause and other physiological stressors are recorded in dental cementum microstructure. Sci Rep. 2020;10: 5381. pmid:32214148
  49. 49. Perrone V, Gocha TP, Randolph-Quinney P, Procopio N. Tooth Cementum Annulation: A Literature Review. Forensic Sci. 2022;2: 516–550.
  50. 50. Wilson N, Koel-Abt K, Schmidt K. The Analysis of Tooth Cementum for the Histological Determination of Age and Season at Death of Active Duty Military Members. Ph. D. Thesis, Uniformed Services University of the Health Sciences, Bethesda, Maryland 20814. 2017.
  51. 51. Obertová Z, Francken M. Tooth Cementum Annulation Method: Accuracy and Applicability. In: Koppe T, Meyer G, Alt KW, Brook A, Dean MC, Kjaer I, et al., editors. Frontiers of Oral Biology. Basel: KARGER; 2009. pp. 184–189. https://doi.org/10.1159/000242415 pmid:19828994
  52. 52. Grosskopf B, McGlynn G. Age diagnosis based on incremental lines in dental cementum: A critical reflection. Anthropol Anz. 2011;68: 275–289. pmid:21905417
  53. 53. Naji S, Colard T, Blondiaux J, Bertrand B, d’Incau E, Bocquet-Appel J-P. Cementochronology, to cut or not to cut? Int J Paleopathol. 2016;15: 113–119. pmid:29539545
  54. 54. Gualdi-Russo E, Saguto I, Frisoni P, Neri M, Rinaldo N. Tooth Cementum Thickness as a Method of Age Estimation in the Forensic Context. Biology. 2022;11: 784. pmid:35625512
  55. 55. Lam YM. Variability in cementum deposition in springbok from the South African Cape. J Archaeol Sci. 2008;35: 1062–1073.
  56. 56. Tang N, Le Cabec A, Antoine D. Dentine and cementum structure and properties. A companion to dental anthropology. John Wiley & Sons; 2015. pp. 204–222.
  57. 57. Cool SM, Forwood MR, Campbell P, Bennett MB. Comparisons between bone and cementum compositions and the possible basis for their layered appearances. Bone. 2002;30: 386–392. pmid:11856646
  58. 58. Ralston CE. Dental cementum increment analysis and estimating season at death in humans. M. Sc. Thesis, Boston University. 2016. Available from: https://www.proquest.com/docview/1767406755?pq-origsite=gscholar&fromopenview=true.
  59. 59. Wall-Scheffler CM, Foley RA. Digital cementum luminance analysis (DCLA): a tool for the analysis of climatic and seasonal signals in dental cementum. Int J Osteoarchaeol. 2008;18: 11–27.
  60. 60. Zazvonilová E, Velemínský P, Černíková A, Danielisová A, Brůžek J. Protocol matters: A need for standardized procedure in cementochronology. Forensic Sci Int. 2022;340: 111439. pmid:36063738
  61. 61. Roksandic M, Vlak D, Schillaci MA, Voicu D. Technical note: Applicability of tooth cementum annulation to an archaeological population. Am J Phys Anthropol. 2009;140: 583–588. pmid:19639640
  62. 62. Sánchez-Hernández C, Gourichon L, Pubert E, Rendu W, Montes R, Rivals F. Combined dental wear and cementum analyses in ungulates reveal the seasonality of Neanderthal occupations in Covalejos Cave (Northern Iberia). Sci Rep. 2019;9: 14335. pmid:31586143
  63. 63. Thomas SC, Corby JH. The Histological Analysis of Tooth Cementum Annulations to Determine Age and Season of Death of Active Duty Military Personnel: Comparative Study Using Brightfield, Polarizing and Phase Contrast Microscopy. Ph. D. Thesis, Uniformed Services University of the Health Sciences, Bethesda, Maryland 20814. 2020. Available from: https://apps.dtic.mil/sti/citations/trecms/AD1182859.
  64. 64. Bertrand B, Robbins Schug G, Polet C, Naji S, Colard T. Age-at-death estimation of pathological individuals: A complementary approach using teeth cementum annulations. Int J Paleopathol. 2016;15: 120–127. pmid:29539546
  65. 65. de Broucker A, Colard T, Penel G, Blondiaux J, Naji S. The impact of periodontal disease on cementochronology age estimation. Int J Paleopathol. 2016; 6.
  66. 66. Kagerer P, Grupe G. Age-at-death diagnosis and determination of life-history parameters by incremental lines in human dental cementum as an identification aid. Forensic Sci Int. 2001;118: 75–82. pmid:11343858
  67. 67. Kershaw I. The great famine and agrarian crisis in England 1315–1322. Past Present. 1973;59: 3–50.
  68. 68. Antonín R, Antonín Malaníková M. Když se jídla nedostává. České středověké reflexe hladových let v evropském kontextu (when one is not fed—the medieval relations about “the years of hunger” in the central European region). Historia naturalna Jedzenia- Miedzy atykem a XIX wiekem. Gdańsk: Muzeum Historyczne Miasta Gdanska; pp. 70–81.
  69. 69. Morgan J. The Invisible Hunger: Is Famine Identifiable from the Archaeological Record? Antrocom Online J Anthropol. 2013;9: 115–129.
  70. 70. McFadden C, Oxenham MF. A paleoepidemiological approach to the osteological paradox: Investigating stress, frailty and resilience through cribra orbitalia. Am J Phys Anthropol. 2020;173: 205–217. pmid:32578874
  71. 71. Wood JW, Milner GR, Harpending HC, Weiss KM, Cohen MN, Eisenberg LE, et al. The Osteological Paradox: Problems of Inferring Prehistoric Health from Skeletal Samples [and Comments and Reply]. Curr Anthropol. 1992;33: 343–370.
  72. 72. Demont P. The Causes of the Athenian Plague and Thucydides. Thucydides between History and Literature. Berlín: De Gruyter; pp. 73–87.
  73. 73. Margerison BJ, Knüsel CJ. Paleodemographic comparison of a catastrophic and an attritional death assemblage: Paleodemography of a Catastrophic Cemetery. Am J Phys Anthropol. 2002;119: 134–143. pmid:12237934
  74. 74. Lubinski PM, O’brien CJ. Observations on Seasonality and Mortality from a Recent Catastrophic Death Assemblage. J Archaeol Sci. 2001;28: 833–842.
  75. 75. Gowland RL, Chamberlain AT. Detecting plague: palaeodemographic characterisation of a catastrophic death assemblage. Antiquity. 2005;79: 146–157.
  76. 76. Keckler CN. Catastrophic mortality in simulations of forager age-at-death: where did all the humans go. Integrating Archaeol Demogr Multidiscip Approaches Prehist Popul Cent Archaeol Investig Occas Pap. 1997;24: 205–228.
  77. 77. Paine RR. If a population crashes in prehistory, and there is no paleodemographer there to hear it, does it make a sound? Am J Phys Anthropol. 2000;112: 181–190. pmid:10813701
  78. 78. DeWitte SN. Age patterns of mortality during the Black Death in London, A.D. 1349–1350. J Archaeol Sci. 2010;37: 3394–3400. pmid:21572598
  79. 79. Gutsmiedl-Schümann D, Päffgen B, Schwarzberg H, Keller M, Rott A, Harbeck M. Digging up the plague: A diachronic comparison of aDNA confirmed plague burials and associated burial customs in Germany. Praehistorische Z. 2018;92: 405–427.