Clarke, A. C., Leary, J., & Kirkland, C. L. (2026). Deliberate prehistoric sourcing of the Devil's Arrows. Proceedings of the Royal Society A. [rspa.2026.0504].
https://royalsocietypublishing.org/rspa/article/482/2345/20260504/483254/Deliberate-prehistoric-sourcing-of-the-Devil-s
I have now had a chance to go through this paper in more detail, and to discuss it with other experts. So here we go. Let's discuss the technical deficiencies in sampling symmetry, statistical sample density, and logical consistency in provenance attribution. It's a bit complicated.......
There are serious doubts about the empirical validity and logical framework of the recently published paper by Clarke et al. (2026) regarding the provenance of the four Devil's Arrows monoliths [rspa.2026.0504]. While the application of high-resolution U-Pb isotopic testing on accessory minerals is a welcome addition to archaeological science, a critical examination of the raw data structure and methodological parameters reveals significant vulnerabilities that undermine the paper’s primary conclusions and indeed its spectacular title. The article has been widely promoted via the media, and the record must be set straight.
The peer review process used by editor of this prestigious Royal Society journal has apparently failed to pick up on six severe methodological and logical flaws:
1. Severe Asymmetry in Sampling and Provenance Controls
The paper compares an undocumented, crushed bulk rock sample from a single, prominent modern tourist location (Brimham Rocks SSSI) (Fig. 1) against non-destructive, surface-adhesive "sticky tape" scrapings taken from the heavily weathered faces of three monoliths and the fresher face of another [rspa.2026.0504] (Fig. 2). The authors fail to present any spatial control sampling across the lateral extent of the target stratigraphic layer. We have no idea where the sample came from within the 30m thick sedimentary sequence, or within the SSI which has a surface area of 1.84 sq km. We are not told whether the rock sample represents the upper or lower boundaries of the sequence, or how the mineral concentrations shift vertically through the 30-metre thickness or laterally across the site. This is an absolute failure of scientific protocol. It means the experiment is completely unreproducible.The Kinderscoutian Brimham Grit is a continuous, sprawling sedimentary formation which occurs across c 140 square kilometres of the North Yorkshire Pennines. Because detrital mineral suites reflect continental river systems, the mineral "barcode" match to this formation must apply to the entire regional sheet, not to a singular geographic coordinate (Morton et al., 2005).
2. Violation of Minimal Grain-Count Statistical Thresholds
To achieve an acceptable 95% to 99% confidence level that minor mineral age populations have not been missed, established sedimentological protocols dictate an absolute minimum of 60 to 117 individual grain analyses per discrete sample (Vermeesch, 2004). While the authors aggregate a total pool of 236 dated zircon grains from the monuments across Figure 3, their comparative breakdowns in Figure 4 expose a fatal lack of data density. The authors split the monument data into two age-segregated sub-panels: a 90-grain subset mapping the deep-time spectrum (500–3500 Ma) and a 146-grain subset mapping the younger spectrum (300–600 Ma) [rspa.2026.0504]. In doing so, they reveal that their primary geographical point-control—the comparison baseline sample for Brimham Rocks—relies on an astonishingly small pool of just 29 grains [rspa.2026.0504]. Attempting to verify an exclusive geographic match to a massive regional sedimentary unit using a sample that fails to reach even half of the industry-accepted minimum threshold is unacceptable.
3. Artificial Variance Dilution via "Phantom" Data Aggregation
Rather than executing a robust point-to-point comparison (individual monolith vs. specific outcrop candidates), the authors combine their low-yield, heterogeneous grain counts into a single, pooled composite "blob” to smooth away internal variance. Furthermore, the Multi-Dimensional Scaling (MDS) diagrams are heavily packed with distal, irrelevant regional "noise" datasets (e.g., Precambrian basement terranes and distinct Triassic strata) [rspa.2026.0504]. This mathematical dilution dramatically compresses the relative variance of the graph, artificially forcing the visually distinct raw histograms of the Arrows and Brimham Rocks into an unrepresentative, close cluster.
4. The Dismissal of the Glacial Transport Option
A fatal flaw of the paper’s internal logic lies in its selective deployment of Late Devensian (LGM) ice-flow models. In the results section, the authors explicitly acknowledge that a literature search yields a lose match between the Devils Arrows samples and the Langsettian sandstones of the Pennine Basin Lower Coal Measures (PLCM) at Binchester Crags [rspa.2026.0504]. They choose to ignore the fact that reconstructed complex ice flow directions (Thorpe and Williams-Thorpe, 1991; Clark et al., 2022) point to a strong possibility of glacial erratic transport from the southern part of County Durham to the Boroughbridge area during several glaciations (Livingstone et al, 2012) (Fig. 3). The Park Farm borehole record, from a site near Boroughbridge, reveals the Vale of York Formation, a glacial diamicton laid down during the Dimlington Stadial and carrying abundant Carboniferous sandstone erratics from the north. But Clarke and his colleagues then use the ice-flow direction map to demonstrate the “impossibility" of glacial erratic transport between Brimham Rocks and the Devils Arrows site, and invoke speculative human transport instead. This demonstrates a powerful and unacceptable bias in their analysis. What is good for the gander must be good for the goose: a model cannot dismiss a statistically superior northern match based on ice direction constraints while simultaneously throwing out those exact same glaciological constraints to protect a highly speculative “human transport" alternative.
5. The Lack of Evidence for the Human Transport Option
Having discarded the ice-flow model, they invent a highly speculative 11-mile human transport route to explain the presence of the Devil’s Arrows monoliths. They present no evidence in support of the human transport hypothesis but instead cite unproven claims of long distance stone transport from Ireland and Scotland — partly based on citations of their own hotly disputed work on the Stonehenge Altar Stone.
6. Inadequate Zircon Age Analyses
Independent analysis of the zircon age data by a specialist geologist suggests the foillowing: (1) the zircon age pattern for stone 1 looks different from stones 2 and 3 — there is a larger group of Archaean zircons and a difference in the structure in the mid-Proterozoic age range — all suggesting at least two different sources for the monoliths; (2) given that problem, it is not appropriate to merge all the data into a single “representative” sample; (3) there is a good match between the Brimham Rocks sample and the samples from stones 2 and 3, but the Early Palaeozoic group appears bimodal, unlike any of the monoliths; (4) the Plumpton Rock ages do look different; (5) the dismissal of the underlying Sherwood Sandstone Group is premature, based on a marked contrast with the Triassic, but the papers quoted by Clarke et al use samples in the Wessex Basin, over 150 miles away; (6) the most important point to emerge from this new analysis is that a number of other published data sets do not look very different from Brimham Rocks. For example, a Marsdenian sample looks very similar to Stone 1, and Yeadonian, Rough Rock N and S samples are not dramatically different to Stones 2 and 3. It must also be borne in mind that there are extensive outcrops of Millstone Grit that have never been analysed by zircon specialists. In summary, the conclusions of Clarke et al are highly questionable in terms of the zircon data.
Conclusion
Ultimately, the central thesis of this paper rests on a suspect, shifting statistical strategy. The authors choose to aggregate data when they need to mask individual sample failures, throwing low-yield, weathered monument surface grains into a collective "phantom blob" to fake a robust, unified monument profile [rspa.2026.0504]. Conversely, they choose to split data when they need to disguise a visual mismatch, slicing the continuous age spectrum into separate sub-panels to hide an unviable control sample size [rspa.2026.0504]. This dynamic manipulation allows the authors to bend a chaotic geological footprint to fit a pre-packaged public relations narrative. By substituting tactical statistical smoothing for rigorous field controls and geomorphological fieldwork, they have published a localized point-source attribution that completely collapses under basic sedimentological review.
Thus Clarke et al. have confirmed that the Devil's Arrows have probably come from a Carboniferous deltaic sandstone sheet [rspa.2026.0504], as previously suspected. However, the data presented is statistically incapable of pinpointing a precise geographic source or ruling out local glacial transport from the extensive northwestern outcroppings of that same layer or from the PLCM formation that provides a better match.
It is vital that published conclusions remain strictly bound to their empirical evidence. The authors of this deeply flawed paper must provide a formal clarification, displaying their raw, unaggregated individual monolith histograms side-by-side against a true, non-diluted matrix of Pennine Basin gritstone candidates. Only when they sample for like-for-like similarities and differences across a large range of localities can they expect their work to be taken seriously.
If the authors cannot provide the data requested, I consider that the paper should be retracted.
References
Clark, C. D., et al, BRITICE-CHRONO Consortium. (2022). Growth and retreat of the last British–Irish Ice Sheet, 31 000 to 15 000 years ago: The BRITICE‐CHRONO reconstruction. Boreas, 51(4), 699–758.
Clarke, A., Leary, J., & Kirkland, C. (2026). Deliberate prehistoric sourcing of the Devil's Arrows, Britain's tallest stone row. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 482(2345), 20260504. [rspa.2026.0504].
Livingstone, S.J., Evans, D.J.A., Cofaigh, C.O. et al. (2012) Glaciodynamics of the central sector of the last British-Irish Ice Sheet in Northern England. Earth Science Reviews, 111 (1-2). 25 - 55. ISSN: 0012-8252
Morton, A., Hallsworth, C., & Moscariello, A. (2005). Stratigraphic evolution of Carboniferous sandstones in the Pennine Basin: Heavy mineral provenance constraints. Proceedings of the Yorkshire Geological Society, 55(3), 163–185.
Thorpe, R.S., Williams-Thorpe, O..(1991) The myth of long-distance megalith transport. Antiquity 65, 64–73.(doi:10.1017/S0003598X00079308)
Vermeesch, P. (2004). How many grains are needed for a detrital zircon age spectrum? Earth and Planetary Science Letters, 224(3-4), 441–451.
Williams-Thorpe, O., Potts, P. J., & Thorpe, R. S. (1999). The provenance of archaeological gritstone monoliths in northern England: A geochemical approach. Journal of Archaeological Science, 26(2), 145–161.