Méthodologie MEGALITHICA

Le moteur MEGALITHICA est une approche computationnelle inédite de l'archéoastronomie qui identifie les alignements astronomiques entre sites mégalithiques en exploitant la lente variation de l'obliquité axiale de la Terre au fil des millénaires. Le pipeline d'analyse comporte six étapes :

1. Calculer le relèvement géodésique entre deux sites
2. Pré-filtrage : la fourchette de relèvement élimine ~50-60% des paires
3. Convertir le relèvement en déclinaison cible
4. Construire une grille temporelle : calculer la déclinaison de l'événement via l'obliquité pour chaque époque
5. Détection de changement de signe pour localiser l'année d'alignement
6. Affiner et calculer le résidu en secondes d'arc

L'axe de rotation de la Terre précesse avec une période d'environ 25 772 ans. Cela fait varier l'obliquité (inclinaison) d'environ 22,1° à 24,5° sur un cycle de ~41 000 ans. À J2000.0 (an 2000), l'obliquité moyenne est de 23°26'21,448" (≈ 23,4393°).

Comme la déclinaison au solstice est égale à l'obliquité, l'azimut où le Soleil se lève et se couche au solstice change graduellement au fil des millénaires — ce qui nous permet de « dater » un alignement en trouvant quelle époque fait correspondre le relèvement.

Le polynôme de Laskar (1986) calcule l'obliquité moyenne avec une précision de ~0,01° sur ±10 000 ans :

ε₀ = 23° 26' 21.448"

U = T / 100    (T = Julian centuries from J2000)

Δε = -4680.93·U - 1.55·U² + 1999.25·U³
     - 51.38·U⁴ - 249.67·U⁵ - 39.05·U⁶
     + 7.12·U⁷ + 27.87·U⁸ + 5.79·U⁹ + 2.45·U¹⁰

ε = ε₀ + Δε / 3600

L'azimut à l'horizon où un corps céleste de déclinaison δ se lève ou se couche à la latitude φ est donné par l'astronomie sphérique :

cos(Az) = sin(δ) / cos(φ)

Rising:  Az (eastern horizon)
Setting: 360° - Az (western horizon)

Cette formule suppose un horizon plat. Les corrections de réfraction atmosphérique (34') et de semi-diamètre solaire (16') peuvent être appliquées pour une plus grande précision.

Le relèvement géodésique du site A au site B utilise la formule d'azimut direct :

θ = atan2(
  sin(Δλ) · cos(φ₂),
  cos(φ₁) · sin(φ₂) - sin(φ₁) · cos(φ₂) · cos(Δλ)
)

La distance utilise la formule de Haversine, précise à ~0,3% pour les distances inférieures à 100 km.

L'algorithme principal parcourt une grille temporelle de 15 000 av. J.-C. à 500 apr. J.-C., calculant l'azimut à l'horizon de l'événement à chaque époque et le comparant au relèvement :

1. Pour chaque année de la grille (pas de 100 ans) :
2. Calculer obliquité → déclinaison de l'événement → azimut à l'horizon
3. Calculer diff = relèvement - azimutÉvénement
4. Quand le signe change entre deux pas → affinage par recherche binaire
5. Calculer le résidu en secondes d'arc
6. Si résidu ≤ seuil → alignement trouvé

Les deux directions (A→B et B→A) sont testées. Les alignements d'équinoxe et cardinaux sont indépendants de l'époque (déclinaison = 0° ou azimut fixe).

L'« Année Parfaite » est le concept central de la méthodologie MEGALITHICA. Pour tout alignement entre deux sites, elle répond à la question : en quelle année historique un événement céleste (lever ou coucher de soleil, lever ou coucher de lune) s'est-il produit exactement dans l'azimut entre les deux sites ?

     Site A                                        Site B
       ●━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━●
       │          bearing = 51.2°
       │
       │    In which year did the Sun rise
       │    at exactly 51.2° azimuth?
       │
       ▼
  ┌─────────────────────────────────────────────┐
  │  Year scan: 15,000 BCE ──────────▶ 500 CE   │
  │                                             │
  │  For each year:                             │
  │    obliquity(year) ─▶ declination ─▶ azimuth│
  │                                             │
  │  Match found: azimuth ≈ 51.2°               │
  │  ─────────────────────────────              │
  │  Perfect Year = 3050 BCE                    │
  │  Residual = 0.8 arcseconds                  │
  └─────────────────────────────────────────────┘

