Stemt de Sumerische Koningenlijst overeen met paleoklimatologische gebeurtenissen?

Context

De Sumerische Koningenlijst begint met acht koningen die regeerden vóór de zondvloed. Hun regeringsperioden zijn enorm en ongebruikelijk regelmatig. Drie voorbeelden zijn 28.800 jaar, 36.000 jaar en 43.200 jaar. De meeste zijn gehele veelvouden van 3.600, en alle acht zijn veelvouden van 600. Het totaal is 241.200 jaar.

De inputsequentie die hier is getest, is de ETCSL-samengestelde antediluviaanse lijst:

Volgorde	Koning	Regeringsduur
1	Alulim	28.800 jaar
2	Alalgar	36.000 jaar
3	Enmenluana	43.200 jaar
4	Enmengalana	28.800 jaar
5	Dumuzid	36.000 jaar
6	Ensipadzidana	28.800 jaar
7	Enmenduranna	21.000 jaar
8	Ubara-Tutu	18.600 jaar
Totaal		241.200 jaar

De transliteraties volgen de hieronder geciteerde samengestelde tekst. Net als de secties na de zondvloed is de antediluviaanse lijst een tekstuele traditie met varianten; de tabel is dus een input voor analyse en geen claim over letterlijke historische regeringsperioden.

Een speculatieve interpretatie beschouwt deze getallen als een vervormde herinnering aan de prehistorie. Volgens deze hypothese coderen de grenzen van de regeringsperioden werkelijke klimaatverschuivingen, uitbarstingen, inslagen of veranderingen in de zeespiegel. Na het schalen en het koppelen van de lijst aan een voorgestelde datum van de zondvloed, zouden de grenzen moeten samenvallen met gedateerde gebeurtenissen in het geologische archief.

In deze verkennende analyse test ik dat idee. De explorer schaalt elke chronologie naar een vaste periode van 241,2 ka (duizend jaar), koppelt één grens aan een vast punt en vergelijkt de resulterende data met een catalogus van Quternaire gebeurtenissen. Ik gebruik 11,6 ka BP (ongeveer 11.600 jaar geleden) als een door de analist gekozen ankerpunt nabij het einde van het Jonger Dryas. De Koningenlijst zelf geeft deze datum niet op en suggereert deze paleoklimatologische interpretatie niet.

Het vinden van overeenkomsten is eenvoudig. Met negen grenzen en de vrijheid om het ankerpunt of de bandbreedte te wijzigen, komen toevallige overeenkomsten vaak voor. De relevante vraag is of de waargenomen volgorde van de Sumerische regeringsperioden ongebruikelijk hoog scoort onder een duidelijk gedefinieerd nulmodel.

Resultaten

In de primaire paleoklimatologische catalogus vertoont de Sumerische sequentie geen statistisch significante uitlijning. Bij het vaste ankerpunt van 11,6 ka en een kernel-bandbreedte van $\sigma = 1,60$ ka is de permutatie p-waarde 0,350. Na correctie voor meervoudige vergelijkingen is $q = 0,622$.

Het uitbreiden van de analyse naar alle 103 bruikbare catalogusvermeldingen verhoogt het aantal schijnbare overeenkomsten, maar verandert de conclusie niet. Bij dezelfde bandbreedte is de p-waarde voor de brede catalogus 0,148 en de gecorrigeerde q-waarde 0,430.

De kleinste ruwe p-waarde voor de Sumerische sequentie komt voor in de catastrofale verkennende catalogus bij $\sigma = 1,60$ ka: $p = 0,021$. Dit was het beste resultaat dat werd gevonden in een grotere verkennende zoektocht. Na correctie voor alle vergelijkingen is $q = 0,222$, waardoor het resultaat de hypothese niet ondersteunt.

Als secundaire controle tel ik ook hoeveel gebeurtenissen binnen een vaste afstand van een grens vallen. De kernel-score is de belangrijkste maatstaf, omdat deze minder gewicht toekent aan gebeurtenissen die verder weg liggen, in plaats van een abrupte afkapgrenzen te gebruiken.

Wat zou gelden als bewijs?

Een geloofwaardige uitlijning zou aan drie voorwaarden moeten voldoen.

Ten eerste moet de tekst bepalen welk ankerpunt wordt gebruikt. De Sumerische Koningenlijst plaatst de zondvloed na de achtste regeringsperiode, dus ik anker het einde van de antediluviaanse sequentie. De tekst kent geen absolute datum toe aan die grens. Ik gebruik 11,6 ka BP als een door de analist gekozen datum nabij het einde van het Jonger Dryas.

