Radiocarbon samples. The analysis is based primarily on radiocarbon dating of charcoal samples collected from archaeological features. These include both previously published samples and newly acquired, unpublished data that are presented and evaluated here for the first time (see Supplementary Table S1). Most of the earlier data had never been subjected to detailed analysis, with the exception of the earliest series of dates (2007–2010) obtained using AMS at the now-defunct Erlangen laboratory in Germany25. Subsequent research was carried out at laboratories in Kraków, Poland (code: MKL), using the LSC technique (years 2010–2016), and in Poznań, Poland (code: POZ), using AMS (since 2017). These methods are considered mutually consistent, and their results statistically comparable34, 35, 36.
A total of 79 radiocarbon determinations were made available for analysis. Five of these—obtained from animal bone samples found at Bojná III—were not included in the integrated analysis and are discussed separately. The integrated analysis encompasses 74 radiocarbon dates, including samples of charcoal (68), slag (2), wood (2), human bone (1), and an animal tooth (1), all originating from the five sites that constitute the core of the Bojná agglomeration.
The dataset was thoroughly evaluated with regard to the origin and method of sample collection, using excavation field documentation provided by the Institute of Archaeology of the Slovak Academy of Sciences (AÚ SAV) in Nitra. This was done to exclude samples that were unsuitable for analysing the chronology of the early medieval agglomeration, or those clearly affected by post-depositional disturbance.
A total of 20 dates were classified as clearly inconsistent with the chronological framework relevant to this study, or demonstrably incompatible with the stratigraphic context of the sampled features. In two cases, the anomalous dating was attributed to contamination by roots or mould; in one instance, to prolonged atmospheric exposure of the sample, as established through excavation records; and in another, an implausibly early date derived from an iron-smelting furnace, which is likely to have been affected by alterations in carbon content due to high-temperature processing. The remaining 16 excluded dates correspond to prehistoric or post-medieval periods, well outside the broadest expected range of the early medieval occupation. Although these fall beyond the scope of the current article, it is important to note that they are likely to be accurate, and to reflect earlier or later phases of human activity within the agglomeration area; this is supported by the presence of associated finds.
During detailed analysis, three further dates from the Bojná VI site were deemed unreliable: two slag dates were rejected following failed agreement tests, and one was excluded based on the dating laboratory’s own recommendation, due to insufficient nitrogen and carbon content.
In total, 51 radiocarbon dates were ultimately included in the final analysis. It is important to emphasise that none of the excluded dates originated from construction elements of early medieval features; rather, they derived exclusively from small charcoal fragments recovered from fill layers, or from mound and rampart deposits.
In the integrated analysis, five radiocarbon dates from animal bones recovered from the fills of two burial mounds at the Bojná III site were excluded, due to concerns regarding their stratigraphic context—and, by extension, their origin and pre-depositional treatment—as well as broader doubts surrounding the reliability of radiocarbon dates from bone material without proper pretreatment37. These concerns arose following the receipt of results that contradicted the expected chronology of the contexts from which the samples were reportedly taken. The decision to exclude these five animal bone samples from the analysis was further reinforced by a similarly problematic radiocarbon result from a human bone (Poz-20138/22/4) from the Bojná VI site, which also proved inconsistent with the archaeological chronology. In both cases, the discrepancies were relatively minor, and the resulting dates still fell within the broader chronological framework of the Bojná agglomeration.
A detailed investigation of the Bojná VI sample revealed insufficient nitrogen and carbon content, thus rendering the date unreliable38, 39. In the case of the Bojná III animal bone samples, however, measurements of carbon and nitrogen content, atomic carbon–nitrogen ratios, and collagen quality provided values that were consistent with reliable radiocarbon dating. This opened up alternative interpretive possibilities, suggesting that the dates themselves may be accurate, and that it was the archaeological context that was misidentified—possibly the bone disposition was unrelated to funerary activity. Given that radiocarbon dating of bone is inherently complex and requires a methodological approach grounded in sound archaeological understanding40, additional analyses using a larger number of samples from Bojná III are planned for the near future.
Analysis and modelling. The calibration, analysis, and modelling of the data were conducted using OxCal version 4.4.4.17528. Bayesian chronological modelling functions provided by the software were applied in interpreting the results. Bayesian statistics offer an explicit, probabilistic framework for integrating different types of evidence in order to estimate the timing of past events and quantify the uncertainty associated with these estimates. This approach—now a standard technique in the construction of chronological models based on radiocarbon datasets41, 42—enables, in short, the refinement and improved precision of dating by combining statistical outputs with archaeological knowledge.
In assessing the evidentiary weight of radiocarbon data for constructing chronological models, a fundamental principle was adopted: dates obtained from charred structural elements are of greater value, as they are likely to reflect—at least approximately—the structure’s period of use. In contrast, isolated charcoal (bone, wood, etc.) fragments retrieved from the fills of archaeological features or stratified layers must be treated with greater caution. As a rule, such samples can only provide a TPQ for the formation of specific features or layers, and more broadly may indicate the general chronology of a site, or highlight phases of intensified human activity44.
In all such cases, there remains the possibility that the charcoal sample was not directly associated with the feature in question, having entered it through post-depositional processes. In the case of hearths, for example, it is possible that fuelwood of considerable age and unknown origin was used. The problem of so-called “old wood” (inbuilt age) can also arise when dating structural wooden elements43. However, in these cases, the uncertainty can be mitigated through appropriate modelling, and outliers can be identified and excluded. This principle has already been successfully applied in the development of a chronological model for the St. George’s Rotunda in Nitrianska Blatnica, part of the broader Bojná agglomeration20, 31. Specific issues related to individual models are discussed in detail below, and in the Supplementary Materials.
Markov chain Monte Carlo (MCMC) analysis was used for almost all multi-parameter Bayesian analyses performed in OxCal. The only exception is the use of the Combine() function, which reduces the problem to a single independent parameter so that it can be solved analytically. In OxCal, a model is considered plausible when its Amodel value exceeds 60. The Amodel parameter represents a measure of agreement between the model’s prior structure and the observational data (likelihood)28.
In our modelling, we adopted the principle of striving for models with high internal consistency, aiming for Amodel values close to or exceeding 100. A value in this range indicates a statistically well-fitting model that is likely to contain no major outliers, and whose structure—such as phasing and boundaries—accurately reflects the underlying data. Amodel values exceeding 150 must, however, be interpreted with caution, as they may indicate overfitting—that is, excessive precision due to overly restrictive assumptions, underestimation of uncertainties, or artificially tight constraints. Such values do not automatically indicate an error, but should be carefully assessed in the broader context of the archaeological and statistical evidence.
Aggregate radiocarbon data for individual sites were also analysed using a kernel density estimation (KDE) model. This method allows a smoothed representation of the temporal distribution of dates while incorporating the full structure of the Bayesian model, including phases, boundaries, and outliers. By applying KDE, it becomes possible to estimate periods of peak activity within a given phase, and to compare these activity patterns across different sets of dates45.
Chronological models
Bojná IV. The only radiocarbon date obtained from the rampart at the Bojná IV site—Poz-16305/20/13 (1,525 ± 30 BP)—is insufficient to construct a standalone model. Its interpretation is therefore only possible in conjunction with the analytical model developed for the Bojná III site (see Supplementary Fig. S5), as described in the Results and Discussion section.
Bojná III. For the Bojná III site, the radiocarbon dataset was grouped into two clusters already at the initial stage of analysis. The first cluster consists of dates obtained from scattered charcoal samples collected from the fortification structures (ramparts) and from the cultural layer located beneath one of the barrows. These dates fall within the range of 1,650 ± 40 to 1,470 ± 30 BP. The second cluster includes dates derived from the construction elements of the barrows, ranging from 1,210 ± 30 to 1,190 ± 30 BP (or 1,208 ± 18 combined BP, to 1,190 ± 30 BP), as well as a single date from the fill of the ditch (1,195 ± 30 BP).
The dates from the first cluster fall within a relatively broad range of 1,650 ± 40 to 1,470 ± 30 BP (Supplementary Fig. S5). These should not be interpreted as indicators of the fortifications’ construction period, as the samples from which they were obtained do not originate from structural elements of the rampart. Most likely, these charcoal fragments entered the rampart fill along with surrounding soil used in its construction. At best, these dates can serve as TPQs for the erection of the fortifications (Supplementary Fig. S1), while nonetheless providing clear confirmation that the defensive structures of Bojná III are of early medieval origin. This latter point is particularly significant, as it effectively rules out the possibility that the rampart represents reused prehistoric fortifications.
The second data cluster consists of four radiocarbon dates obtained from two burial mounds. These derive from samples collected from the cremation pyre and structural elements of the mounds, and thus provide an approximate indication of the time of their construction. The dates are tightly grouped within a narrow range of 1,210 ± 30 to 1,190 ± 30 BP (or 1,208 ± 18 combined BP, to 1,190 ± 30 BP). Although this cluster demonstrates remarkable consistency and dating precision, its utility for establishing an exact absolute chronology is limited. The dates fall within the so-called “early medieval plateau” on the calibration curve, which constrains resolution and permits only a broad estimate for the construction of these features—somewhere between the late 8th and late 9th centuries.
A third, separate date—Poz-186685 (1,195 ± 30 BP)—was obtained from the fill of the hillfort’s inner ditch. This date is nearly identical to those from the younger of the two barrows. While it does not directly contribute to narrowing the chronology, it may tentatively be interpreted as confirmation of some form of human activity at the site during the period of barrow construction; this possibility is also supported by the associated archaeological material.
In the basic chronological model (Supplementary Fig. S2), developed on the basis of many years of archaeological research, it was assumed that the fortifications of the hillfort are contemporaneous with the initial phase of settlement, while the barrows were constructed in the final stage of site use—either during or after the decline of occupation. Since no traces of earlier cultural layers or associated artefacts were recorded beneath or within the rampart body, it is reasonable to presume that the rampart was among the earliest structures at the site, erected on terrain that was previously unoccupied or only lightly utilised. Over time, continued use of the hillfort led to the formation of a cultural layer. In this context, dates obtained from unstratified charcoal of uncertain origin found within the archaeologically sterile rampart fill must be treated as TPQs for the fortifications’ construction. The most reliable and informative of these is sample Poz-186694 (1,470 ± 30 BP), taken from the outer face of the rampart and found in direct association with ceramic material. This date may therefore serve as a boundary marker for the onset of occupation at the site. Similarly, sample Poz-186596 (1,545 ± 30 BP)—derived from charcoal found in the cultural layer beneath Barrow 2—should also be considered a TPQ for the burial mound’s construction. The stratigraphic sequence is supported by the presence of mixed cultural material from the occupation layer within the barrow fill, and by its absence in the surrounding area. The lack of settlement-related deposits in the immediate vicinity of the mounds indirectly suggests that intensive human activity at the site had ceased by the time the barrows were erected.
The model made it possible to establish a TPQ for the fortifications’ construction and the onset of early medieval settlement at the Bojná III site (Supplementary Fig. S2), which corresponds approximately to the point defined by date Poz-186694 (1,470 ± 30 BP). Due to the lack of further data, the duration of the settlement cannot be precisely determined; however, the approximate end of site use is marked by the barrows’ construction dates.
Between these chronological boundaries lies a “temporal gap” in the radiocarbon record, which most likely corresponds to the actual period during which the Bojná III hillfort was in use. On this basis, the maximum span of settlement activity may be estimated as extending from the mid-7th to the early 9th centuries.
The absence of radiocarbon dates for the period between the mid-7th and early 9th centuries does not reflect a true archaeological hiatus. This is clearly contradicted by the excavated material, which is consistently dated within that range21. The archaeological assemblage indicates uninterrupted activity from the mid-7th century to the early 9th century—roughly matching the period for which radiocarbon data are lacking. The proposed model was subsequently tested through simulation (Supplementary Fig. S3), which confirmed its validity.
Test. The preliminary test models simulate radiocarbon dates for the period between the early 7th and early 9th centuries. Several simulations were conducted using calendar years from AD 600 to 800, at 25-year intervals and with uncertainties of 50, 40, and 30 years—thus reflecting the precision typical of modern radiocarbon laboratories. One such example is presented in Supplementary Fig. S3.
In the second test, for which two models were constructed, it was assumed that the charcoal fragments found within the rampart represent traces of early medieval human activity at the site—activity that directly preceded the fortifications’ construction and the establishment of permanent settlement. The modified simulation of the settlement phase was incorporated into the broader chronological sequence. The test models (Supplementary Fig. S4) demonstrate that simulated dates corresponding to the period before AD 625 (± 30–40) should be excluded, while calibrated values approximating AD 625 appear to represent a threshold. These models support the chronological observation derived from the simple model: that the settlement phase is likely to fall between the mid-7th century and the turn of the 8th and 9th centuries (ca. AD 650–800).
Bojná II. The radiocarbon dates that define the chronological framework of the Bojná II hillfort form two clusters (Supplementary Fig. S10). The first cluster comprises charcoal samples collected from the fortifications and ditch in the north-eastern part of the site. As these were not taken from structural elements, the dates—ranging from 1,250 ± 30 to 1,190 ± 30 BP—should be treated with caution; and, following the adopted methodology, they should be interpreted primarily as TPQs for construction of the fortifications. Notably, no older material was identified within this context, thereby supporting the interpretation that the ramparts were constructed during the early medieval period.
The second cluster of radiocarbon dates consists of samples obtained from burnt structural elements of the fortifications, recovered from the bottom of the ditch in the north-western part of the site (Supplementary Fig. S9). These dates indicate that this section of the defences is likely to have been constructed between the late 9th and late 10th centuries. Importantly, they build upon the TPQ established through the analysis of the first cluster, and thus reinforce the overall chronological framework for the Bojná II hillfort’s development.
In the chronological model (Fig. 4), the charcoal samples from the north-eastern rampart (which were not structural elements) were treated as TPQs, as they are likely to have entered the rampart fill from the surrounding area at the time of construction—analogous to the situation at Bojná III. Within this model, the TPQ for the fortifications’ construction falls at the end of the 9th century. This aligns perfectly with the earliest dates obtained from the burnt structural elements found at the bottom of the north-western ditch, and hence reinforces the interpretation of a construction phase in the late 9th to early 10th centuries.
This model is both safe to interpret and methodologically sound, as it rests on the reasonable assumption that the charcoal within the rampart originated from a phase of human activity predating its construction. Archaeological finds indicate the possibility of occupation at the site during the first half of the 9th century—a scenario supported by the model. However, these finds cannot confirm the existence of fortifications during that early phase.
Test. In the test, two models were developed (Supplementary Fig. S11; Supplementary Fig. S12) based on the assumption that the Bojná II site was occupied prior to the rampart’s construction. Since the radiocarbon data from the NE and NW sections of the rampart clearly cluster into two distinct chronological horizons, the test models applied a combination of data from these sections as indicators of two hypothetical phases—an earlier and a later one.
In the models, it was assumed that the radiocarbon data from the NE rampart reflect the chronology of a hypothetical event—either the construction of two phases of the rampart (earlier and later), or of two distinct sections (older and newer). The modelling results in both cases closely resemble those of the preliminary model, but they do not clearly confirm the existence of a discrete “event”. Instead, they require us to assume the presence of an earlier settlement phase, which the models constrain to the period between the late 8th and late 9th centuries. This interpretation is consistent with the archaeological evidence. However, it remains unclear whether the rampart was already constructed during this earlier phase—with the later phase reflecting a repair or modification—or whether it was built only at the end of the 9th century. From a methodological perspective, the latter hypothesis is more cautious and conservative interpretation.
Bojná I. The radiocarbon dates used to construct the chronological models for the Bojná I hillfort were grouped into three clusters. The first cluster comprises data obtained from the fortifications of the main enclosure, while the other two clusters correspond to archaeological features identified within the hillfort and the adjacent suburbium. For analytical purposes, the first cluster was further subdivided into two groups corresponding to samples from the eastern and western sections of the rampart. All samples in this first cluster were taken from structural elements of the fortifications and, in line with the adopted methodology, can be treated not merely as TPQs but as approximate indicators of the construction or use of the defences. Samples were collected from three locations along the rampart (the eastern section, the adjacent eastern gate, and the western section), and are interpreted as belonging to a single construction phase, with the exception of sample MKL-1420 (1,290 ± 60 BP). According to the field documentation, this sample may derive from an earlier (stratigraphically distinct) structure; and indeed, modelling consistently identifies it as an outlier across all scenarios. Although individual radiocarbon determinations from the remaining samples may suggest a relatively broad chronological spread—ranging from 1,290 ± 80 to 1,020 ± 35 BP—this variation falls within the 2σ confidence intervals for each sample, and is not inconsistent with the dendrochronological dating framework (see Supplementary Fig. S13)
Since dendrochronology has provided precise dating for structural elements associated with the final phase of the hillfort’s fortifications (after AD 866 – after AD 908), and given that the broader chronological framework for the Bojná I site (not necessarily limited to the hillfort) has been outlined through archaeological investigation, these findings were incorporated into several test models as a posteriori modifiers. Based on this combined evidence, the following assumptions can be made: (a) the rampart was constructed around AD 890; (b) the last repairs to the rampart were carried out after AD 908; and (c) the overall chronology of the Bojná I site, based on the archaeological material, spans approximately AD 800–950, and is likely to include more than one phase of use. The purpose of the test models was to evaluate whether, and to what extent, the radiocarbon data support these assumptions, and whether Bayesian modelling could reveal additional insights into the site’s chronological development.
The basic single-phase chronological model did not provide answers to detailed questions concerning the chronology of the site. Its main outcome was the clear indication that the radiocarbon data needed to be subdivided into distinct phases. Moreover, it was not possible to apply the TPQ-based approach used at the other sites, since no loose charcoal samples were collected from the Bojná I fortifications; only structural elements were dated. It was only through the construction and evaluation of multiple test models that a more refined and reliable chronological model for the development of the Bojná I hillfort could be established (Fig. S5).
Test. Repeated modelling of the dataset revealed a high degree of statistical consistency among the radiocarbon dates. In the first model, the samples were arranged in a depositional sequence reflecting the stratigraphic positions of structural elements: load-bearing components → inner wall elements → outer wall elements → outer basketwork. In the second model, all elements were assumed to belong to a single construction phase. The third model combined the assumptions of the previous two. These models (Supplementary Fig. S15) enabled the identification of outliers and elements that weakened the model’s robustness. A subsequent test model assumed that the rampart was in use between AD 800 and 1000 (Supplementary Fig. S16).
The models are weakened only by samples with chronologically extreme dates; however, these do not invalidate the models, and it is known with certainty that they originated from structural elements of the rampart. Moreover, the 1σ ranges of the modelled chronological phase in all test models align—either fully or partially—with the range of absolute dendrochronological dates, thus supporting the argument that the fortifications were constructed around AD 900. It is likely that older elements or long-lived timber (e.g., sample MKL-1422) were incorporated into the construction. It is also possible that the rampart replaced an earlier structure that was dismantled. During its use, the outer part of the structure appears to have undergone maintenance and repairs, as evidenced by the youngest radiocarbon dates.
The radiocarbon dates obtained from residential and economic features discovered within the hillfort and its suburbium exhibit greater chronological variation than those from the fortifications. These dates suggest the possibility of an earlier settlement phase predating the preserved fortifications (Supplementary Fig. S13). This, in turn, raises the possibility that the ramparts themselves may have had an earlier construction phase—a hypothesis already proposed on the basis of archaeological observations18, 19. Radiocarbon dating of samples from these features, combined with targeted chronological modelling, appears to support this interpretation (Supplementary Fig. S17; Supplementary Fig. S18; Supplementary Fig. S19; Supplementary Fig. S20).
At present, two distinct clusters of radiocarbon dates can be observed within the area of the hillfort: the first ranging from 1,290 ± 40 to 1,260 ± 30 BP preceded by two markedly earlier dates: 1400 ± 30 and 1370 ± 90 BP; and the second from 1,225 ± 30 to 1,090 ± 35 BP (or 1,197 ± 19 BP combined, to 1,090 ± 35 BP). Interestingly, neither this older cluster nor the very early individual dates have yet been identified within the suburbium (Supplementary Fig. S13 B; Supplementary Fig. S14), where only a small number of features have been investigated. The samples from those contexts show striking consistency, clustering tightly within the range of 1,185 ± 30 to 1,165 ± 30 BP (or 1,178 ± 18 BP combined, to 1,165 ± 30 BP); this suggests a later and relatively short-lived phase of occupation in the suburbium—an interpretation that is supported by archaeological observations.
The only earlier date from the suburbium—Poz-16305/20/6: 1,245 ± 30 BP—was obtained from a large posthole pit, possibly associated with fortification construction. Notably, this date is surprisingly close to the outlier from the rampart, MKL-1420: 1,290 ± 60 BP, which perhaps indicates the existence of an earlier phase of fortification. However, these results cannot yet be considered definitive for determining the absolute chronology of Bojná I, or for clearly distinguishing individual chronological phases. The data are derived from charcoal samples retrieved from feature fills; they must therefore be treated primarily as indicators of periods of human activity, and as TPQs. These dates cannot be confidently correlated with specific phases of fortification construction. That said, their overall range is strikingly similar to that of the radiocarbon chronology for the ramparts, which suggests a degree of temporal alignment between settlement activity and fortification use.
In the final test models (Supplementary Fig. S21; Supplementary Fig. S22; Supplementary Fig. S23), the phasing of the Bojná I site was simulated based on archaeological data and the results of previous radiocarbon analyses of the fortifications. The modelling and interpretation of samples from archaeological features indicated several equally plausible ways to group the data. However, all options require acknowledging the existence of at least one settlement phase predating the construction of the preserved rampart. This rampart was built at the end of the 9th century and underwent repairs during the first decade of the 10th century. This chronological transition is clearly visible in every radiocarbon model.
The earlier settlement phase most likely dates generally to the 9th century, although it is currently not possible to definitively determine its chronological boundaries. Based on the adopted methodology, it is possible to hypothesise that the duration of this earlier phase may be defined by the TPQs of the respective radiocarbon clusters. According to models for features 1–2 (Supplementary Fig. S17; Supplementary Fig. S18) and the overall test model 2 (Supplementary Fig. S22), the settlement phase preceding the rampart’s construction can be framed within ca. AD 840–900. In contrast, models 3 and 4 for individual features (Fig. S19; Fig. S20) and overall test model 1 (Supplementary Fig. S21) extend this phase much further into the 8th century. While not an implausible hypothesis, this extension does not correlate with the archaeological evidence gathered during excavations at the site. This discrepancy remains unexplained, and it is suspected that the radiocarbon calibration curve may be partly responsible. The modelling process helps to highlight this issue. The two oldest dates from Bojná I (1,400 ± 30 and 1,370 ± 90 BP) are so far removed from the mean values of other samples that they can only serve as very general TPQs, roughly corresponding to the initial stages of regional settlement (see Supplementary Fig. S5). Among all the tested scenarios, test model 3 (Supplementary Fig. S23), which assumes that intensive settlement activity ended before AD 1000, appears to be the most probable.
The test chronological models for the complete dataset (Supplementary Fig. S24) indicate that the period of intensive use of the site is likely to fall between ca. AD 840 and 940. This finding allows the terminus ante quem (TAQ) for the radiocarbon dates to be further narrowed in the final chronological model (Fig. 6), thus potentially setting it at AD 950.
All presented models exhibit a high degree of statistical plausibility; however, the most robust and informative are those that assume the existence of more than one occupational phase at the Bojná I site. These results clearly indicate that future archaeological investigations must seriously consider the possibility of a multi-phased settlement development. The final chronological model (Fig. 6) may serve as a testable hypothesis, capable of guiding subsequent archaeological investigations. Particular attention should be given to the earliest phase of occupation, which extends significantly earlier than the chronology inferred from the artefactual assemblage. Additionally, the refinement of the terminal date for settlement activity would further enhance the precision of the modelling process.
Bojná VI. The dataset obtained from the Bojná VI site demonstrates a high degree of incoherence—most likely due to the limited number of samples—which hinders the identification and potential rejection of outliers. This issue is compounded by the selection of materials that are inherently problematic for radiocarbon dating (e.g., slag and human bone). As a result, the two dates obtained from slag samples taken from a furnace (1,165 ± 30 and 1,015 ± 30 BP) fail the χ² agreement test when combined (T = 12.498, critical value for 5% = 3.8), despite appearing individually plausible. The situation is further complicated by an entirely unreliable radiocarbon date of 10,150 ± 60 BP derived from charcoal found in the same furnace, which strongly reinforces doubts regarding the slag dates46. Another questionable result is the date of 1,120 ± 30 BP obtained from a human bone burial beneath Mound 3. While this would imply an interment in the 10th century AD—which is not inherently inconsistent with the overall chronology of the Bojná agglomeration—it does not align with the archaeological evidence recovered from Bojná VI. Chemical analysis of sample Poz-20138/22/4 revealed 0% nitrogen and only 0.8% carbon, which confirms its unsuitability for radiocarbon dating39, and the result must therefore be deemed unreliable. Consequently, the dating of Bojná VI must rely exclusively on two samples—POZ-4 and Poz-172285—whose calibrated dates are consistent with the archaeological record (see Supplementary Fig. S28). These two dates were used in the KDE modelling.
General model. For all sites, a chronological modelling test was also conducted using the KDE method (Supplementary Fig. S26; S27). However, in this case, models constructed for individual sites reflect relatively short chronological spans, and do not yield any conclusive or transformative results beyond indicating peaks of activity (example on Supplementary Fig. S25). Potential occupational phases are not revealed within the models45.
Bayesian modelling was demonstrably more effective, as confirmed by comparative analysis. KDE modelling of the entire Bojná agglomeration dataset did, however, confirm the general chronological succession of settlement activity between sites (Fig. 2), though it failed to disclose the finer detail evident in the Bayesian models. This limitation arises from the relatively narrow time window under investigation, and the uneven sample representation, which is heavily skewed by data from a single site.