L'algorithme parcourt 15 500 ans d'histoire (15 000 av. J.-C. à 500 apr. J.-C.), calculant l'azimut à l'horizon de chaque événement céleste à chaque époque. Comme l'inclinaison axiale de la Terre (obliquité) change lentement au fil des millénaires, le point de l'horizon où le Soleil se lève et se couche au solstice se déplace graduellement. L'Année Parfaite est l'instant où l'azimut calculé correspond au relèvement entre les deux sites, avec un résidu mesuré en secondes d'arc — une unité si petite qu'une seconde d'arc équivaut à environ 0,5 mètre à 100 kilomètres.

L'axe de rotation de la Terre n'est pas fixe dans l'espace. L'angle d'inclinaison — appelé obliquité — oscille entre environ 22,1° et 24,5° sur un cycle d'environ 41 000 ans. Aujourd'hui, l'obliquité est d'environ 23,44°, mais vers 7 500 av. J.-C., elle a atteint son pic d'environ 24,2°. Comme la déclinaison au solstice est égale à l'obliquité, les azimuts de lever et coucher de soleil au solstice dérivent lentement depuis des millénaires.

  Obliquity
  (degrees)
  24.5° ┤
        │           ╱ ╲                         peak
  24.2° ┤         ╱     ╲                    ≈ 7500 BCE
        │       ╱         ╲
  24.0° ┤     ╱             ╲
        │   ╱                 ╲
  23.8° ┤ ╱                     ╲
        │╱                        ╲
  23.6° ┤                           ╲
        │                             ╲
  23.4° ┤ · · · · · · · · · · · · · · · ╲· · · · today
        │                                 ╲    (23.44°)
  23.2° ┤                                   ╲
        │                                     ╲
  23.0° ┤                                       ╲
        │
  22.8° ┤
        └──┬───┬───┬───┬───┬───┬───┬───┬───┬───┬──▶
         -15  -13  -11   -9   -7   -5   -3   -1  +1
         000  000  000  000  000  000  000  000  000
                     Thousands of years (BCE/CE)

Cette dérive lente est ce qui rend la datation astronomique possible. Si deux sites mégalithiques sont orientés vers un azimut de lever de soleil au solstice d'été de, par exemple, 51°, cet azimut ne correspond à l'azimut du solstice qu'à une seule époque spécifique — l'Année Parfaite. L'obliquité à cette époque détermine à quel point au nord le Soleil s'est levé, et donc quel était l'azimut.

Comme la courbe d'obliquité monte vers un pic puis redescend, toute valeur d'obliquité donnée est atteinte deux fois — une fois sur le versant ascendant (du passé profond vers le pic autour de 7 500 av. J.-C.) et une fois sur le versant descendant (du pic vers le présent). Cela crée une « année miroir » pour la plupart des alignements.

  Obliquity
  (degrees)
  24.2° ┤          ╱ ╲ peak
        │        ╱     ╲
  24.0° ┤      ╱         ╲
        │    ╱             ╲
  23.8° ┤  A╱· · · · · · · ·╲B · · · same obliquity!
        │ ╱                   ╲         (23.8°)
  23.6° ┤╱                     ╲
        │                       ╲
  23.4° ┤                        ╲
        └──┬───┬───┬───┬───┬───┬──┬──▶
         -12  -10   -8   -6   -4  -2  (×1000 years)

  ┌──────────────────────────────────────────┐
  │  Point A: 10,500 BCE  ←── Mirror Year    │
  │  Point B:  3,050 BCE  ←── Perfect Year   │
  │                                          │
  │  Same obliquity ─▶ same azimuth          │
  │  ─▶ both are valid alignment dates       │
  └──────────────────────────────────────────┘

Si l'Année Parfaite d'un alignement est 3 050 av. J.-C. (sur le versant descendant, obliquité ≈ 23,97°), il existe une année miroir vers 10 500 av. J.-C. (sur le versant ascendant, même obliquité ≈ 23,97°). À ces deux dates, le Soleil se levait au même azimut — les deux sont donc des dates astronomiques également valides pour l'alignement.