Ten tweede zou de Sumerische sequentie er ongebruikelijk uit moeten zien naast andere oude chronologieën. Ik vergelijk deze met een sequentie afgeleid van Bijbelse patriarchale leeftijden, zeven regeringsperioden van goden en halfgoden uit het Manethonische fragment bewaard in de Excerpta Latina Barbari, en de eerste acht vermelde heersers van Kish I. Dit zijn voorbeelden ter vergelijking, geen onafhankelijke statistische controles; ik test niet of de ene chronologie beter presteert dan de andere. Ze laten zien hoe gemakkelijk niet-gerelateerde oude sequenties schijnbare overeenkomsten kunnen produceren onder dezelfde procedure.

Ten derde zou een geloofwaardig resultaat significant moeten blijven na correctie voor de geteste chronologieën, catalogusniveaus en bandbreedtes. Geen van de resultaten voldoet hieraan.

Interactieve Explorer

De interactieve explorer maakt het mogelijk om te schakelen tussen chronologieën, catalogusniveaus en bandbreedtes. De primaire paleoklimatologische catalogus definieert de hoofdanalyse. De brede en catastrofale catalogussen dienen voor verkennende gevoeligheidsanalyses. De sectie met het vaste ankerpunt gebruikt p-waarden die vooraf zijn berekend in Rust. De sectie met het schuivende ankerpunt berekent de sweep van het ankerpunt live in de browser.

Hoe de cijfers te lezen

De belangrijkste statistiek is een Gaussiaanse kernel-nabijheidsscore:

$S = \sumi \sumj \exp(-(ti - bj)^2 / 2\sigma^2)$

Hier is $ti$ de datum van een gebeurtenis en $bj$ een grens van een regeringsperiode. Een nabije gebeurtenis draagt bijna één punt bij, terwijl de bijdrage van een verre gebeurtenis naar nul nadert. Dit voorkomt de abrupte discontinuïteit van een harde afkapgrens, waarbij een gebeurtenis net binnen het venster telt en een gebeurtenis net buiten het venster niet.

Er zijn twee verschillende statistische vragen. Het eenvoudigste onderscheid is wat er mag bewegen:

Methode	Vastgehouden	Mag bewegen	Wat wordt gerapporteerd	Rol
Uitputtende permutatie	11,6 ka anker, catalogusniveau, bandbreedte en dezelfde set regeringslengtes	De volgorde van de regeringslengtes	Fractie van alle gelabelde volgordes die minstens zo hoog scoren als de waargenomen volgorde	Primaire p-waarde
Random-anker Monte Carlo	Volgorde regeringslengtes, catalogusniveau en bandbreedte	Het ankerpunt, of het lokale anker-venster	Fractie van bemonsterde ankers of vensters die minstens zo hoog scoren	Secundaire gevoeligheidscontrole

De twee bandbreedtes zijn gekalibreerd om hetzelfde totale gewicht te hebben als harde vensters die 1 ka en 2 ka aan weerszijden van een grens reiken. Het gebruik van $\sigma = \tau \cdot \sqrt{2/\pi}$ geeft $\sigma \approx 0,80$ ka en $\sigma \approx 1,60$ ka. De browser kapt de kernel af bij $4\sigma$ voor de snelheid, wat overeenkomt met de Rust-preberekening.

De primaire analyse stelt het anker vast op 11,6 ka BP en rapporteert de waargenomen kernel-score, een ruwe permutatie p-waarde en een q-waarde gecorrigeerd voor meervoudige vergelijkingen. De p-waarde is afkomstig van een uitputtende permutatietest, niet van Monte Carlo-bemonstering. Het programma controleert elke mogelijke volgorde van regeringsperioden. Met het anker, catalogusniveau, bandbreedte en de set regeringslengtes vastgehouden, is de p-waarde de fractie van die volgordes waarvan de kernel-score minstens zo groot is als de waargenomen score. Onder dit nulmodel wordt elke volgorde als even plausibel beschouwd. De test houdt geen rekening met het kiezen van het anker, de catalogus, de bandbreedte of de score na inspectie van de data.

De primaire catalogus bevat het einde van het Jonger Dryas, dat werd gebruikt om het ankerpunt van 11,6 ka te motiveren, waardoor de terminale match in de opzet is ingebouwd. De bijdrage hiervan is constant over de permutaties van de volgorde en is geen bewijs dat de interne sequentie ongebruikelijk is.

De secundaire ankerzoektocht onderzoekt hoe het resultaat verandert wanneer het anker wordt geoptimaliseerd na inspectie van de data. Het vergelijkt het maximum in het vaste venster van 10 tot 13 ka met maxima uit 3.000 vensters waarvan de centra uniform zijn bemonsterd van 0 tot 100 ka. Zowel de vaste als de willekeurige vensters gebruiken een breedte van 3 ka en een anker-afstand van 0,1 ka. De weergegeven Monte Carlo-schatting voegt één toe aan zowel het aantal overschrijdingen als het aantal pogingen, zodat een eindige simulatie nooit een waarschijnlijkheid van precies nul rapporteert. Dit is niet de primaire p-waarde van het artikel.