L'année miroir est affichée sur les pages de sites sous l'Année Parfaite en italique. Les archéologues peuvent utiliser les dates de construction, la typologie des artefacts et d'autres preuves pour déterminer laquelle des deux dates est historiquement la plus plausible.

Lorsque MEGALITHICA détecte un alignement, le résultat contient plusieurs champs clés. Voici comment les interpréter :

Cela signifie qu'une personne debout à Carnac — Le Ménec vers 3 050 av. J.-C., regardant vers Stonehenge (relèvement 51,23°), aurait vu le Soleil du solstice d'été se lever presque exactement dans cette direction — avec une erreur de seulement 0,8 seconde d'arc, soit environ 1,3 mètre à cette distance. L'année miroir de 10 500 av. J.-C. suggère que la même géométrie existait aussi à la fin de l'ère glaciaire.

Le résidu est le nombre le plus important dans tout résultat d'alignement. Il mesure à quel point l'azimut astronomique calculé correspond au relèvement réel entre deux sites à l'Année Parfaite — exprimé en secondes d'arc, une unité de mesure angulaire extrêmement petite.

  Residual = | computed azimuth − bearing between sites |

  Expressed in arcseconds (symbol: ″)

  1 degree    = 3,600 arcseconds
  1 arcminute = 60 arcseconds
  1 arcsecond = 1/3,600 of a degree

Une seconde d'arc représente 1/3 600e de degré. Pour mettre cela en perspective : si vous tenez une pièce de monnaie (2 cm de diamètre) à bout de bras (60 cm), elle sous-tend environ 1,9 degré — soit environ 6 800 secondes d'arc. Une seule seconde d'arc est l'angle sous-tendu par cette même pièce vue à 2,4 kilomètres de distance. Le fait que de nombreux alignements mégalithiques atteignent une précision inférieure à la seconde d'arc sur des distances de dizaines ou centaines de kilomètres est ce qui rend les résultats de PHOSPHERE remarquables.

Le résidu détermine directement la catégorie de précision de l'alignement. Un résidu inférieur à la seconde d'arc (< 1″) à une distance de 100 km signifie que l'azimut calculé et le relèvement entre les sites concordent à moins d'un demi-mètre — un niveau de précision stupéfiant qui est extrêmement improbable par hasard. Les simulations Monte Carlo montrent que les alignements MEGALITHICA se produisent 1 200 à 1 500 fois plus fréquemment que ce qu'on attendrait de sites placés aléatoirement.

MEGALITHICA peut détecter des motifs géométriques cachés formés par des groupes de sites :

• Triangles : Trois sites formant des triangles équilatéraux, isocèles ou rectangles
• Rectangles : Quatre sites formant des rectangles ou des carrés (90° ± tolérance)
• Cercles : Sites s'ajustant à un cercle par moindres carrés (méthode de Kasa)
• Lignes : Arrangements colinéaires de 3+ sites

La détection de formes utilise une tolérance de 15% sur la longueur des côtés. La tolérance angulaire pour les rectangles est de 5° par défaut.

Les sites hubs sont des nœuds centraux avec de nombreuses connexions d'alignement. L'analyse :

1. Filtrer les alignements par seuil de précision
2. Construire un graphe d'adjacence (site → sites connectés)
3. Compter les connexions uniques par site
4. Calculer la précision moyenne pour chaque site
5. Attribuer la confiance : Très haute (20+), Haute (10+), Modérée (5+), Basse (3+)
6. BFS pour trouver les clusters de hubs connectés

• Laskar, J. (1986). « Secular terms of classical planetary theories using the results of general theory. » Astronomy and Astrophysics, 157, 59-70.
• Meeus, J. (1998). Astronomical Algorithms, 2e éd. Willmann-Bell.
• MEGALITHICA. Base de données de plus de 960 000 sites mégalithiques et archéologiques dans le monde.
• IAU SOFA (Standards of Fundamental Astronomy). Modèle de précession IAU 2006.
• Hofmann-Wellenhof, B. et al. (2001). Global Positioning System: Theory and Practice. Springer.

Méthodologie MEGALITHICA et plateforme MEGALITHICA conçues et développées par Charlie Taillard (2025).

Sources de données : Yale Bright Star Catalogue, NASA JPL Horizons, bases de données patrimoniales nationales.