De resultatenkaart voor het vaste anker rapporteert ook een aparte gevoeligheidswaarde voor een willekeurig anker. Voor elke chronologie, catalogusniveau en bandbreedte vergelijkt het Rust-programma de score bij 11,6 ka met scores van zes miljoen ankers die uniform zijn bemonsterd van 0 tot 100 ka. Deze berekening is verschillend van de Monte Carlo-schatting van 3.000 vensters in de browser. Geen van beide waarden is de primaire p-waarde.

Data en vergelijkingen

De broncatalogus bevat 104 gedateerde gebeurtenissen uit het Quaternair, verdeeld over 11 categorieën. De primaire analyse gebruikt een beperktere paleoklimatologische catalogus met data van Heinrich-gebeurtenissen, het begin van Greenland Interstadials en geselecteerde Holocene klimaatgebeurtenissen. De resulterende primaire catalogus bevat 39 gebeurtenissen binnen het analysevenster van 0 tot 260 ka BP.

De brede verkennende catalogus bevat alle 103 bruikbare vermeldingen, inclusief de 39 primaire vermeldingen. Deze bevat daarnaast smeltwaterpulsen, grenzen van Marine Isotope Stages, grote vulkaanuitbarstingen, inslagstructuren en betwiste inslaghypotheses, geomagnetische excursies, extreme zonne-protongebeurtenissen, extinctieknopen van megafauna en belangrijke culturele overgangen. De catastrofale en "deep-time" niveaus bieden overlappende gevoeligheidsanalyses. Lonar is opgenomen in de broncatalogus, maar uitgesloten van het dashboard omdat het buiten het analysevenster valt.

Elke vermelding registreert een geselecteerde datum, een geschatte onzekerheid, broninformatie en of de gebeurtenis betwist is. De onzekerheidswaarden werden op verschillende manieren geschat en zijn daarom niet direct vergelijkbaar; ze moeten niet allemaal als standaardfouten worden gelezen. De huidige analyse gebruikt alleen de geselecteerde data. De brede catalogus is nog steeds verkennend, ook al heeft elke vermelding een gedocumenteerde bron.

Voor de Bijbelse vergelijking zijn de eerste acht waarden de levensduur van de patriarchten en de laatste 600 de leeftijd van Noach bij de zondvloed. Voor de vergelijker Excerpta Latina Barbari gebruik ik de zeven genoemde regeringsperioden van goden en halfgoden in dat fragment: Hephaestus 680, Helios 77, Sosinosiris 320, Orus 28, Typhon 45, Anubes 83 en Amusis 67 jaar. Voor Kish I gebruik ik de eerste acht heersers in de ETCSL-samengestelde tekst; deze traditie heeft manuscriptvarianten, dus de vergelijker moet beschrijvend worden gelezen en niet als een vaste historische chronologie.

Elke chronologie wordt op dezelfde manier getransformeerd. Regeringslengtes worden geschaald zodat de totale span gelijk is aan het Sumerische antediluviaanse totaal van 241,2 ka, en vervolgens omgezet in absolute data door één grens te ankeren en de duurstappen terug in de tijd op te tellen. Deze normalisatie verwijdert de totale duur als factor en vergelijkt alleen de relatieve spreiding van de grenzen.

Het Rust-programma leest dezelfde events.json-catalogus die bij dit artikel is gepubliceerd. Het voert zes miljoen random-anker-tests uit per chronologie, catalogusniveau en bandbreedte, voert uitputtende permutatietests uit, past de Benjamini-Hochberg-correctie toe op de permutatie p-waarden en schrijft de resultatentabel die door de pagina wordt gebruikt. De TypeScript-explorer laadt de publieke catalogus en resultatentabel en berekent en toont de live anker-sweep.

Implementatienota

Het rekenintensieve werk wordt offline in Rust uitgevoerd in plaats van in de browser. De primaire test is klein genoeg om exact te enumereren: de Sumerische en Kish-sequenties hebben elk $8! = 40.320$ gelabelde volgordes, de Excerpta Latina Barbari-sequentie $7! = 5.040$, en de Bijbelse vergelijker $9! = 362.880$. Herhaalde regeringslengtes worden nog steeds behandeld als gelabelde posities, zodat het nulmodel vraagt of deze specifieke volgorde ongebruikelijk is tussen alle herordeningen van dezelfde waarden.

De kernberekeningen zijn expliciet in de broncode: Gaussiaanse kernel-scoring, Heap's algoritme voor permutaties, de random-anker gevoeligheidscontrole, de Benjamini-Hochberg step-up aanpassing en de deterministische JSON-schrijver. De random-anker controles gebruiken een vaste seed, zodat het opnieuw uitvoeren van cargo run --release het gepubliceerde resultatenbestand regenereert in plaats van elke keer een nieuwe simulatie te produceren.

Kanttekeningen

Een Benjamini-Hochberg-aanpassing voor meervoudige vergelijkingen dekt de 32 vaste-anker permutatietests. Het corrigeert niet voor eerdere experimenten met de inclusie van gebeurtenissen, ankerkeuzes, scoredefinities of de constructie van de vergelijkers. Omdat dit een verkennende analyse is en de volledige set tests niet vooraf was gespecificeerd, vatten de q-waarden deze zoektocht samen in plaats van een bevinding te bevestigen. Geen enkele vergelijking heeft een $q < 0,05$.

De catalogus is redactioneel samengesteld in plaats van via een systematische review of een vooraf gedefinieerd inclusieprotocol. Meerdere vermeldingen beschrijven gerelateerde delen van dezelfde klimaatsequentie, waardoor de 39 primaire data niet moeten worden geïnterpreteerd als 39 onafhankelijke waarnemingen. De permutatietest houdt deze catalogus vast; het kan de catalogus niet onafhankelijk of volledig maken.

De dichtheid van de catalogus varieert aanzienlijk over de tijd. Recente gebeurtenissen zijn talrijker en over het algemeen beter gedateerd, waardoor chronologieën met Holocene grenzen meer kansen hebben om overeen te komen met een gebeurtenis. De permutatietest pakt deze onbalans gedeeltelijk aan door het anker en de gebeurteniscatalogus vast te houden terwijl alleen de volgorde van de regeringsperioden verandert.

Het "deep-time" gevoeligheidsniveau verwijdert het drukke recente record door alleen data ouder dan 60 ka te gebruiken. De Sumerische sequentie scoort bij beide bandbreedtes op of onder het niveau van toeval.

De analyse neemt de dateringsonzekerheid niet mee in de berekeningen. Meerdere vermeldingen hebben onzekerheden die breder zijn dan de vaste match-vensters, en de impact-hypothese van het Jonger Dryas blijft betwist. Elke vermelding behandelen als een exacte datum overschat de precisie van het resultaat.

De gerapporteerde p-waarden en q-waarden houden geen volledig rekening met de keuzes van de analist, inclusief beslissingen die zijn genomen na inspectie van de overeenkomsten. De interface legt sommige van deze keuzes bloot, maar de analyse blijft verkennend.

Data en Code

De gebeurteniscatalogus en de gegenereerde resultatentabel zijn publiek en machineleesbaar:

• Event catalog with source metadata
• Precomputed dashboard result JSON
• Fixed-anchor Rust precompute, browser script, catalog, and results release v0.2.1

De gebeurteniscatalogus registreert het analyseniveau, de datumschatting, onzekerheid, bron-URL, betwistingsvlag en redactionele opmerkingen voor elke vermelding. Het resultatenbestand bevat de waarden die door het Rust-programma zijn gegenereerd en door de browser-explorer worden weergegeven.

Geselecteerde Bronnen

• The Electronic Text Corpus of Sumerian Literature, Sumerian King List
• CDLI record for the Sumerian King List tradition
• Hemming 2004, Heinrich events
• Rasmussen et al. 2014, INTIMATE event stratigraphy for Greenland ice cores
• Lisiecki and Raymo 2005, LR04 benthic stack
• Earth Impact Database
• The Fragments of Manetho, Book I, including Excerpta Latina Barbari

Volgende Stappen

Een toekomstige versie zou de dateringsonzekerheden op een consistentere basis moeten plaatsen en deze door de berekeningen moeten voeren. Het herhaaldelijk bemonsteren van plausibele data voor elke gebeurtenis zou een reeks scores en "hit counts" opleveren in plaats van een enkele waarde per item.

Een sterkere analyse zou gebruikmaken van een catalogus die onafhankelijk is samengesteld en gedefinieerd voordat enige chronologie wordt gescoord. Het herhalen van de procedure met verschillende van dergelijke catalogussen zou testen of het resultaat afhangt van de onderliggende gebeurtenisdata of van de redactionele curatie.

Deze analyse onderzoekt een specifieke claim over uitlijning. Het vindt geen statistisch significant bewijs dat de Sumerische volgorde van regeringsperioden de hier gebruikte paleoklimatologische catalogus codeert.