THE LEF MODEL · A SERIES OF TWELVE PAPERS
Lunar Elevation & Fission

7. Les pôles

Pourquoi les pôles sont plus denses
Parties A–E

Zohar Duran
Chercheur indépendant, Ramla, Israël
ORCID https://orcid.org/0009-0005-3194-9387
Résumé

Analyse de la densité des cratères : calottes polaires contre ceinture équatoriale

D’après les données du catalogue Robbins, on observe un écart significatif de densité des cratères de diamètre égal ou supérieur à 20 km entre les régions polaires et la ceinture équatoriale de la Lune. Dans le comptage brut, les calottes polaires (dans un rayon de 30° autour de chaque pôle) présentent une densité de 219 cratères par million de km², tandis que la ceinture équatoriale présente une densité de 166 cratères par million de km², ce qui donne un rapport de densité de 1.32. Pour les cratères de diamètre égal ou supérieur à 40 km, ce rapport monte à 1.51.

Le comptage brut est affecté par un biais structurel, car la plupart des maria et des plaines basaltiques – caractérisés par des fonds jeunes et pauvres en cratères – sont concentrés dans la région de l’équateur. Afin de réaliser une comparaison neutre entre surfaces anciennes uniquement, on applique une procédure de « comptage propre », dans laquelle les surfaces des maria et du bassin Pôle Sud–Aitken (South Pole–Aitken - SPA) sont retirées du calcul dans les deux régions.

Dans le comptage propre, on établit un rapport de densité de en faveur des pôles pour l’ensemble des cratères de diamètre supérieur à 20 km. Ce rapport présente une tendance croissante en fonction du diamètre du cratère, comme le détaille le tableau suivant :

Plage de diamètre des cratères (km)Rapport de densité (pôles contre équateur)
20–301.09
60–1001.53

Mécanisme dynamique : modèle d’aplatissement, de corotation et d’arrondissement de la jeune Lune

L’hypothèse géophysique présentée rejette l’idée que les cratères se soient formés par une chute uniforme d’astéroïdes ou de matériau terrestre de retour (comme proposé dans l’Article 6), car une chute uniforme ne peut pas produire une concentration de densité bipolaire le long d’un axe de symétrie unique. Le schéma actuel s’explique par l’évolution morphologique de la jeune Lune au cours des étapes de sa séparation d’avec la Terre :

1. Phase de rassemblement et d’aplatissement initial : une masse sialique de faible densité, d’environ 1,300 kg/m³, s’est formée à l’extrémité d’une structure en forme de poire, a commencé à se déplacer vers l’extérieur et s’est aplatie sur une base initiale large de 6,000 km.

2. Expansion vers la ligne de corotation : à mesure que la masse se déplaçait vers l’extérieur, sa base s’est réduite à une largeur de 3,000 km. Le corps est resté verrouillé en rotation avec la Terre, avec une période d’environ 5.5 heures. À proximité de la ligne de corotation – la distance dynamique à laquelle la vitesse de rotation égalait la vitesse orbitale libre – le corps était soumis à deux vecteurs de contrainte opposés : la gravité terrestre tirait vers l’intérieur (côté proche) et la force centrifuge tirait vers l’extérieur.

3. Géométrie triaxiale et volume étiré : sous ces contraintes, le corps a pris une forme triaxiale aplatie et étirée (rapport d’axes d’environ 1.8 : 1.4 : 1), définie comme longue en direction de la Terre, moyenne le long de la trajectoire de rotation et courte selon l’axe des pôles. Dans cette phase gonflée, la densité du corps s’élevait à environ 1,700 kg/m³. La plupart des grands cratères se sont formés près du sommet de cette phase, sous l’effet de matériau éjecté puis retombé.

4. Libération dynamique et arrondissement d’une croûte rigide : au moment où le centre de masse de la Lune a franchi la ligne de corotation, la force de compression sur la base s’est annulée. La Lune s’est libérée de la Terre à environ 210 km au-delà de la ligne de corotation. La libération de l’étirement a conduit à une diminution spontanée du volume, à l’arrondissement du corps et à la hausse de la densité jusqu’à 2,600 kg/m³, alors que l’enveloppe extérieure était déjà rigide.

Ce passage géométrique d’un corps très aplati à un corps sphérique a provoqué une contraction de la surface, ce qui a entraîné un resserrement mécanique des cratères existants. Le taux de resserrement était maximal aux pôles et minimal à l’équateur. Les cratères anciens, formés près du pic d’aplatissement, ont subi le maximum de contraction et de resserrement (), tandis que les cratères plus tardifs, formés lorsque le corps était déjà proche de sa forme sphérique, n’ont presque pas été resserrés ().

Tests empiriques et résultats statistiques (Partie D)

Voici la synthèse des quatre tests principaux effectués dans la Partie D pour examiner la validité du modèle de contraction et d’arrondissement :

N° du testDéfinition du test scientifiqueRésultats empiriques et valeurs numériquesStatut
1Rapport de densité pôles/équateur en comptage propreRapport de pour l’ensemble des cratères ( km). Hausse de 1.09 à 1.53 en fonction du diamètre.Réussi (conforme à la prédiction de l’ordre de formation).
2Variation du diamètre des cratèresLes cratères ont conservé leur morphologie et leur taille d’origine sur une croûte rigide. La médiane polaire est 1.06 fois plus grande que la médiane équatoriale (intervalle de confiance : 1.01–1.12).Réussi (confirme une richesse en grands cratères aux pôles).
3Distribution des cratères jeunes (post-arrondissement)Des rapports de 0.92 et 1.12 ont été mesurés dans deux échantillons indépendants.Cohérent (en attente d’un élargissement de l’échantillon daté ; aucune preuve significative d’excédent).
4Symétrie entre les deux pôlesAprès retrait du bassin SPA, on a compté 300 cratères par million de km² à un pôle contre 305 cratères par million de km² au pôle opposé.Réussi (écart à l’intérieur de la marge d’erreur de comptage).

En outre, cinq tests morphométriques complémentaires ont été effectués (examen des plus grands cratères à surface égale, forme de la distribution des tailles, comparaison des cratères anciens et frais, degré d’ellipticité, et transition structurelle du cratère simple au cratère complexe), destinés à vérifier si les cratères eux-mêmes se sont physiquement contractés avec la surface. Quatre tests ont établi de manière univoque que le diamètre des cratères a été conservé, et aucun test n’a montré de signe de contraction des cratères eux-mêmes, résultat qui confirme que les cratères ont été transportés comme des éléments rigides et intacts sur la croûte qui s’est contractée.

Ce que propose l’article

Synthèse de la proposition scientifique : modèle de déformation structurelle et d’arrondissement de la jeune Lune

L’article propose un mécanisme géophysique et dynamique expliquant l’excédent de cratères aux pôles de la Lune par rapport à son équateur. Selon le modèle, la densité variable des cratères n’est pas le résultat d’impacts non uniformes, mais d’une contraction mécanique de la croûte lunaire lors de la phase de passage d’un corps triaxial très aplati à un corps sphérique rigide.

1. L’affirmation centrale : l’évolution morphologique de la Lune

2. Les indicateurs numériques et les prédictions scientifiques

Diamètre du cratère (km)Rapport de densité empirique (pôles contre équateur)
20–301.09
30–401.15
40–601.26
60–1001.53
100–300

3. Les exigences structurelles et mécaniques du modèle

Pour que le mécanisme proposé soit valide, le système géophysique doit satisfaire aux conditions mécaniques suivantes :

Domaine d’exigenceDescription mécanique de l’exigenceRéférence dans l’article
Radial et structurelUn intérieur suffisamment mou pour permettre l’écoulement et l’arrondissement pendant la phase attachée, sous une croûte extérieure visqueuse/rigide.C.7
Conservation morphologiqueUne croûte (enveloppe) suffisamment rigide pour préserver l’intégrité des cratères et leur diamètre d’origine pendant que la surface sous-jacente se raccourcissait aux pôles.C.7
Soutien mécanique lors de la libérationUne croûte dotée d’une résistance mécanique capable de porter et de stabiliser le corps lors du passage géométrique d’une forme de 1.20 à un corps rond à la libération.C.1
Résistance aux contraintes de tractionLa croûte doit résister aux contraintes de traction résiduelles subsistant au-delà de la ligne de corotation, jusqu’à une distance dynamique d’environ 20,100 km.C.6
Absorption de surface (fermeture)La fermeture de la surface pendant l’arrondissement – estimée à environ un quart de la surface moyenne et à environ la moitié de la surface des pôles – a été absorbée au sein de la structure crustale.C.7, E.3
Non-contraction des extrémités libresL’extrémité éloignée de la Terre (le sommet du dôme) a agi mécaniquement comme une extrémité libre et n’a pas subi de contraction ; elle est donc cohérente avec les indicateurs de la face cachée.A.9

L’environnement de calcul (a) : données du corps et dynamique orbitale

Voici la définition des paramètres physiques, des distances et des états dynamiques qui servent de base au calcul du modèle de la Lune en formation (modèle LEF) pendant la phase attachée :

Grandeur ou condition scientifiqueValeur numériqueSignification physique et justification dynamiqueSource / référence
Période de rotation pendant la phase attachée (aplatissement maximal)Environ 5.5 heures ( à la libération)La Terre et le corps attaché tournent ensemble en raison d’une contrainte mécanique. La poire allongée exerce un couple qui ralentit la Terre, tandis que le noyau qui se concentre l’accélère. Cette vitesse de rotation élevée est le mécanisme physique qui a aplati le corps jeune et mou.Article 2
Période de rotation après verrouillage sur l’orbite heures à la libération ; 14 heures à 30,000 km ; 41 heures à 60,000 kmAprès la libération : la rotation propre suit l’orbite et ralentit à mesure que le corps s’éloigne. À ce stade, le corps est déjà rond et son enveloppe rigide ; les forces de marée et de rotation ne soulèvent qu’un petit renflement (rapport d’axes de 1.03 à 30,000 km et de 1.00 à 60,000 km). L’aplatissement précoce a entièrement disparu avant la libération.Article 4
Rayon de la ligne de corotation16,331 km (période de 5.5 heures) ; 16,399 km ( km avec correction à deux corps pour une période de 5.77 heures)Point d’équilibre où l’attraction terrestre égale la contribution centrifuge dans le référentiel en rotation. Au-delà de cette distance ( rayons terrestres), le vecteur de contrainte nette est dirigé vers l’extérieur et le corps se déplace à la vitesse d’une orbite libre.Article 2
Densité du corps au fil des phases kg/m³ ; rayons de la sphère équivalente : kmPassage continu d’une masse sialique légère sur l’extrémité de la poire à un corps étiré, gonflé et volumineux au pic d’aplatissement , qui se réduit et durcit à la libération jusqu’à la densité actuelle. Un corps léger est plus sensible à l’aplatissement pour une même vitesse de rotation.Article 2 ; Article 3 §B.2 ; Annexes C.1, C.3
Limite de Roche rigide km à une densité de 1,300 kg/m³ (10,314 km à une densité de 2,600 ; variable en fonction de )Distance au-delà de laquelle les forces de marée sont plus faibles que la gravité propre du corps, ce qui permet le maintien d’une forme propre par gravité et résistance mécanique. Le corps s’aplatit sur une base qui se réduit de 6,000 à 3,000 km. Un corps fluide libre d’une densité de 2,600 kg/m³ aurait nécessité une distance de 19,971 km ( d’environ 20,100 km) pour se maintenir – distance qui n’a pas été atteinte avant la libération, au cours de laquelle le corps était maintenu par son enveloppe rigide.Article 2 ; Annexes C.2, C.5
Rayon de libération mécaniqueLigne de corotation à km sur l’équateur (17,287 km à ) ; libération environ 210 km au-delàLe début de la séparation survient lorsque le centre de masse atteint la ligne de corotation. La libération finale s’effectue environ 210 km au-delà, lorsque environ 58% de sa masse a franchi la ligne. Le contact mécanique offre une légère résistance et le corps est entraîné vers l’extérieur. À ce stade, le corps est rigide et rond (rapport d’axes ).Articles 2 et 4
Plage d’aplatissement et d’arrondissementAplatissement : de l’extrémité jusqu’à la réduction de la base à 3,000 km près de la corotation. Arrondissement : du franchissement de la ligne par le centre de masse jusqu’à la libération.Le corps s’aplatit au cours de son déplacement vers l’extérieur et atteint le pic d’aplatissement près de la ligne de corotation. Au franchissement de la ligne, la force de compression sur la base s’annule, l’étirement se libère, et le corps s’arrondit, diminue de volume et se détache. Tard dans l’arrondissement, environ 75% de l’aplatissement avait disparu, et à la libération le corps est rond.Article 3 ; la présente discussion
Pression source et élévation du matériauPression source : 150 GPa ; pression attribuée à l’élévation : 140 GPaPression interne de la Terre, générée par la rotation et le tri des densités, qui alimente et pousse le matériau à travers la structure de la poire jusqu’à la Lune. Le détachement mécanique interrompt le flux de pression, expose les évents de la poire, libère l’énergie et empêche la poursuite de son allongement.Article 3 §A ; Annexe C.6

L’environnement de calcul (b) : données de l’enregistrement géométrique et du comptage empirique

Voici la définition des indicateurs, des bases de données et des hypothèses géométriques de base utilisés pour calculer les surfaces et convertir les comptages de cratères en densités normalisées :

Grandeur ou condition scientifiqueValeur numérique / géométriqueSignification physique et justification dynamiqueSource / référence
Masse statique de la Terre kgMasse cible de la Terre utilisée dans les équations de gravité, d’orbite et de marée, une fois achevé le transfert de masse vers le corps lunaire.Article 2
Fixation de la forme de la Lune (aplatissement fossile)Plage de distance : avant 16 rayons terrestres ; aplatissement fossile : Définition de la distance orbitale à laquelle la forme globale de la Lune s’est figée. Cette ligne n’impose pas une solidification simultanée de tout l’intérieur lunaire. L’aplatissement fossile représente l’écart actuel de la Lune par rapport à une sphère parfaite, conservé en conséquence de la rotation rapide précoce.Annexe 4 §C.3
Base de données de l’enregistrement des cratèresCatalogue Robbins (2019) : 6,973 cratères ( km) ; 2,329 cratères ( km)Les deux groupes de contrôle morphométriques servant à comparer les densités par unité de surface. Le comptage de base est effectué sans masques géologiques ; la procédure de « comptage propre » (retrait des maria et des bassins) est définie séparément à la section A.11.Article 11 §B.8
Système d’axes et détermination des pôlesPoint subterrestre : ( N, E) ; pôles géographiques : Les ceintures de distance angulaire et la géométrie spatiale sont calculées directement à partir du point tourné en permanence vers la Terre. Les calottes polaires sont définies mathématiquement autour des pôles Nord et Sud afin de distinguer absolument les zones de comptage.Article 11 §B.8 ; carte de l’Article 3
Rayon lunaire pour les calculs de surface1,737.4 kmLe rayon moyen officiel (IAU), utilisé comme constante géométrique fixe pour convertir les nombres bruts de cratères en unités de densité standard (cratères par million de km²).Calculs de surface de l’article

Dynamique des phases de détachement : évolution morphologique et bilan des contraintes

Voici un déroulé chronologique et mécanique des phases de développement de la jeune Lune – depuis la formation de la masse initiale sur l’extrémité de la structure en poire jusqu’à son état moderne – avec la définition du bilan des forces, des variations de densité et de la transformation géométrique :

Phase dynamiqueDistance orbitale estimée (km)Densité ()Période de rotation propre (heures)Géométrie et rapport d’axes ( ou )État mécanique et bilan des contraintes du système
Masse rassemblée à l’extrémité de la poire1,3005.50—Matériau sialique reposant sur l’extrémité d’un long prolongement de la Terre (base large de 6,000 km) ; début de la sortie et développement d’un aplatissement initial.
Corps autonome attaché (limite de Roche rigide)5.50Tendance à l’augmentation de l’aplatissementLe corps atteint la distance à laquelle sa gravité propre et sa résistance l’emportent sur les forces de marée (). Le corps maintient sa propre forme mais reste attaché ; il s’aplatit et s’étire vers l’extérieur tandis que la base de la poire se réduit de 6,000 à 3,000 km.
Pic d’aplatissementPrès de la ligne de corotation5.50Le corps est étiré, gonflé et atteint un maximum de déformation géométrique. Il est soumis à un équilibre de contraintes extrême : la gravité terrestre tire vers l’intérieur (sur le côté proche) contre l’entraînement centrifuge vers l’extérieur. Mécaniquement comprimé sur une base de 3,000 km.
Début de l’arrondissement (franchissement de la ligne)Au-delà de la ligne de corotation1,7005.551.40Le centre de masse franchit la ligne de corotation. La force de compression nette sur la base s’annule ; la partie œuf (axe ) se libère, et le volume global commence à diminuer.
Phase d’arrondissement intermédiaire aPhase intermédiaire2,0005.621.30Libération continue des contraintes de traction ; le corps s’épaissit et se resserre progressivement.
Phase d’arrondissement intermédiaire bPhase intermédiaire2,3005.701.24Poursuite de l’effondrement mécanique vers l’intérieur des axes étirés et hausse de la densité du matériau du noyau.
Tard dans l’arrondissementPrès de la ligne de libération2,6005.771.20Environ 75% de l’excédent maximal d’aplatissement du grand axe par rapport à l’axe polaire a disparu (contraction d’environ 0.6 du rayon). L’enveloppe extérieure se fige, se rigidifie et fixe l’enregistrement des cratères.
Libération mécanique (détachement) (environ 210 au-delà de la ligne)2,6005.77Ouverture de la couche de contact et séparation finale d’avec la Terre. Le corps se libère dans un état presque entièrement rond, protégé par une croûte rigide qui porte les cratères intacts.
Phase d’éloignement a30,0002,60014.001.03La rotation propre se verrouille sur l’orbite et ralentit progressivement. Le 1.03 n’est pas un aplatissement fossile mais un petit renflement élastique soulevé sur le corps rigide par des forces de marée et de rotation encore rapides.
Phase d’éloignement b60,0002,60041.001.00Atténuation du renflement dynamique jusqu’à une sphère parfaite à mesure que la distance augmente et que la rotation propre ralentit.
L’état moderne (aujourd’hui)3,34427.3 jours1.00Le corps est rond (il porte un aplatissement fossile mineur de , figé avant la distance de 16 rayons terrestres).

Chronologie des événements et fenêtre de formation des cratères

Le modèle géodynamique définit une échelle de temps totale d’environ 600 millions d’années depuis le début de la formation du renflement originel jusqu’à la phase de libération finale (datée de 3.873 milliards d’années avant le présent), divisée en trois sous-périodes qui constituent la fenêtre morphogénétique exclusive des cratères de l’article :

Définitions de base et indicateurs mécaniques pour la lecture du tableau

Critère de stabilité fluide libre (indicateur )

L’évaluation de l’étirement mécanique est effectuée au moyen de l’indicateur du système :

Introduction — la question et la manière dont elle est examinée

1. Méthodologie de recherche : définition des zones de comptage et des indicateurs empiriques

Afin d’effectuer la comparaison morphométrique et d’éviter les biais dus à l’incomplétude de l’identification des petits cratères, un seuil minimal strict a été fixé pour le diamètre des cratères comptés. Les données reposent sur le catalogue Robbins (2019) et sont analysées selon les unités géométriques suivantes :

2. Résultats du comptage brut comparés à la procédure de comptage propre

Le comptage empirique initial montre un excédent de densité significatif dans les régions polaires, qui s’accentue à mesure que le diamètre des cratères examinés augmente. Toutefois, l’analyse géologique montre que le comptage brut souffre de biais structurels sévères qui exigent une normalisation mathématique :

Indicateur de cratères examinéDensité dans la calotte nordDensité dans la calotte sudDensité dans la ceinture équatorialeRapport de densité (pôles contre équateur)Statut méthodologique
Diamètre km (brut)227 par million de km² (577 cratères sur ~2.54M km²)211 par million de km²166 par million de km²1.32 (densité combinée : 219)Biaisé (inclut maria et SPA)
Diamètre km (brut)———1.51Biaisé (tendance à l’accentuation chez les grands)
Ensemble des cratères (comptage propre)———Le comptage normalisé officiel

3. L’hypothèse dynamique de départ et le mécanisme de contraction mécanique

1. Prédictions du modèle géométrique et mécanisme de resserrement différentiel

Le degré de resserrement mécanique des cratères n’est pas une fonction directe du seul âge chronologique, mais un enregistrement géométrique de la forme du corps et de son aplatissement dynamique au moment de la formation du cratère. Comme du matériau a été éjecté et est retombé tout au long des phases de transition, le niveau de déformation structurelle a été conservé dans la morphologie de la croûte selon les principes suivants :

2. Corrélation empirique : prédiction du modèle contre indicateurs du comptage propre

Voici une comparaison stricte et normalisée entre les rapports de densité théoriques prédits (dérivés de la forme de la Lune au moment de la formation) et les indicateurs empiriques effectivement mesurés dans le cadre du comptage propre (retrait des maria et du bassin SPA), par classe de diamètre :

Phase dynamique attribuéeForme théorique du corps ()Rapport de densité prédit (calottes/ceinture)Classe de diamètre des cratères (km)Rapport de densité empirique mesuréForme reconstituée directement (lecture inverse)
Pic d’aplatissement structurel1.49100–300
Franchissement de la ligne de corotation (libération de l’œuf)1.401.3460–100
Phase d’arrondissement intermédiaire1.301.2640–60
Phase d’arrondissement tardive1.201.1730–40
Après la libération (corps rond rigide)1.01.020–30
Ensemble de la population ( km)1.21Ensemble des cratères —

3. Conservation du diamètre des cratères sur une croûte en contraction

L’une des questions centrales dans l’élaboration du modèle géophysique est la suivante : si la surface polaire s’est contractée mécaniquement pendant l’arrondissement, pourquoi les cratères eux-mêmes ne se sont-ils pas contractés avec elle (ce qui les aurait fait passer sous le seuil de comptage de 20 km et aurait fait disparaître l’effet) ?

4. Structure et ordre des chapitres de l’article

La suite de l’analyse scientifique de l’article est organisée selon la structure suivante :

PART A

Partie A — L’observation empirique

La Partie A présente les résultats empiriques bruts tels qu’ils sont comptés à la surface, avant toute interprétation physique ou tout modèle de formation. L’observation repose sur un comptage normalisé : le nombre de cratères dans chaque cellule de surface est divisé par la surface géométrique de cette cellule.

Dans cette partie, quatre comptages indépendants sont effectués sur la même base de données et selon la même méthodologie :

A.1 Densité des cratères par latitude (données brutes)

L’analyse spatiale de cette section vise à déterminer s’il existe une non-uniformité de la densité des cratères le long des axes de latitude géographique.

Critères méthodologiques et calibrage du comptage :

1. Complétude de la base (catalogue Robbins) : le comptage est limité aux cratères d’un diamètre ajusté égal ou supérieur à 20 km. Cette limite est stricte et garantit une complétude totale du comptage à toutes les latitudes : elle neutralise les biais optiques (conditions de faible éclairage et ombres longues aux pôles qui masquent les petits cratères) et filtre les cratères secondaires (amas et chaînes de débris de diamètre inférieur à 20 km, qui produiraient une densité locale trompeuse).

2. Groupe de contrôle et confirmation de la tendance : une seconde série est effectuée pour les cratères de diamètre égal ou supérieur à 40 km afin de vérifier que le schéma spatial ne repose pas sur la population des plus petits cratères limites.

3. Absence initiale de masques : le premier comptage présente la surface telle quelle, sans omission artificielle des maria, considérant qu’ils font partie de la prédiction structurelle (l’effet de leur retrait est analysé séparément dans la Partie A.3).

4. Calcul des surfaces géométriques différentielles : la Lune est divisée en 12 bandes de latitude d’une largeur constante de . La surface de chaque bande est calculée selon l’équation spatiale :

où le rayon moyen retenu est (norme IAU), et représentent les limites de la bande.

D’après cette équation, la surface d’une bande sur une sphère n’est pas déterminée par sa largeur angulaire mais par son étendue le long de l’axe des pôles (comme la hauteur d’une étiquette sur un cylindre). Près de l’équateur, la bande est haute le long de l’axe, tandis que près des pôles la surface s’incurve fortement et la circonférence de l’anneau est petite. En conséquence, les surfaces des bandes varient considérablement : la bande équatoriale () s’étend sur 4.909 millions de km², tandis que la bande polaire () ne s’étend que sur 0.646 million de km² (rapport d’environ 7:1). Cette différence géométrique impose l’usage de l’indicateur de densité (cratères par 10⁶ km²) et non de nombres absolus.

Synthèse des données de comptage des cratères par bandes de latitude (A.1)

Voici les données numériques complètes et relatives mesurées dans les 12 bandes de latitude globales, sans application de masques géologiques :

SecteurDescription géographique de la zoneBande de latitude ()Surface calculée ()Cratères ()Densité par Rapport à l’équateur ()Cratères ()Densité par Rapport à l’équateur ()
◐Sud — adjacente au pôle75–90° S0.6461392151.2953821.54
◐Sud — extérieure à la calotte polaire60–75° S1.8953972101.26132701.31
◐Sud — voisine de la calotte polaire45–60° S3.0146932301.38230761.44
◐Sud — latitudes moyennes30–45° S3.9288062051.23272691.30
◐Sud — basses latitudes15–30° S4.5748171791.07287631.18
◐Équateur — demi-bande sud0–15° S4.9098731781.07283581.08
◐Équateur — demi-bande nord0–15° N4.9097611550.93239490.92
◐Nord — basses latitudes15–30° N4.5746931510.91233510.96
◐Nord — latitudes moyennes30–45° N3.9286921761.06207530.99
◐Nord — voisine de la calotte polaire45–60° N3.0145251741.05171571.07
◐Nord — extérieure à la calotte polaire60–75° N1.8953832021.21157831.56
◐Nord — adjacente au pôle75–90° N0.6461943001.80651011.89

Analyse et conclusions du profil spatial

La lecture des colonnes de densité et du profil graphique fait ressortir les résultats principaux suivants :

Échelle de couleurs : à gauche, « low density », densité faible (jaune doux) ; à droite, « high density », densité élevée (rouge doux).

איור 1

Figure 1. Densité des cratères par latitude, cratères de 20 km et plus, dans les douze bandes de 15° (A.1) ; le catalogue entier, sans masque. Les deux pôles sont les plus élevés. (Sur le graphique : axe vertical — cratères ≥ 20 km par 106 km², 100–300 ; axe horizontal — bande de latitude, le pôle Sud à gauche et le pôle Nord à droite, de 75–90S à 75–90N ; la ligne bleu foncé — la densité dans chaque bande.)

A.2 Densité des cratères par longitude (données brutes)

Alors que la section A.1 a trouvé un excédent de densité significatif dans les régions polaires, cette section examine l’uniformité de la répartition spatiale dans la direction perpendiculaire – le long de la circonférence globale de la Lune. L’objectif de l’analyse est d’isoler et d’éviter un mélange statistique entre l’effet de latitude (le sujet central de la recherche) et d’éventuelles asymétries entre hémisphères le long des axes de longitude.

Critères méthodologiques et géométrie des fuseaux :

1. Structure des cellules de comptage : la surface de la Lune a été divisée en 12 bandes de longitude (fuseaux) d’une largeur constante de . Chaque bande s’étend de façon continue d’un pôle à l’autre et inclut toutes les latitudes (comme un quartier d’orange).

2. Égalité des surfaces intrinsèque : contrairement aux bandes de latitude, tous les fuseaux de longitude sur une sphère sont de surface parfaitement égale. Cette forme géométrique rend inutile l’application de facteurs de correction ou de normalisation à la taille de la cellule, et permet de comparer directement le nombre absolu de cratères (les comptages bruts) et les rapports de densité.

3. Dimensions de surface : la surface de chaque bande est calculée exactement selon :

Pour un rayon , la surface de chaque fuseau de longitude est constante et égale à 3.161 millions de km² (1/12 de la surface totale de la Lune).

4. Le système de coordonnées dynamique : la bande correspondant au secteur Est contient le point subterrestre () – le centre de la face tournée en permanence vers la Terre en raison du verrouillage de marée.

5. Définition des composantes orbitales (« avant » contre « arrière ») : les termes sont appliqués selon les définitions mécaniques usuelles :

6. Base de normalisation fixe : le point de référence (1.00) reste strictement identique à celui fixé à la section A.1 – la densité de la bande équatoriale brute complète (). Sa valeur est de 166.436 cratères par million de km² pour le groupe , et de 53.170 pour le groupe . Le maintien de cette base fixe permet une comparaison directe entre les quatre tableaux d’analyse de la Partie A (de A.1 à A.4).

Synthèse des données de comptage des cratères par bandes de longitude (A.2)

Voici les données numériques complètes mesurées dans les 12 fuseaux de longitude globaux, sans application de masques géologiques :

SecteurDéfinition géographique et morphologique du secteurBande de longitude ( est)Surface calculée ()Cratères ()Densité par Rapport à l’équateur ()Cratères ()Densité par Rapport à l’équateur ()
◐Autour du méridien du point subterrestre15–3453.1614791520.91168531.00
◐Face proche — à l’est du point subterrestre45–153.1615051600.96135430.80
◐Face proche — vers le méridien arrière75–453.1616111931.16170541.01
◐Autour du méridien arrière ( est)105–753.1617502371.43269851.60
◐Face cachée — au-delà du méridien arrière135–1053.1617932511.51285901.70
◐Face cachée — vers le point anti-terrestre165–1353.1616672111.27245781.46
◐Autour du méridien du point opposé à la Terre195–1653.1617112251.35233741.39
◐Face cachée — au-delà du point opposé à la Terre225–1953.1617362331.40239761.42
◐Face cachée — vers le méridien avant255–2253.1616432031.22205651.22
◐Autour du méridien avant ( ouest)285–2553.1614871540.93170541.01
◐Face proche — au-delà du méridien avant315–2853.1613241020.62111350.66
◐Face proche — à l’ouest du point subterrestre345–3153.161267840.5199310.59

Analyse et conclusions du profil longitudinal global

L’examen des indicateurs et du profil graphique fait ressortir les résultats principaux suivants :

Échelle de couleurs : à gauche, « low density », densité faible (jaune doux) ; à droite, « high density », densité élevée (rouge doux).

איור 2

Figure 2. Densité des cratères par longitude, cratères de 20 km et plus, dans les douze bandes de 30° (A.2) ; le catalogue entier, sans masque. Le maximum se trouve un peu à l’est du méridien arrière, et le minimum un peu à l’ouest du point subterrestre. (Sur le graphique : axe vertical — cratères ≥ 20 km par 106 km², 50–300 ; axe horizontal — bande de longitude (° est), du point subterrestre vers l’est en passant par le méridien arrière, le point opposé à la Terre et le méridien avant, de 345–15 à 315–345 ; la ligne bleu foncé — la densité dans chaque bande.)

A.3 Densité des cratères par latitude après retrait des bassins et des maria

Alors que les sections A.1 et A.2 ont compté l’intégralité de la surface globale, cette section applique un masque géologique strict pour retirer les unités de surface dont l’enregistrement des cratères a été affecté par des facteurs endogènes (activité thermique et volcanique et effacement mécanique) sans lien avec la latitude géographique.

Justification du retrait des unités géologiques :

1. Le système des maria et des plaines basaltiques (masque LROC) : le fond des maria représente de jeunes ouvertures d’alimentation au-dessus de dômes magmatiques (Article 5), où les cratères anciens ont été fondus. Les plaines basaltiques (telles que Procellarum, Tranquillitatis, Fecunditatis, Insularum, Frigoris) représentent une croûte fondue sous l’effet de la chaleur globale rayonnée par la Terre, sans dôme d’alimentation local. Dans les deux cas, le déficit de cratères provient d’une remise à zéro thermique de l’enregistrement. Comme ces unités sont concentrées sur la face proche et aux basses latitudes, elles abaissent artificiellement la densité équatoriale.

2. Les grands bassins () : ces structures d’impact s’étendent sur de larges portions des bandes de latitude (surtout dans les régions polaires). Leur fond et leur bordure ont une histoire morphologique propre qui influe de manière indépendante sur la statistique de la bande.

3. Objectif de la procédure : isoler les portions crustales historiquement les plus stables (Highlands), où l’enregistrement est resté uniforme, afin d’obtenir un indicateur fiable de la répartition originelle des cratères.

Méthodologie d’application du masque géologique

Une valeur positive indique que la surface retirée était pauvre en cratères par rapport au reste de la bande (et que son retrait a donc augmenté la densité résiduelle).

Synthèse des données de densité des cratères masquée par latitude (A.3)

Voici les valeurs numériques exactes obtenues dans les 12 bandes de latitude après application du masque de filtrage :

SecteurDescription géographique de la zoneBande de latitude ()Surface résiduelle ()Cratères ()Densité par Rapport à la référence ()Variation de densité ()Cratères ()Densité par Rapport à la référence ()Variation de densité ()
◐Sud — adjacente au pôle75–90° S0.275843051.83+42.0%341242.32+50.7%
◐Sud — extérieure à la calotte polaire60–75° S1.0332382301.38+10.0%71691.29-1.3%
◐Sud — voisine de la calotte polaire45–60° S1.5443802461.48+7.1%109711.33-7.5%
◐Sud — latitudes moyennes30–45° S2.0674762301.38+12.2%150731.37+4.8%
◐Sud — basses latitudes15–30° S2.1944552071.25+16.1%156711.34+13.3%
◐Équateur — demi-bande sud0–15° S2.9766342131.28+19.8%215721.36+25.3%
◐Équateur — demi-bande nord0–15° N2.8805822021.21+30.3%186651.21+32.6%
◐Nord — basses latitudes15–30° N2.1534842251.35+48.4%171791.49+55.9%
◐Nord — latitudes moyennes30–45° N2.4115802411.45+36.6%171711.33+34.6%
◐Nord — voisine de la calotte polaire45–60° N1.6393752291.37+31.4%123751.41+32.3%
◐Nord — extérieure à la calotte polaire60–75° N1.7523662091.25+3.3%148841.59+1.9%
◐Nord — adjacente au pôle75–90° N0.6461943001.800.0%651011.890.0%

Analyse des effets spatiaux et du biais intrinsèque de la base statique

La comparaison du profil masqué avec le profil brut fait ressortir les enseignements scientifiques suivants :

1. Comportement des cas limites de la colonne de variation :

2. Robustesse de l’excédent polaire endogène : l’application du masque relève significativement les basses et moyennes latitudes (surtout au nord), aplanissant ainsi une partie de l’écart spatial. Toutefois, les deux bandes adjacentes aux pôles restent les zones les plus denses de la surface lunaire. Dans le groupe , le pôle Sud atteint 305 et le pôle Nord 300 par million de km², tandis que les 10 bandes intermédiaires sont plus basses (valeurs de 202 à 246). Dans le groupe , le pôle Sud domine tout le système avec une valeur maximale de 124 (rapport de 2.32). Ce résultat prouve que l’excédent polaire n’est pas un artefact de la présence des maria à l’équateur ; il est inscrit dans le sol solide lui-même.

3. Précision méthodologique sur la comparaison à l’identique : la colonne des rapports de ce tableau divise une bande masquée (numérateur ayant perdu ses maria) par un équateur brut (dénominateur les ayant conservés), produisant ainsi une exagération mathématique de l’excédent (rapports de 1.40–1.48). Si l’on compare les deux calottes polaires masquées (882 cratères sur 3.706 millions de km², soit une densité de 238.0) avec la ceinture équatoriale masquée par exactement le même masque (1,216 cratères sur 5.856 millions de km², soit une densité de 207.7), le rapport réel s’établit à 1.15. Par conséquent, pour le test scientifique final de l’article, on applique le comptage propre des chapitres A.11 et D.1 (retrait simultané des maria et du bassin SPA des deux côtés), qui fournit la valeur normalisée et stricte du rapport : .

Échelle de couleurs : à gauche, « low density », densité faible (jaune doux) ; à droite, « high density », densité élevée (rouge doux).

איור 3

Figure 3. Les mêmes douze bandes, après retrait des vingt bassins de 500 km et plus et des maria, et déduction de leur surface (A.3). Le masque relève l’équateur et laisse les pôles en place. (Sur le graphique : axe vertical — cratères ≥ 20 km par 106 km², 100–300 ; axe horizontal — bande de latitude, le pôle Sud à gauche et le pôle Nord à droite ; la ligne rouge — le profil masqué.)

A.4 Densité des cratères par longitude après retrait des bassins et des maria

Cette section applique le masque géologique strict (les 20 grands bassins de Neumann et les polygones de maria du LROC) aux 12 bandes de longitude de la section A.2. Son objectif principal est de déterminer quelle part de l’asymétrie extrême entre l’est et l’ouest provient directement de la présence des jeunes maria et des bassins, et quelle part est inscrite dans la croûte montagneuse ancienne elle-même.

Méthodologie de calcul des surfaces et calibrage du système :

1. Perte de l’égalité des surfaces des bandes de longitude : contrairement au comptage brut, où tous les fuseaux de longitude sont de taille égale, l’application du masque modifie les surfaces des bandes de manière différentielle. Une bande qui contenait une large portion de maria ou de bassins perd davantage de surface. Les surfaces résiduelles effectives s’étendent désormais sur une large plage, de 1.073 à 2.729 millions de km², ce qui impose de s’appuyer exclusivement sur l’indicateur de densité et non sur des comptages absolus.

2. Cohérence et fermeture mathématique du masque : en sommant les 12 bandes, on obtient un accord parfait avec les données de la section A.3 : les deux méthodes de découpage géométrique (en latitude et en longitude) partagent exactement la même surface crustale résiduelle, qui représente un peu plus de la moitié de la Lune (21.570 millions de km²) et contient exactement le même nombre de cratères (4,848 dans le groupe km, et 1,599 dans le groupe km).

3. Conservation de la base de référence fixe : le point de référence (valeur 1.00) reste parfaitement constant – la densité de la ceinture équatoriale brute complète sans masque (166.436 pour km, et 53.170 pour km).

Synthèse des données de densité des cratères masquée par longitude (A.4)

Voici les valeurs numériques exactes obtenues dans les 12 fuseaux de longitude après application du masque de filtrage :

SecteurDéfinition géographique et morphologique du secteurBande de longitude ( est)Surface restante ()Densité Rapport Variation ()Densité Rapport Variation ()
◐Autour du méridien du point subterrestre15–3451.8794202241.34+47.5%149791.49+49.2%
◐Face proche — à l’est du point subterrestre45–151.5103512321.40+45.5%92611.15+42.6%
◐Face proche — vers le méridien arrière75–451.8374382381.43+23.3%121661.24+22.5%
◐Autour du méridien arrière ( est)105–751.6923702191.31-7.8%131771.46-9.0%
◐Face cachée — au-delà du méridien arrière135–1052.7297252661.60+5.9%251921.73+2.0%
◐Face cachée — vers le point anti-terrestre165–1351.9954572291.38+8.6%166831.57+7.4%
◐Autour du méridien du point opposé à la Terre195–1651.3453562651.59+17.7%96711.34-3.1%
◐Face cachée — au-delà du point opposé à la Terre225–1951.9044802521.51+8.3%165871.63+14.6%
◐Face cachée — vers le méridien avant255–2251.8133692041.22+0.1%121671.26+2.9%
◐Autour du méridien avant ( ouest)285–2552.0444172041.23+32.4%152741.40+38.3%
◐Face proche — au-delà du méridien avant315–2851.7482761580.95+54.0%90510.97+46.6%
◐Face proche — à l’ouest du point subterrestre345–3151.0731891761.06+108.5%65611.14+93.4%

Analyse dynamique de l’effet du masque sur l’asymétrie longitudinale

L’examen des indicateurs résiduels et de la colonne de variation géologique fait ressortir les enseignements suivants :

1. Comblement du creux occidental extrême : la bande de longitude est – qui abrite le bassin de Mare Imbrium et la plupart des maria occidentaux de la face proche – présente la hausse la plus forte du système : un bond de densité de +108.5% (la densité fait plus que doubler, de 84 à 176 par million de km²). Ce résultat prouve que le minimum extrême mesuré à la section A.2 provenait directement de la remise à zéro de l’enregistrement par les jeunes fonds basaltiques.

2. Comportement atypique dans la zone arrière : le secteur entourant le méridien arrière ( est) est la seule bande du système où la densité résiduelle baisse après application du masque (chute de -7.8% dans le groupe des 20 km, et de -9.0% dans le groupe des 40 km). Cela signifie que les portions retirées dans ce secteur (bassins locaux) étaient plus riches et plus denses en cratères que la croûte restante autour d’elles.

3. Réduction de l’ampleur de l’écart du schéma longitudinal : le masque de filtrage réduit de façon spectaculaire l’asymétrie entre l’est et l’ouest et resserre le profil dans une plage plus étroite. Le rapport extrême entre le fuseau maximal et le fuseau minimal chute de 2.97 à 1.68 pour les cratères km (maximum de 266 dans le secteur contre minimum de 158 dans le secteur ), et de 2.88 à 1.79 pour les cratères km (92 contre 51).

4. Robustesse du schéma résiduel de circonférence : malgré la réduction mathématique, le schéma spatial de base ne disparaît pas. La face cachée au-delà du méridien arrière conserve fermement son statut de zone la plus dense dans les deux groupes de contrôle, tandis que la face proche au-delà du méridien avant reste la plus pauvre. Ce schéma résiduel est analysé plus loin dans l’article au moyen du modèle des six axes (A.9). Pour résoudre l’énigme polaire, la mesure isolée elle-même et la séparation des effets longitudinaux fournissent l’infrastructure méthodologique nécessaire.

Échelle de couleurs : à gauche, « low density », densité faible (jaune doux) ; à droite, « high density », densité élevée (rouge doux).

איור 4

Figure 4. Les mêmes douze bandes de longitude, après retrait des vingt bassins de 500 km et plus et des maria, et déduction de leur surface (A.4). Le minimum occidental remonte ; le maximum de la face cachée demeure. (Sur le graphique : axe vertical — cratères ≥ 20 km par 106 km², 50–300 ; axe horizontal — bande de longitude (° est), du point subterrestre vers l’est en passant par le méridien arrière, le point opposé à la Terre et le méridien avant ; la ligne rouge — le profil masqué.)

A.5 Synthèse et conclusions des quatre examens par bandes

Cette section rassemble et analyse les résultats combinés des quatre examens spatiaux par bandes effectués aux chapitres A.1–A.4. L’analyse vise à répondre à une question stricte et structurelle : l’excédent de densité des cratères dans les régions polaires est-il une propriété géophysique réelle de la croûte lunaire, ou s’agit-il d’un biais statistique artificiel lié au mode de découpage de la surface ou à la présence de grandes formations géologiques locales ?

La méthodologie repose sur le croisement de quatre examens complémentaires :

1. Distribution par latitude sans masque (A.1) : examen de la tendance équateur contre pôles sur la surface brute.

2. Distribution par longitude sans masque (A.2) : examen de la tendance sur la circonférence globale (est contre ouest) sur la surface brute.

3. Distribution par latitude masquée (A.3) : nettoyage des effets endogènes (maria, bassins km) le long des axes de latitude.

4. Distribution par longitude masquée (A.4) : nettoyage des effets endogènes le long des axes de longitude.

La combinaison des quatre examens permet d’isoler deux directions de découpage perpendiculaires (latitude contre longitude) tout en contrôlant le degré de dépendance de l’enregistrement vis-à-vis des jeunes surfaces basaltiques et des bassins géants.

A.5.1 Séparation des schémas : deux configurations spatiales indépendantes

Le premier enseignement scientifique de la synthèse est qu’aucun profil spatial unique ne peut expliquer la répartition globale des cratères. Les axes de latitude et de longitude révèlent deux schémas géographiques entièrement distincts et séparés :

1. Caractéristiques du profil de latitude (tendance bipolaire) :

2. Caractéristiques du profil de longitude (asymétrie de circonférence) :

Analyse méthodologique : exclusion du mélange des schémas

Le fait que l’asymétrie le long des méridiens atteigne un facteur trois pourrait conduire à affirmer que l’excédent polaire (plus faible, d’un facteur deux) n’est qu’un sous-produit ou une distorsion statistique causée par le schéma plus important de la longitude. La recherche rejette cette affirmation au moyen de l’explication géométrique suivante :

A.5.2 Le masque comme test de contrôle

Les sections A.3 et A.4 ne constituent pas un nouveau comptage d’une nouvelle population de cratères, mais une répétition méthodologique exacte des deux examens par bandes (en latitude et en longitude) sous des conditions de filtrage par un masque géologique uniforme. Ce masque retire du calcul les 20 bassins de diamètre km et les polygones de maria du LROC [A.3, A.4]. Un cratère dont le centre géométrique tombe à l’intérieur des limites du masque est omis du comptage, et la surface du secteur est réduite en conséquence dans le dénominateur de densité [A.3, A.4].

Justification scientifique du choix des surfaces filtrées :

La dynamique mathématique de l’indicateur de densité :

Une baisse du nombre absolu de cratères () à la suite de l’application du masque n’entraîne pas nécessairement une baisse de la densité résiduelle [A.3]. L’équation de densité s’effectue selon le rapport variable entre le numérateur et le dénominateur sur le même sol [A.3] :

Si le masque retire une part relative de surface supérieure à la part de cratères qu’il omet (situation typique des surfaces de maria dépourvues de cratères), la densité globale de la surface résiduelle présentera une tendance à la hausse [A.3]. Le mécanisme fonctionne comme le calcul d’un indicateur d’arbres par hectare dans une forêt contenant un lac : retirer la surface du lac (dépourvue d’arbres) du dénominateur de surface totale augmente mathématiquement la valeur de densité de la terre restante, bien qu’aucun arbre physique n’ait été ajouté au système [A.5.2].

L’application de la même opération analytique aux deux axes de découpage perpendiculaires produit un puissant outil de contrôle : si les deux effets spatiaux (en latitude et en longitude) étaient exclusivement le produit de la présence des maria et des bassins, l’application du masque aurait dû modifier les deux profils de manière cohérente et semblable. Les résultats empiriques l’excluent et montrent un comportement différentiel net [A.5.2].

A.5.3 Réponse du profil de longitude au masque de filtrage

Contrairement au profil de latitude, l’effet du masque sur le profil de longitude est spectaculaire et modifie fondamentalement l’amplitude (Mesa) brute du système [A.5.3] :

1. Convergence des valeurs extrêmes et réduction de l’amplitude le long de la circonférence :

2. Effondrement de l’asymétrie entre le méridien arrière et le méridien avant :

L’analyse statistique des bandes de longitude complètes (s’étendant d’un pôle à l’autre) révèle l’effondrement de la significativité de l’asymétrie orbitale au niveau de la bande [A.5.3] :

3. L’effet géométrique : moyennage par bandes contre concentration en calottes :

Malgré l’effondrement de la significativité au niveau de la bande longue, le schéma longitudinal ne disparaît pas entièrement de la croûte solide [A.4]. La raison tient à la limite géométrique des bandes de longitude complètes : une bande s’étendant d’un pôle à l’autre tend à diluer et à brouiller les excédents locaux en les moyennant avec des surfaces ordinaires tout le long du méridien [A.5.3].

Lorsque l’on change de méthode géométrique et que l’on passe à l’examen de calottes ciblées d’un rayon de autour des deux points extrêmes orbitaux (section A.9), la population masquée conserve des valeurs de densité de 252 par million de km² sur la face arrière contre 174 sur la face avant [A.5.3]. Cet écart donne un rapport de 1.45, représentant 5 écarts-types (), et reste fortement significatif. Ce résultat prouve que l’excédent de cratères qui survit au masque est concentré localement autour des coordonnées centrales du méridien est, et n’est pas réparti de façon homogène sur toute la bande [A.5.3].

Synthèse méthodologique :

L’application du masque de contrôle produit un double résultat : elle modifie et réduit une part considérable de l’amplitude longitudinale brute (et confirme le rôle des jeunes maria), mais n’élimine pas entièrement la non-uniformité résiduelle le long de la circonférence [A.5.3]. Toutefois, pour résoudre l’énigme polaire et confirmer le modèle d’arrondissement et de contraction mécanique de l’article, il suffit que le schéma longitudinal ait été mesuré, défini et séparé statistiquement ; l’explication physique de l’asymétrie entre l’est et l’ouest reste donc en dehors de l’argument central de la présente recherche [A.5.3].

A.5.4 Réponse du profil de latitude au masque de filtrage : stabilité de l’excédent polaire

Contrairement au profil de longitude, où le masque efface une part considérable de l’amplitude asymétrique [A.5.3], le profil de latitude répond à la procédure de filtrage d’une manière qui confirme la robustesse de l’effet polaire dans la croûte rocheuse ancienne [A.5.4]. Le comportement des bandes extrêmes reflète la mécanique de l’équation de densité résiduelle [A.3, A.5.4] :

A.5.5 Les enseignements scientifiques de la comparaison des pôles et du test d’égalité des surfaces

La convergence des valeurs entre le nord et le sud constitue l’un des principaux résultats méthodologiques du test de contrôle, et fournit une leçon ferme pour le traitement des données planétaires [A.5.5] :

A.5.6 Synthèse de la première phase des observations spatiales et limites de la géométrie des bandes

Les quatre examens par bandes (A.1–A.4) règlent entièrement la question de la validité du résultat morphologique et permettent de réduire le problème scientifique à quatre conclusions fermes et à une limite structurelle [A.5.6] :

Conclusions établies :

1. L’effet est global et bipolaire : la densité polaire élevée n’est pas une caractéristique locale du seul hémisphère nord ; elle est commune aux deux bandes polaires extrêmes après nettoyage des biais [A.5.6].

2. L’excédent ne dépend pas du masque : l’effet ne disparaît pas et ne s’estompe pas à la suite du retrait des grands bassins et des maria. Au contraire, au pôle Sud, le nettoyage de la surface a fait ressortir et augmenté la densité réelle de la croûte restante [A.5.6].

3. Isolement complet vis-à-vis de l’asymétrie longitudinale : la densité polaire n’est pas une expression indirecte de l’asymétrie entre l’est et l’ouest. Le profil de longitude présente une amplitude différente, un comportement géométrique distinct et une réponse entièrement opposée au masque (effondrement de la majeure partie de l’amplitude longitudinale au niveau de la bande contre conservation et renforcement de la densité aux extrémités polaires) [A.5.3, A.5.6].

4. Aucune exigence de profil monotone régulier : les données empiriques de terrain ne soutiennent pas l’affirmation d’une gradation régulière montant continûment de l’équateur vers les pôles. Le résultat scientifique précis – et physiquement plus solide – établit que les bandes adjacentes aux pôles se fixent nettement à l’extrémité supérieure absolue de la distribution globale en latitude [A.5.6].

La limite géométrique intrinsèque des bandes (The Ring Averaging Limitation) :

Malgré le règlement de la question initiale et la confirmation de la robustesse du résultat, l’usage de bandes de latitude complètes (75°-90°) se heurte à une limite analytique claire :

A.5.7 — Fin de la phase des bandes et passage aux calottes

À ce stade, le travail analytique que les bandes de latitude et de longitude peuvent accomplir dans le système est entièrement achevé [A.5.7].

Synthèse des résultats en une ligne :

On a établi l’existence d’une densité de cratères élevée et exceptionnelle aux deux extrémités polaires du profil de latitude, entièrement isolée du schéma spatial est-ouest, et qui reste stable et ferme même après nettoyage des effets des bassins et des maria [A.5.7].

Dans un détail méthodologique plus large, l’analyse des bandes a atteint quatre objectifs de contrôle principaux [A.5.7] :

1. Localisation spatiale : définition de la position exacte de l’anomalie de densité par rapport aux axes de latitude.

2. Séparation des schémas : isolement complet de l’effet polaire vis-à-vis de l’asymétrie orbitale le long des méridiens.

3. Réponse géologique : examen quantitatif de la variation des deux schémas sous la contrainte d’un masque géologique fixe et uniforme.

4. Équilibre bipolaire : mise en évidence du fait qu’après nettoyage des biais géologiques, on obtient deux bandes polaires (nord et sud) aux valeurs de densité semblables, élevées et équilibrées dans le système.

Changement de méthode analytique : le passage à la question géométrique

Malgré la robustesse des résultats, la mécanique des bandes se heurte à une limite structurelle et ne permet pas de répondre à la question scientifique suivante [A.5.7]. Jusqu’ici, les questions examinées étaient : les latitudes polaires se caractérisent-elles par une densité de cratères plus élevée, et ce résultat survit-il aux procédures de contrôle spatial [A.5.7] ?

Désormais, l’argument scientifique change de nature et devient un problème géométrique et morphologique strict [A.5.7] :

Comment la densité des cratères est-elle disposée et organisée à l’intérieur de l’espace polaire lui-même ?

Pour déchiffrer la structure interne de l’anomalie, la recherche modifie fondamentalement la forme de la mesure et du découpage géométrique de la croûte [A.5.7] :

Par conséquent, le chapitre suivant – la section A.6 (les calottes polaires) – ne constitue pas une répétition statistique ni un comptage supplémentaire de la même population [A.5.7]. Il représente l’étape suivante du développement de l’argument, découlant directement de la limite mécanique de l’examen par bandes : un passage net de la preuve de l’existence même d’une densité polaire élevée au déchiffrement de la géométrie spatiale et de la disposition interne des cratères à l’intérieur de la calotte polaire [A.5.7].

A.6 Les calottes polaires, sans masque

Dans cette section, la recherche modifie formellement la géométrie de la mesure et passe de l’usage de bandes circonférentielles à l’analyse de six calottes sphériques concentriques (emboîtées) autour de chaque pôle, de rayons angulaires : et [A.6]. L’objectif du test est de déchiffrer la structure interne de l’anomalie et d’examiner comment la densité cumulée varie à mesure que l’espace analytique converge vers l’intérieur, vers le point extrême du pôle de rotation [A.6].

Principes de la structure géométrique des calottes :

1. Définition du rayon angulaire () : une calotte de rayon englobe tous les cratères dont la distance d’arc à l’apex du pôle satisfait (par exemple : une calotte de s’étend de la latitude jusqu’à du pôle) [A.6].

2. L’équation de surface différentielle : la surface de chaque cellule sphérique est calculée exactement selon :

Pour un rayon , le facteur de surface croît de façon non linéaire : la calotte centrale restreinte () s’étend sur 0.288 million de km², tandis que la calotte large () s’étend sur 9.483 millions de km² (une surface plus de 30 fois plus grande) [A.6].

3. Structure cumulative (emboîtée) : les lignes de ce tableau ne représentent pas des zones annulaires isolées et indépendantes, mais un profil de sommation cumulative [A.6]. Le mécanisme ressemble aux anneaux d’une cible de fléchettes : la calotte de n’est pas seulement l’anneau extérieur, mais la cible entière depuis le centre () jusqu’à [A.6].

4. Conservation de la base de référence statique : le point de référence (valeur 1.00) reste parfaitement constant – la densité de la ceinture équatoriale brute complète sans masque (). Sa valeur est de 166.436 (pour ) et de 53.170 (pour ), ce qui permet une comparaison directe avec les données des parties A.1–A.4 [A.6].

5. Sensibilité de l’échantillon et facteur de taille : les calottes intérieures () couvrent de petites cellules de surface et contiennent un nombre restreint de cratères ; leurs intervalles de confiance sont donc larges et portent l’incertitude statistique la plus élevée (l’analyse complète est présentée en D.4) [A.6]. Le groupe de contrôle () représente une population de cratères plus ancienne ; un rapport plus élevé que celui du groupe des 20 km indique une accentuation de l’effet chez les cratères anciens (discuté en A.8.5) [A.6].

Synthèse des données de densité des cratères dans six calottes emboîtées (A.6)

Voici les valeurs numériques et relatives mesurées dans les calottes concentriques autour des deux pôles, sans application de masques géologiques :

Pôle géographiqueRayon angulaire de la calotte ()Surface de la calotte ()Cratères ()Densité par Rapport à la référence ()Cratères ()Densité par Rapport à la référence ()
Pôle Nord60°9.4831,7941891.14600631.19
Pôle Nord45°5.5551,1021981.19393711.33
Pôle Nord30°2.5415772271.36222871.64
Pôle Nord20°1.1443092701.621191041.96
Pôle Nord15°0.6461943001.80651011.89
Pôle Nord10°0.288862981.79311082.02
Pôle Sud10°0.288722501.5026901.70
Pôle Sud15°0.6461392151.2953821.54
Pôle Sud20°1.1442372071.2483731.36
Pôle Sud30°2.5415362111.27185731.37
Pôle Sud45°5.5551,2292211.33415751.41
Pôle Sud60°9.4832,0352151.29687721.36

Analyse morphologique du profil radial brut

La lecture des indicateurs et du comportement graphique des deux hémisphères fait ressortir les résultats suivants :

1. Prépondérance de l’excédent polaire global : aucune calotte polaire, quel que soit son diamètre et quel que soit l’hémisphère, ne présente de déficit par rapport à la ceinture équatoriale [A.6]. Le rapport minimal du système s’établit à 1.14 (la calotte nord la plus large, de ), ce qui prouve que l’intérieur de toutes les calottes se trouve en excédent de densité intrinsèque par rapport à l’équateur brut [A.6].

2. Le comportement différentiel du pôle Nord : le profil nord (la ligne bleu foncé) présente une convergence et une concentration mécanique classiques vers l’intérieur [A.6]. À mesure que le rayon de mesure se réduit, la densité monte continûment : de 189 dans la calotte large (), en passant par 227 dans la calotte de test moyenne (), jusqu’à un maximum de 300 cratères par million de km² dans la calotte de [A.6]. Dans la calotte centrale extrême (), le profil se stabilise sur une valeur presque identique (298), soit un rapport de 1.79 [A.6]. Dans le groupe de contrôle des grands cratères (), la tendance de concentration au nord s’accentue et atteint un rapport maximal de 2.02 dans la calotte de [A.6].

3. Le comportement plat du pôle Sud brut : le profil sud (la ligne rouge) présente à première vue une tendance entièrement différente [A.6]. Sur la majeure partie de la plage de mesure ( à ), la courbe de densité reste plate, horizontale et presque indépendante de la distance au pôle, avec des valeurs constantes comprises dans la plage étroite de 207–221 par million de km² [A.6]. Ce n’est que dans la calotte centrale la plus intérieure () qu’on enregistre une forte hausse vers une valeur de 250 (rapport de 1.50) [A.6].

Limite de l’interprétation géophysique à ce stade :

Les résultats empiriques bruts montrent en apparence un écart structurel profond entre les deux pôles : le nord se caractérise par une concentration intérieure nette, tandis que le sud présente une distribution plate et homogène [A.6]. Toutefois, la physique du modèle identifie que le profil sud est découpé et moyenné sur le fond du bassin géant South Pole–Aitken, qui a effacé une large part de l’enregistrement [A.6, A.11]. Pour déterminer quelle part de la variation provient de la déformation de la forme du corps et quelle part constitue la déformation d’un bassin d’impact isolé, il faut appliquer le masque de filtrage géologique [A.6].

Légende des couleurs : transition continue du jaune doux (« densité faible ») au rouge doux (« densité élevée »).

איור 5

Figure 5. Densité des cratères de diamètre égal ou supérieur à 20 km dans des calottes de rayon croissant autour des deux pôles (A.6), le catalogue entier. Le nord monte vers son pôle ; le sud est plat — le fond de South Pole–Aitken, dont la calotte sud de 30° est en majeure partie constituée (A.3, D.4). (Sur le graphique : axe horizontal — « rayon de la calotte autour du pôle », 10°, 15°, 20°, 30°, 45°, 60° ; axe vertical — « cratères ≥ 20 km par 106 km² », 150–300 ; la ligne bleu foncé — « pôle Nord » ; la ligne rouge — « pôle Sud ».)

A.7 Les calottes polaires après retrait des bassins et des maria

Cette section examine si le comportement différentiel entre les deux hémisphères observé à la section A.6 (profil nord concentré contre profil sud plat) est une propriété de la population originelle de cratères, ou s’il s’agit d’une déformation causée par le type de sol et par les fonds géologiques renouvelés [A.7]. Pour cela, le masque de filtrage strict (20 bassins de Neumann de diamètre km et les polygones de maria du LROC) est appliqué aux six calottes concentriques [A.7]. Cette procédure garantit une comparaison précise et à l’identique sur une surface ancienne et uniforme, en dehors des sites d’effacement thermique et mécanique [A.7].

Indicateurs et calibrage méthodologique :

1. Calcul de la surface et du numérateur résiduels : les cratères dont le point central tombe dans les limites du masque sont omis, et la surface du masque est déduite exactement du dénominateur [A.7]. La grille de mesure repose sur des cellules sphériques exactes d’une résolution de , les chevauchements n’étant retirés qu’une seule fois [A.7]. Là où les géométries de découpage se recouvrent (les calottes de contre les bandes de latitude en A.3), on obtient une cohérence et une fermeture mathématique complètes des données de surface et des numérateurs [A.7].

2. Fixation de la règle de référence fixe : le point de référence (valeur 1.00) reste inchangé – la densité de la ceinture équatoriale brute complète sans masque ( pour 20 km, et pour 40 km) [A.7]. Ce choix est délibéré et sert de règle statique fixe pour suivre quantitativement le degré de variation géologique de chaque calotte, et non d’indicateur de l’excédent polaire final [A.7].

3. Indicateur de variation de densité () : compare chaque calotte concentrique à elle-même dans son état brut en A.6 :

Une valeur positive indique une hausse de la densité de la surface restante à la suite du retrait d’unités pauvres en cratères [A.7].

Synthèse des données de densité des cratères masquée dans six calottes emboîtées (A.7)

Voici les indicateurs numériques et les valeurs relatives obtenus dans les calottes concentriques masquées :

Pôle géographiqueRayon angulaire de la calotte ()Surface restante ()Densité Rapport Variation ()Densité Rapport Variation ()
Pôle Nord60°6.4481,5152351.41+24.2%507791.48+24.3%
Pôle Nord45°4.0379352321.39+16.7%336831.57+17.6%
Pôle Nord30°2.3995602331.40+2.8%213891.67+1.6%
Pôle Nord20°1.1443092701.620.0%1191041.960.0%
Pôle Nord15°0.6461943001.800.0%651011.890.0%
Pôle Nord10°0.288862981.790.0%311082.020.0%
Pôle Sud10°0.104363482.09+39.1%141352.54+49.8%
Pôle Sud15°0.275843051.83+42.0%341242.32+50.7%
Pôle Sud20°0.5321492801.68+35.2%52981.84+34.7%
Pôle Sud30°1.3083222461.48+16.7%105801.51+10.3%
Pôle Sud45°2.8527022461.48+11.3%214751.41+0.5%
Pôle Sud60°4.9181,1782401.44+11.6%364741.39+2.2%

Analyse géodynamique du profil radial masqué

L’examen des indicateurs résiduels et du comportement graphique commun (Figure 6) fait ressortir les enseignements suivants :

איור 6

Figure 6. Les mêmes six calottes autour des deux pôles après retrait des bassins et des maria (A.7), cratères de 20 km et plus. La ligne nord ne change pas dans la plage de 20° ; la ligne sud monte dans ses calottes intérieures, et les deux se rapprochent. (Sur le graphique : axe horizontal — « rayon de la calotte autour du pôle », 10°, 15°, 20°, 30°, 45°, 60° ; axe vertical — « cratères ≥ 20 km par 106 km² », 150–350 ; la ligne bleu foncé — « pôle Nord » ; la ligne rouge — « pôle Sud ».)

1. Effet géologique nul dans le cœur nord : les trois calottes intérieures au nord () présentent des valeurs de variation constantes de [A.7]. Cette donnée prouve mathématiquement qu’il n’y a ni bassins de Neumann ni maria basaltiques dans un rayon de autour de l’apex du pôle Nord ; par conséquent, le profil central élevé qui y est mesuré (300 par million de km²) reflète précisément l’état vierge de la croûte [A.7].

2. L’effet extrême dans le cœur sud et l’effacement par le bassin SPA : les plus petites calottes au sud présentent l’effet de l’effacement mécanique extrême du bassin SPA [A.7]. Dans la calotte sud de , la surface de mesure se réduit de plus de moitié (chute de 0.288 à 0.104 million de km²), ce qui indique que la majeure partie du cœur sud intérieur était couverte par le fond du bassin [A.7]. Le nettoyage de ce fond fait bondir la densité réelle de la croûte montagneuse résiduelle de dans le groupe des 20 km (valeur de 348, rapport de 2.09) et de dans le groupe des 40 km (valeur de 135, rapport de 2.54) [A.7].

3. Convergence morphologique vers une courbe de concentration uniforme : la comparaison de la Figure 6 avec la Figure 5 révèle un renversement morphologique [A.7]. Alors qu’à l’état brut la ligne sud était entièrement plate, l’application du masque de contrôle a égalisé le type de croûte et relevé les calottes intérieures du sud [A.7]. Désormais, la courbe sud (rouge) change entièrement de forme et monte en pente raide vers l’intérieur en direction du pôle, comme la courbe nord (bleu foncé) [A.7]. Dès que les deux pôles sont comptés sur le même type de surface, tous deux montrent une seule corrélation ferme : des cratères de plus en plus denses à mesure qu’on s’approche de l’apex du pôle, résultat qui confirme le modèle géodynamique de contraction de surface [A.7].

Précision méthodologique sur les valeurs d’excédent propres :

La colonne des rapports de ce tableau compare un numérateur masqué à un dénominateur brut (qui conserve les maria vides de l’équateur), et sert donc uniquement de règle de suivi [A.7]. Pour dériver l’excédent polaire propre réel, l’article applique une comparaison à l’identique (retrait simultané des maria et du bassin SPA des deux côtés) sur les calottes standard de contre la ceinture équatoriale () [A.7] :

L’écart-type () reflète la fluctuation aléatoire du comptage ; l’intervalle de confiance de est deux fois plus large, et établit avec une pleine confiance statistique que la valeur réelle du rapport se situe dans la plage de 1.10 à 1.30 – valeur nettement supérieure à la base de distribution uniforme (1.00), sur laquelle repose la conclusion scientifique de la recherche [A.7].

A.8 Synthèse et conclusions de l’examen des calottes

Les sections A.6 et A.7 ont examiné les régions polaires selon une géométrie concentrique (emboîtée), fondamentalement différente du moyennage spatial des examens par bandes [A.8]. Au lieu de mesurer la densité moyenne dans des anneaux de latitude globaux faisant le tour du corps sur 360°, chaque pôle de rotation a été défini comme le centre structurel de six calottes sphériques concentriques de rayons angulaires et [A.5.6, A.8]. Le test a été effectué dans un premier temps sur le catalogue Robbins brut sans masques, puis sous la contrainte stricte d’un masque géologique (retrait des 20 grands bassins de Neumann et des polygones de maria du LROC) [A.8].

Le croisement des deux procédures sépare méthodologiquement deux questions centrales [A.8] :

1. Quelle est la géométrie spatiale observée de la densité des cratères à mesure que le vecteur de distance converge vers le point extrême du pôle de rotation [A.8] ?

2. Quelle part de cette géométrie radiale est inscrite dans la croûte vierge ancienne, et quelle part constitue une déformation résiduelle causée par une activité endogène tardive (bassins et maria) [A.8] ?

Le résultat principal établit que l’application du masque de contrôle modifie fondamentalement la morphologie statistique et la manière dont il faut définir physiquement la différence structurelle entre les deux pôles [A.8].

A.8.1 L’état brut (sans masque) : asymétrie de morphologie interne

Dans le comptage empirique initial sans masque (A.6), les deux hémisphères présentent deux profils radiaux différents et clairement distincts [A.8.1] :

Si l’analyse scientifique s’était arrêtée à ce stade (A.6), la seule conclusion physique aurait été que les deux pôles diffèrent fondamentalement à la fois par l’intensité de leur densité absolue et par leur géométrie structurelle interne (nord concentré contre sud homogène et plat) [A.8.1]. Les données de la section A.7 prouvent que cette description est incomplète et résulte d’un biais géologique sévère [A.8.1].

A.8.2 Comportement du profil nord sous la contrainte du masque

L’un des résultats les plus clairs et les plus importants de l’application du masque de contrôle géologique (A.7) est la robustesse absolue du cœur structurel du pôle Nord [A.8.2] :

Le masque modifie donc l’arrière-plan spatial qui entoure la région polaire nord, mais laisse le cœur intact [A.8.2]. Avant le masque, le contraste structurel entre le cœur () et l’environnement () était extrême (300 contre 189) [A.8.2]. Après le masque, alors que le cœur reste stable à 300, l’environnement de est monté mathématiquement à 235 en raison du nettoyage des maria qui s’y trouvent [A.8.2]. Par conséquent, la concentration vers l’intérieur survit fermement, mais son amplitude (Contrast) se réduit [A.8.2]. Le masque prouve qu’une partie du contraste brut d’origine provenait du biais de distribution des surfaces jeunes hors du cœur, et non d’une asymétrie d’impact [A.8.2].

A.8.3 Le renversement géodynamique du profil sud

Contrairement au nord, où le cœur intérieur était insensible au masque, dans l’hémisphère sud le masque de filtrage modifie fondamentalement l’image morphologique, et ce dès les calottes les plus restreintes [A.8.3] :

Le mécanisme géophysique :

Dans un rayon de autour de l’apex du pôle Sud, la seule composante structurelle du masque est le bassin géant South Pole–Aitken (SPA) [A.8.3]. Le fond de ce bassin se caractérise par une pauvreté extrême en cratères en raison de l’effacement mécanique complet de l’enregistrement précoce [A.8.3]. Le retrait des surfaces du bassin soustrait mathématiquement une part énorme de surface (le dénominateur) et très peu de cratères (le numérateur), et révèle donc la densité élevée réelle de la croûte montagneuse vierge résiduelle qui l’entoure [A.8.3]. Ce résultat ferme exclut entièrement la possibilité de lire le profil sud plat d’origine comme l’enregistrement d’une population de cratères vierge ; il s’agit d’une déformation locale due à un bassin d’impact isolé, neutralisée avec succès [A.8.3].

A.8.4 Convergence bipolaire et égalité statistique après le masque

Dès que les deux pôles sont mesurés et calculés sur le même type de surface crustale (Highlands vierges hors des sites de remise à zéro thermique et mécanique), la gradation morphologique change et présente une image uniforme et cohérente [A.8.4] :

Rayon angulaire de la calotteDensité nord masquée ()Densité sud masquée ()Rapport empirique (sud/nord)Intervalle de confiance normalisé ()Signification statistique et structurelle
Calotte de 1.020.78 – 1.32Égalité mécanique absolue (l’intervalle contient la valeur de symétrie 1.00) [A.8.4].
Calotte de ——Proximité structurelle étroite des valeurs de densité crustale [A.8.4].
Calotte de —0.77 – 1.74Le sud est plus élevé que le nord, mais sur la base d’un échantillon restreint de seulement 36 cratères (incertitude large) [A.8.4].

Conclusion méthodologique finale :

L’analyse statistique complète des intervalles de confiance () détaillée au chapitre D.4 établit que, comme la valeur de symétrie parfaite (1.00) est contenue dans les intervalles d’incertitude des échantillons pour toutes les calottes et pour les deux seuils de diamètre, les deux pôles sont définis comme égaux et parfaitement équilibrés dans les limites de l’erreur de comptage aléatoire [A.8.4]. Par conséquent, la formulation scientifique précise n’affirme pas que « le masque rend les pôles entièrement identiques », mais que : le masque de filtrage géologique neutralise et retire la majeure partie du contraste brut apparent qui existait entre les profils intérieurs du nord et du sud, et les amène à converger sur des plages de densité crustale semblable, élevée et parfaitement équilibrée [A.8.4].

A.8.5 Comparaison des seuils de diamètre linéaires et facteur de dimension du cratère (horloge chronologique approximative)

L’usage différentiel de deux seuils de diamètre stricts ( et ) constitue un outil de contrôle méthodologique essentiel, qui vérifie que les tendances géodynamiques observées ne sont pas le produit d’une anomalie statistique propre à une population de taille spécifique [A.8.5] :

Cette cohérence dans les deux hémisphères prouve que la gradation morphologique générale ne dépend pas du choix du seuil de diamètre [A.8.5]. Toutefois, on enregistre un écart constant d’une signification physique profonde : les rapports à la ceinture équatoriale sont systématiquement plus élevés dans le groupe que dans le groupe [A.8.5]. Dans la calotte nord brute de , le rapport s’établit à 1.64 (à 40 km) contre 1.36 (à 20 km) ; dans la calotte sud masquée de , il s’accentue à 2.32 contre 1.83 [A.6, A.7, A.8.5].

Le mécanisme de l’horloge chronologique du diamètre des cratères :

Dans le cadre du modèle physique, le diamètre d’un cratère fonctionne comme une horloge chronologique approximative – les cratères géants () sont statistiquement attribués aux phases les plus précoces de l’enregistrement des impacts, tandis que les cratères plus petits se sont accumulés tout au long de la séquence ultérieure [A.8.5].

Plus un cratère s’est formé tôt (près du pic d’aplatissement géométrique de la jeune Lune), plus il a subi une part étendue de la séquence d’arrondissement et de diminution de volume qui a suivi, et il a donc subi le maximum de contraction mécanique de surface convergeant vers les pôles [A.8.5].

A.8.6 Intégration structurelle : croisement des données des bandes et des calottes

La combinaison des deux niveaux de découpage géométrique de la Partie A permet de constituer une image physique cohérente et isolée de la croûte du corps [A.8.6] :

Le croisement révèle une distinction critique entre deux types d’information : l’excédent polaire aux hautes latitudes est une propriété physique stable qui survit au masque, tandis que la variation géométrique initiale entre le nord et le sud est une déformation locale passagère qui ne survit pas au contrôle [A.8.6]. Cette distinction évite une erreur méthodologique grave, et permet de séparer l’état actuel de la croûte conservée aujourd’hui des propriétés initiales du phénomène géodynamique [A.8.6].

A.8.7 Synthèse des conclusions empiriques de la première phase

Après le traitement des sections A.1–A.7, on établit cinq conclusions observationnelles fermes, qui délimitent le phénomène que le mécanisme physique doit expliquer [A.8.7] :

1. Non-uniformité latitudinale significative : la population globale de cratères n’est pas distribuée de façon homogène selon la latitude ; les bandes polaires se fixent absolument à l’extrémité supérieure du profil de densité [A.8.7].

2. Robustesse endogène : l’excédent de densité polaire n’est pas fonction de la présence des 20 grands bassins ou des surfaces de maria, et il survit entièrement à l’application du masque géologique [A.8.7].

3. Symétrie bipolaire résiduelle : le nettoyage des biais géologiques révèle des valeurs de densité semblables et équilibrées entre les deux pôles, tant au niveau de la moyenne globale des bandes qu’à l’intérieur de la structure interne des calottes intérieures () [A.8.7].

4. Instabilité de l’asymétrie nord-sud : l’écart morphologique entre la courbe nord concentrée et la courbe sud plate s’effondre sous le contrôle, et se révèle être le produit d’un impact de bassin local au sud et non une propriété initiale [A.8.3, A.8.7].

5. Cohérence multi-seuils : la tendance générale est reproduite et visible simultanément au seuil de comme au seuil de , et ne repose pas sur le choix d’un diamètre unique [A.8.7].

A.8.8 Limites du test de contrôle et frontières de la connaissance au stade actuel

Afin de préserver la précision scientifique, la recherche définit trois limites analytiques strictes à l’interprétation des données à ce stade [A.8.8] :

A.8.9 Fin de l’examen des calottes et définition du phénomène scientifique

Cette section résume et clôt formellement la question morphologique qui a conduit de l’analyse par bandes à la phase des calottes emboîtées, et constitue la base de données vierge de la croûte résiduelle [A.8.9].

Synthèse de la dynamique analytique :

Effondrement de l’asymétrie et convergence des données résiduelles

Le nettoyage des enregistrements géologiques a conduit à une réduction spectaculaire des écarts bipolaires et à une convergence structurelle des valeurs de densité à l’intérieur des cœurs intérieurs [A.8.9] :

Rayon angulaire de la calotteRapport brut d’origine (nord contre sud)Densité nord masquéeDensité sud masquéeÉtat résiduel après filtrage
Calotte de 300 contre 215 par million de km²Convergence complète [A.8.9]
Calotte de 270 contre 207 par million de km²Convergence complète [A.8.9]
Calotte de —Le sud dépasse le nord (échantillon de 36 cratères) [A.8.9]

Le résultat restreint et précis établit que la majeure partie de l’asymétrie géométrique observée initialement entre les calottes intérieures n’est pas stable et s’estompe dans des conditions de contrôle uniformes [A.8.9]. Ce résultat complète et confirme les résultats des bandes : au niveau des bandes, le masque a égalisé les secteurs à des valeurs de 300 et 305 ; au niveau des calottes, la même convergence ferme est reproduite et ne disparaît pas même lorsque l’espace analytique se réduit progressivement vers l’apex polaire [A.3, A.7, A.8.9]. Les deux géométries de découpage différentes indiquent précisément le même vecteur physique : l’excédent de densité élevé n’est pas une caractéristique locale d’un seul pôle [A.8.9].

Définition des limites du contrôle mécanique du bassin SPA

La recherche maintient une séparation scientifique stricte entre la description empirique et le mécanisme mécanique [A.8.9] :

1. Identification des composantes du masque : le masque de filtrage est défini comme l’union des polygones de maria et des 20 grands bassins de Neumann [A.7, A.8.9].

2. Le rôle géodynamique du bassin SPA : le bassin Pôle Sud–Aitken (South Pole–Aitken) constitue la composante géante principale qui pénètre profondément dans l’espace polaire sud, et sa géométrie morphologique fournit la justification physique du fait que le masque agit et fait bondir les valeurs au sud de manière exceptionnelle par rapport au nord vierge [A.8.3, A.8.9].

3. Isolement des composantes étape par étape : le calcul masqué total ne constitue pas un test du bassin SPA seul [A.8.9]. Par conséquent, il est méthodologiquement interdit d’attribuer tout le renversement morphologique au sud au seul SPA, sans isoler sa contribution résiduelle de celle des autres composantes du masque [A.8.9]. Cette séparation est effectuée et confirmée sur deux plans stricts :

La base de données consolidée : définition de l’anomalie crustale

À l’issue de l’examen des calottes, on définit le phénomène physique vierge que les autres parties de l’article doivent expliquer au moyen des équations de la géodynamique et du modèle de contraction mécanique [A.8.9] :

Sur la base de ces quatre indicateurs centraux, la recherche passe maintenant à la phase d’intégration géométrique et de test physique des axes du corps [A.8.9].

A.9 Les six axes et la définition de la géométrie triaxiale de la jeune Lune

Cette section élargit les limites de l’analyse et pose un test géométrique combiné. Alors que les sections A.1–A.8 examinaient la surface selon une seule direction de comptage à la fois (latitude ou longitude séparément), cette section combine les données sur la base des trois axes naturels d’un corps verrouillé par marée [A.9] :

1. L’axe de rotation (polaire) : il passe par les deux pôles. Cet axe gouverne tous les phénomènes influencés par les contraintes centrifuges, par l’aplatissement rotationnel et par son relâchement géodynamique (le mécanisme principal proposé pour la réduction et le resserrement de la surface aux pôles) [A.9].

2. La ligne de la Terre (radiale) : elle va du point subterrestre (la face proche, tournée en permanence vers la Terre) au point anti-terrestre (le centre de la face cachée). Cet axe gouverne les contraintes de marée extérieures et l’interaction mécanique à deux corps (comme la zone de contact ancienne avec la poire du corps parent) [A.9].

3. La direction de l’orbite (tangentielle) : elle va du point avant (l’apex du mouvement dans l’espace, ) au point arrière (l’arrière orbital, ). Cet axe reflète des asymétries dynamiques découlant du mouvement du corps à travers le vol des projections ou des matériaux d’éjection [A.9].

La méthodologie repose sur la comparaison de la densité à l’intérieur de six calottes de surface égale, d’un rayon angulaire de (rayon d’environ 910 km à la surface), autour des six extrémités de ces axes, après application du masque de filtrage géologique strict (retrait des 20 bassins géants de Neumann et des polygones de maria du LROC) [A.9]. L’égalité intrinsèque des surfaces permet une comparaison quantitative directe et exempte de biais de dimension [A.9].

Synthèse des données de densité des cratères dans les six calottes principales (A.9)

Voici les indicateurs numériques exacts, classés strictement selon la gradation de la densité résiduelle – de la calotte la moins dense à la plus dense du système [A.9] :

Rang dans la gradationCalotte géométrique examinéeCoordonnées du centreSurface restante ()Densité Rapport à la référence ()Variation de densité ()Densité Rapport à la référence ()Variation de densité ()
1Point subterrestre1.0421541480.89+57.8%50480.90+50.6%
2Point avant1.7653081741.05+34.7%107611.14+29.4%
3Point anti-terrestre1.6683682211.33+2.3%105631.18-8.6%
4Pôle Nord2.3995602331.40+2.8%213891.67+1.6%
5Pôle Sud1.3083222461.48+16.7%105801.51+10.3%
6Point arrière1.4683702521.51-0.3%135921.73+9.2%
איור 7

Figure 7. Les six axes (A.9) : densité des cratères de 20 km et plus dans des calottes de 30° autour des points subterrestre, avant, pôle Sud, anti-terrestre, pôle Nord et arrière. Les pôles se situent au niveau de la face cachée ; le point subterrestre est le moins dense. Les barres représentent le catalogue complet, avant le masque ; lues de gauche à droite, elles montent du point tourné vers la Terre vers la face arrière. (Sur le graphique : axe vertical — « cratères ≥ 20 km par million de km² » ; les barres de gauche à droite — subterrestre 94, avant 129, pôle Sud 211, anti-terrestre 216, pôle Nord 227, arrière 253 ; légende inférieure — « calotte de 30° autour de chacun des six points principaux, le catalogue complet ».)

Analyse dynamique de l’effet du masque et confirmation du minimum subterrestre

La lecture de la colonne de variation de densité () et de la gradation statistique fournit une réponse nette à une objection méthodologique attendue [A.9] :

Lecture de la gradation empirique comme enregistrement d’un corps triaxial aplati

La gradation numérique conservée avant comme après le masque (subterrestre avant anti-terrestre / pôles arrière) prouve que la distribution des cratères n’est pas gouvernée par un axe unique, mais constitue le reflet d’une géométrie triaxiale d’un corps déformé [A.9]. Le modèle géophysique interprète le schéma spatial de la manière suivante [A.9] :

1. Équilibre de l’axe de rotation (pôles contre équateur) : les pôles sont fixés dans un excédent de densité prononcé au-dessus des points équatoriaux (subterrestre et avant), comme l’impose le modèle d’arrondissement et de contraction de l’aplatissement rotationnel [A.9]. Après le masque, la moyenne des deux pôles est 1.12 fois supérieure à la moyenne des deux faces équatoriales (avant et arrière) [A.9].

2. Asymétrie de l’axe de la ligne de la Terre (tête libre contre disque de contact) : les deux extrémités de l’axe de la ligne de la Terre présentent un comportement entièrement différent [A.9]. L’extrémité anti-terrestre (la face cachée) se trouve en excédent de densité élevé () et rejoint entièrement le niveau des pôles [A.9]. À l’inverse, l’extrémité subterrestre (la face proche) constitue la zone la moins dense du système () [A.9].

La combinaison de ces contraintes dessine un corps soumis à un aplatissement rotationnel symétrique autour de l’axe polaire (), et en même temps étiré à l’extrême le long de l’axe radial () entre la gravité terrestre vers l’intérieur et l’entraînement centrifuge vers l’extérieur [A.9]. Cette structure définit une forme triaxiale précise qui satisfait la gradation physique : un axe long en direction de la Terre (), un axe moyen le long du vecteur orbital (), et un axe court passant par les pôles () [A.9]. La forme théorique d’un corps de marée libre (C.2) restitue précisément une gradation de rapports identique (rapports de 1.66, 1.56 et 1.06 pour les pôles contre le subterrestre, les faces contre le subterrestre, et les pôles contre les faces, respectivement) [A.9].

Détermination des indicateurs de forme géométrique à partir de l’enregistrement des cratères

Au lieu de postuler une forme hypothétique, les équations d’arrondissement et de diminution de volume (C.3) permettent de lire les dimensions originelles du corps directement à partir de l’enregistrement empirique de l’excédent de cratères [A.9] :

A.9.1 L’asymétrie orbitale : la face arrière contre la face avant

Les deux extrémités de l’axe de la direction orbitale présentent un écart statistique cohérent et prononcé dans tous les indicateurs d’analyse spatiale appliqués dans la Partie A, la direction de l’anomalie restant ferme sans inversion de signe, tandis que l’amplitude quantitative varie selon la géométrie de découpage [A.9.1] :

Indicateur de cratères examiné (≥ 20 km par 10⁶ km²)La face arrière (90° E)La face avant (90° W)Rapport empirique (arrière/avant)Significativité et structure spatiale
Calotte de 30° (le catalogue brut complet)2531291.96Excédent extrême (presque un facteur deux) [A.9.1].
Calotte de 30° (après application du masque)2521741.45Survit avec une force d’environ 5σ (pleinement significatif) [A.5.3, A.9.1].
Bande de 30° (le catalogue brut complet)2371541.54Asymétrie significative au niveau de la bande brute [A.9.1].
Bande de 30° (après application du masque)2192041.07S’atténue au niveau de 1σ (dans le domaine du hasard) [A.5.3, A.9.1].

A.9.2 La corrélation géologique : distribution des surfaces de maria et de basalte

L’isolement de l’anomalie locale autour de 90° est soulève la question de la variation morphologique : qu’existe-t-il dans cette cellule de surface qui n’existe pas sur la face opposée [A.9.2] ? Le découpage précis des polygones de maria du LROC selon la géométrie des calottes de 30° donne les indicateurs quantitatifs suivants [A.9.2] :

Indicateur géologique à l’intérieur de la calotte de 30°Autour de 90° est (arrière)Autour de 90° ouest (avant)Rapport empirique
Nombre d’unités de mare isolées (polygones)132632.10
Surface totale de maria (km²)205,830279,4700.74
Dont : surfaces du bassin Oceanus Procellarum—192,134—
Surface de maria propre à la face (sans Procellarum)205,83087,3362.36

A.9.3 L’équation du bilan morphologique : les maria comme ouvertures d’éjection et d’alimentation

Le paradoxe scientifique suivant exige une solution : pourquoi le secteur le plus riche en cratères est-il celui qui abrite l’excédent de maria, alors que les maria eux-mêmes sont des fonds jeunes dépourvus de cratères [A.9.3] ? La réponse repose sur la mécanique du modèle LEF [A.9.3] :

1. Mécanisme d’éjection et de retour (Article 6) : dans le cadre du modèle, chaque mare constitue une ouverture crustale au-dessus d’un dôme d’alimentation magmatique actif [A.9.3]. Ces ouvertures ont éjecté et projeté vers le haut d’énormes quantités de matériau [A.9.3]. Le matériau projeté a conservé son énergie cinétique et son moment cinétique, et a effectué des vols balistiques qui sont retombés sur la croûte à des vitesses de 2 à 5 km/s (exactement la vitesse à laquelle il avait quitté l’ouverture) [A.9.3]. Ces retombées se sont produites dans un rayon de 50 à 1,000 km autour de l’ouverture source lors des phases précoces, et se sont réduites plus tard avec l’arrondissement et le resserrement du corps [A.9.3].

2. L’effet bidirectionnel d’une ouverture d’éjection : la mare agit sur l’enregistrement crustal selon deux vecteurs opposés et simultanés : elle produit des cratères d’impact autour d’elle au moyen du matériau de retour, tandis que son propre fond reste chaud, fond et remet à zéro les cratères qui s’y sont imprimés (le fond ne s’est figé et fixé qu’après le détachement) [A.9.3]. L’enregistrement géométrique sur le terrain s’exprime précisément au moyen de l’équation du bilan morphologique :

En dehors des limites des maria, le terme d’effacement thermique s’annule tandis que le terme du matériau de retour est maximal ; à l’intérieur, le terme d’effacement domine et annule le numérateur [A.9.3].

Asymétrie de l’action du masque :

L’équation explique pourquoi l’application du masque (A.3) modifie les valeurs de manière non symétrique [A.9.3] :

A.9.4 Intégration des mécanismes de densification de la face arrière

L’excédent de cratères conservé sur la face arrière s’explique par la combinaison de deux mécanismes physiques autonomes agissant dans la même direction dynamique [A.9.4] :

1. Mécanisme de dispersion locale du matériau de retour :

L’échelle des vols balistiques du matériau (jusqu’à environ 1,000 km) correspond précisément aux dimensions géométriques des calottes de 30° (≈ 910 km), et explique pourquoi l’effet est conservé dans la calotte et s’effondre dans les bandes longues (qui moyennent la surface avec un sol vierge que les projections n’ont pas atteint) [A.9.1, A.9.4]. Le lien mécanique est défini par la chaîne logique stricte suivante [A.9.4] :

Concentration accrue d’ouvertures (maria) dans la structure → volume accru de matériau projeté et retombé → densité de cratères exceptionnelle dans la croûte qui les entoure

Le calcul quantitatif complet de l’histoire des éjections, du rendement en cratères et des géométries de vol est développé et détaillé dans la Partie E.3 de l’article et dans la Partie K de l’Article 6 [A.9.4].

2. Mécanisme de retard de rotation dynamique (Lag Effect) :

Pendant la phase attachée, la Lune est restée verrouillée en rotation et a tourné simultanément avec la Terre [A.9.4]. Le matériau projeté avec force vers le haut à partir d’ouvertures d’éjection situées dans d’autres régions du corps a subi un retard physique par rapport à la vitesse de rotation rapide du système (période de 5.5 heures) [A.9.4]. La dynamique fonctionne comme le lancer d’une balle verticalement depuis le plancher d’un manège en rotation : le corps avance, et le matériau qui retombe atterrit, figurativement, derrière son point de départ [A.9.4]. Par conséquent, l’arrière orbital (la face arrière) a fonctionné comme un piège dynamique de quantité de mouvement, qui a recueilli et absorbé des quantités résiduelles de matériau éjecté de toutes les parties du corps et retombé à l’arrière, en plus du matériau de retour de ses propres ouvertures [A.9.4].

A.10 Le bassin Pôle Sud–Aitken (SPA) comme objet spatial autonome

Au cours du traitement des données bipolaires, on a établi que le bassin South Pole–Aitken ne constitue pas seulement une composante de filtrage passive ; il possède une structure interne particulière et une non-uniformité significative de la distribution de la population en son sein [A.10]. Le bassin s’étend sur une largeur d’environ 2,400 km sur la face cachée et pénètre profondément dans la calotte polaire sud (il couvre 24% de la surface de l’ensemble des calottes polaires) [A.10].

La recherche rejette l’usage d’une moyenne globale unique pour le bassin, car l’application d’un seul nombre monolithique produirait une déformation méthodologique grave : une moyenne combinée d’un centre entièrement vide et de marges montagneuses ordinaires présenterait l’illusion crustale d’un bassin modérément et uniformément pauvre en tout point [A.10]. Par conséquent, afin de préserver une rigueur académique stricte, la section A.10 isole la structure du bassin et effectue une série de mesures radiales continues, déployées par paliers depuis son point central géométrique vers l’extérieur jusqu’aux marges, dans le but de clore les questions ouvertes de l’enregistrement des cratères de l’Article 7 [A.10].

A.10.1 Méthodologie de l’analyse radiale du bassin South Pole–Aitken (SPA)

Afin de déchiffrer la non-uniformité structurelle interne du bassin South Pole–Aitken, les coordonnées du centre géométrique du bassin ont été fixées selon les mêmes paramètres que ceux utilisés pour construire le masque de filtrage : sud, est, avec une définition stricte du diamètre de bordure extérieure de 2,400 km [A.10.1].

Le facteur neutralise la courbure de la sphère, et impose qu’un anneau de largeur éloigné du centre soit plus grand en surface qu’un anneau de même largeur proche de celui-ci [A.10.1].

Synthèse des données de densité des cratères dans les anneaux autour du centre de SPA (A.10)

Voici les indicateurs numériques exacts obtenus dans les six anneaux indépendants, calculés pour le groupe et le groupe [A.10.1] :

Secteur de distance au centre de SPAPosition géologique dans la structureSurface de l’anneau ()Densité Rapport local ()Densité Rapport local ()
0 – 10°À l’intérieur de SPA (le cœur central)0.288381320.5915520.69
10 – 20°À l’intérieur de SPA0.8561651930.8758680.90
20 – 30°À l’intérieur de SPA1.3972792000.9098700.93
30 – 39.57°À l’intérieur de SPA, jusqu’à la limite de la bordure1.8063642020.91129710.95
Bordure à 39.57°· · Limite mécanique de la bordure · ·———————
39.57 – 50°Immédiatement à l’extérieur de la bordure du bassin2.4285402221.00183751.00
50 – 60°En dehors des limites du système SPA2.7085752120.95196720.96

A.10.2 La zone centrale comme foyer du minimum (structure en marche)

L’analyse des colonnes de densité et du profil spatial (Figure 8) fait ressortir les résultats scientifiques principaux suivants :

Légende des couleurs : du jaune (« faible ») au rouge (« élevé »).

איור 8

Figure 8. Densité des cratères de 20 km et plus dans six anneaux séparés autour du centre de South Pole–Aitken (A.10). La zone centrale est le minimum du profil ; la densité monte dès l’anneau 10–20°, et se rapproche de celle de l’environnement à l’approche de la bordure. La partie ombrée est l’intérieur du bassin ; la ligne pointillée est sa bordure. (Sur le graphique : axe vertical — « cratères ≥ 20 km par million de km² » ; axe horizontal — « distance au centre de SPA (points médians des anneaux) », avec des repères à 0°, 10°, 20°, 30°, 50°, 60° et « bordure » ; légendes — « à l’intérieur de SPA », « à l’extérieur de SPA » ; sur la courbe rouge — « 132 — le minimum » et 222.)

A.10.3 Calibrage du groupe de contrôle () et intervalles d’incertitude

Le groupe des grands cratères () reproduit exactement la tendance morphologique et confirme la cohérence de la direction du résultat [A.10.3] :

A.10.4 Conséquences dynamiques pour la procédure de comptage propre (A.11)

Les résultats empiriques prouvent que définir le bassin SPA comme un bloc unique sous le terme monolithique de « déficit » nuit à la précision scientifique, car la pauvreté en cratères est limitée et physiquement concentrée à l’intérieur de la seule zone centrale [A.10.4]. Les trois anneaux extérieurs de l’intérieur du bassin ( jusqu’à la bordure) présentent des valeurs de densité nettement supérieures à celles du centre (193–202), qui se rapprochent du niveau de l’arrière-plan montagneux ordinaire de l’environnement [A.10.2, A.10.4].

Cet enseignement structurel est essentiel pour comprendre la procédure de comptage propre des chapitres A.11 et D.1 [A.10.4] :

A.10.5 Interface de la recherche : le point de contact avec l’Article 8

Le résultat morphologique concernant l’existence d’une zone centrale pauvre en cratères à l’intérieur du bassin South Pole–Aitken est apparu de manière indépendante lors du test de contrôle des cratères des pôles dans le présent article (Article 7) [A.10.5]. Par conséquent, l’identification de l’anomalie et sa quantification numérique font partie intégrante de la documentation observationnelle empirique de cette recherche [A.10.5].

Toutefois, l’établissement de l’existence de ce minimum central pose une nouvelle question physique et géodynamique au système : pourquoi le cœur géographique intérieur du bassin SPA porte-t-il une densité de cratères aussi exceptionnellement basse par rapport aux anneaux périphériques et aux surfaces des marges extérieures du bassin [A.10.5] ?

La répartition structurelle des rôles entre les parties de la série scientifique

Afin de préserver une structure académique sobre et retenue, une répartition des rôles définie et stricte a été établie entre les parties de la série, qui délimite le champ de l’Article 7 [A.10.5] :

Précision méthodologique stricte

Toutes les valeurs et données numériques présentées dans ce chapitre ont été dérivées et traitées de manière cohérente à partir du catalogue Robbins (2019), qui sert de base de données fixe tout au long de la Partie A [A.10.5]. On souligne fermement que la définition de la limite de du cœur central constitue uniquement un élément de découpage analytique destiné à définir l’anneau intérieur, et n’est ni présentée ni interprétée comme une limite physique ou géologique démontrée de la zone de déficit en cratères [A.10.5].

A.11 Méthodologie du comptage propre et limites de la mesure empirique

Cette section rassemble les règles mécaniques et les critères structurels qui sous-tendent tous les comptages de la Partie A, et définit les termes de l’équation du rapport principal et officiel de l’article [A.11]. L’analyse vise à normaliser la surface, à séparer les biais géologiques et à encadrer les limites méthodologiques de la base de données [A.11].

1. Définition des zones de comparaison (la géométrie spatiale) :

איור 9

Figure 9. Où sont effectués les comptages : les calottes polaires (orange) contre la ceinture équatoriale de ±15° (bleu foncé ; la ceinture de ±30° en bleu clair pour comparaison). Les maria (gris) sont retirés des deux côtés ; South Pole–Aitken (contour rouge) est retiré du comptage propre. Points : cratères de 60 km et plus.

La procédure de comptage propre équilibré (retrait des biais)

Pour effectuer une comparaison équitable et neutre, la physique impose de comparer des unités crustales d’histoire évolutive semblable (sol ancien contre sol ancien) [A.11]. Pour cela, on applique un nettoyage simultané et équilibré de deux surfaces géologiques qui biaisent les enregistrements dans des directions inverses et opposées [A.11] :

1. Le système des maria et des plaines basaltiques (gris) : des fonds jeunes où l’enregistrement originel a été effacé et remis à zéro thermiquement [A.11]. Ces surfaces sont concentrées aux basses latitudes et représentent 22% de la ceinture équatoriale mais seulement 1% des calottes polaires [A.11]. Leur inclusion dans le comptage brut abaisse artificiellement la densité du dénominateur équatorial et gonfle à tort le rapport des pôles (rapport brut de 1.32) [A.7, A.11].

2. Le bassin Pôle Sud–Aitken (contour rouge) : une structure d’impact géante qui a effacé mécaniquement la population précoce de cratères au cœur de sa zone centrale (A.10) [A.11]. Le bassin se trouve sur la face cachée et couvre 24% de la surface de l’ensemble des calottes polaires (surtout au sud) [A.10, A.11]. Son inclusion dans le comptage brut abaisse artificiellement le numérateur polaire et biaise le rapport en faveur de l’équateur [A.11].

3. Le bilan normalisé : le retrait simultané de ces deux éléments des deux côtés de l’équation neutralise les deux biais opposés, et laisse une croûte montagneuse vierge et d’histoire stable (Highlands) aux pôles comme à l’équateur [A.11]. C’est sur ce comptage propre et équilibré que repose fermement la conclusion de l’article (D.1–D.2), qui fixe la valeur du rapport du système à [A.7, A.11].

Les cinq limites de la connaissance et les contraintes analytiques du comptage

Afin de préserver un standard académique sobre et retenu, la recherche définit explicitement cinq règles de base qui limitent et délimitent ce qu’un comptage de cratères peut prouver [A.11] :

Clôture de la phase descriptive et passage aux tests d’explication

Avec l’achèvement des sections A.1–A.11, la partie descriptive de l’article est entièrement close [A.11]. Les cinq étapes spatiales ont établi la base de données suivante :

Désormais, l’argument scientifique change de visage et passe de la question de la description empirique (ce qui se trouve sur le terrain et où) à la question du mécanisme physique (ce qui est capable de produire un tel schéma) [A.11]. Tout modèle ou hypothèse alternative examiné dans les chapitres suivants de l’article devra satisfaire à six contraintes limites strictes établies par l’observation, sans possibilité d’écart [A.11] :

PART B

Partie B — Ce que nous avons essayé avant la réponse

La Partie A a défini quantitativement et spatialement l’anomalie de densité : les cratères de la Lune sont nettement plus denses dans les régions des deux pôles que dans la ceinture équatoriale, l’effet étant organisé selon une configuration globale autour d’un anneau complet de la surface et non comme un point singulier isolé [B]. La Partie C développera et présentera les équations géodynamiques du mécanisme d’arrondissement et de contraction mécanique proposé comme solution scientifique [B].

Ce chapitre (Partie B) constitue la phase d’examen mécanistique préliminaire, et présente les six hypothèses et interprétations alternatives examinées et rejetées au cours de la recherche [B]. C’est leur rejet mécanique qui a conduit la recherche à la solution de l’arrondissement, car l’analyse combinée des six échecs prouve que la physique du système ne laisse place qu’à un seul type de cause génératrice [B].

La matrice des quatre questions structurelles :

Chaque hypothèse alternative est examinée de manière cohérente et sobre au regard de quatre questions strictes, dont les valeurs découlent directement des données d’observation de la Partie A [B] :

1. Le test du signe (Sign Test) : le mécanisme physique proposé crée-t-il effectivement un excédent de densité aux pôles par rapport à l’équateur, ou agit-il en sens inverse en les appauvrissant [B] ? (Un mécanisme qui appauvrit les pôles est rejeté immédiatement d’emblée) [B].

2. Le test de l’ampleur (Magnitude Test) : si le signe dynamique est correct, le mécanisme est-il capable de produire les rapports empiriques mesurés – dans le comptage brut, dans le comptage propre, ainsi que la tendance nette à la hausse à mesure que le diamètre du cratère augmente (D.1–D.2) [B] ?

3. Le test de la forme (Morphology Test) : le mécanisme place-t-il et organise-t-il l’excédent de densité dans l’espace géographique précis où il est observé – tenu aux deux pôles dans une mesure semblable, et disposé selon la gradation triaxiale des six axes (A.9) [B] ?

4. Le test du masque (Mask Test) : l’excédent de densité produit par le mécanisme est-il capable de survivre à l’application du masque géologique (retrait des maria, des grands bassins et du bassin SPA des deux côtés de l’équation), comme il survit effectivement dans les comptages résiduels de la Partie A [B] ?

Les tests sont ordonnés selon une gradation méthodologique, du test le plus grossier (signe) au test le plus fin (masque), si bien que la plupart des hypothèses sont tranchées et rejetées dès les premières étapes de l’analyse [B].

B.1 Conservation sur des terres hautes anciennes (Highlands Preservation Hypothesis)

Définition de l’interprétation :

Cette hypothèse soutient que l’excédent de densité aux pôles n’est que le reflet du degré de conservation de la surface et d’un âge géologique variable [B.1]. Les régions polaires constituent une croûte montagneuse ancienne et vierge (Highlands) qui n’a pas subi l’activité d’inondation volcanique tardive des plaines basaltiques et des maria sombres, lesquels ont effacé et recouvert l’enregistrement précoce aux basses latitudes de la face proche [B.1]. Par conséquent, les pôles sont denses parce qu’ils sont plus anciens et n’ont pas été affectés par des altérations endogènes, et non parce qu’un mécanisme dynamique quelconque y a activement concentré ou rassemblé des cratères [B.1]. C’est l’explication initiale admise dans la communauté planétaire, car la densité des cratères sert d’indicateur classique pour déterminer l’âge relatif des surfaces [B.1].

La prédiction morphologique de l’hypothèse (ce qu’elle prédit) :

Le modèle impose que la courbe de densité suive linéairement l’âge chronologique et le degré de conservation de la cellule de surface [B.1]. La densité doit être parfaitement uniforme et homogène sur toutes les unités de croûte montagneuse non inondées, et ne présenter des valeurs basses (creux) qu’aux endroits où s’étend une mare ou un bassin jeune [B.1]. L’hypothèse exclut l’existence d’un centre ou d’un apex anomal, car une terre haute ancienne du même âge est uniforme dans l’espace [B.1]. Si la conservation thermique était le seul facteur, le pôle présenterait des valeurs identiques à celles de tout autre espace montagneux du même âge à la surface du corps, et la carte de densité globale serait composée de plateaux plats et de creux locaux, sans courbe de pic concentrique [B.1].

Analyse des résultats empiriques des comptages (ce que montre le terrain) :

L’effet de conservation géologique et l’impact thermique des maria sont parfaitement réels, et façonnent effectivement la surface lunaire à la plus grande échelle (comme le reflètent les zones d’exposition de l’Article 11, où la face cachée présente une densité brute de 219 cratères par 10⁶ km², la ceinture des marges 164, et le disque de contact subterrestre seulement 99) [B.1]. Toutefois, trois observations fermes de la Partie A excluent l’interprétation et conduisent à son rejet [B.1] :

1. Robustesse de l’excédent sous le test décisif du masque : le masque de filtrage de la section A.3 omet du calcul précisément les éléments sur lesquels repose l’hypothèse – les polygones de maria et les bassins géants [B.3, B.1]. Cette procédure normalise le dénominateur et garantit une comparaison propre de terre haute contre terre haute des deux côtés de l’équation [B.1]. Si les différences de conservation étaient la seule cause du phénomène, le rapport de densité résiduel tomberait précisément à la valeur de distribution uniforme () avec le retrait de la croûte remplacée [B.1]. En réalité, le rapport brut chute de à (nettoyage des maria), et remonte à une valeur ferme de dès que le bassin SPA est lui aussi neutralisé [B.1]. Le rapport ne s’annule pas à 1, et continue de présenter une tendance à la hausse en fonction du diamètre du cratère jusqu’à une valeur de dans les grandes classes [B.1]. Toute différence de densité qui survit à cette procédure à l’identique ne peut être attribuée à l’écart entre sol inondé et sol vierge [B.1].

2. Le point maximal bipolaire au-dessus de la terre haute la plus ancienne du système : la couronne de la face cachée (autour des coordonnées du point anti-terrestre) représente dans le modèle dynamique l’unité crustale la plus ancienne et la plus sûre, que les surfaces de maria n’ont presque pas touchée géographiquement [B.1]. Après application du masque, ce secteur porte une densité de 221 cratères par 10⁶ km² [B.1]. À l’inverse, les deux calottes polaires résiduelles présentent des valeurs nettement plus élevées : 233 au nord et 246 au sud [B.9, B.1]. Les pôles se situent morphologiquement au-dessus de la terre haute la plus ancienne, et non au même niveau qu’elle [B.1]. Si l’âge chronologique de la croûte était le facteur déterminant, la couronne de la face cachée aurait dû mener la gradation globale ; les résultats de terrain l’excluent [B.1].

3. Existence d’une concentration radiale (un centre) dans le cœur nord vierge : dans l’espace du pôle Nord, où aucun polygone de mare ni bassin de Neumann ne pénètre dans un rayon de autour de l’apex, la densité n’est pas homogène et plate mais présente une courbe de concentration nette [B.1]. L’indicateur de densité monte continûment de 189 dans la calotte large () jusqu’à une valeur maximale de 300 cratères par million de km² dans la calotte intérieure de [B.6, B.1]. Une croûte montagneuse uniforme et statique ne peut pas produire une convergence radiale de densité vers un point géographique défini [B.1]. Le cœur nord constitue le site d’observation le plus propre pour réfuter l’hypothèse, car il ne contient aucun élément géologique étranger nécessitant un filtrage par le masque [B.1].

Verdict :

L’hypothèse de conservation sur des terres hautes anciennes ne réussit que le test du signe, mais échoue entièrement au test du masque et au test de la forme [B.1]. La large variation spatiale à l’échelle globale (comme la gradation face proche contre face cachée) s’explique par la géométrie des zones de contact de l’extrémité de la poire (Article 3) et par la distribution des vecteurs de retombées balistiques les plus violentes, concentrées sur la face cachée (Article 6), et non en fonction d’une conservation statique [B.1].

La recherche ne retient le mécanisme de conservation que comme la cause physique du fait que les surfaces de maria et les jeunes bassins d’impact présentent des valeurs de densité basses (remise à zéro de l’enregistrement), et non comme la source ou l’explication de l’excédent de densité polaire [B.1]. La présence des terres hautes anciennes de la face proche au sud du cratère Tycho, qui présentent une densité exceptionnelle de 302 cratères dans la plage , constitue un rappel mécanique supplémentaire que la densité d’une surface ancienne n’est pas fonction du seul âge chronologique [B.1]. Ce que l’hypothèse de conservation laisse sans réponse physique est précisément le cœur de la recherche de l’article : un excédent de densité significatif conservé au niveau de la comparaison terre haute contre terre haute, et qui possède un centre géométrique ancré sur l’axe de rotation du corps [B.1].

B.2 Aplatissement rotationnel pendant l’attachement (Rotational Flattening Hypothesis)

Définition de l’interprétation :

Cette hypothèse repose sur le principe mécanique selon lequel un corps visqueux/mou tournant rapidement autour de son axe subit une déformation d’aplatissement – il s’élargit le long de la ceinture équatoriale sous l’effet des contraintes centrifuges et se raccourcit le long de l’axe polaire, comme une pâte à pizza qui s’étire et s’aplatit quand on la fait tourner en l’air [B.2]. Comme la jeune Lune, pendant la phase attachée, tournait en synchronisation rapide avec la Terre (période d’environ 5.5 heures), elle a subi un aplatissement global prononcé (C.1) [B.2]. L’idée initiale soutient que ces contraintes de compression ont physiquement pressé les régions polaires vers l’intérieur, rapprochant ainsi mécaniquement les cratères qu’elles portaient et augmentant leur densité [B.2].

La prédiction morphologique (ce qu’elle prédit) :

La géométrie physique d’une déformation d’aplatissement à volume constant (situation dans laquelle aucun matériau n’est ajouté ni retiré du système, seule sa configuration changeant) présente un comportement entièrement opposé à la prédiction de resserrement [B.2] :

En conséquence, le nombre de cratères par unité de surface () diminue et s’appauvrit dans la région des pôles, et l’équation de répartition différentielle de la surface est définie précisément selon le rapport d’axes du sphéroïde (, où est le rayon équatorial et le rayon polaire) [B.2] :

Cette équation impose que plus un corps est aplati, plus ses pôles deviennent clairsemés et vides de cratères par unité de surface (l’inverse exact de l’équation d’arrondissement développée en C.3) [B.2]. Pour les rapports d’axes du sphéroïde de Maclaurin au pic d’aplatissement (, densité 1,300 , rotation de 5.5 heures) et dans les phases tardives d’arrondissement (, densité 2,600 , période de 5.77 heures), le mécanisme d’aplatissement laisserait les pôles dans une position d’infériorité théorique de seulement 0.59 et 0.83 par rapport à l’équateur (déficit structurel de et respectivement) [B.2].

Ce que disent les comptages empiriques :

Les données physiques de terrain montrent un excédent de densité prononcé aux pôles (rapport de dans le comptage brut et de dans le comptage propre) [B.2]. Par conséquent, le signe mécanique de l’hypothèse est fondamentalement erroné ; l’aplatissement rotationnel appauvrit les cratères à l’apex, alors que les observations de terrain les montrent à leur densité maximale [B.2].

Verdict :

L’hypothèse de l’aplatissement rotationnel échoue entièrement au test du signe [B.2]. Si l’on presse une balle élastique par le haut et par le bas, la surface de l’extrémité supérieure ne se contracte pas – elle s’étale et s’élargit, et les points qui y sont imprimés s’éloignent les uns des autres [B.2].

Toutefois, l’hypothèse saisit correctement la géométrie générale du corps lorsqu’on la lit en sens inverse : une surface crustale étirée malgré elle pendant la phase d’aplatissement, et libérée plus tard de la contrainte, se contractera fortement dans les régions polaires [B.2]. Cette compréhension inverse constitue le cœur mécanique du chapitre C.3 [B.2]. Ce constat a entraîné un changement méthodologique fondamental au cours de la recherche : au lieu de chercher un mécanisme dynamique qui produit l’aplatissement, on a commencé à chercher le processus géodynamique qui annule et libère l’aplatissement (le processus d’arrondissement de retour à la sphère) [B.2].

B.3 L’axe polaire raccourci par la marée — l’effet de Roche (Roche Tidal Stretching Hypothesis)

Définition de l’interprétation :

Pendant la phase attachée, la jeune Lune était soumise à l’effet de Roche – des contraintes de marée extrêmes qui ont étiré son volume le long de la ligne de la Terre (l’axe radial ) en raison de la différence dynamique entre la gravité terrestre, qui tirait la face proche vers l’intérieur, et l’entraînement centrifuge de la corotation, qui tirait la masse vers l’extérieur, alors que l’intérieur du corps était tenu sur l’extrémité de la structure en poire [B.3]. Comme le volume du corps restait presque constant, l’allongement de l’axe radial s’est fait aux dépens de la réduction des deux axes perpendiculaires [B.3]. Dans le système en rotation, le terme centrifuge a agi pour stabiliser et conserver l’axe de la direction du mouvement sur l’orbite () ; l’essentiel de la contrainte de raccourcissement et d’effondrement mécanique a donc porté sur l’axe polaire () [B.3]. Les pôles ont été pressés et rapprochés l’un de l’autre, et l’idée alternative soutient que les cratères qu’ils portaient ont été entraînés et resserrés avec eux [B.3].

La prédiction morphologique (ce qu’elle prédit) :

La géométrie du corps fluide libre d’une densité de 2,600 sous ces contraintes de marée se calcule précisément avec un rapport d’axes triaxial de (pour les axes ) [B.3]. Cette structure produit un corps ellipsoïdal très long en direction de la Terre et très court aux pôles (un œuf posé sur l’axe de marée) [B.3]. Lorsqu’on applique cette déformation d’étirement de marée à une croûte portant des cratères distribués uniformément, l’indicateur de densité de chaque cellule de surface varie en proportion inverse des facteurs de contraction (l’inverse exact des valeurs de C.3) [B.3] :

Le secteur spatial examinéLe vecteur de variation sous l’étirement de RocheLe vecteur de variation dans l’arrondissement de retour (C.3)
Les régions polaires (l’axe court )0.83× (dilution et éloignement)1.21× (resserrement mécanique)
Les points avant et arrière (axe )0.88× (dilution modérée)1.14× (resserrement modéré)
Le point subterrestre (l’axe long )1.38× (concentration et resserrement)0.73× (élargissement et dilution)

La lecture de la matrice dans sa largeur prouve qu’un axe polaire raccourci mécaniquement sous des contraintes d’étirement de marée étale et élargit la surface située au-dessus de lui, exactement comme dans le mécanisme d’aplatissement rotationnel [B.3]. Le pôle constitue le secteur physiquement le plus étiré de la forme ; les cratères y sont donc éloignés les uns des autres et dilués () [B.3]. À l’inverse, les extrémités de la ligne de la Terre (le point subterrestre) sont les zones où la surface converge et se resserre jusqu’au maximum du système () [B.3].

Ce que disent les comptages empiriques :

Les observations réelles de terrain présentent un état entièrement inverse : les pôles sont les zones les plus denses du système, tandis que le point subterrestre est le plus clairsemé et le plus vide des six calottes principales (valeur brute de 94, et valeur masquée de seulement 148) [A.9, B.3]. Le mécanisme d’étirement de marée prédit précisément les inversions opposées dans ces deux secteurs [B.3].

Verdict :

L’hypothèse de l’étirement de Roche échoue au test du signe en tant que mécanisme générateur de resserrement [B.3].

Toutefois, l’effet physique lui-même est réel et son ampleur est fermement conservée dans les équations du modèle : c’est lui qui définit la forme de départ triaxiale de la section C.2, et l’arrondissement géodynamique de retour à partir de celle-ci vers une sphère parfaite est la réponse scientifique de la recherche [B.3]. L’effet de Roche est responsable de la fixation du pic de déformation et d’étirement du corps attaché ; c’est l’arrondissement qui suit (la libération de l’étirement et le retour à la sphère) qui déplace et rapproche physiquement les cratères [B.3]. La lecture de la colonne de droite du tableau (le vecteur d’arrondissement) restitue en accord complet la gradation des comptages : pôles maximaux, faces intermédiaires, et point subterrestre à la base du minimum [B.3].

B.4 Gravité plus forte au pôle (Polar Gravity Excess Hypothesis)

Définition de l’interprétation :

Dans un corps planétaire de forme aplatie, l’apex du pôle se situe géométriquement plus près du centre de masse du corps que la ceinture équatoriale [B.4]. En outre, au pôle géographique (situé directement sur l’axe de rotation), la vitesse linéaire s’annule, et la composante de poussée centrifuge vers l’extérieur n’existe donc pas [B.4]. La combinaison de ces deux facteurs produit aux pôles une gravité effective (différence entre l’attraction propre et la contribution centrifuge) supérieure à celle de l’équateur [B.4]. L’interprétation alternative propose que ce champ de gravité accru aux pôles ait agi comme une pompe mécanique qui aurait attiré et déplacé le matériau de surface élastique et meuble le long de la croûte vers les régions polaires, resserrant ainsi l’enregistrement autour d’elles [B.4].

La prédiction morphologique (ce qu’elle prédit) :

La loi physique impose que le vecteur de gravité agisse vers le bas (radialement vers le centre) et non latéralement (tangentiellement le long de la surface) [B.4]. Pour que du matériau crustal subisse un mouvement latéral actif le long de la surface vers le pôle, le pôle doit fonctionner physiquement comme le point le plus bas du potentiel effectif du corps (le fond d’une vallée de potentiel, vers lequel le matériau glisse et s’écoule comme l’eau sur une pente) [B.4].

Dans le système dynamique du corps attaché tournant rapidement avec la Terre, la situation physique est entièrement inverse [B.4]. Dans le système en rotation, le potentiel effectif dans la ceinture équatoriale est inférieur à celui du pôle d’un facteur strict de (où est la vitesse de rotation et le rayon), valeur qui représente l’énergie cinétique de mouvement excédentaire que possède la ceinture équatoriale et qui fait entièrement défaut sur l’axe polaire [B.4].

Dans les phases tardives de l’arrondissement (densité 2,600 , période de 5.77 heures, rayon moyen ), pour une gravité de surface de (C.7), l’équation de hauteur théorique est définie par [B.4] :

Cette équation prouve que, sur un corps aplati qui n’a pas encore achevé son effondrement vers sa forme d’équilibre hydrostatique, les régions polaires se dressent comme des collines topographiques géantes d’une hauteur d’environ 120 km « au-dessus » de la ceinture équatoriale au sens énergétique [B.4]. Elles constituent les points de potentiel le plus élevé (sommet de colline), et tout matériau meuble ou toute croûte élastique subira donc un vecteur de poussée inverse et glissera à partir d’elles vers l’extérieur, vers la vallée de potentiel équatoriale [B.4]. C’est le mécanisme qui produit l’aplatissement (B.2), et non un resserrement polaire [B.4].

Ce que disent les comptages empiriques :

L’excédent se situe aux pôles [B.4]. Pendant la rotation du corps, la combinaison de la gravité et de l’accélération centrifuge a agi pour éloigner et diluer le matériau des pôles, et non pour l’y concentrer [B.4].

Verdict :

L’hypothèse de la gravité accrue au pôle échoue entièrement comme mécanisme de déplacement latéral tant que le corps se trouve dans la phase attachée et en rotation [B.4].

L’effet ne trouve sa place physique correcte qu’au moment de la transition dynamique : c’est la gravité propre de la Lune qui l’a maintenue comme un corps entier et a empêché sa dislocation au-delà de la limite de Roche rigide ( à une densité de 1,300) [B.4]. Au moment où le centre de masse de la Lune a franchi la ligne de corotation et où la compression mécanique sur la base s’est annulée, le corps s’est libéré de la contrainte de rotation synchrone excédentaire [B.4]. Dans cet état, l’aplatissement que la rotation maintenait auparavant devient une charge excessive que le corps ne peut pas porter [B.4]. Ce n’est qu’à ce moment que le pôle du corps aplati s’effondre et passe d’un sommet de colline énergétique au fond d’une profonde vallée de potentiel [B.4]. Désormais, la gravité propre agit comme le moteur principal qui entraîne l’écoulement et la contraction de la surface vers l’intérieur, vers les pôles – c’est le mécanisme d’arrondissement de retour (C.4) [B.4]. L’écart géométrique de 120 km représente l’échelle précise de l’énergie que la rotation maintenait en hauteur, et sa libération mécanique lors de l’arrondissement est ce qui resserre les cratères [B.4].

B.5 Puits de drainage au pôle (Polar Drain/Vortex Hypothesis)

Définition de l’interprétation :

Cette hypothèse repose sur l’analogie mécanique intuitive d’une « aspiration vers l’intérieur » : l’existence d’un mécanisme de drainage, d’un affaissement crustal ou d’un tourbillon dynamique interne sous chaque pôle géographique [B.5]. Selon ce modèle, le système magmatique ou géodynamique agit comme l’eau sur le fond d’un évier qui s’écoule et se rassemble vers la bonde, attirant, concentrant et rassemblant ainsi activement la surface et les cratères qu’elle porte vers un point extrême unique [B.5].

La prédiction morphologique (ce qu’elle prédit) :

La dynamique d’un puits de drainage impose une concentration dirigée exclusivement vers un point singulier (The Vertex Point) [B.5]. Par conséquent, la carte de densité doit présenter une tendance à la hausse raide et illimitée (singularité statistique) à mesure qu’on s’approche physiquement du point du pôle lui-même, la topographie géophysique devant présenter un creux ou un puits profond au point central [B.5]. Le modèle exige l’action d’un vecteur de force mécanique latéral dirigé directement vers les coordonnées du pôle [B.5].

Ce que disent les comptages empiriques :

Les données mesurées sur le terrain excluent la structure de singularité ponctuelle [B.5] :

1. Arrêt et stabilisation du cœur nord vierge : dans l’espace propre du pôle Nord, la densité des cratères () monte effectivement de manière radiale : de 189 dans la calotte de , à 227 (), et à 270 () [A.6, B.5]. Toutefois, en pénétrant dans le cœur extrême, la courbe de densité s’arrête complètement et se stabilise : valeur de 300 dans la calotte de contre 298 dans la calotte de [A.6, B.5]. Les données révèlent une concentration globale vers une large région polaire (Highlands Plateau) et non une convergence ponctuelle dans un puits de drainage [B.5].

2. Sensibilité de l’échantillon sud résiduel : au sud, après application du masque, les calottes intérieures présentent des valeurs de 305 (à 15°) et 348 (à 10°) [A.7, B.5]. Toutefois, l’analyse du système montre que, dans un rayon de autour de l’apex du pôle Sud, le masque se compose uniquement d’une union coïncidant avec le bassin SPA, et la calotte de repose sur un échantillon restreint de seulement 36 cratères [B.5]. Dans un échantillon aussi petit, une fluctuation aléatoire de quelques cratères seulement déplace la valeur numérique dans de larges proportions ; la hausse au sud porte donc la géométrie locale du bassin d’impact et les marges d’erreur du comptage, et non la preuve d’un puits de drainage [B.5]. Le poids scientifique du test repose sur la séquence nord propre [B.5].

3. Absence de preuves topographiques : les données topographiques planétaires réelles ne montrent aucune paire de dépressions gravitationnelles ou structurelles correspondant aux pôles géographiques [B.5].

Verdict :

L’hypothèse du puits de drainage réussit le test du signe (positif), mais échoue entièrement au test de la forme [B.5].

Dans un corps planétaire en rotation, tout vecteur de force physique doit conserver une symétrie axiale autour de l’axe de rotation, et une telle symétrie produit une structure géométrique en anneau (Ring) et non en point unique ; il n’existe aucune force physique capable de converger vers le point du pôle lui-même [B.5]. L’excédent de densité réel règne le long d’un anneau complet – l’anneau situé à une distance de du point subterrestre, qui passe de façon continue par les deux pôles et les deux faces (A.9, C.3) – et les pôles mènent la gradation uniquement parce que l’axe de rotation fixe leur centre de convergence géométrique [B.5]. L’interprétation du puits de drainage possède le bon signe (excédent) mais un schéma spatial qui produit une déformation : elle exige un pic ponctuel, alors que les observations montrent un niveau constant sur une large cellule de surface et une bande de resserrement globale qui entoure le corps entier [B.5].

B.6 Changement uniforme (Isotropic Volume/Scale Change Hypothesis)

Définition de l’interprétation :

Cette hypothèse propose d’expliquer l’excédent de densité par des processus isotropes (uniformes) globaux, tels que : un gonflement de volume uniforme, une contraction thermique uniforme de tout le corps au cours de ses phases de refroidissement, ou l’action d’une forte gravité propre attirant tous les matériaux de la croûte de manière égale vers le centre de masse – c’est-à-dire un changement mécanique de la dimension et de la taille du corps () sans changement de sa forme géométrique (structure) [B.6]. Chacun de ces événements constitue un processus physique réel et documenté dans l’histoire précoce d’un corps planétaire jeune et chaud, et chacun d’eux déplace et modifie la position de l’ensemble de la surface [B.6].

La prédiction morphologique (ce qu’elle prédit) :

Un changement uniforme de la dimension du corps modifie chaque élément et chaque unité de surface crustale exactement du même facteur et dans la même proportion [B.6]. L’intervalle géographique entre les cratères augmente ou diminue de manière homogène en chaque coordonnée de la surface du corps, et l’indicateur de densité varie donc d’un facteur constant et uniforme aux pôles comme dans la ceinture équatoriale [B.6] :

Pour une surface ayant subi une cratérisation uniforme, un changement global de taille multiplie le numérateur (le nombre de cratères) et le dénominateur (la surface) par exactement le même facteur mathématique ; la valeur du rapport statistique entre eux ne peut donc pas changer [B.6]. Le mécanisme fonctionne comme le gonflement uniforme d’un ballon : l’opération rend tous les points qui y sont imprimés plus clairsemés, mais elle préserve l’égalité et ne peut pas produire une zone plus dense ou plus clairsemée qu’une autre [B.6].

Ce que disent les comptages empiriques :

Les observations de terrain établissent un écart morphologique prononcé : rapport de dans le comptage brut, de dans le comptage propre, et tendance à la hausse régulière en fonction du diamètre du cratère jusqu’à atteindre une valeur de dans les grandes classes [B.1, B.6].

Verdict :

L’hypothèse du changement uniforme échoue absolument tant au test du signe qu’au test de l’ampleur, car elle ne peut produire ni excédent de densité ni variation spatiale [B.6].

Toutefois, c’est l’interprétation méthodologique qui a tracé la voie analytique de la recherche [B.6]. Une fois exclu physiquement un processus de changement uniforme de taille, la physique a imposé une conclusion ferme : un changement de taille ne peut pas produire un excédent polaire ; seul un changement de forme du corps (Anisotropic Distortion) peut produire une telle anomalie [B.6]. Ce changement de forme doit être non uniforme dans l’espace de manière définie : il doit contracter mécaniquement la surface dans les régions polaires au taux maximal, et il doit conserver une symétrie axiale complète autour de l’axe de rotation [B.6].

Comme les hypothèses B.2 et B.3 représentent des changements de forme réels mais avec le mauvais signe mécanique (dilution des pôles), la solution scientifique de l’énigme exige leur inversion géodynamique exacte [B.6]. Cette étape a resserré la question scientifique : la recherche s’est orientée exclusivement vers le changement de forme spécifique qui contracte la surface polaire plus que toute autre zone, et il n’existe qu’un seul processus géodynamique dans l’histoire du modèle qui satisfasse cette contrainte [B.6].

B.7 La réponse, et pourquoi c’est elle qui demeure (The Flattening Release / Spherical Rebound)

Un seul changement de forme géodynamique dans le système fournit une réponse précise et complète à toutes les contraintes limites fixées par les observations et auxquelles les interprétations précédentes ont échoué : le passage mécanique d’un corps triaxial très aplati qui se libère de ses contraintes et s’arrondit pour revenir à la forme d’une sphère parfaite avec une diminution spontanée de volume [B.7].

Clôture des contraintes et accord avec le dispositif physique :

1. Le signe mécanique (positif) : le mécanisme constitue l’inversion géométrique et mécanique exacte des hypothèses B.2 et B.3 [B.7]. La libération de l’étirement et la contraction qui suit agissent avec le bon signe et resserrent la croûte [B.7].

2. Le moteur physique : il est actionné par les forces de gravité propre de la section B.4 à leur place physique exacte – au moment où le centre de masse de la Lune franchit la ligne de corotation et où la force de compression nette du corps sur la base de la poire s’annule [B.4, B.7]. Le corps se libère de la contrainte centrifuge excédentaire tout en restant attaché et en tournant en synchronisation avec la Terre [B.7].

3. La géométrie spatiale : elle produit précisément la structure en anneau symétrique autour de l’axe de rotation (B.5), les régions polaires menant le pic de resserrement et de contraction mécanique parce que l’axe de rotation passe par elles et que tous les méridiens se rétractent vers elles (les pôles sont les seuls points que la rotation n’a jamais déplacés) [B.5, B.7].

4. Le test de l’ampleur et la cohérence chronologique : les dimensions physiques découlent directement de la forme que les cratères eux-mêmes enregistrent autour des six axes (A.9) et des tests d’accord avec la forme d’équilibre hydrostatique du sphéroïde de Maclaurin (C.1) [B.7]. Le modèle d’arrondissement ne prédit pas une valeur statique, mais une gradation chronologique dynamique qui concorde parfaitement avec les résultats empiriques du comptage propre [B.7] :

5. Le test du masque : l’effet est robuste et survit entièrement au nettoyage des maria, des bassins et du système SPA, car il est inscrit dans la mécanique de la croûte solide elle-même [B.7].

6. Moment et fenêtre de formation définis : le processus se situe chronologiquement dans la phase d’attachement précoce (avant le détachement mécanique final, lors duquel la Lune s’est libérée dans un état rond et rigide avec un rapport d’axes de ) [B.7]. L’intérieur du corps était suffisamment mou et visqueux pour permettre l’écoulement et l’arrondissement, tandis que l’enveloppe extérieure était suffisamment rigide pour porter les cratères comme des éléments intacts et isolés (la contraction a refermé les intervalles élastiques entre eux et a accumulé la contrainte de C.4) [B.7].

Signification de la hausse avec le diamètre du cratère :

La tendance à la hausse du rapport de densité à mesure que le diamètre du cratère augmente constitue la preuve empirique la plus claire et la plus directe de l’existence d’un processus d’arrondissement se déroulant au fil du temps [B.7]. Un cratère crustal ne peut conserver et enregistrer que la part de contraction de surface survenue après sa formation [B.7]. Les plus grands cratères (les plus anciens) ont subi la transformation entière, tandis que les petits cratères (les plus tardifs) n’en ont enregistré que l’accord final géométrique [B.7]. Aucune hypothèse alternative (de B.1 à B.6) ne peut produire une dépendance structurelle entre l’excédent de densité et le diamètre du cratère ; le modèle d’arrondissement l’exige physiquement [B.7].

La matrice de synthèse des parties de l’examen mécanistique

Voici un regroupement quantitatif et sobre des six interprétations rejetées face à l’hypothèse d’arrondissement retenue, qui définit la contribution résiduelle de chaque mécanisme rejeté à la formulation scientifique finale [B.7] :

N° de sectionL’hypothèse examinée dans le systèmeTest du signeTest de l’ampleurTest de la formeTest du masqueLa composante mécanique conservée dans la description finale
B.1Conservation sur des terres hautes anciennesOui—NonNonDéfinit le biais des maria et des bassins bas (remise à zéro de l’enregistrement) [B.7].
B.2Aplatissement rotationnelNon———Définit la contrainte d’étirement initiale que l’arrondissement libère [B.7].
B.3Effet de Roche (étirement de marée)Non———Définit la géométrie de départ triaxiale de la section C.2 [B.7].
B.4Gravité polaire accrueNon———Agit comme le moteur principal qui entraîne l’effondrement d’arrondissement (C.4) [B.7].
B.5Puits de drainage / tourbillonOui—Non—Confirme la géométrie axiale de l’anneau global [B.7].
B.6Changement uniforme (dimension/taille)NonNon——Prouve la nécessité physique exclusive d’un changement de forme non uniforme [B.7].
CLe modèle d’arrondissement et de diminution de volumeOui (propre)1.20 / 1.5–1.2Anneau polaireOuiLa réponse scientifique confirmée de la recherche [B.7]

Lisez la matrice de haut en bas : chaque interprétation rejetée a laissé derrière elle une composante physique ferme, et la dernière ligne les réunit en une seule structure cohérente [B.7]. L’étirement de B.2 et B.3 fournit les indicateurs de déformation à annuler, la gravité propre de B.4 fournit le vecteur moteur, la symétrie axiale de B.5 définit la géométrie, B.6 établit la loi d’un changement de forme nécessaire, et B.1 isole l’effet des surfaces jeunes [B.7].

Nous disposons maintenant d’une cause physique exclusive qui réussit les quatre tests empiriques, et de six hypothèses rigoureusement rejetées [B.7].

PART C

Partie C — Le mécanisme géodynamique et ses chiffres

La Partie B de la recherche a rejeté les six explications alternatives [B.7]. Le rejet analytique a laissé une seule cause mécanique possible : le passage géométrique d’un corps triaxial très aplati qui se libère de ses contraintes et s’arrondit pour revenir à la forme d’une sphère parfaite, avec une diminution spontanée de volume et une hausse de densité [B.7, C]. Cette partie transforme cette réponse en mécanisme physique quantitatif : elle définit les dimensions de la forme de la jeune Lune, les vecteurs de contrainte qui l’ont modifiée, l’effet spatial sur la croûte, et l’énergie qui a entraîné le processus [C].

La séquence des forces de déformation et l’évolution morphologique :

La géométrie de la jeune Lune a été déterminée par une équation d’équilibre dynamique entre trois vecteurs de force principaux [C] :

1. L’attraction de la gravité propre : elle tendait à presser uniformément toutes les parties du corps vers le centre, et à elle seule elle aurait produit une sphère parfaite [C].

2. L’accélération centrifuge de la rotation rapide : pendant la phase d’attachement synchrone avec la Terre (période d’environ 5.5 heures), la rotation a poussé la ceinture équatoriale vers l’extérieur et abaissé les régions polaires [C].

3. Les contraintes d’étirement de marée (effet de Roche) : elles ont étiré le corps à l’extrême le long de l’axe radial entre la Terre et la Lune [B.3, C]. La gravité de la Terre a tiré la face proche vers l’intérieur, tandis que la force centrifuge de la corotation a tiré la masse vers l’extérieur, la base du corps étant physiquement tenue sur l’extrémité de la structure en poire [C]. Cet étirement a raccourci les deux axes perpendiculaires, et a contraint l’axe polaire à absorber le maximum de raccourcissement mécanique [C].

La chronologie de la séquence de forme se déploie au fil des phases mécaniques de séparation [C] :

C.1 La forme initiale au pic d’aplatissement : rotation et densité

La question de cette section quantifie mathématiquement à quel point la Lune était aplatie et déformée au pic de la phase d’aplatissement, et comment sa forme a évolué au fil des phases d’effondrement vers l’intérieur [C.1]. Le calcul de la composante de rotation ne requiert que deux paramètres physiques : la vitesse angulaire de rotation () et la densité moyenne du matériau () [C.1].

Pour un fluide homogène autogravitant (un corps de densité uniforme maintenu uniquement par sa propre gravité, dépourvu de résistance mécanique initiale), la forme d’équilibre hydrostatique exacte est définie comme le sphéroïde de Maclaurin (Maclaurin Spheroid) [C.1]. L’écart du sphéroïde par rapport à une sphère parfaite se mesure par son excentricité, (où est le rayon équatorial et le rayon polaire ; donne une sphère parfaite, et un disque aplati) [C.1]. Le corps reste stable au repos lorsque la vitesse angulaire (exprimée par , où est la période de rotation) satisfait l’équation structurelle [C.1] :

Matrice de solution du rapport d’axes dynamique () :

La résolution de l’équation de Maclaurin pour les valeurs des paramètres des phases de détachement donne la matrice suivante des rapports du rayon équatorial au rayon polaire [C.1] :

Période de rotation propre ()
5.50 heures (pic d’aplatissement structurel)Pas d’équilibre1.691.461.411.321.221.16
5.77 heures (tard dans l’arrondissement de la croûte)2.221.561.391.351.281.201.14

La limite inférieure de stabilité pour la densité du matériau :

Pour une période de rotation de 5.5 heures et une densité de 1,000 , l’équation ne fournit pas de solution d’équilibre, car la contrainte centrifuge dépasse le point maximal admissible du membre de droite [C.1]. Dans cette situation hypothétique, le corps se disloquerait et perdrait de la masse aux marges de son équateur [C.1]. Ce point définit une limite inférieure stricte pour la densité du matériau nécessaire au maintien de l’intégrité du corps à toute vitesse de rotation :

Pour une période de rotation de 5.5 heures, la limite de densité minimale est fixée mathématiquement à 1,070 ; un corps planétaire plus léger que cette valeur ne peut pas se maintenir comme fluide sous une telle vitesse de rotation [C.1].

Position du pic de déformation et géométrie triaxiale

L’enregistrement empirique de l’excédent de cratères restitue et confirme précisément ces dimensions (sections A.9, C.3) : la forme d’origine au pic présentait un rapport d’environ , dans une plage de dispersion morphométrique comprise entre et [C.1]. La population globale () converge vers une moyenne combinée de [C.1]. Il s’agit d’une déformation modérée et physiquement claire : l’axe long est 1.6 à 1.95 fois plus grand que l’axe polaire, et l’axe moyen 1.33 à 1.45 fois plus grand que celui-ci [C.1].

Confirmation des indicateurs de densité au pic () :

La détermination de la valeur de densité du corps au pic d’aplatissement repose sur le croisement de deux vérifications dynamiques indépendantes, qui concordent parfaitement [C.1] :

1. Confirmation de la partie pita selon le rapport de Maclaurin : la partie pita de la forme () est déterminée exclusivement par les contraintes de rotation [C.1]. En première approximation, la différence entre l’axe moyen et l’axe court d’un corps étiré par la marée dépend de la vitesse de rotation propre, et peut donc être lue directement dans l’équation de Maclaurin [C.1]. Pour une période de rotation de 5.5 heures, le rapport de Maclaurin donne une valeur exacte de pour un corps d’une densité de 1,700 (la valeur en gras du tableau) [C.1]. La plage de dispersion empirique () encadre la densité entre 1,615 et 1,960 , le cœur se concentrant à 1,600–1,750 [C.1].

2. La limite géométrique de fermeture de la surface : le passage d’un corps étiré à une sphère parfaite entraîne une libération de l’étirement et une diminution globale de volume [C.1]. Les équations géométriques (C.3) prouvent que, si la densité du corps au pic avait été supérieure à 1,800 , le volume d’étirement disponible aurait été trop petit ; dans ce cas, les extrémités de l’axe de marée auraient dû s’élargir et croître pendant l’arrondissement au lieu de se contracter et de se refermer [C.1].

Les deux tests fixent la valeur du corps au pic comme un corps étiré, gonflé et d’une densité d’environ 1,700 , ce qui donne un rayon de sphère équivalente (rayon d’une sphère de même volume) d’environ 2,180 km [C.1]. Cette valeur contraste avec le rayon de 1,890 km à la libération et le rayon moderne de 1,737 km [C.1]. La masse sialique d’origine sur l’extrémité de la poire (, rayon équivalent d’environ 2,380 km) constitue la base de départ ancienne, avant le déclenchement de la phase d’aplatissement dynamique [C.1].

Les phases de l’arrondissement mécanique

Au franchissement de la ligne de corotation par le centre de masse, le corps a commencé le processus d’arrondissement de retour à la sphère, qui s’est effectué en deux phases mécaniques distinctes [C.1] :

1. Libération de la partie œuf (l’axe de marée ) : elle s’est libérée en premier et rapidement, car sur la ligne de corotation le vecteur d’attraction nette qui pressait le corps sur la base de la poire est tombé à zéro [C.1, C.6]. Cette part de l’étirement le long de la ligne de la Terre s’est effondrée vers l’intérieur [C.1].

2. Libération de la partie pita (l’aplatissement rotationnel ) : elle a progressé parallèlement à l’épaississement et au resserrement du matériau sous une période de rotation presque constante (C.4) [C.1]. Cette phase se définit par la hausse des indicateurs de densité () et un léger allongement du jour terrestre (jusqu’à atteindre heures au moment de la libération) [C.1]. Sous cette dynamique, le rapport d’axes a chuté et s’est fixé à une valeur d’environ 1.20 dans les phases tardives de l’arrondissement [C.1].

Un corps fluide qui effectue une rotation propre toutes les quelques heures doit physiquement présenter une déformation aplatie et ne peut exister comme sphère parfaite à l’état visqueux [C.1]. Le passage final et strict d’un rapport d’environ 1.20 à un corps sphérique parfait au moment de la libération est stabilisé et renforcé par la rigidification de l’enveloppe extérieure [C.1]. C’est une contrainte structurelle du modèle : une croûte supérieure dotée d’une résistance mécanique et d’une rigidité de est capable d’absorber les contraintes de rotation résiduelles et de stabiliser un corps rond d’un rapport d’axes de 1.04 au plus [B.7, C.1]. Au moment de la libération, la Lune est ronde et son enveloppe rigide [C.1]. L’aplatissement résiduel mineur qui subsiste après la libération ne constitue qu’un petit renflement élastique, qui s’atténue et disparaît entièrement avec l’éloignement du corps et le ralentissement de sa vitesse de rotation sur l’orbite (C.5) [C.1].

L’équation de Maclaurin est une loi physique fondamentale, et les solutions de la matrice ne dépendent de rien d’autre que des valeurs de la vitesse angulaire () et de la densité du matériau () [C.1].

C.2 La composante de marée : la forme libre et l’effet de Roche

Alors que la section C.1 a calculé l’effet des seules contraintes de rotation, cette section intègre le second vecteur de force qui a façonné morphologiquement le corps : les forces de marée de la Terre (effet de Roche) et le bilan complexe des potentiels dans le système [C.2]. L’analyse sert d’outil géométrique pour illustrer la contribution de l’étirement radial, et ne constitue pas la base du test statistique des densités [C.2].

La dynamique à deux corps dans le système en rotation :

Un corps planétaire situé à proximité immédiate de la Terre et en verrouillage de rotation synchrone (période orbitale égale à la période de rotation propre ) est étiré à l’extrême le long de l’axe radial entre lui et la Terre [C.2]. Pendant la phase attachée, la masse est tenue à l’extrémité de la structure en poire, sous un équilibre de contraintes opposées : la gravité de la Terre tire le côté proche vers l’intérieur, et la force centrifuge le tire vers l’extérieur [C.2].

Le modèle physique classique restitue cette structure par le calcul de la forme d’un fluide homogène libre (sans contrainte de contact), qui prend la configuration d’un ellipsoïde de Roche (Roche Ellipsoid) à trois axes distincts () : l’axe de marée en direction de la Terre (), l’axe de la direction du mouvement orbital (), et l’axe court reliant les pôles () [C.2]. Contrairement à la forme biaxiale de Maclaurin (qui apparaît comme un cercle parfait vue du dessus du pôle), l’ellipsoïde de Roche est déformé et ovoïde sous tous les angles de vue [C.2].

L’équation du bilan des potentiels hydrostatiques :

La surface d’un corps fluide au repos doit physiquement constituer une surface équipotentielle (Equipotential Surface), afin d’empêcher l’écoulement du matériau le long de celle-ci [C.2]. À l’intérieur d’un ellipsoïde homogène, le potentiel total se compose de la somme mathématique de trois composantes d’énergie [C.2] :

1. La gravité propre du corps : l’énergie croît comme le carré de la distance au centre le long des axes (, parallèles aux axes ), et les coefficients statiques (qui satisfont l’équation d’indice de Chandrasekhar : ) déterminent la raideur de la hausse [C.2]. L’axe le plus court () reçoit le coefficient maximal (), car la masse est concentrée plus près de lui [C.2].

2. Le potentiel de marée de la Terre : exprimé par , où le facteur d’étirement est ( est la masse de la Terre et la distance orbitale) [C.2]. Cette composante abaisse le potentiel le long de l’axe de marée () et l’élève transversalement (), action qui exerce un étirement radial et une compression tangentielle [C.2].

3. Le potentiel centrifuge orbital : exprimé par (selon la troisième loi de Kepler dans le système : ) [C.2]. Cette composante abaisse le potentiel vers l’extérieur, loin de l’axe de rotation, et laisse la composante polaire () inchangée [C.2].

Pour que la surface s’équilibre en une forme fermée, l’attraction nette vers l’intérieur (la gravité moins les facteurs de rotation et de marée), multipliée par le carré de la longueur de l’axe, doit être exactement égale aux trois extrémités [C.2] :

Cette équation combinée fournit deux systèmes de conditions pour les deux rapports d’axes indépendants : et [C.2].

Matrice de solution des équations de l’ellipsoïde de Roche

Voici les valeurs numériques exactes obtenues par la résolution du système d’équations, calculées à une distance de 20,595 km (un peu au-delà de la limite de Roche fluide d’un corps d’une densité de 2,600 ), qui sert d’outil d’illustration morphologique [C.2] :

Distance orbitale ()Densité du matériau ()Rapport d’axes triaxial ()Facteur de soulèvement polaire ()Facteur d’enfoncement radial ()État mécanique et contraintes de la limite de Roche
19,660 km2,600Pas d’équilibre——À l’intérieur de la limite de Roche fluide ; le corps se disloque par effet de marée [C.2].
20,595 km2,6001.2080.726Convergence de la forme étirée de Roche [C.2].
20,595 km3,3441.1310.821Corps plus dense ; la déformation de marée s’atténue [C.2].

Conséquences géodynamiques et limites de la limite de Roche rigide

Le calcul mathématique établit deux lois physiques fondamentales qui contraignent les phases de détachement [C.2] :

1. Rejet de la libération fluide légère : à la densité source de la poire (), la limite de Roche fluide se situe très loin, à une distance d’environ 25,200 km [C.2]. Cela signifie que le corps ne peut pas physiquement atteindre sa phase de libération dynamique tout en conservant cette faible densité, car il se disloquerait par effet de marée ; cette donnée confirme l’hypothèse de base de l’Article 2, qui a établi que la densité doit monter à 2,600 pour ramener les limites de Roche vers l’intérieur [C.2].

2. Effondrement de la limite de Roche fluide exacte : à une densité de 2,600, la résolution des équations montre qu’à 19,660 km la forme se disloque, tandis qu’à 20,595 elle est stable ; la limite de Roche fluide exacte se fixe entre les deux – à 19,971 km selon le coefficient classique, et à 20,100 km () selon la solution exacte des potentiels [C.2].

La limite de la croûte rigide dans la séquence dynamique :

Dans le modèle LEF mis à jour, la Lune se libère et se détache mécaniquement bien plus tôt et bien plus près de la Terre – à une distance d’environ 16,400 km sur l’équateur ( pour une période de 5.77 heures) et de 17,287 km aux latitudes de [C.2]. À ces distances proches, la physique impose qu’aucune forme fluide libre stable n’existe ni ne peut exister à une densité de 2,600 (la zone se trouve profondément à l’intérieur de la limite de Roche fluide) [C.2].

La signification physique est nette : la jeune Lune a tenu et conservé son intégrité structurelle dans cette zone grâce à son enveloppe extérieure, rigidifiée et figée tôt, et non en vertu d’un équilibre hydrostatique fluide [C.2]. Par conséquent, le modèle triaxial de ce chapitre sert uniquement d’indicateur et d’outil d’illustration montrant comment les contraintes de marée extérieures ont imposé un étirement radial prononcé le long de l’axe en plus de l’aplatissement rotationnel de l’axe (C.1) [C.2]. L’excédent de densité et l’anomalie statistique mesurée sur le terrain sont exclusivement le produit de cette libération combinée des contraintes au cours de la séquence d’arrondissement et de diminution de volume (C.1, C.5, D.2) [C.2].

Alors que les sections C.1 et C.2 ont défini les dimensions de la déformation initiale, cette section calcule la transformation de la surface déformée au cours de son effondrement mécanique vers une forme sphérique parfaite, sous condition de diminution de volume (resserrement et durcissement du noyau) [C.3]. L’objectif est de déterminer mathématiquement quelles cellules de surface crustale se referment et se contractent au taux maximal, lesquelles au taux minimal, et où les cratères existants se resserrent [C.3]. Le résultat du calcul est un indicateur de rapport de densités différentiel entre les pôles et l’équateur, robuste et indépendant du taux de contraction isotrope général du corps [C.3].

L’analogie du ballon écrasé :

Pour rendre la mécanique accessible, imaginons un ballon rond couvert de points (cratères) répartis uniformément par unité de surface, pressé et écrasé fermement entre deux paumes plates (les contraintes à deux corps) [B.5, C.3]. Sous l’effet de la pression, ses parties supérieure et inférieure deviennent larges, aplaties et plates, tandis que la bordure périphérique (l’équateur) s’étire fortement [B.5, C.3].

Au moment où les contraintes extérieures se relâchent (franchissement de la ligne de corotation), le ballon libère l’énergie élastique résiduelle et reprend d’un bond sa forme sphérique [C.3]. Dans ce processus géométrique, les surfaces larges et plates des extrémités se rassemblent rapidement vers l’intérieur et convergent vers une zone restreinte : les points qui y étaient imprimés se rapprochent physiquement les uns des autres au taux maximal, tandis que la bordure médiane subit la variation de surface minimale [B.5, C.3].

La mécanique géophysique de la fixation des points extrêmes aux pôles

Pourquoi sont-ce précisément les régions des pôles géographiques qui fonctionnent comme les principaux centres de convergence et de resserrement de la croûte [C.3] ?

Distribution des contraintes de traction différentielles sur la surface

Pendant la phase de déformation, la croûte des régions polaires était la plus plate, la plus étirée et la plus élargie du système [C.3]. À l’inverse, l’étirement de marée de la Terre (effet de Roche) a étiré le corps à l’extrême le long de l’axe (la ligne radiale), sa partie proche étant tenue à l’extrémité de la structure en poire [C.3]. En conséquence, la surface aux deux extrémités de la ligne de la Terre (les points subterrestre et anti-terrestre) a été resserrée et était moins étirée que tout autre secteur [C.3].

Par conséquent, au moment de la libération des contraintes et du relâchement du corps :

L’équation du bilan géométrique : le facteur de densité combiné

Le resserrement mécanique ne se limite pas aux extrémités des pôles ; il affecte un anneau global complet situé à une distance angulaire de de l’axe de la ligne de la Terre (l’anneau polaire-tangentiel qui passe par les deux pôles et les deux faces orbitales) [C.3]. La raison en est que les deux extrémités de l’axe de marée, en s’effondrant vers l’intérieur pendant l’arrondissement, poussent et envoient leur part de surface latéralement, vers l’anneau médian [C.3].

Voici un calibrage mathématique des facteurs de variation de surface bidirectionnels, en utilisant l’ellipsoïde de Roche (C.2) comme modèle d’illustration [C.2, C.3] :

La position spatiale sur la croûteFacteur de variation selon l’axe de la ligne de la Terre ()Facteur de variation selon l’axe perpendiculaire (/)Le facteur de densité unifié ()La signification géodynamique dans la croûte
Les régions polaires (apex de l’axe)1.38× (contraction de marée)0.88× (étirement latéral)1.21×Pic de resserrement général [C.3].
Les faces (avant et arrière)1.38× (contraction de marée)0.83× (étirement polaire)1.14×Resserrement modéré le long de l’anneau [C.3].
Le point subterrestre (l’axe de marée)0.88× (élargissement de marée)0.83× (élargissement polaire)0.73×Dilution maximale (la vallée de densité) [C.3].

Confirmation de la gradation à partir des données propres des six axes

Voici les données empiriques restituées à partir du catalogue Robbins dans les calottes de , après application de la procédure de « comptage propre » (retrait complet des maria et du bassin SPA des deux côtés de l’équation), qui confirment précisément les prédictions de l’équation du bilan géométrique [A.9, C.3] :

Extrémité de l’axe géométriqueDensité des cratères ( par 10⁶ km²)Densité des cratères ( par 10⁶ km²)Le statut mécanique du secteur dans la croûte lunaire
Le point subterrestre9432Base du minimum ; site du disque de contact de la poire [C.3].
Le point opposé à la Terre (anti-terrestre)20460Se situe sur l’anneau ; fonctionne comme le sommet du dôme libre [C.3].
La face avant ()16258Se situe sur l’anneau ; neutralisée du biais du bassin Orientale [C.3].
Le pôle Sud ()24982Se situe sur l’anneau ; en tête morphométriquement (sans le bassin SPA) [C.3].
Le pôle Nord ()22787Se situe sur l’anneau ; en tête morphométriquement [C.3].
La face traînante (arrière, )25384Se situe sur l’anneau ; enrichie par la géométrie du matériau de retour et le retard de rotation [C.3].

Les résultats empiriques sont fermes : chaque point situé le long de l’anneau polaire présente un excédent de densité significatif par rapport à la vallée de densité du point subterrestre, qui contient moins de la moitié du nombre de cratères des pôles [C.3]. L’asymétrie interne entre la face arrière et la face avant s’explique par la distribution des ouvertures d’éjection locales et par le flux de quantité de mouvement du matériau de retour qui a pris du retard sur la vitesse de rotation et a atterri à l’arrière (A.9) [C.3]. Les régions polaires constituent le site géométrique où se produit le maximum de resserrement et où le phénomène est mesuré sous sa forme la plus vierge et la plus propre, d’où le nom des chapitres de l’article [C.3].

C.3 Le calcul mathématique de la contraction de surface différentielle et les prédictions morphométriques

Cette section constitue le cœur analytique strict de l’article [C.3]. Elle traduit les changements de forme géométriques du corps en équations de transformation de la surface, et calcule précisément les rapports de resserrement différentiels entre les différentes zones de la croûte [C.3]. Le calcul repose sur une séparation physique absolue entre deux composantes d’effondrement autonomes :

1. Le changement de forme (Anisotropic Deformation) : calculé sous la contrainte stricte d’un volume constant vers une sphère équivalente de rayon [C.3].

2. Le changement de taille isotrope (Isotropic Volume Contraction) : la Lune a subi une contraction thermique globale et une libération de l’étirement qui ont réduit son volume d’un facteur d’environ 1.5 pendant l’arrondissement [C.3]. Comme cette contraction commune multiplie et réduit chaque cellule de surface du corps par le même facteur constant, elle s’annule mathématiquement et ne modifie aucun rapport de densité [B.6, C.3]. Par conséquent, l’indicateur de resserrement des pôles par rapport à l’équateur est parfaitement insensible au taux de contraction général du corps [C.3].

L’équation de transformation ellipsoïdale homogène :

Le passage géométrique est représenté comme une déformation élastique homogène qui transforme la forme de l’ellipsoïde en sphère, comme si la surface avait été dessinée sur une feuille de caoutchouc uniforme [C.3]. Chaque coordonnée le long des axes long, moyen et court est multipliée par les facteurs d’étirement ciblés : et respectivement [C.3]. Une petite cellule de surface dont le vecteur normal extérieur (perpendiculaire à la surface et de longueur unité) est défini par modifie sa surface selon le facteur géométrique combiné [C.3] :

Area\_Factor}\mathbf{=J}^{\mathbf{-1}}\mathbf{=}\left[\mathbf{n}_{\mathbf{x}}^{\mathbf{2}}\left(\frac{\mathbf{R}_{\mathbf{s}}}{\mathbf{a}}\right)^{\mathbf{2}}\mathbf{+n}_{\mathbf{y}}^{\mathbf{2}}\left(\frac{\mathbf{R}_{\mathbf{s}}}{\mathbf{b}}\right)^{\mathbf{2}}\mathbf{+n}_{\mathbf{z}}^{\mathbf{2}}\left(\frac{\mathbf{R}_{\mathbf{s}}}{\mathbf{c}}\right)^{\mathbf{2}}\right]^{\mathbf{-1/2}

Les cratères imprimés uniformément sur la surface d’origine modifient leurs densités en proportion inverse de ce facteur de surface (), car la même quantité de cratères comprimée dans une cellule de surface plus petite présentera une densité plus élevée [C.3].

Signification de la formule dans l’espace de la croûte :

L’équation établit que la surface d’une tache orientée le long d’un axe spatial est gouvernée par le produit des étirements des deux axes qui lui sont perpendiculaires [C.3]. Comme le produit des trois étirements est exactement égal à 1 (conservation du volume dans la phase de forme), la variation de surface est exactement égale à l’inverse de l’étirement le long de l’axe normal de la tache [C.3] :

Calibrage du modèle biaxial (le sphéroïde aplati)

Pour un corps aplati symétriquement, dont l’axe long et l’axe moyen sont égaux (), la condition de conservation du volume impose [C.3].

Le rapport statistique entre ces deux zones définit l’équation de resserrement absolue du pôle par rapport à l’équateur, dans laquelle le rayon de la sphère () s’annule entièrement [C.3] :

On établit une loi stricte : le resserrement théorique maximal au point du pôle par rapport à l’équateur est exactement égal au rapport d’axes équatorial/polaire () de la forme aplatie d’origine [C.3]. Un corps qui était 1.69 fois plus large que haut au pic de déformation laissera les cratères qui y étaient imprimés à une densité exactement 1.69 fois plus élevée à l’apex du pôle qu’à l’équateur [C.1, C.3]. Toutefois, lorsque les moyennes sont calculées sur les calottes finies empiriques, le rapport est plus petit, car la contraction s’atténue avec la distance angulaire au pôle [C.3].

Calibrage du modèle triaxial et accord avec les six axes

Pour un corps triaxial aplati et étiré (), la même équation différentielle établit une gradation de densités à trois niveaux : la densité monte d’un facteur aux pôles, d’un facteur sur les faces équatoriales, et d’un facteur à la base de la ligne de la Terre [C.3]. Cette gradation (pôles faces ligne de la Terre) restitue et concorde parfaitement avec la gradation empirique des six axes mesurée à la section A.9 [C.3].

Pour tester la validité du modèle, on reconstitue la forme d’origine du corps à partir des données des six axes, et celle-ci est transférée, sans modification, à la prédiction du rapport propre des calottes de contre la ceinture équatoriale () sur la croûte vierge [A.9, C.3] :

La base de données empirique examinéeRapport pôles/facesRapport pôles/extrémitésForme triaxiale reconstituée ()Rapport calottes/ceinture préditRapport calottes/ceinture mesuré (comptage propre)
Ensemble des cratères ()1.121.301.21 (intervalle : en D.1)
Cratères 1.111.531.281.39
Cratères de 60 – 100 km1.191.781.421.53
Cratères de 60 – 300 km1.281.731.481.61
Cratères de 100 – 300 km1.581.611.641.84

Déchiffrement de la forme au pic à partir des grands cratères

Pour les classes de taille anciennes (60–100 km et 60–300 km), les intervalles d’incertitude des six axes sont plus larges en raison d’échantillons restreints dans chaque calotte (seulement 13–68 cratères, qui dispersent les valeurs sur de larges plages de 1.63–2.87 pour et de 1.02–1.46 pour ) [C.3]. Par conséquent, la reconstitution de la forme physique au pic s’effectue au moyen d’une analyse analytique reposant sur une famille de formes définie selon l’équation de déformation dynamique [C.3] :

où l’indicateur représente la part de l’excédent de l’axe long que porte aussi l’axe moyen ( représente un œuf pur, et une pita pure) [C.3]. Pour chaque valeur de , on choisit la valeur du rapport d’axes qui donne le rapport propre exact de la classe 60–100 km (valeur de 1.44, le point subterrestre étant neutralisé de la ceinture des deux côtés en raison de la dilution du disque de contact en A.9.2) [C.3] :

Facteur de forme ()Configuration de la forme triaxiale ()La prédiction pour les six calottes (pôles/faces , pôles/extrémités)
0.191.95 / 1.22
0.351.72 / 1.28
0.501.61 / 1.32
0.601.55 / 1.34
0.701.51 / 1.36
1.001.42 / 1.41

Le modèle mécanique de l’horloge chronologique géométrique

L’analyse mathématique de la séquence de transformation est définie comme une horloge chronologique géométrique [C.3]. Chaque étape géométrique le long du chemin d’effondrement vers l’intérieur, du pic d’aplatissement jusqu’à la sphère, représente un moment dynamique de l’histoire, et le resserrement résiduel fixé dans la croûte dépend exclusivement du reste de l’arrondissement survenu après la formation du cratère [C.3] :

Pic d’aplatissement (1.8:1.4:1) → prédit 1.49 arrondissement intermédiaire (1.30) → prédit 1.26 tardif (1.20) → prédit 1.17

Un cratère formé au pic a absorbé l’intégralité de la séquence de contraction de surface, tandis qu’un cratère formé tardivement n’en a enregistré que l’accord final [C.3]. Comme les éjections de matériau se sont aussi produites pendant la croissance, une forme intermédiaire a été traversée deux fois (à la montée et à la descente) ; par conséquent, un grand cratère dense reflète une formation près du point maximal, et pas nécessairement une formation chronologiquement plus précoce [C.3].

Le calcul s’effectue par intégration directe sur une grille de de la fonction du facteur de densité spatiale [C.3] :

où sont les coordonnées angulaires à la surface ; la précision atteinte entre les deux méthodes de calcul analytiques est de , et le code est entièrement intégré dans le paquet de reproduction du système [C.3].

Le constat morphométrique est ferme : les cratères eux-mêmes ont conservé leur diamètre d’origine et ont été transportés comme des éléments rigides (des tuiles) sur une croûte supérieure solide ; ce qui s’est contracté et refermé de manière élastique, c’est la surface et les intervalles des joints entre eux (C.7, D.2) [C.3]. Les plus grands cratères () présentent une valeur empirique de , qui suggère une lecture de forme au-delà du pic théorique, mais qui tombe dans les limites de l’erreur de comptage large de l’échantillon [C.3].

C.3.1 Calibrage mathématique de la diminution de volume et des facteurs de surface absolus

Pour compléter le développement du mécanisme, le calcul mathématique combine les deux composantes de la transformation de la croûte : le changement de forme géométrique (à volume constant) et le changement de taille isotrope (la contraction thermique et le compactage du noyau) [C.3]. Dans la phase de transition dynamique du pic d’aplatissement () vers la libération sphérique (), le volume global du corps chute d’un facteur strict de 1.53 [C.3.1].

Cette chute réduit le rayon équivalent de la jeune Lune de 13% (de 2,180 km au pic à 1,890 km à la libération), et réduit mathématiquement la surface globale moyenne d’un quart [C.3.1]. Voici le déroulé des facteurs de variation de surface absolus (qui combinent la diminution de volume et l’arrondissement) calculés pour le centre de la famille de formes au pic () à une densité de 1,700 [C.3, C.3.1] :

Le secteur spatial examiné dans la croûteLe vecteur de variation de surface absolueSignification dynamique et mécanique de la surface
La surface moyenne du corps0.75×Contraction commune (isotrope) découlant de la diminution générale de volume [C.3.1].
La surface dans les régions polaires (axe )0.55×Maximum de contraction de surface du système (cellule foncée) [C.3.1].
La surface des faces (avant et arrière, axe )0.80×Contraction de surface modérée le long de l’anneau [C.3.1].
La surface aux extrémités de la ligne de la Terre (axe )0.97×Minimum de contraction de surface du système [C.3.1].
La surface des calottes polaires par rapport à la ceinture équatoriale0.67×La moyenne combinée des surfaces des calottes (0.58) contre la ceinture (0.86) [C.3.1].
Indicateur de densité des cratères (calottes par rapport à la ceinture)1.49×Le resserrement prédit que le test statistique lit sur le terrain [C.3.1].

Confirmation de la limite théorique de densité au pic

Le calcul absolu fournit une seconde vérification physique de la limite de densité au pic () [C.1, C.3.1] :

On souligne fermement que le test scientifique de l’article ne repose pas sur ces valeurs de facteurs absolus, qui sont sensibles à la densité théorique au pic [C.3.1]. Le test repose exclusivement sur le rapport statistique entre les deux dernières lignes (), qui est robuste et reste parfaitement constant pour toute valeur de densité au pic, car le terme de contraction isotrope commun s’en annule entièrement [C.3.1]. L’affirmation scientifique de la recherche repose sur ce facteur de rapport : les pôles se resserrent jusqu’au maximum du système et l’équateur jusqu’au minimum, sans dépendance ni besoin d’une détermination numérique exacte du taux de contraction de l’équateur en lui-même [C.3.1]. (La présence d’une résistance mécanique dans la croûte aurait partiellement freiné la déformation, produisant des facteurs de contraction plus modérés, mécanisme développé à la section E.3) [C.3.1].

L’évolution de la forme de la Lune en trois phases dynamiques

La Figure 10 résume tout le mécanisme géodynamique en une séquence de trois phases d’évolution successives, qui illustrent comment les points des cratères, imprimés uniformément sur le corps aplati, convergent et se resserrent aux deux pôles lorsque le corps tend vers une forme sphérique [C.3.1] :

איור 10

Figure 10. La forme de la Lune en trois phases : une masse rassemblée à l’extrémité de la poire, tenue sur sa base large et incurvée ; au pic d’aplatissement, en approchant de la base de 3,000 km près de la ligne de corotation — un corps à part entière, encore attaché, étiré entre la gravité de la Terre sur sa face proche et l’entraînement centrifuge de la corotation — triaxial, 𝑎:𝑏:𝑐 ≈ 1.8:1.4:1, vu ici de côté, la ligne de la Terre étant horizontale — avec les cratères répartis uniformément sur lui ; et rond, plus petit et plus dense, après que son centre de masse a franchi la ligne, les cratères formés sur le corps aplati se resserrant aux deux pôles. Les points sont illustratifs ; les chiffres sont ceux des sections C.1 et C.3. Les diamètres du corps au pic (le long de la ligne de la Terre, le long de l’orbite et à travers les pôles) sont calculés pour la forme 1.75 : 1.45 : 1 et pour la masse de la Lune à environ 1,700 kg/m³ (rayon équivalent d’environ 2,180 km), et le corps rond à 2,600 kg/m³ ; l’affirmation de l’article repose sur le rapport entre le pôle et l’équateur, qui ne dépend pas de la taille globale. (Sur le schéma : axe horizontal — la distance à la Terre, pas à l’échelle ; rouge — les pôles ; vert — la ligne entre la Terre et la Lune ; points — cratères. Phase 1, « extrémité de la poire » : masse rassemblée à l’extrémité de la poire, tenue sur la base large et incurvée. Phase 2, « pic d’aplatissement » : un corps à part entière, attaché, près de la ligne ; cratères uniformes, 𝑎:𝑏:𝑐 ≈ 1.8:1.4:1 ; étiré et gonflé, environ 1,700 kg/m³ ; 5,585 × 4,628 × 3,191 km. Phase 3, « arrondi » : au-delà de la ligne, plus petit et plus dense, pôles denses, rond au moment de la libération (enveloppe rigide), 3,778 km ; surface polaire 0.55×, calottes/ceinture 1.49.)

1. Phase 1 (« extrémité de la poire ») : une masse sialique initiale se forme à l’extrémité de la structure allongée de la Terre, et est tenue sur sa base mécanique large et incurvée (largeur de base initiale d’environ 6,000 km à une distance orbitale de ) [C.1, C.3.1].

2. Phase 2 (« pic d’aplatissement structurel ») : la masse se sépare du matériau de la poire et devient un corps autonome attaché, qui se déplace au-delà de la limite de Roche rigide () vers la ligne de corotation [C.1, C.3.1]. La base se réduit à une largeur de 3,000 km, et le corps est étiré à l’extrême entre la gravité de la Terre vers l’intérieur et la force centrifuge de la corotation vers l’extérieur [C.1, C.3.1]. Le corps prend une géométrie triaxiale précise de (calculée pour le centre de la famille ) [C.1, C.3.1]. Dans cet état étiré et gonflé, la densité du matériau s’établit à environ 1,700 et ses dimensions physiques sont fixées à [C.3.1]. La population de cratères est déposée et répartie sur cette morphologie aplatie selon une distribution parfaitement uniforme par unité de surface [C.3.1].

3. Phase 3 (« le corps arrondi à la libération ») : le centre de masse de la Lune franchit la ligne de corotation ; l’étirement se libère, et le corps s’effondre vers l’intérieur avec une diminution spontanée de volume vers une forme ronde, plus petite et plus dense [C.1, C.3.1]. La densité monte à 2,600 et la dimension du corps converge vers un diamètre sphérique de 3,778 km au moment de la libération mécanique, son enveloppe extérieure étant déjà rigide [C.1, C.3.1]. Le passage géométrique contracte la surface locale aux pôles d’un facteur 0.55×, et engendre l’excédent de densité bipolaire officiel des calottes de par rapport à la ceinture équatoriale, fixé à une valeur prédite de 1.49× [C.3.1].

Synthèse du mécanisme combiné :

Il a été démontré géodynamiquement qu’un corps aplati qui s’arrondit sous une condition de diminution de volume ultime referme et réduit sa surface polaire au taux maximal, et la ceinture équatoriale au taux minimal [C.3.1]. Les cratères crustaux imprimés sur la croûte avant la phase d’arrondissement achèvent la séquence concentrés et resserrés dans les régions polaires selon un facteur précis déterminé exclusivement par la forme géométrique d’origine (valeurs d’excédent d’environ 1.2 à 1.5 à l’intérieur des calottes comptées, selon le moment de leur formation par rapport au point maximal) [C.3.1]. Le modèle des six axes restitue et confirme la validité du phénomène empirique dans toutes les classes de taille, sans nécessiter de paramètres libres [C.3.1].

C.4 La gravité comme moteur physique, le profil de pente et le bilan d’énergie et de pression

Cette section définit le moteur dynamique qui a entraîné le processus d’arrondissement, quantifie le rendement énergétique libéré dans le système, et déchiffre ses effets sur l’intérieur du corps lunaire [C.4]. Le mécanisme prouve qu’il n’est pas nécessaire de postuler une nouvelle force extérieure : au franchissement de la ligne de corotation, le corps se libère de la compression sur la base de la poire, et c’est le champ de gravité propre de la Lune, qui s’intensifie à mesure que le volume diminue et que le matériau se comprime, qui contraint le corps à s’arrondir [C.4]. Comme l’enveloppe extérieure est rigide, cette déformation accumule une pression interne prononcée qui trouve une issue par les ouvertures géologiques clés déjà présentes dans la croûte (les maria) [C.4].

1. Le mécanisme de libération à deux composantes pendant le verrouillage de rotation

La déformation d’origine du corps était maintenue par un équilibre de forces temporaire : les contraintes de la rotation synchrone rapide, ainsi que l’étirement de marée à deux corps le long de la ligne de la Terre (effet de Roche), alors que le corps était mécaniquement pressé sur la base de la poire [C.2, C.4]. Le relâchement de la déformation et l’effondrement vers l’intérieur ont progressé par la libération séparée des deux composantes de la forme [C.4] :

a chuté de 0.28 au pic à 0.17 dans les phases tardives [C.1, C.4]. Cette chute paramétrique a contraint la forme d’équilibre théorique à descendre en conséquence de 1.41 à 1.20 [C.1, C.4]. La gravité est le moteur exclusif de cette phase : une vitesse de rotation donnée ne peut plus maintenir une configuration aplatie face au champ de gravité interne d’un corps devenu plus dense et attirant plus fortement vers l’intérieur [C.1, C.4]. La dynamique est comparable à une patineuse qui tourne bras écartés et dont le poids physique augmente soudain sans que sa vitesse de rotation augmente – elle ne pourra plus maintenir ses bras tendus au même rayon [C.4].

2. Le profil de pente du potentiel et l’écoulement crustal

Sur un corps très aplati qui s’est écarté de son champ d’équilibre hydrostatique, le vecteur du potentiel effectif change et inverse sa direction dynamique [C.4]. L’aplatissement de la structure que la rotation maintenait auparavant à faible densité () devient une surcharge excessive que la gravité accrue à densité élevée () ne peut pas porter [C.4].

En conséquence, les régions polaires deviennent géométriquement les points minimaux absolus du potentiel effectif (le fond d’une vallée de potentiel), tandis que la ceinture équatoriale se trouve en tête de la gradation [C.4]. Les matériaux crustaux élastiques et l’ensemble de la surface subissent un vecteur de poussée latéral et s’écoulent le long de cette pente énergétique – de la ceinture équatoriale directement vers l’intérieur, vers les pôles, le long de tous les méridiens globaux [C.4]. Le processus mécanique fonctionne comme un bol de miel épais qu’on a incliné puis remis à l’horizontale : le fluide visqueux s’écoulera de lui-même sous son seul poids jusqu’à retrouver une surface plane et un nouvel équilibre, sans qu’il soit nécessaire d’appliquer une pression extérieure supplémentaire [C.4].

3. Les équations du rendement énergétique et les valeurs de compression thermique

Le processus d’arrondissement et d’effondrement vers l’intérieur libère dans le système deux termes d’énergie spatiaux indépendants, qui se transforment en chaleur et alimentent thermiquement l’intérieur du corps [C.4] :

a. Le terme de changement de forme (l’énergie d’excentricité) :

L’arrondissement vers l’intérieur abaisse l’énergie gravitationnelle globale du corps du niveau de déformation de l’ellipsoïde au niveau minimal d’une sphère équivalente de même volume [C.4]. Pour un sphéroïde aplati d’excentricité (par exemple un rapport d’axes ), l’énergie de liaison gravitationnelle (le travail nécessaire au démantèlement mécanique complet de la structure, défini comme négatif) se calcule exactement selon [C.4] :

où est la masse de la Lune et le rayon de la sphère équivalente [C.4]. Pour une sphère parfaite, le facteur géométrique s’établit à () ; la forme aplatie porte un coefficient supplémentaire légèrement inférieur à 1 et dépendant exclusivement de son excentricité [C.4]. Comme la sphère est la forme physiquement la plus liée, le passage géométrique libère de l’énergie élastique, comme une pierre qu’on descend dans un puits [C.4]. Pour la forme au pic triaxiale du système (), la différence d’énergie libère exactement , valeur qui représente de l’énergie de liaison totale du corps (valeur d’environ ) [C.4]. Cette énergie, libérée sous forme de chaleur interne dans un noyau visqueux (qui se comporte mécaniquement comme un goudron épais), élève la température globale de l’intérieur d’environ 30 kelvins [C.4]. Ce facteur thermique fournit de l’énergie pour maintenir un intérieur mobile et visqueux, mais ne peut pas faire fondre une croûte solide [C.4].

b. Le terme de diminution de volume (l’énergie de compression) :

Le second terme d’énergie, dominant dans le système, découle de la convergence radiale de la masse [C.4]. L’énergie de liaison gravitationnelle d’une sphère, , s’approfondit fortement à mesure que le rayon du corps () diminue [C.4]. Avec la chute du rayon équivalent d’environ 2,180 km à environ 1,890 km, l’énergie libérée dans le système bondit à une valeur de (valeur de ) [C.4].

Ce terme de compression est sept fois supérieur au terme de changement de forme géométrique, et il libère une puissance thermique capable de chauffer l’intérieur du corps d’environ 200 degrés Celsius en moyenne globale (valeur exacte de 206 kelvins, calculée avec une capacité thermique fixe de ) [C.4]. L’énergie agit comme un ressort bandé libéré de ses contraintes qui restitue le travail d’étirement qui le maintenait ; le travail est converti directement en valeurs de pression hydrostatique et en chaleur interne élevée, qui maintiennent le noyau dans un état mobile et alimentent les ouvertures d’éjection crustales (les maria) [C.4].

4. Les contraintes de résistance de l’enveloppe et le mécanisme de l’horloge chronologique de la pression

Alors que l’intérieur visqueux du corps tend à s’arrondir conformément aux lignes de gravité, l’enveloppe extérieure solide (la croûte de la jeune Lune) s’est figée tôt et possède une rigidité élevée, nécessaire à la préservation de l’intégrité des cratères qui y sont imprimés (C.7) [C.4].

La séquence d’arrondissement et la variation de la pression dans la croûtePhase a : pic d’aplatissement structurelPhase b : arrondissement intermédiairePhase c : arrondissement tardif (libération)
Rapport d’axes théorique attribué ()1.401.20
Désaccord polaire résiduel ()Maximum du système (1.00)Indicateur décroissant le long de la courbeValeur résiduelle basse
Le profil de pression interne dans les ouvertures de la croûtePic d’intensité hydrostatiquePression intermédiaire décroissantePression minimale (tendant vers zéro)
Morphologie des cratères du matériau de retourCratères géants ()Cratères moyens ()Petits cratères ()
Taux de resserrement des cratères aux pôlesContraction maximale (rapport prédit ~1.5)Contraction intermédiaire (rapport prédit ~1.26)Contraction minimale (rapport prédit ~1.17)

L’enveloppe n’a pas besoin de se briser ou de se disloquer entièrement pour libérer ces contraintes [C.4]. L’énorme pression hydrostatique émanant de l’intérieur trouve une issue et un échappement rapides à travers le système de conduits et d’ouvertures géologiques clés déjà présents et répartis dans la croûte – le dispositif des ouvertures des maria (Article 5) [C.4]. Le mécanisme fonctionne comme une bouteille au contenu mou et à la peau rigide percée en plusieurs points : l’application d’une pression manuelle éjectera et projettera les matériaux intérieurs vers l’extérieur par les trous avec une force élevée, qui ira en diminuant à mesure que la pression se relâchera [C.4].

Les matériaux magmatiques et la roche fondue ont été projetés et éjectés avec force à travers les ouvertures des maria, ont effectué des vols balistiques (Article 6), et sont retombés sur la croûte solide en imprimant la population ancienne de cratères (le matériau de retour) [C.4].

La mécanique physique de l’équation du désaccord polaire

Le profil dynamique de la pression interne est calculé précisément (Figure 12) en fonction du vecteur de désaccord géométrique restant entre la surface polaire de l’enveloppe rigide et la forme sphérique vers laquelle tend l’intérieur :

איור 12

Figure 12. La pression interne au fil des phases d’arrondissement, par rapport à sa valeur au pic d’aplatissement (elle suit le désaccord polaire restant, 1− (𝑐/𝑎) 2/3), avec les tailles de cratères dont le resserrement polaire propre désigne chaque phase (D.2). Les étiquettes de l’axe de la figure conservent les anciens noms et l’échelle aplatie des phases : « 𝑎/𝑐) »start of rounding 1.69) est le pic d’aplatissement, environ 1.8 : 1.4 : 1, près de la ligne de corotation ; « on the pear, rising », « on the pear, densifying » et « far side crosses line » sont les phases d’arrondissement, à lire 1.40, 1.30 et 1.24.

Pressure\_Profile}\mathbf{\propto 1-}\left(\frac{\mathbf{c}}{\mathbf{a}}\right)^{\mathbf{2/3}

C’est la base physique stricte qui définit le diamètre du cratère comme horloge chronologique du système (C.5, D.2) [C.4]. Lorsque la forme sphérique parfaite est atteinte, la surface devient équipotentielle hydrostatique, le vecteur de pente s’annule, le désaccord structurel disparaît entièrement et avec lui s’effondre la pression interne dans les ouvertures, et l’écoulement crustal résiduel s’arrête entièrement dans le système [C.4].

איור 11

Figure 11. La Lune s’arrondit à l’intérieur d’une enveloppe rigide. La pression de l’intérieur mou s’échappe par les ouvertures existantes (les maria), projette du matériau vers l’extérieur, et le matériau retombe sous forme de cratères — éjection forte et grands cratères dans la phase la plus aplatie, plus faible à mesure qu’il s’arrondit, faible tard dans l’arrondissement, près de la libération. Les panneaux portent les anciennes étiquettes « Rounding on the pear » et « At release: almost round » pour ces deux phases de l’arrondissement, ainsi que l’échelle aplatie d’une version antérieure (𝑎/𝑐 1.69, 1.44, 1.20) : il faut les lire comme la forme au pic (environ 1.8 : 1.4 : 1), l’arrondissement après la libération de la partie œuf (1.40), et tard dans l’arrondissement (1.20). Schématique.

C.5 L’étirement dans le temps : dynamique des charges, atténuation de l’aplatissement et échelle de temps de l’arrondissement

Cette section ancre l’évolution morphologique et les forces géodynamiques des chapitres C.1–C.4 le long d’une échelle de temps chronologique stricte [C.5]. Elle définit, à chaque phase du système attaché, l’intensité des vecteurs d’étirement de la Terre, le mécanisme mécanique qui a préservé l’intégrité du corps, et le niveau exact de déformation triaxiale [C.5]. Ces valeurs déterminent sur quelle forme géométrique chaque cratère a été imprimé, et c’est de là que découle mathématiquement le niveau d’excédent de densité polaire résiduel qui y est mesuré aujourd’hui (l’horloge chronologique géométrique) [C.5].

1. Dynamique de l’indicateur de charge de marée () face au seuil de stabilité

Les forces de marée de la Terre agissent pour étirer le volume de la Lune le long de l’axe radial (axe ), tandis que le champ de gravité propre de la Lune tend à resserrer le matériau en sphère [C.5]. Le facteur du bilan des contraintes du système est calculé le long de la séquence au moyen de l’indicateur physique strict [C.5] :

Matrice d’évolution de l’indicateur et des états mécaniques au fil des phases d’attachement :

Voici l’analyse numérique exacte des contraintes mécaniques le long de la trajectoire de séparation, classée chronologiquement de haut en bas [C.5] :

Phase dynamique dans la séquence d’évolutionDistance orbitale (, km)Valeur calculée de l’indicateur ()Densité du matériau ()Le mécanisme mécanique de préservation et l’état de stabilité de la croûte
Phase 2 : formation de la masse à l’extrémité de la poire8,9912.011,300Pas d’équilibre libre ; le corps est maintenu et stabilisé mécaniquement par la base de la poire [C.5].
Phase du corps autonome (limite de Roche rigide)12,9940.671,300Pas d’équilibre libre ; le corps est maintenu par la base de la poire. C’est la charge mécanique la plus extrême depuis que le corps maintient une forme propre ( 7 fois supérieur au seuil de Roche fluide) [C.5].
Franchissement de la ligne de corotation par la face cachée0.45Pas d’équilibre libre ; le corps se trouve en pleine phase d’aplatissement dynamique et est stabilisé par la base [C.5].
Phase d’approche de la ligne (base de 3,000 km)Pas d’équilibre libre ; le corps atteint son pic d’aplatissement structurel, tenu et stabilisé sur sa base [C.5].
Phase de franchissement par le centre de masse (début de la séparation)16,4000.172,600Pas d’équilibre libre ; la pression sous la base s’annule, le corps s’arrondit (rapport d’axes ) et est stabilisé par son enveloppe rigide [C.5].
Phase de libération orbitale à la latitude 17,2870.142,600Pas d’équilibre libre ; libération mécanique finale d’avec la Terre ; le corps est stabilisé et maintenu par son enveloppe rigide [C.5].
Seuil de stabilité fluide théorique (Roche)20,5950.0842,600Obtention d’un équilibre libre fluide () ; le corps peut pour la première fois se maintenir sans l’aide d’une croûte rigide (la limite de Roche fluide se situe à 19,971–20,100 km) [C.2, C.5].
La branche de contrôle rejetée de l’Article 419,6600.0962,600Pas d’équilibre libre ; la distance se trouve à l’intérieur de la limite de Roche fluide pour cette densité, et la forme est donc rejetée mécaniquement [C.2, C.5].

2. La mécanique géophysique du soutien par la base et par la croûte rigide

L’analyse numérique révèle une anomalie structurelle prononcée : tout au long des phases d’attachement, la valeur de la charge de marée ( comprise entre 2.01 et 0.17) n’atteint jamais le seuil de stabilité fluide libre (0.09) avant le moment de la libération [C.5]. La jeune Lune se déplace vers l’extérieur et se libère à une distance d’environ 16,400 km (au-delà de la ligne de corotation équatoriale), profondément à l’intérieur de la limite disloquante de Roche fluide (qui exige une distance minimale d’environ 20,000 km pour cette densité) [C.2, C.5].

La survie de la structure s’explique par une chaîne de soutien mécanique en deux étapes [C.5] :

1. Phase de soutien par la base de la poire (phases d’aplatissement) : le corps ne se déchire pas et ne se disloque pas par effet de marée à la limite de Roche rigide (), malgré une charge extrême sept fois supérieure au seuil de rupture, car il repose et est fixé sur une base d’attachement large de 6,000 km reliée au matériau parent [C.5]. La largeur de la base porte et équilibre les vecteurs de contrainte radiaux et de marée que la gravité propre de la masse ne peut pas porter seule (Article 3) [C.5]. À mesure que le corps se déplace vers l’extérieur, la base se réduit de 6,000 à 3,000 km, et la charge diminue progressivement avec la distance d’arc [C.5].

2. Phase de soutien par la croûte rigide (phase de libération et d’arrondissement) : lorsque le corps se resserre et fixe une densité de 2,600 , la base de la poire se réduit et la pression sur elle s’annule [C.5]. À ce moment, l’enveloppe extérieure rigidifiée et figée porte seule les contraintes résiduelles et stabilise le corps rond face à la dislocation par effet de marée [C.2, C.5]. Le calcul dynamique des marées et les équations mécaniques de libération de l’Article 2 convergent vers la même contrainte physique d’une croûte rigide au moment du détachement [C.5].

3. Chronologie des quatre moments de transition géométriques

L’évolution morphologique, de l’extrémité de la poire jusqu’à l’état moderne, se déroule selon quatre segments dynamiques stricts [C.5] :

Analyse des courbes de diamètres et de l’horloge chronologique de la pression

L’examen des indicateurs issus des courbes de diamètres (Figure 13) fait ressortir les enseignements scientifiques suivants :

איור 13

Figure 13. Le diamètre long et le diamètre court (polaire) de la Lune, du pic d’aplatissement jusqu’à aujourd’hui, en fonction du temps (l’axe est déformé, les phases précoces étant dilatées). L’aplatissement a augmenté avant le premier point, lorsque le corps est sorti sur l’extrémité de la poire et s’est déplacé vers l’extérieur ; il n’est pas tracé. L’écart entre les deux courbes est l’aplatissement : le plus large à son pic, lorsque le corps s’est approché de la base de 3,000 km près de la ligne de corotation (1.8 : 1.4 : 1, à environ 1,700 kg/m³ ; le diamètre long est la ligne de la Terre, et le moyen, le long de l’orbite, 4,628 km, n’est pas tracé), se refermant pendant l’arrondissement après que le centre de masse a franchi la ligne — d’abord lorsque la partie œuf se libère (1.40), puis avec l’épaississement — jusqu’à environ 1.20, et fermé au moment de la libération, environ 210 km au-delà de la ligne de corotation (la ligne à 16,399 km à l’équateur, 17,287 km à 22.5°), lorsque le corps est rond à l’intérieur de son enveloppe rigide (𝑎/𝑐 ≤ 1.04). Après la libération, il ne reste que le petit renflement de marée du corps rigide (1.03 à 30,000 km, rotation de 14 heures ; 1.00 à 60,000 km, rotation de 41 heures). Les flèches indiquent les tailles de cratères dont le resserrement polaire propre (D.2) désigne chaque phase : les grands près du pic, les petits tard dans l’arrondissement, avant la libération. Les diamètres sont calculés à partir du volume de la masse de la Lune à la densité de la séquence attribuée à chaque phase et de la forme de la phase (C.1 ; au pic, 1.75 : 1.45 : 1) ; les deux courbes descendent parce que le volume diminue lorsque l’étirement se libère — le diamètre équivalent descend d’environ 4,350 km au pic à 3,778 km à la libération ; le rapport entre le pôle et l’équateur, que l’article examine, ne dépend pas de la taille globale. (Une version antérieure de la figure dessinait le corps encore aplati, 𝑎/𝑐 1.20, au moment de la libération.) (Sur le graphique : axe vertical — diamètre de la Lune en km ; axe horizontal — les phases : pic d’aplatissement près de la ligne, phases d’arrondissement au-delà de la ligne, tard dans l’arrondissement, la libération environ 210 km au-delà de la ligne, l’éloignement et aujourd’hui ; rouge — le diamètre long (équatorial), bleu foncé — le diamètre polaire ; au pic 5,585 et 3,191 km, à 1.40 4,871 et 3,479 km.)

1. Déchiffrement de l’écart d’origine entre les courbes : la distance et l’intervalle spatial entre la courbe rouge (le diamètre long équatorial le long de l’axe ) et la courbe bleue (le diamètre court polaire le long de l’axe ) représentent précisément la valeur de l’aplatissement géométrique [C.5]. Cet écart atteint son maximum du système au point du pic d’aplatissement (valeurs de diamètres de en longueur, contre à travers le pôle ; l’axe moyen de n’est pas tracé) [C.3.1, C.5]. Les courbes se referment et convergent rapidement au fil des phases d’arrondissement : valeurs de 4,871 km contre 3,479 km à la phase de 1.40, jusqu’à l’égalité et la fermeture complète en un diamètre sphérique uniforme de 3,778 km au moment de la libération mécanique [C.5].

2. Chute combinée de la ligne de base (effondrement du volume) : les deux courbes présentent une tendance générale à la baisse le long de la séquence [C.5]. Cette tendance provient du fait que le volume du corps s’est réduit fortement et que le rayon équivalent est passé d’environ 2,180 km au pic étiré et gonflé à 1,890 km à la libération dense et rigide [C.1, C.5]. Le rapport morphométrique entre le pôle et l’équateur est robuste et n’est pas affecté par la taille globale du corps [C.5].

3. L’épiphanie de l’horloge de la taille des cratères : les flèches dynamiques qui relient les groupes de diamètres empiriques aux phases de la courbe (D.2) représentent la loi de la pression hydrostatique interne dans les ouvertures [C.4, C.5]. Les cratères géants se sont formés au moment le plus aplati, où la pression dans les ouvertures des maria était à son pic extrême ; ils ont donc subi le maximum mécanique de contraction de surface pendant l’arrondissement et enregistrent aujourd’hui l’excédent de densité maximal aux pôles [C.4, C.5]. Les petits cratères ont été imprimés lors de phases tardives, lorsque la pression dans les ouvertures avait déjà décru ; ils ont subi moins de contraction résiduelle de surface, et présentent donc un excédent modéré sur le terrain [C.4, C.5]. L’ensemble du dispositif tombe précisément dans les limites de la plage prédite du comptage propre équilibré () [C.5].

1. Dynamique du resserrement différentiel : pourquoi les grands cratères sont-ils les plus denses ?

Le degré de resserrement mécanique d’un cratère à la surface de la croûte n’est pas déterminé par son âge chronologique absolu, mais constitue un enregistrement géométrique de la forme du corps et de son aplatissement dynamique au moment de l’impact et de la formation du cratère [B.7]. Les cratères ont conservé leur diamètre d’origine et ont été transportés comme des éléments intacts et rigides sur la croûte ; ce qui s’est contracté et refermé de manière élastique, c’est la surface et les intervalles des joints entre eux (C.7, D.5) [C.3]. La mécanique géométrique est définie selon l’équation du reste d’arrondissement [C.3] :

Le mécanisme de la pression dans les ouvertures et le cratère comme horloge de la forme :

Le lien causal entre le diamètre du cratère et le moment de sa formation est gouverné par le profil de pression interne dans les ouvertures des maria [C.4] :

1. Le moment du pic d’aplatissement : il représentait le moment du maximum de la pression hydrostatique interne, découlant du désaccord morphologique extrême [C.4]. La pression a poussé et projeté d’énormes quantités de matériau de retour avec une force maximale (Article 6), et c’est pourquoi la plupart des plus grands cratères du système se sont formés à ce moment [C.4].

2. Les phases d’arrondissement : à mesure que le corps convergeait vers l’intérieur, le désaccord diminuait, et la pression dans les ouvertures décroissait continûment, produisant ainsi des cratères de diamètres de plus en plus petits [C.4].

Du matériau a été projeté et est retombé tout au long des phases de transition (à la montée vers le pic comme à la descente), si bien que le rapport de densité accru des cratères enregistre précisément à quel point le corps était aplati au moment de leur impression [C.3]. La taille du cratère fonctionne comme une horloge géométrique stricte de la forme : plus il est grand, plus il s’est formé près du pic d’aplatissement, plus il a absorbé d’arrondissement, et plus il finit dense aux pôles [C.3, C.4]. Cette prédiction mécanique est confirmée précisément dans le comptage propre (D.2), où le rapport de densité monte de 1.09 dans la plus petite classe de taille (20–30 km) jusqu’à 1.53 dans la classe des cratères géants de 60–100 km.

2. L’évolution en trois phases des régions polaires

Les régions des pôles géographiques sont plus anciennes que la formation de l’aplatissement, et leur évolution dynamique se déploie selon trois phases structurelles strictes :

[Phase 1 : arrondissement vers l’extérieur] ──> [Phase 2 : phase d’aplatissement] ──> [Phase 3 : arrondissement vers l’intérieur]

(fixation des pôles et de l’axe de rotation) (croissance de la déformation, formation des grands) (contraction de la surface, resserrement des cratères)

3. La naissance chronologique du mois et l’orbite de la jeune Lune

Au moment de la libération et du détachement mécanique final d’avec la Terre, la période de rotation propre et la période orbitale de la Lune coïncident près d’une valeur d’un peu plus de 5.77 heures :

La courbe dynamique (Figure 14) commence par un profil entièrement plat (un mois court de quelques heures à proximité de la Terre), et ne monte fortement que loin à l’extérieur : 14 heures à une distance de 30,000 km 3.7 jours à une distance de 16 rayons terrestres jusqu’à se stabiliser à 27.3 jours dans son état moderne aujourd’hui [C.5]. Une fois fixé le point de libération mécanique, la durée du mois à chaque distance géographique découle directement de la loi de la gravitation newtonienne [C.5].

איור 14

Figure 14. Où est né le mois : la période orbitale de la Lune — le mois — en fonction de sa distance à la Terre, de 5.77 heures au moment de la libération (égale au jour terrestre de l’époque) jusqu’à 27.3 jours aujourd’hui, selon la loi 𝜇 = 𝐺(𝑀⊕ +𝑚) ,𝑇 = 2𝜋√𝑟3/𝜇) de Kepler, où 𝑟 est la distance entre les centres, 𝑀⊕ la masse de la Terre et 𝑚 la masse de la Lune). La courbe commence presque plate — un mois de quelques heures — et ne monte fortement que loin à l’extérieur. La courbe commence à la ligne de corotation, 16,400 km, et la libération se trouve environ 210 km au-delà. (Sur le graphique : axe vertical — le mois en jours ; axe horizontal — la distance à la Terre en km, de 16,400 (la ligne) en passant par 𝑅⊕ 16) 101,900) jusqu’à 384,400 (aujourd’hui).)

4. Chronologie de la phase d’arrondissement attaché et confirmation de la sphère rigide

Le processus d’arrondissement vers l’intérieur a progressé comme un événement long et graduel au cours de la phase attachée, et n’est pas contraint dans le temps ; le rythme de l’effondrement vers l’intérieur a été dicté par la libération de la compression mécanique sur la base et par l’épaississement du noyau, et non en fonction du relâchement des matériaux de la croûte (C.7) [C.5].

Le calendrier chronologique du modèle (la fenêtre de 600 millions d’années) :

Le modèle géodynamique aligne le déroulé temporel total d’environ 600 millions d’années depuis le début du renflement jusqu’à la libération finale (datée de 3.873 milliards d’années avant le présent, soit environ 4.5 milliards d’années avant aujourd’hui) [C.5] :

Quand la structure de la sphère a-t-elle été achevée ?

La Lune s’est détachée de la Terre entièrement ronde (rapport d’axes de 1.04 au plus), son enveloppe extérieure étant suffisamment rigide pour absorber les contraintes de rotation résiduelles (Article 2 §B.7) [C.5].

À une distance d’environ 60,000 km (moins de 10 rayons terrestres, et profondément à l’intérieur de la limite des 16 rayons à laquelle la forme se fige selon l’Article 4), le renflement dynamique s’est entièrement annulé et la Lune est devenue une sphère parfaite selon tous les indicateurs [C.5]. Le minuscule renflement fossile qu’elle conserve aujourd’hui (, comparé à une déformation au pic d’environ ) constitue une preuve géophysique univoque : l’aplatissement extrême de la phase attachée a entièrement disparu, et il n’en reste qu’un petit vestige 400 fois plus petit que son état d’origine [C.5].

La base de données combinée à l’issue de la Partie C

Le modèle géodynamique a achevé la construction de l’infrastructure théorique et mécanistique, et fixe les indicateurs centraux suivants [C.5] :

C.6 La base large : mécanisme de soutien mécanique, dynamique des accélérations et point de détachement

Cette section déchiffre l’équilibre mécanique et géodynamique de l’intérieur de la Lune pendant la phase attachée, et explique comment le corps a été préservé d’une rupture par effet de marée [C.6]. Comme établi à la section C.5, la charge de marée extérieure sur le corps a constamment dépassé la limite de stabilité hydrostatique d’un fluide libre () [C.5, C.6]. L’analyse suivante quantifie ce qui a porté cet écart énergétique le long de la trajectoire de séparation, et pourquoi le soutien mécanique a cessé au moment de la libération [C.6].

a. Le bilan global : le rôle structurel de la largeur de la base

Pendant les phases de croissance et d’aplatissement, la jeune Lune n’a pas agi comme un corps libre, mais était physiquement attachée et ancrée au matériau parent [C.6]. Cet attachement s’effectuait sur une base crustale très large, physiquement nécessaire pour absorber et porter l’excédent de contraintes de marée que le champ de gravité propre de la masse ne pouvait pas équilibrer seul [C.6].

La mécanique s’intègre le long de la séquence des phases de détachement [C.6] :

Matrice d’évolution de la charge de marée () et de la largeur de la base d’appui :

Voici les valeurs numériques exactes qui gouvernent les phases de détachement de la structure, intégrées dans la Figure 15 [C.6] :

איור 15

Figure 15. La base large : la charge de marée qu’elle a portée (𝑘, en bleu foncé, multipliée par 1000) et sa largeur (en rouge) en fonction de la distance à la Terre. La charge est la plus grande lorsque le corps maintient une forme propre (12,994 km) et diminue avec la distance ; la base se réduit de 6,000 à 3,000 km lorsque le corps s’aplatit et se déplace vers l’extérieur, et se relâche au moment de la libération : la séparation commence lorsque le centre de la Lune atteint la ligne de corotation (16,400 km), et le corps se libère environ 210 km au-delà. Le seuil fluide 𝑘 = 0.09 (pointillé) n’est atteint que près de 20,100 km (une version antérieure donnait 20,600), au-delà de la libération : le corps libéré était maintenu par son enveloppe rigide. Les points extrêmes proviennent de l’Article 3 (base de 6,000 km sur l’extrémité de la poire, 3,000 km au dernier contact) et de l’Article 2 (densité de 2,600 au moment de la libération) ; les points intermédiaires sont interpolés. (Sur le graphique : axe horizontal — la distance à la Terre en km : 8,991 extrémité de la poire, 12,994 corps à part entière, 14,100 face cachée sur la ligne, 16,400 la ligne, 20,100 la limite fluide (non atteinte) ; axe vertical — 𝑘×1000 | largeur de la base en km.)

Distance orbitale (, km)La phase dynamique dans la séquence d’évolutionDensité du matériau ()Indicateur de charge de marée () [C.5]Largeur physique de la base d’appui (km)
8,991Masse rassemblée à l’extrémité de la poire1,3002.016,000 [A.11]
12,994Corps à part entière, aplati (limite de Roche rigide)1,300Pic de charge du système — 0.676,000
La face cachée franchit la ligne de corotation0.45 (interpolé)
À l’approche de la ligne (base de 3,000 km)3,000
16,400Centre de masse sur la ligne de corotation (l’arrondissement)2,6000.17 (se relâche au moment de la libération)
La limite fluide théorique après la libération2,6000.09— (stabilisé par une croûte rigide)

Au moment où le centre de masse atteint la ligne de corotation (), la base se relâche entièrement [C.6]. Le corps détaché se libère environ 210 km au-delà de la ligne, profondément à l’intérieur de la limite fluide (qui n’est atteinte qu’à 20,100 km), l’enveloppe extérieure rigidifiée à portant désormais seule la charge des contraintes résiduelles [C.5, C.6].

b. Analyse dynamique du profil des accélérations et de la gravité propre

Pour déchiffrer de manière isolée la physique du bilan à deux corps, on peut traduire les indicateurs statiques en un système d’accélérations physiques [C.6]. L’équation d’équilibre est gouvernée par deux vecteurs opposés dont la valeur change au fil du temps [C.6] :

1. Le vecteur de gravité propre à la surface () : il subit une hausse régulière [C.6]. Comme la masse de la Lune est constante, le compactage du noyau et la diminution de volume (réduction du rayon ) rapprochent les matériaux de surface du centre et accroissent fortement l’attraction vers l’intérieur : de 0.87 dans les phases de départ de la poire, à 0.95 dans les phases intermédiaires, jusqu’à 1.37 au moment de la libération sphérique ( aujourd’hui) [C.6].

2. Le vecteur d’accélération de marée à travers le corps () : il représente la différence mécanique pure entre l’attraction de la Terre sur la surface la plus proche et son attraction sur le centre de masse de la Lune () [C.6]. Cette valeur subit une chute et une décroissance prononcées le long de la trajectoire d’éloignement [C.6].

Synthèse du profil des accélérations et du rapport direct entre elles :

Voici les indicateurs quantitatifs qui reflètent le bilan des accélérations, représentés graphiquement dans la Figure 16 [C.6] :

איור 16

Figure 16. La gravité propre contre la marée, de la masse à l’extrémité de la poire jusqu’à aujourd’hui. La gravité du corps lui-même (en bleu foncé) augmente à mesure qu’il se resserre à masse constante — 0.87 m/s² à l’état d’écume, 1.37 à 2,600 kg/m³ au moment de la libération, 1.62 aujourd’hui à 3,344 ; l’accélération de marée à travers lui (en rouge) diminue avec la distance, de trois fois la gravité propre à l’extrémité de la poire jusqu’au quart de celle-ci au moment de la libération (environ 210 km au-delà de la ligne située à 16,400 km), et ne passe sous le seuil fluide (0.135, le seuil 𝑘 = 0.09 de C.5 exprimé dans ce rapport) que près de 20,100 km, au-delà de la libération (une version antérieure donnait 20,600). Les deux courbes se croisent à la limite de Roche rigide, là où le corps maintient pour la première fois une forme propre : 1,300 kg/m³, 12,994 km (10,314 km à 2,600, Article 2). Les densités comprises entre 1,300 et 2,600 sont des étapes de la séquence. (Sur le graphique : axe vertical — accélération en m/s² ; axe horizontal — la distance à la Terre en km, de 8,991 (extrémité de la poire) en passant par 12,994 (corps à part entière), 16,400 (la ligne) et 20,100 (la limite fluide) jusqu’à 384,400 (aujourd’hui).)

Le secteur dynamique dans la séquenceDistance orbitale (, km)Densité ()Rayon du corps (, km)Gravité propre vers l’intérieur ()Accélération de marée à travers le corps ()Rapport direct des accélérations ()
L’extrémité ancienne de la poire8,9911,3002,3810.872.613.0 (la marée domine d’un facteur 3)
Limite de Roche rigide12,9941,3002,3810.870.871.0 (le point de croisement mécanique)
La face cachée sur la ligne14,1001,5002,2690.950.640.67
Le moment de la libération (le détachement)16,4002,6001,8891.370.340.25 (la gravité domine d’un facteur 4)
Secteur intermédiaire ()101,9002,6001,8891.370.0010.001 (la marée est négligeable)
L’état moderne aujourd’hui384,4003,3441,7381.620.000020

Les deux courbes d’accélération (Figure 16) se croisent et se verrouillent précisément à la limite de Roche rigide (), où le corps stabilise pour la première fois une forme autonome et où l’accélération de marée à travers lui égale précisément sa gravité propre (rapport de 1.0) [C.6]. Avant le point de croisement, la marée extérieure domine et tente de disloquer le corps ; après le croisement, la gravité propre de la Lune l’emporte et domine le système, et à partir d’une distance de 16 rayons terrestres la composante de marée devient statistiquement négligeable [C.6].

c. La précision méthodologique : différences de définition des indicateurs

Afin de préserver une cohérence mathématique stricte, la recherche explique l’écart numérique entre les deux indicateurs de marée utilisés aux chapitres C.5 et C.6, et prouve qu’il s’agit d’une différence de définition analytique et non d’un écart physique [C.6] :

1. L’indicateur de charge théorique () : défini à la section C.5 par [C.5, C.6]. C’est par rapport à cet indicateur qu’est calibré et formulé le seuil de Roche fluide classique, strictement égal à 0.09 [C.5, C.6].

2. L’indicateur de rapport direct des accélérations : défini dans la dernière colonne de ce chapitre comme la division directe de l’équation de l’accélération de marée à travers le corps () par l’équation de la gravité propre à la surface () [C.6].

Après simplification mathématique, le rayon du corps () s’annule entièrement, et l’équation donne la valeur du rapport simplifié [C.6] :

Acceleration\_Ratio}\mathbf{=}\frac{\mathbf{3}}{\mathbf{2}}\left[\frac{\mathbf{M}_{\mathbf{E}}}{\mathbf{\pi \rho r}^{\mathbf{3}}}\right]\mathbf{=1.5\times k

On établit une loi mathématique : le rapport direct entre les accélérations est constamment fois plus grand que l’indicateur de charge de la section C.5 [C.6].

Cette dépendance constante impose que les deux systèmes évoluent en synchronisation complète à toutes les étapes de la trajectoire [C.6] :

Cette cohérence (tracée dans les graphiques des Figures 15 et 16) prouve une fermeture mathématique complète du système géodynamique [C.6].

C.6.1 Chronologie de la gravité effective locale et mécanique dynamique du détachement

Alors que la première partie de la section C.6 traitait du bilan global des accélérations de l’ensemble du corps, cette partie décompose les données à une résolution locale [C.6]. La gravité effective – la force réelle que subit un grain libre posé à la surface (combinaison de la gravité propre vers l’intérieur, de l’étirement de marée et de la projection centrifuge vers l’extérieur) – n’est pas une valeur monolithique uniforme, mais varie fortement selon la coordonnée géographique où se trouve le grain [C.6] :

1. Dans les régions polaires (l’axe de rotation) : elle s’exprime par (où est la vitesse angulaire de rotation et le rayon du corps) [C.6]. Perpendiculairement à l’axe de marée, les forces de marée compriment et pressent physiquement le corps vers l’intérieur, et s’ajoutent positivement à l’attraction de la gravité propre [C.6].

2. Dans les régions des faces (l’axe du vecteur de mouvement orbital) : elle est exactement égale à la valeur de la gravité propre pure (), car dans ce secteur la composante centrifuge compense et neutralise entièrement le facteur de marée [C.6].

3. Aux deux extrémités de l’axe de marée (les points subterrestre et anti-terrestre) : elle s’exprime par [C.6]. À ces coordonnées, l’étirement de marée extérieur et la projection centrifuge de la rotation agissent simultanément et conjointement dans une direction inverse et opposée à l’attraction de la gravité propre du corps [C.6].

Matrice de la distribution spatiale de la gravité effective locale :

Voici les valeurs numériques exactes (en ) qui montrent comment les contraintes locales ont convergé le long de la trajectoire de séparation [C.6] :

Phase dynamique dans la séquence d’évolutionDistance orbitale (, km)Apex des pôles ()Côtés de l’axe de mouvement ()Extrémités de l’axe de marée ()
Corps à part entière (limite de Roche rigide)12,9941.300.870.43− (valeur nettement négative)
Franchissement de la ligne par la face cachée14,1001.270.950.01− (équilibre proche de zéro)
Le moment de la libération (le détachement mécanique)16,4001.541.370.86
L’état moderne aujourd’hui384,4001.621.621.62 (uniformité homogène)

C.6.2 La mécanique de l’effondrement de la pression inférieure sur la ligne de corotation

Le passage dynamique d’une base large à un détachement complet est gouverné par l’intersection simultanée de trois contraintes physiques sur la ligne de corotation [C.6] :

1. Équilibre des vitesses (à l’intérieur de la ligne) : à l’intérieur du rayon de la ligne de corotation, un corps tournant en synchronisation avec la Terre se déplace à une vitesse linéaire inférieure à celle d’une orbite képlérienne libre [C.6]. En conséquence, la gravité de la Terre attire la masse vers l’intérieur avec une force supérieure à la force de projection centrifuge vers l’extérieur [C.6]. Le corps est pressé et fermement fixé contre sa base, et la base presse en retour par le bas (c’est le mécanisme d’équilibre radial de l’Article 3) [C.6].

2. Le point zéro (sur la ligne) : au point de corotation, les deux vecteurs de force (la gravité vers l’intérieur et la projection vers l’extérieur) s’égalisent exactement [C.6]. La force de compression verticale sur la base tombe mathématiquement à un zéro absolu [C.6].

3. Le vecteur d’attraction inversé (au-delà de la ligne) : avec le franchissement de la ligne et l’éloignement au-delà, la vitesse du corps est supérieure à la vitesse orbitale libre, et la compression devient une attraction active dirigée vers l’extérieur, loin de la base [C.6]. La charge de marée décroît à – une charge entièrement gérée par l’enveloppe rigidifiée, tandis que la surface de contact s’est réduite à 3,000 km [C.5, C.6].

Comme la base n’a plus à porter de charge de compression, la Lune est lancée et libérée mécaniquement environ 210 km au-delà de la ligne, lorsque environ 58% de sa masse a déjà franchi la limite du secteur (Article 4) [C.6]. Avec le détachement, le flux de pression hydrostatique de la Terre qui passait par la poire jusqu’à la Lune cesse entièrement, et le chapitre des puissantes éjections magmatiques se clôt ainsi [C.6].

C.6.3 Analyse de l’équation de la compression face à l’étirement de marée résiduel

Pour quantifier le moteur mécanique qui a conduit à la libération de la partie œuf (l’axe de marée ), le tableau suivant suit l’attraction nette au centre (la gravité vers l’intérieur moins l’accélération centrifuge vers l’extérieur, mesurée au centre de masse de la Lune) le long des 1,000 derniers kilomètres de la trajectoire avant la ligne jusqu’au moment de la libération [C.6] :

Base\_Load\ (GPa)}=\frac{\text{Net\_Attraction}\times \text{Mass}_{\text{Moon}}}{\text{Area}_{\text{Base}}

Pour une base fixe large de 3,000 km (surface de ), on calcule la gradation de la compression face à l’étirement de marée à travers la surface () et la charge secondaire qu’il exerce à travers le plan médian [C.6] :

Position du centre du corps par rapport à la ligne (km)Attraction nette au centre ()Charge mécanique sur la base ()Étirement à la surface ()Charge d’étirement dans le plan médian ()Attraction nette (extrémité proche / éloignée, )
1,000−0.2933.04 (compression positive)0.5851.140.23− / 0.97
750−0.2162.240.5661.100.29− / 0.87
500−0.1421.470.5481.070.35− / 0.77
250−0.0700.720.5301.030.41− / 0.68
100−0.0280.290.5211.010.44− / 0.62
La ligne (0)00 (équilibre parfait)0.5141.000.47− / 0.59
100+0.027−0.28−0.5080.990.49− / 0.55
Libération (+210)0.057−0.59− (attraction vers l’extérieur)0.5020.980.51− / 0.52

L’échelle de comparaison des contraintes dans le système :

Pour comprendre l’intensité des forces géodynamiques agissant dans la croûte, on peut confronter les données aux valeurs limites mécaniques du corps [C.6] :

La conclusion géophysique :

Le croisement quantitatif prouve qu’à l’intérieur des limites de la ligne, la compression verticale de la Terre produit exactement l’ordre de grandeur des contraintes nécessaire () pour maintenir l’intégrité et la déformation de la forme aplatie [C.6]. Au moment du passage au-delà de la ligne, l’effondrement de la pression par le bas et sa transformation en attraction extérieure fournissent la libération mécanique et énergétique complète de l’effondrement vers l’intérieur de la partie œuf de la forme, tandis que le vecteur d’étirement statique ne change pas et ne peut pas produire cette libération [C.6].

Comme la charge de compression est parfaitement symétrique autour de la ligne de la Terre, elle a gouverné et conservé la partie œuf (l’axe de marée ), tandis que la libération de la partie pita (l’axe de rotation ) a été entraînée exclusivement par la hausse de densité et le resserrement du noyau (C.4) [C.6]. L’étirement de marée qui subsiste au-delà de la ligne dépasse la limite de stabilité fluide jusqu’à une distance de 20,100 km, et la Lune conserve son intégrité et stabilise une forme sphérique rigide () grâce à sa croûte élastique résistante [C.2, C.6].

C.6.4 Le mécanisme d’extraction du bouchon et l’arrêt de l’activité de la poire

Le modèle géodynamique établit un lien causal direct entre le détachement mécanique de la Lune et l’arrêt de l’évolution de la structure en poire terrestre [C.6] :

Avec la disparition de la pression interne motrice, le vecteur qui allongeait mécaniquement la structure de la poire a entièrement cessé, ce qui explique géophysiquement trois énigmes centrales :

La Terre a libéré sa pression au moment de l’extraction du bouchon, et le mécanisme s’est ainsi achevé [C.6].

איור 17

Figure 17. La gravité effective selon la position sur le corps qui tourne avec la Terre — pôles, faces et les deux extrémités de l’axe de marée — depuis le moment où le corps maintient pour la première fois une forme propre jusqu’à aujourd’hui. À 12,994 km, les extrémités sont négatives (la barre rouge pend sous la ligne du zéro) : le matériau y aurait été tiré vers l’extérieur sans la base ; les pôles portent une fois et demie plus. Là où la face cachée franchit la ligne (environ 14,100 km, à une densité de phase d’environ 1,500 kg/m³), les extrémités se situent près de zéro. À la libération, la dispersion s’est réduite, et aujourd’hui, à 384,400 km, elle a disparu : 1.62 m/s² partout. (Sur le graphique : bleu foncé — pôles (𝑔+𝜔2𝑅) ; vert — faces, le long de l’orbite (𝑔) ; rouge — extrémités de l’axe de marée (𝑔−3𝜔2𝑅) ; axe vertical — gravité effective en m/s² ; axe horizontal — la distance à la Terre en km : 12,994 corps à part entière, 14,100 face cachée sur la ligne, 16,400 la ligne, 384,400 aujourd’hui.)

C.7 La composition mécanique de la croûte de la jeune Lune : le modèle à deux couches

Le mécanisme géodynamique détaillé aux chapitres C.1–C.6 pose deux exigences physiques opposées qui tirent dans des directions contraires [C.7] :

1. L’exigence élastique/d’écoulement (l’intérieur) : le corps de la jeune Lune devait être suffisamment mou et visqueux pour s’arrondir et annuler une déformation triaxiale extrême d’environ 1.8 : 1.4 : 1 tout en étant attaché et en tournant rapidement [C.7].

2. L’exigence mécanique rigide (la croûte) : la surface devait être suffisamment rigide et solide pour conserver les cratères d’impact entièrement intacts – leur bordure, leurs marges et leurs diamètres d’origine – tout au long du processus de transformation [C.7].

La solution scientifique à ce paradoxe mécanique est mise en œuvre au moyen du modèle à deux couches (Two-Layer Rheological Model) : la division de la structure du corps en deux couches distinctes aux propriétés de viscosité entièrement différentes, l’argument géodynamique de l’arrondissement étant l’œuvre exclusive de la couche intérieure [C.7].

Le rythme de relaxation géodynamique et la fenêtre de temps disponible :

Le cycle de vie de la phase d’attachement global s’étend sur une large fenêtre de temps d’environ 600 millions d’années (depuis le début du renflement il y a environ 4.5 milliards d’années jusqu’au moment de la libération mécanique finale il y a 3.873 milliards d’années) [C.5, C.7]. Dans cette séquence, la phase d’aplatissement maximal a duré environ 100 millions d’années, et la phase d’arrondissement vers l’intérieur (de l’effondrement de l’aplatissement jusqu’au détachement) a duré encore environ 100 millions d’années [C.5, C.7].

La vitesse physique à laquelle un corps planétaire peut s’effondrer et ajuster la morphologie de sa surface est gouvernée par son indicateur de viscosité dynamique (), exprimé en pascals-secondes () [C.7]. Le temps de relaxation mécanique () – la durée nécessaire à un corps maintenu par sa propre gravité pour converger vers sa nouvelle forme d’équilibre hydrostatique après un changement de contraintes – se calcule exactement selon l’équation limite [C.7] :

Pour les paramètres du corps au moment de la libération (densité , rayon , gravité de surface ), on obtient la matrice suivante des temps de relaxation [C.1, C.7] :

Viscosité du matériau intérieur ()Temps de relaxation mécanique estiméRapport à la fenêtre d’arrondissement vers l’intérieur ( ans)
(niveau du manteau terrestre moderne)Environ 4,700 ansSuivi dynamique immédiat ; le corps s’équilibre immédiatement en parallèle au resserrement [C.7].
Environ 47,000 ansSuivi dynamique immédiat sans retard structurel [C.7].
(indicateur de l’enveloppe de l’Article 3)Environ 470,000 ansSuivi cohérent avec un retard chronologique court et négligeable [C.7].
Environ 4.7 millions d’annéesSuivi continu avec un retard modéré de quelques millions d’années [C.7].
Environ 47 millions d’annéesLimite stricte ; la forme géométrique suit à peine le rythme de l’effondrement [C.7].

La lecture de la matrice prouve que chaque multiplication par 10 du facteur de viscosité ralentit le rythme de relaxation d’un facteur 10 [C.7]. Comme la fenêtre d’arrondissement vers l’intérieur a duré environ 100 millions d’années, toute valeur de viscosité du noyau inférieure à la limite stricte de (et jusqu’à la plage de ) peut achever l’arrondissement en une infime fraction du temps disponible, et suivre étroitement les variations des contraintes gravitationnelles [C.7].

La mécanique rhéologique des deux couches

La séparation des propriétés mécaniques du corps est définie précisément selon les caractéristiques suivantes [C.7] :

a. La couche intérieure (le noyau visqueux) :

Elle se caractérise comme une zone chaude, mobile et partiellement fondue, aux valeurs de viscosité basses de [C.7]. Cette faible viscosité répond immédiatement au champ de gravité propre qui se renforce (C.4), et fait s’écouler l’intérieur et le fait se rétracter en plein accord vers la forme d’une sphère parfaite en 1,000 à 10,000 ans après toute variation de contrainte [C.4, C.7].

b. La couche extérieure (l’enveloppe rigide) :

La croûte de la jeune Lune se caractérise par une viscosité élevée et exceptionnelle de , physiquement nécessaire à la préservation de l’intégrité des cratères qui y sont gravés (Article 3) [C.7]. Comme l’enveloppe est trop rigide pour s’écouler librement, elle ne peut pas suivre l’effondrement de l’intérieur fluide, et est portée passivement par celui-ci [C.7].

En conséquence, la croûte est contrainte de s’adapter physiquement à la surface qui se réduit et se referme sous elle, de manière non uniforme dans l’espace [C.7] :

La conservation de l’enregistrement d’origine des cratères

Le modèle à deux couches résout l’énigme de la survie de la morphologie des cratères et explique pourquoi l’enregistrement est lisible et normalisé aujourd’hui [C.7] :

L’effacement ou le brouillage morphologique exige un écoulement local rapide ou une fusion de la croûte, éléments qui n’existent pas au cours de la lente séquence d’arrondissement [C.7]. Les cratères polaires sont restés à la fois intacts et denses parce que l’enveloppe était entièrement rigide tandis que le corps sous-jacent s’écoulait librement [C.7].

Prédictions tectoniques et test de l’épaisseur de la croûte par les satellites GRAIL

L’image mécanique des deux couches pose des contraintes claires qui laissent des marques géologiques vérifiables empiriquement et de manière indépendante sur le terrain [C.7] :

1. Le réseau des structures tectoniques : le modèle prédit l’existence de structures de raccourcissement et de compression (telles que des rides de chevauchement et des plis crustaux) concentrées avec une densité élevée vers les pôles et de plus en plus absentes vers la ceinture équatoriale, ainsi que l’absence totale d’une ceinture de structures d’extension (rifting) dans les régions équatoriales des terres hautes anciennes de la face cachée [C.7].

2. La résolution de la contrainte d’épaisseur de la croûte (GRAIL) : la diminution de volume et l’annulation de l’aplatissement imposent une contrainte stricte : la surface globale s’est réduite d’un quart, et aux pôles on a enregistré une chute de presque la moitié de la surface d’origine [C.3.1, C.7]. Si l’enveloppe avait conservé son volume de manière passive, la compression de la surface aux pôles aurait contraint les matériaux de la croûte à s’épaissir d’environ un tiers au total, produisant une croûte polaire 1.5 fois plus épaisse que celle de la ceinture équatoriale [C.7].

Les mesures des satellites GRAIL ont réfuté ce scénario, et ont établi que la croûte polaire n’est pas plus épaisse que celle de la ceinture de latitude où elle se trouve [C.7]. Par conséquent, le modèle pose une exigence géodynamique fermée : l’enveloppe a absorbé et traité la part de surface refermée aux pôles par une combinaison de raccourcissement mécanique, de chevauchement tectonique et de transport résiduel de matériaux crustaux fondus vers l’intérieur, dans le noyau mou sous-jacent [C.7]. La solution du calcul quantitatif du dispositif de compression et du transfert de masse est détaillée et entièrement exposée dans la Partie E.3 de l’article [C.7].

PART D

Partie D — Les quatre tests et le contrôle de l’hypothèse de base

Alors que la Partie C a établi mathématiquement et quantitativement le mécanisme d’arrondissement et de diminution de volume du corps triaxial, la Partie D confronte ce modèle géodynamique à l’enregistrement empirique des cratères sur le terrain au moyen d’une série de tests stricts [D].

La méthodologie du test (les règles derrière les équations) :

Chaque test est conduit par étapes selon un protocole académique uniforme et sobre [D] :

1. Définition des éléments de comparaison : identification des zones comparées, avec une confirmation rigoureuse de leur similitude morphologique complète pour chaque facteur, à l’exception de la seule variable testée [D].

2. La prédiction théorique et la fixation du seuil de réussite : définition de la valeur prédite par le modèle, accompagnée d’un seuil supérieur et inférieur (Pass/Fail Threshold) fixé à l’avance, avant l’exécution des comptages sur le catalogue et les échantillons datés (ce qui a entièrement évité tout ajustement artificiel des seuils aux réponses du terrain) [D].

3. La mesure empirique sur le terrain et la dérivation du résultat [D].

La procédure « des deux côtés » du comptage propre :

Le passage initial (Partie A) a traité l’intégralité du catalogue brut ; dès qu’il est apparu que les fonds basaltiques des maria abaissent la densité dans la ceinture équatoriale, l’ensemble des tests a été refondé selon le format du comptage propre [D]. Selon cette méthode, les surfaces des maria et du bassin SPA sont activement et équitablement retirées des deux côtés de l’équation simultanément (règle « des deux côtés »), ce qui garantit que les deux zones de mesure (les pôles et l’équateur) ne diffèrent que par le paramètre de latitude géographique, et sont normalisées sur le même type de croûte vierge (Highlands) [D].

Calibrage mathématique des indicateurs d’incertitude :

Chaque rapport empirique est rapporté avec un intervalle de confiance de , à l’exception des rapports de densité propres (D.1, D.2), écrits avec la notation où la valeur représente un écart-type de comptage ( ; l’intervalle de est environ deux fois plus large) [D]. L’intervalle de est défini mathématiquement de sorte que, si l’ensemble du dispositif de mesure était reproduit de nombreuses fois, il contiendrait la valeur réelle dans 95 cas sur 100 ; par conséquent, tout rapport dont l’intervalle de ne contient ni ne franchit la valeur de symétrie uniforme (1.00) est défini comme une différence physique significative située au-delà du domaine du hasard [D].

La mesure est effectuée au moyen de deux méthodes statistiques distinctes :

Synthèse de la matrice des cinq tests de la Partie D

Voici la synthèse quantitative et sobre des résultats des tests face aux seuils théoriques définis à l’avance [D] :

Code du testLe sujet statistique examinéLa prédiction théorique du modèle LEFLe seuil fixé à l’avanceLa mesure empirique effective (intervalle de confiance de )Le résultat
D.1Le rapport de densité bipolaireEnviron 1.2–1.5 (croissant en fonction du diamètre du cratère) ; valeur combinée prédite de 1.21 pour l’ensemble de la population () [C.3, C.5].1.10 – 1.60 pour les cratères (propre ; intervalle : 1.10–1.30) [A.11].Réussi entièrement [D]
D.2La gradation des diamètres (l’horloge géométrique)Les cratères conservent leur diamètre ; la population polaire est plus riche en grands cratères [C.3, C.4].Médiane polaire 1.06 (intervalle : 1.01–1.12) pour les cratères (comptage propre).Réussi entièrement [D]
D.3La population des cratères jeunesAbsence d’excédent polaire ; la population a été imprimée sur un corps rond et rigide après l’arrondissement [C.3].0.92 (intervalle : 0.39–2.04) ; et 1.12 (intervalle : 0.60–2.03).Cohérent ; échantillon trop restreint pour trancher.
D.4L’équilibre de la symétrie (nord/sud)Les deux hémisphères présentent des valeurs égales et équilibrées après nettoyage des biais géologiques [A.8.4].Absence de différence inexpliquéeLe rapport sud/nord s’établit à 1.02 (intervalle : 0.78–1.32) dans la calotte de [A.8.4].Réussi entièrement [D]
D.5Contrôle de la conservation de la taille des cratèresLes cratères ont conservé leur diamètre d’origine dans la croûte (tuiles rigides) [C.3, C.7].Conservation de la taille dans les testsLes plus grands par surface : 1.16 (1.04–1.29) ; population anciens/frais : 1.02 (0.89–1.16).Réussi entièrement (4 sur 5) [D]

D.1 Le test du rapport de densité bipolaire normalisé

La procédure des étapes de nettoyage du comptage :

Le tableau suivant quantifie le mouvement des valeurs en fonction de l’application des masques de contrôle spatiaux [D.1] :

Étape de nettoyage du comptage appliquée au système (cratères )La valeur du rapport statistique obtenueLa signification géodynamique et méthodologique de l’étape
Étape a : le catalogue brut entier (maria inclus)1.32Les maria jeunes et vides concentrés à l’équateur abaissent artificiellement la densité du dénominateur et gonflent le rapport [D.1].
Étape b : retrait des seuls polygones de mariaLe nettoyage des fonds basaltiques résiduels (retrait de surface vide de l’équateur) abaisse le rapport, et révèle la croûte montagneuse réelle [D.1].
Étape c : retrait des maria et du système du bassin SPA (le comptage propre)La comparaison propre et finale à l’identique ; le nettoyage du fond pauvre du bassin SPA au sud relève le numérateur et fixe la valeur de l’anomalie crustale [D.1].

Le résultat résiduel final (, avec un intervalle de confiance de compris entre 1.10 et 1.30) tombe avec une précision absolue au centre du domaine théorique prédit, et coïncide de manière pleinement cohérente avec la valeur de 1.21 du modèle des six axes, sans nécessiter de paramètres libres [D.1]. Le test est pleinement réussi [D.1].

Rejet des prédictions des modèles de pluie d’impacteurs extérieurs (Impactor Flux Models)

Une objection classique parmi les spécialistes planétaires propose que l’excédent de densité aux pôles pourrait provenir d’un écart de flux probabiliste, les hautes latitudes ayant reçu un volume plus important d’impacteurs (Bolides) venus de l’espace [D.1]. Les calculs physiques et dynamiques admis dans la communauté pour la répartition du flux des corps extérieurs provenant de l’extérieur du système Terre-Lune réfutent ce scénario, et donnent un signe mécanique entièrement inverse [D.1] :

Le comptage propre empirique () se situe au-dessus de ces prédictions de flux spatiaux avec un écart énorme de 6 à 8 écarts-types de comptage () [D.1]. Comme l’erreur de mesure maximale s’établit à , un écart de exclut entièrement la probabilité d’une fluctuation aléatoire ou d’un accord avec des flux spatiaux extérieurs [D.1]. Dans le cadre du modèle LEF, la densité initiale des impacts était parfaitement homogène sur la morphologie plate du renflement attaché (les pôles anciens ont reçu un volume de matériau de retour identique à tout autre secteur de la croûte), et c’est l’effondrement dynamique de l’arrondissement vers l’intérieur qui a mécaniquement resserré, contracté et concentré l’enregistrement ; aucune pluie d’impacteurs extérieurs ne peut produire une telle anomalie bipolaire [D.1].

La limite du secteur des 20 km (raison de l’omission des petits cratères) :

Sous un seuil de diamètre strict de 20 km, la procédure de comparaison spatiale est entièrement arrêtée et n’est pas appliquée [D.1]. La raison en est documentée dans la mécanique de la surface : à ces petits diamètres, l’enregistrement des comptages est entièrement dominé et faussé par trois facteurs étrangers – des amas de cratères secondaires (impacts de débris locaux produits par la projection de matériau d’un grand cratère), des champs de projections fines, et un état de saturation crustale (équilibre destructeur dans lequel chaque nouveau petit cratère efface physiquement un ancien cratère sous-jacent, et le comptage global s’arrête et ne peut plus croître) [D.1].

Aucun de ces éléments (secondaires, projections, saturation) ne conserve ni n’enregistre la position géométrique et l’état de la surface au moment de l’impact planétaire d’origine ; ils ne peuvent donc pas porter la signature mécanique du processus d’arrondissement [D.1]. Une comparaison de champs de projections secondaires avec des cratères primaires nuit à la précision scientifique, et a donc été omise [D.1].

D.2 Le test de la gradation des diamètres et l’horloge géométrique

Cette section soumet les tailles des cratères (les diamètres) à deux questions structurelles strictes [D.2] :

1. Le test de l’horloge : l’excédent de densité polaire augmente-t-il en fonction de la taille et du diamètre du cratère, comme l’impose le modèle d’arrondissement dans lequel les plus grands cratères se sont formés près du pic d’aplatissement [C.4, D.2] ?

2. Le test de la rigidité de la taille : les cratères ont-ils conservé leurs diamètres d’origine pendant que la surface et les intervalles des joints entre eux se contractaient et se refermaient élastiquement [C.7, D.2] ?

a. Le rapport de densité selon le diamètre du cratère (calibrage de l’horloge géométrique) :

Dans le cadre du modèle géodynamique, le diamètre d’un cratère reflète la phase d’évolution au cours de laquelle il a été imprimé sur la croûte [C.4, D.2]. Les cratères géants se sont formés au pic de l’activité d’éjection des ouvertures des maria, lorsque le corps était le plus aplati (), et ils ont donc absorbé l’intégralité de la séquence de contraction de surface [B.7, C.4]. À mesure que le profil de pression dans les ouvertures décroissait, la population s’est accumulée vers des tailles de plus en plus petites, qui n’ont enregistré que l’accord final de l’arrondissement [C.4, D.2].

Pour tester cette prédiction, le comptage propre (retrait des maria et du système SPA des deux côtés) est divisé en classes de diamètre isolées, et exécuté à rebours à travers l’équation de transformation ellipsoïdale (C.3) afin de reconstituer le rapport d’axes (, l’axe polaire étant fixé à 1) et la phase dynamique au cours de laquelle la population a été imprimée [C.3, D.2] :

La classe de diamètre des cratères examinéeRapport de densité propre mesuré ()Reconstitution de la forme d’origine du corps ()La phase géodynamique estimée dans la séquence (Figure 13)
20 – 30 kmTard dans l’arrondissement, juste dans les phases de libération mécanique finale [D.2].
30 – 40 kmTard dans l’arrondissement, près du point de détachement [D.2].
40 – 60 kmEn pleine phase d’effondrement de l’arrondissement vers l’intérieur [D.2].
60 – 100 kmAu pic d’aplatissement structurel (conforme à la forme ) [C.3, D.2].
100 – 300 kmAu pic de déformation extrême, à l’intérieur des marges d’erreur [D.2].

L’ordre et le domaine global (entre 1.00 après la libération et 1.50 au pic) sont ce que le modèle d’arrondissement prédit et exige physiquement, et ils sont restitués sur le terrain en plein accord [D.2]. La classe 100–300 km présente des valeurs légèrement au-delà du pic théorique (), mais en raison d’un échantillon restreint (seulement 46 cratères aux pôles et 50 à l’équateur), son incertitude statistique est large () et englobe la forme au pic [D.2].

b. Preuve de la rigidité de la taille des cratères (non-contraction du diamètre) :

Le modèle établit que les cratères ont été transportés comme des éléments intacts et rigides sur une croûte supérieure solide (C.7) ; ils ont conservé leur diamètre d’origine, et la contraction de la surface polaire a été absorbée et refermée dans les intervalles des joints entre eux [C.7, D.2]. La dynamique est comparable à une foule se tenant sur une place urbaine dont la surface se réduit : les gens se rapprochent et se serrent les uns contre les autres, mais personne ne devient plus petit ni plus menu dans ses dimensions physiques [D.2].

Cette prédiction géométrique est prouvée au moyen de deux résultats fermes [D.2] :

1. Confirmation de la significativité de l’excédent à seuils de diamètre fixes (rejet de la contraction du cratère) : la distribution cumulée des tailles des cratères lunaires dans la plage de 20 à 100 km décroît selon une loi de puissance constante : , le dispositif de mesure actuel fixant un exposant de dans les calottes et de dans la ceinture équatoriale [D.2]. Cela signifie que diviser par deux le diamètre du cratère multiplie la population cumulée par 4 [D.2]. Si les cratères s’étaient contractés physiquement avec la surface d’un facteur linéaire , la cellule de surface aurait gagné un excédent de densité de , mais aurait en même temps subi une perte d’un facteur dans le numérateur, car les cratères proches de la limite du seuil se seraient contractés en dessous et seraient sortis du calcul [D.2]. Ces deux effets opposés se seraient compensés et annulés exactement, et aucun excédent de densité bipolaire n’aurait été mesuré sur le terrain [D.2]. L’existence même d’un excédent de densité significatif mesuré à seuils de diamètre fixes constitue la preuve mécanique exclusive que les cratères ont conservé leur diamètre d’origine et que c’est la surface entre eux qui s’est contractée [D.2].

2. Le test de la médiane du diamètre linéaire : comme les grands cratères sont concentrés au pic et sont les plus resserrés, la population polaire doit présenter un enrichissement structurel en grands cratères, qui rend sa médiane de diamètre plus grande que celle de la ceinture équatoriale au même seuil de mesure [D.2].

Voici les valeurs empiriques des médianes de diamètre, calculées sur le comptage propre au moyen d’une matrice Bootstrap de 20,000 rééchantillonnages [D.2] :

Le seuil strict de diamètre de mesureLes calottes polaires : échantillon () et médianeLa ceinture équatoriale () : échantillon et médianeLe rapport des médianes mesuré ()Intervalle de confiance normalisé ()
Seuil , médiane = 32.5 km, médiane = 30.5 km1.061.01 – 1.12 (exclut l’égalité absolue)
Seuil , médiane = 59.2 km, médiane = 56.0 km1.060.98 – 1.15

Aux deux seuils de mesure distincts, la médiane du diamètre polaire présente une valeur nettement supérieure à la médiane de la ceinture équatoriale (valeurs de 32.5 contre 30.5 au seuil de 20), comme l’exige le mécanisme de l’horloge de la taille [D.2]. Au seuil de 20 km, l’intervalle de confiance () ne contient pas la valeur de symétrie 1.00, et exclut ainsi entièrement la probabilité d’une égalité aléatoire [D.2]. Au seuil de 40 km, le rapport reste fermement à 1.06, mais en raison de la réduction de l’échantillon d’un facteur trois, l’intervalle de confiance s’élargit et franchit légèrement le 1 [D.2]. Le test est pleinement réussi [D.2].

D.3 Le test de la population des cratères jeunes (The Post-Release Control Test)

Analyse des deux échantillons indépendants :

La datation des cratères jeunes est complexe, et pour éviter les biais, la recherche combine deux échantillons distincts, construits selon des méthodologies entièrement différentes, qui se vérifient mutuellement [D.3] :

1. Échantillon a — la base de données copernicienne (LPI) : fondé sur la Lunar Impact Crater Database (Losiak et al., 2015) [D.3]. L’échantillon isole 80 cratères nommés attribués géologiquement à l’ère copernicienne (la période la plus jeune du calendrier lunaire, caractérisée par des systèmes de rayons brillants et nets comme le cratère Copernicus), dont 65 cratères de diamètre [D.3]. Cet échantillon compilé porte un biais intrinsèque en faveur de la face proche en raison de la nature de l’identification historique [D.3].

2. Échantillon b — l’échantillon thermique des blocs (Mazrouei et al., 2019) : il isole 111 cratères de diamètre dont l’âge physique est strictement inférieur à un milliard d’années () [D.3]. La datation est effectuée au moyen d’une méthode physique exempte de biais d’observation : la mesure de l’abondance des blocs et des roches exposés autour de la bordure au moyen de l’instrument thermique Diviner à bord de la sonde LRO [D.3]. Un cratère frais est parsemé de blocs rocheux qui émettent un rayonnement thermique distinct, et qui s’érodent en poussière fine au fil du temps à un rythme connu [D.3]. L’échantillon garantit une couverture homogène et uniforme sur tout le corps dans la plage de latitudes comprise entre et (les surfaces des calottes et de la ceinture sont coupées en conséquence à la ligne de afin de rester compatibles avec les limites de l’instrument) [D.3].

Le résultat et la décision statistique :

Les deux échantillons indépendants convergent de manière cohérente autour de la valeur de symétrie uniforme ( contre ), et s’inscrivent précisément dans le seuil de réussite prédit [D.3]. Les résultats confirment les données de contrôle : la population des cratères jeunes ne présente pas l’excédent de densité bipolaire extrême présent dans la population ancienne [D.3].

Toutefois, en raison de la taille restreinte des échantillons (quelques dizaines de cratères pour tout le globe), les intervalles de confiance de sont très larges ( et ) et franchissent le domaine de dispersion aléatoire [D.3]. Par conséquent, les données sont pleinement cohérentes avec le modèle et confirment sa direction, mais l’échantillon existant est encore trop petit pour conférer une décision statistique autonome absolue, et la recherche détaille au chapitre D.3 les dimensions exactes d’échantillon qui seront nécessaires à l’avenir pour atteindre une significativité fermée [D.3].

Synthèse des données des cratères jeunes et analyse de la probabilité statistique (D.3)

Voici les valeurs numériques exactes, les comptages () et les densités calculées (par 10⁶ km²) pour les deux échantillons de contrôle jeunes [D.3] :

Le secteur spatial examinéÉchantillon LPI copernicien () : / densitéÉchantillon Mazrouei (l’ensemble des 111 cratères) : / densitéÉchantillon Mazrouei (plus jeunes que 1 , 103 cratères) : / densité
Calotte nord 30°7 / 2.7512 / 5.3311 / 4.88
Calotte sud 30°3 / 1.186 / 2.666 / 2.66
Les deux calottes (les extrémités)10 / 1.9718 / 3.9917 / 3.77
La bande équatoriale ()21 / 2.1435 / 3.5732 / 3.26
Découpage face proche / face cachée—3.11 / 2.843.00 / 2.52
Indicateur de rapport pôle / équateur0.921.121.16

Justification de la neutralisation du filtre des maria dans la population jeune :

Contrairement à la population ancienne, pour laquelle les maria exigeaient un filtrage mécanique (parce que les fonds basaltiques ont effacé les cratères anciens sous-jacents), la population des cratères jeunes est comptée sur toute la surface globale sans masque des maria [D.3]. Pour un cratère jeune formé après la solidification du basalte, les maria constituent un sol ordinaire ; un cratère copernicien s’imprimera sur une mare avec la même facilité et la même probabilité que sur les terres hautes anciennes [D.3].

Matrice des intervalles de confiance et estimation des probabilités du hasard

Comme les échantillons jeunes sont très restreints à la résolution globale (seulement 10 et 18 cratères dans les calottes polaires), les intervalles de confiance s’élargissent considérablement [D.3]. Le tableau suivant quantifie les intervalles de ainsi que la probabilité statistique (valeur ) d’obtenir le résultat de mesure dans un scénario hypothétique où la population jeune serait réellement resserrée comme la population ancienne [D.3] :

Code de l’échantillon statistiqueRapport pôle / équateur (intervalle )Probabilité d’obtenir le rapport si le rapport réel était (propre)Probabilité d’obtenir le rapport si le rapport réel était (brut)
LPI copernicien ()0.92 (2.04 – 0.39)0.310.22
Mazrouei (l’ensemble des 111 cratères)1.12 (2.03 – 0.60)0.470.34
Mazrouei (plus jeunes que 1 )1.16 (2.15 – 0.60)0.520.39

Preuve de l’uniformité longitudinale et réfutation de l’hypothèse de conservation statique

Les données de l’échantillon thermique des blocs de Mazrouei ajoutent un élément géodynamique critique, et prouvent que la population jeune se caractérise par une uniformité complète également le long des méridiens (est contre ouest) [D.3] :

Le constat physique est ferme : le moteur qui a mécaniquement resserré les cratères anciens aux pôles – et sur la face arrière des maria (A.9) – s’est arrêté, figé et a entièrement cessé d’agir avant l’impression de la population jeune [D.3].

C’est la confirmation précise de la prédiction de la forme figée, et une réfutation absolue de l’hypothèse de conservation statique (B.1) [B.1, D.3] : si l’excédent de densité polaire ancien était le produit de conditions environnementales polaires conservant les cratères ou empêchant leur érosion plus efficacement que dans la ceinture équatoriale, ce facteur de conservation aurait dû agir simultanément et faire ressortir aussi le nombre de cratères jeunes dans ces mêmes régions [D.3]. L’absence d’excédent chez les jeunes prouve qu’il s’agit d’un effet passager de contraction de surface [D.3].

Rejet des cartes géologiques classiques (biais de l’USGS) :

La recherche rejette activement l’usage des unités de cratères coperniciens de la carte géologique unifiée (USGS Unified Geologic Map, Fortezzo et al. 2020 ; 222 polygones) [D.3]. Le traitement des données de ces polygones présente une anomalie artificielle prononcée, dans laquelle la densité de la face proche est quatre fois supérieure à celle de la face cachée (9.2 contre 2.4 par 10⁶ km²), ainsi qu’un déficit polaire apparent de 0.4–0.5 par rapport à l’équateur [D.3].

Cette variation extrême n’est pas une propriété physique de la Lune, mais une déformation méthodologique grave de la cartographie historique : la face tournée vers la Terre a été cartographiée en premier, à haute résolution et avec un niveau de détail extrême par rapport à la face cachée, si bien que la carte unifiée compilée à partir de ces feuilles englobe un excédent d’unités jeunes là où les cartographes humains ont observé avec le plus de soin [D.3]. Par conséquent, il est interdit d’utiliser des cartes géologiques comme jeux de comptage statistiques, et le catalogue géométrique uniforme de Robbins (2019) est la seule base de données valide pour la recherche [A.11, D.3].

D.4 Le test de la symétrie bipolaire (nord contre sud)

Synthèse des données de la comparaison bipolaire dans six calottes emboîtées (D.4)

Voici les indicateurs numériques exacts, les comptages (N entre parenthèses) et les densités calculées (par 10⁶ km²) obtenus dans les deux hémisphères après application du masque géologique, ainsi que les intervalles de confiance exacts (95%) des classes de taille [D.4] :

Rayon angulaire de la calotteDensité nord (≥ 20 km)Densité sud (≥ 20 km)Rapport empirique (sud/nord, ≥ 20)Intervalle de confiance (95%, ≥ 20)Rapport empirique (sud/nord, ≥ 40)Intervalle de confiance (95%, ≥ 40)
Calotte de 10°(86) 298(36) 3481.160.77 – 1.741.260.62 – 2.43
Calotte de 15°(194) 300(84) 3051.020.78 – 1.321.230.79 – 1.89
Calotte de 20°(309) 270(149) 2801.040.85 – 1.260.940.66 – 1.31
Calotte de 30°(560) 233(322) 2461.050.92 – 1.210.900.71 – 1.15
Calotte de 45°(935) 232(702) 2461.060.96 – 1.170.900.76 – 1.07
Calotte de 60°(1,515) 235(1,178) 2401.020.94 – 1.100.940.82 – 1.08

Analyse statistique et morphologique de la convergence bipolaire

Une lecture transversale des indicateurs résiduels et des courbes de densité fait ressortir les enseignements suivants :

1. Neutralisation de l’asymétrie artificielle brute : dans le comptage brut d’origine (A.6), les deux hémisphères présentaient l’apparence trompeuse d’une variation extrême : dans le cœur intérieur de 15°, le nord était nettement plus dense que le sud (300 contre 215), et la courbe nord montait fortement vers l’intérieur tandis que la courbe sud restait entièrement plate [A.6, D.4]. Cette différence avait une cause géologique claire : le bassin géant South Pole–Aitken (SPA) couvre une part décisive de la surface des calottes sud, et dans la plage étroite de 20° autour de l’apex du pôle Sud, il constitue la seule composante du masque [A.8.3, D.4]. Dès que le fond pauvre du bassin SPA est nettoyé et retiré du dénominateur, les deux densités résiduelles convergent vers des valeurs serrées et proches, et le rapport entre elles se stabilise près de la valeur de symétrie uniforme (1.00) sur toute la plage de 10° à 60° [D.4].

2. Conservation de l’équilibre statistique dans tous les intervalles de confiance : dans les six calottes emboîtées et pour les deux seuils de diamètre isolés (≥ 20 km et ≥ 40 km), les intervalles de confiance à 95% contiennent et ancrent la valeur d’égalité mécanique absolue (1.00) [D.4].

Verdict :

Le test est pleinement et fermement réussi [D]. Les résultats empiriques prouvent une répartition morphologique nette des rôles : ce que les deux pôles ont en commun (des valeurs de densité élevées et une convergence radiale continue vers l’intérieur après filtrage) appartient entièrement au mécanisme géodynamique de l’arrondissement ; à l’inverse, ce qui les séparait à l’état brut appartient exclusivement au biais de l’impact du seul bassin SPA [D.4]. Dès que le type de croûte est normalisé à l’identique des deux côtés de la comparaison, les deux pôles portent précisément le même facteur de resserrement mécanique [D.4]. (La structure interne particulière du déficit en cratères au cœur du bassin South Pole–Aitken lui-même est renvoyée à l’analyse ciblée de l’Article 8) [A.10, D.4].

D.5 Contrôle de la conservation de la taille des cratères dans la croûte : cinq tests morphométriques stricts

Cinq tests de taille et de morphométrie interrogent directement l’enregistrement des cratères sur ce qui s’est produit dans la croûte, en utilisant les données du comptage propre (retrait des maria et du système SPA des deux côtés simultanément), afin de garantir un isolement complet du paramètre de latitude [D.5] :

Matrice des résultats des cinq tests de taille dans la croûte (D.5)

Voici la synthèse des indicateurs statistiques, des seuils et des valeurs empiriques mesurées sur le terrain pour les cinq tests morphométriques [D.5] :

Code du testL’élément morphométrique mesuré dans la croûteLa prédiction : taille conservée (tuiles)La prédiction : cratère contracté ()La valeur empirique mesurée sur le terrain (intervalle à 95%)Le verdict et le résultat
(a)Les plus grands cratères à surface égale1.16 (intervalle : 1.04 – 1.29)Taille conservée [D.5]
(b)La forme de la famille globale (R-plot)Position de la bosse identique/agrandieLa bosse est décalée vers la gauche ()Échelle optimale : Taille conservée [D.5]
(c)Rapport anciens contre frais (Old/Fresh)1.000.791.02 (intervalle : 0.89 – 1.16)Taille conservée [D.5]
(d)Le profil d’ellipticité et d’orientationUniformité par rapport au reste du corpsEllipticité accrue ( à )Excédent mineur de ; uniformité complète.Taille conservée [D.5]
(e)La règle de transition morphologique interne km km14.0 – 16.4 kmNon concluant

Décomposition et analyse des cinq tests de taille et de morphométrie

(a) Les plus grands cratères à surface égale (le test du rang) :

(b) Le test de la forme de la famille (The R-plot Landmark) :

(c) Le test anciens contre frais (The Degradation Ratio) :

(d) Le test de l’ellipticité et de l’orientation tectonique :

(e) La règle à l’intérieur du cratère (Simple-to-Complex Transition Benchmark)

La mesure empirique et la contrainte méthodologique :

L’ajustement du point de rupture géométrique dans le rapport profondeur-diamètre de la population la plus fraîche a donné deux résultats distincts selon la méthode de calcul [D.5] :

Les deux méthodes de mesure analytiques n’atteignent pas de fermeture numérique entre elles dans le secteur des calottes, car les cratères de ces régions se caractérisent par une faible profondeur morphologique générale [D.5]. En outre, une limite chronologique est apparue : les cratères les plus vierges et les plus frais, qui déterminent le point de rupture dans le système, sont pour la plupart plus jeunes que le processus d’arrondissement vers l’intérieur [D.5]. Comme ils ont été imprimés après l’effondrement, ils n’ont jamais été portés par une surface en contraction, et ne peuvent donc porter la signature de la contraction ni dans un sens ni dans l’autre [D.3, D.5].

Définition de la contrainte quantitative pour le travail futur :

En raison de ces limites, ce cinquième test spécifique est défini comme non concluant au niveau de résolution actuel, et il est formulé comme une contrainte spécifique pour la prochaine recherche : effectuer un comptage ciblé des cratères présentant un pic central sous un diamètre de 15 km dans les calottes par rapport à la ceinture équatoriale (Article 12, sujets de recherche) [D.5]. Un excédent quantitatif significatif de cratères à pic à ces petits diamètres à l’intérieur des calottes polaires constituerait le signe physique absolu d’une contraction du cratère [D.5].

Verdict final du chapitre des tests empiriques

Le test morphométrique global est pleinement et fermement réussi [D].

Quatre des cinq tests statistiques de la Partie D (le test du rang, de la bosse, de l’indicateur d’érosion et de l’ellipticité) établissent à l’unanimité et avec une large significativité que les dimensions des cratères ont été entièrement conservées, et aucun test ne soutient le scénario de contraction du cratère [D.5]. Les cratères ont été imprimés et transportés sur une enveloppe solide qui a conservé sa rigidité, et le processus d’arrondissement vers l’intérieur n’a refermé élastiquement que la surface et les intervalles des joints entre eux [C.7, D.5].

La base de données combinée et la clôture de la Partie D

Avec l’achèvement des cinq sections, l’article consolide la base de données empirique suivante, qui prouve la validité du mécanisme [D] :

1. La valeur normalisée de l’anomalie : le comptage propre (retrait des biais des maria et du bassin SPA des deux côtés simultanément) fixe un excédent de densité bipolaire officiel de 1.20 ± 0.05, valeur qui tombe au centre du domaine théorique et qui est entièrement opposée au signe que donnent les modèles de flux d’impacts venus de l’espace (0.80-0.91) [D.1].

2. La cohérence de l’horloge géométrique : l’excédent de densité monte en fonction directe de la taille du cratère (de 1.09 dans la classe 20–30 km jusqu’à 1.53 dans la classe 60–100 km) [D.2]. Les plus grands cratères bien comptés restituent et projettent à rebours précisément les dimensions de la forme triaxiale au pic d’aplatissement (1.8 : 1.4 : 1) [C.3, D.2].

3. La conservation des dimensions dans la croûte : les cratères polaires sont morphométriquement plus grands (médiane de diamètre de 32.5 km contre 30.5 km dans la ceinture) et non plus petits, donnée qui exclut un effondrement du diamètre des cratères [D.2].

4. L’équilibre bipolaire résiduel : les deux hémisphères présentent des valeurs de densité semblables et parfaitement équilibrées (rapport sud/nord serré de 1.02 dans la calotte de 15°), dès que le fond du bassin SPA est retiré du dénominateur [A.8.4, D.4].

5. L’enregistrement de contrôle jeune figé : la population des cratères jeunes (cratères coperniciens et blocs thermiques), imprimée après l’arrondissement, présente une uniformité complète sans excédent polaire significatif (0.92 et 1.12), dans les limites des marges d’erreur de leurs échantillons restreints [D.3].

Les seuils de réussite et les seuils statistiques ont tous été fixés à l’aveugle avant l’exécution des comptages, et la phase empirique de l’article s’achève ainsi [D].

PART E

Partie E — Ce qui reste

Cette partie clôt et résume l’article [E]. Les Parties A–D ont cartographié la distribution des cratères, établi le mécanisme géodynamique de l’arrondissement bipolaire avec diminution de volume, et l’ont soumis à des tests de significativité stricts [E]. Le présent chapitre rassemble les résultats en trois volets : les éléments expliqués physiquement (E.1), les anomalies conservées qui ne découlent pas du modèle actuel (E.2), et les tâches de travail futures (E.3) [E]. Aucune nouvelle mesure n’est présentée dans cette partie ; tous les chiffres et données sont tirés des chapitres précédents détaillés [E].

E.1 Ce qui est expliqué : le mécanisme combiné et l’intégration des prédictions

Un seul processus géodynamique exclusif – l’arrondissement vers l’intérieur de la jeune Lune aplatie, accompagné d’un effondrement du volume et d’un resserrement du noyau, alors que le corps était attaché et tournait en synchronisation avec la Terre – fournit une réponse complète, quantitative et systémique à six observations planétaires centrales [C.4, E.1]. Toutes les prédictions ont été calculées comme un modèle aveugle avant l’exécution des comptages empiriques [D].

La séquence d’évolution :

La masse sialique s’est formée à l’extrémité de la poire terrestre (), et s’est déplacée vers l’extérieur sur une base large de 6,000 km [C.1, C.5, E.1]. Avec son éloignement, la base s’est réduite à 3,000 km et la masse a pris forme comme un corps autonome étiré à l’extrême entre la gravité de la Terre vers l’intérieur et la force centrifuge de la corotation vers l’extérieur [C.5, E.1]. Le pic de déformation a été atteint près du franchissement de la ligne de corotation (distance ), où le corps a pris une configuration triaxiale de , dans un état gonflé et à une faible densité de [C.1, C.5, E.1].

Au moment où le centre de masse a franchi la ligne de corotation (), la compression verticale sur la base s’est annulée, et la partie œuf (axe ) s’est effondrée vers l’intérieur [C.5, C.6, E.1]. Le corps s’est arrondi avec une forte diminution de volume (hausse de la densité à 2,600 sans perte ni ajout de masse), qui a réduit le rayon équivalent de 13% et ramené le rapport d’axes à 1.20 dans les phases tardives de l’arrondissement [C.1, C.3.1, E.1]. Le détachement mécanique final (la libération) s’est produit environ 210 km au-delà de la ligne de corotation, alors que la Lune était déjà passée à une forme sphérique parfaite () et que son enveloppe extérieure était rigide [C.5, E.1].

L’arrondissement qui a resserré les cratères s’est entièrement produit alors que la Lune était attachée [E.1]. Le détachement dynamique a agi comme l’extraction d’un bouchon : la Lune a quitté l’extrémité de la poire et a exposé ses ouvertures à l’espace, ce qui a provoqué l’effondrement de la pression dans la poire et l’arrêt de son allongement, a empêché l’éjection d’une deuxième lune, et a dispersé le matériau d’éjection dans toutes les directions (C.6) [C.6, E.1]. Après la libération, le ralentissement de la rotation verrouillée n’a modifié que le petit renflement élastique du corps rigide (1.03 à 1.00), sans lien avec l’aplatissement précoce qui avait disparu auparavant [C.5, E.1].

La matrice de synthèse des six observations expliquées (E.1)

Voici un regroupement sobre des valeurs des prédictions théoriques (formulées à l’aveugle avant le comptage) face aux indicateurs empiriques effectivement mesurés sur le comptage propre [E.1] :

Code de l’observationLe phénomène géologique examinéLa valeur prédite par le modèle (avant le comptage)La valeur empirique mesurée (dans le comptage propre)Section source (calcul et test)
1L’excédent de densité bipolaireRapport de 1.20–1.50 (croissant en fonction du diamètre) ; valeur combinée de 1.21 pour toute la classe . pour ; gradation continue de 1.09 1.53 avec la taille.C.3, D.1, D.2 [E.1]
2Gradation de la médiane du diamètre linéaireLa médiane polaire doit être plus grande ; les grands cratères ont absorbé l’essentiel de la contraction de surface.Rapport des médianes de 1.06 (intervalle : 1.01 – 1.12) au même seuil de mesure.D.2 [E.1]
3L’équilibre de la symétrie bipolaireÉgalité mécanique absolue (valeur 1.00) ; symétrie axiale complète autour de l’axe de rotation.Rapport sud/nord de 1.02 (intervalle : 0.78 – 1.32) dans la calotte de ; égalité dans chaque calotte.D.4 [E.1]
4La gradation globale des six axesLes pôles et les faces se situent au-dessus de la courbe de la ligne de marée ; forme combinée de qui prédit un rapport propre de 1.21.Le point subterrestre est le moins dense du système (148) ; le rapport propre à la ceinture s’établit à 1.19–1.20.A.9, C.3 [E.1]
5L’enregistrement de contrôle jeuneAbsence d’excédent de densité (valeur ) ; la population d’impacts a été imprimée sur un corps entièrement rond.Valeurs de 0.92 et 1.12 dans deux échantillons distincts ; cohérent avec le modèle, échantillon restreint.D.3 [E.1]
6Géométries du point maximalDéformation triaxiale d’environ à une densité de matériau fondu d’environ 1,700 .Reconstituée à rebours à partir de la classe 60–100 km à une valeur de (composante de rotation conforme à Maclaurin).C.1, C.3, D.2 [E.1]

Décomposition et analyse des conclusions fondamentales de la recherche

1. L’anomalie de l’excédent polaire et la variation du diamètre des cratères (horloge chronologique géométrique) :

2. Le contrôle de la rigidité des dimensions des cratères (le modèle des tuiles) :

3. L’équilibre de la symétrie axiale bipolaire :

4. La gradation des six axes et la géométrie triaxiale :

5. L’enregistrement de contrôle des cratères jeunes :

6. Reconstitution de la géométrie du point maximal à partir des cratères géants :

Synthèse de la méthodologie et fermeture rhéologique

L’intégration des parties de l’article repose sur des lois physiques strictes et sur une fermeture cohérente entre les parties de la série scientifique [E.1] :

1. L’analyse de la forme globale : la forme de l’ensemble de la population combinée est dérivée directement des données des six axes, et prédit le rapport du comptage propre sans aucun réglage ni ajustement artificiel [C.3, E.1]. La forme d’origine au pic est isolée à partir des cratères géants, et sa composante de rotation mécanique est confirmée directement par le rapport de Maclaurin (Maclaurin Spheroid) pour la vitesse de rotation de l’Article 2 et les limites de densité imposées par la diminution de volume (C.1) [C.1, E.1]. L’équation de Maclaurin constitue une loi physique classique qui définit le rapport d’axes d’un corps fluide homogène en rotation sans paramètres libres ; les vitesses de rotation et les densités ont été tirées de manière cohérente des articles précédents, et la traduction des dimensions géométriques en facteurs de resserrement de surface relève des mathématiques pures [C.1, E.1].

2. L’élément empirique appris (l’horloge de la taille) : la seule variable apprise directement des données de terrain, et non fixée comme prédiction à l’avance, est l’attribution spécifique de chaque groupe de diamètre de cratères à sa phase dynamique (D.2) [D.2, E.1]. Le mécanisme exige physiquement que les plus grands cratères se forment plus près du point maximal et absorbent le maximum de contraction de surface, mais la fixation de la gradation numérique exacte (quel diamètre reflète quel moment) est apprise des comptages empiriques sur le terrain [D.2, E.1].

3. Unification des contraintes de la croûte (l’enveloppe de C.7) : la croûte mécanique jugée essentielle à la section C.7 pour la conservation des dimensions des cratères est exactement la même croûte que celle définie dans l’Article 3 à partir d’une justification entièrement différente [C.7, E.1]. Dans l’Article 3, une croûte d’une viscosité élevée de était nécessaire pour garantir la survie morphologique des cratères imprimés alors que le corps se trouvait à proximité de la Terre ; dans le présent article, la même enveloppe rigide fonctionne comme une tuile mécanique qui empêche l’effondrement du diamètre des cratères, et contraint la contraction de surface à être absorbée dans les intervalles des joints entre eux [C.7, E.1]. Une seule exigence structurelle sert et relie deux articles indépendants de la série [E.1].

Conclusion scientifique finale :

Nous disposons maintenant d’un processus géodynamique unifié, gouverné et fixé par des valeurs numériques tirées des articles précédents et par la forme que les cratères eux-mêmes enregistrent autour des six axes [E.1]. Ce processus fournit une réponse quantitative complète et explique précisément l’excédent de densité, la gradation des tailles (l’horloge), l’équilibre de la symétrie axiale et la gradation spatiale des axes des cratères anciens, tout en restant pleinement cohérent avec l’enregistrement de contrôle des cratères jeunes uniformes [E.1].

E.2 Anomalies et résultats non expliqués par le modèle actuel

Au cours du développement de l’argument scientifique et du traitement des données empiriques, cinq observations principales ont été enregistrées qui ne découlent pas du mécanisme d’arrondissement sous sa forme actuelle et ne sont pas expliquées par lui [E.2]. Ces résultats sont isolés et détaillés ci-dessous avec les outils qui permettront de les trancher à l’avenir, afin de bloquer tout biais méthodologique et d’éviter qu’ils ne soient pris à tort pour des résultats centraux de l’article [E.2].

1. L’anomalie de l’excédent biaxial de la face arrière

2. L’anomalie de densité de la couronne de la face cachée (le point anti-terrestre)

3. L’anomalie de faible profondeur morphologique des cratères polaires

4. L’anomalie d’orientation méridienne aux latitudes moyennes

5. La structure interne des cratères du bassin South Pole–Aitken

E.3 Travail supplémentaire et tâches de recherche futures

Cette section rassemble la série d’étapes, d’analyses analytiques et d’extensions statistiques qui transformeront les composantes restantes du modèle LEF d’exigences géodynamiques strictes en calculs physiques précis et fermés [E.3]. Chaque point est défini sobrement et se termine par l’indication de l’objectif scientifique précis qu’il est censé démontrer [E.3] :

1. Le modèle de la forme mécanique dotée de résistance (The Rigid Shell Shaping Effect) :

2. Vérification élargie de la liste des cratères jeunes :

3. Résolution de l’équation de relaxation rhéologique à deux couches :

4. Isolement et traitement des deux sujets d’arrière-plan crustaux :

5. Dérivation de la contribution des contraintes de marée (la partie œuf) à partir des vecteurs de force :

6. Le test de l’espacement tangentiel absolu des cratères (Nearest-Neighbor Distance) :

7. Cartographie des structures tectoniques et des empreintes géologiques dans l’enveloppe :

8. Développement de l’équation du bilan des ouvertures d’éjection et du matériau de retour :

9. Le test de résolution extrême des calottes les plus intérieures :

10. Solution dynamique des contraintes de compression sur la base et de la partie œuf :

11. Calcul des contraintes dans l’enveloppe au-delà de la ligne de corotation :

12. Résolution du bilan de masse de la croûte disparue (GRAIL Coherence) :

13. Précision des dimensions de la forme au pic au moyen d’un catalogue de cratères élargi :

14. Reprise des cinq tests de taille sur la morphologie triaxiale :

Disponibilité du code numérique et de la base de données de la recherche

L’ensemble des calculs, des traitements de comptages, des calculs de surface des calottes et des ceintures, des solutions des équations d’équilibre du sphéroïde de Maclaurin et de l’ellipsoïde de Roche, des matrices des facteurs de relaxation, des fonctions d’ajustement triaxial aux six calottes principales, des données de diminution de volume et des courbes d’attraction nette du chapitre C.6, est entièrement reproduit de manière fermée et propre au moyen des scripts du paquet de reproduction officiel du modèle [E.3]. Le paquet est développé en langage Python 3, à l’aide des bibliothèques numpy et scipy, et est entièrement déposé avec la série scientifique (un identifiant d’enregistrement Zenodo dédié sera attribué et publié ultérieurement) [E.3].

Références

[R1] Robbins, S. J. (2019). A new global database of lunar impact craters >1–2 km: 1. Crater locations and sizes, comparisons with published databases, and global analysis. Journal of Geophysical Research: Planets, 124, 871–892. doi:10.1029/2018JE005592.

[R2] Chandrasekhar, S. (1969). Ellipsoidal Figures of Equilibrium. Yale University Press. (Maclaurin spheroids; Roche ellipsoids; the index symbols used in C.2.)

[R3] Neumann, G. A., et al. (2015). Lunar impact basins revealed by Gravity Recovery and Interior Laboratory measurements. Science Advances, 1, e1500852. doi:10.1126/sciadv.1500852. (Table 1, basins of 500 km and more, used in A.3.)

[R4] Losiak, A., et al.; revised by T. Öhman (2015). Lunar Impact Crater Database, v08Sep2015. Lunar and Planetary Institute. (Copernican-age craters used in D.3.)

[R5] Nelson, D. M., et al. (2014). Mapping lunar maria extents and lobate scarps using LROC image products. LPSC 45, abstract 2861. (LROC global mare shapefile, used in A.3.)

[R6] Mazrouei, S., Ghent, R. R., Bottke, W. F., Parker, A. H., & Gernon, T. M. (2019). Earth and Moon impact flux increased at the end of the Paleozoic. Science, 363, 253–257. (Second young-crater sample of D.3; the authors’ input file on GitHub was used, not Table S1.)

[R7] Grier, J. A., McEwen, A. S., Lucey, P. G., Milazzo, M., & Strom, R. G. (2001). Optical maturity of ejecta from large rayed lunar craters. JGR, 106, 32847–32862; as compiled in Tai Udovicic, C. J., et al. (2021), Zenodo 4345058. (Rayed craters, D.3.)

[R8] Fortezzo, C. M., Spudis, P. D., & Harrel, S. L. (2020). Release of the digital Unified Global Geologic Map of the Moon at 1:5,000,000 scale. LPSC 51, abstract 2760. (Copernican units examined and not used, D.3.)

[R9] Le Feuvre, M., & Wieczorek, M. A. (2008). Nonuniform cratering of the terrestrial planets. Icarus, 197, 291–306. (Pole-to-equator ratio about 0.90, D.1.)

[R10] Le Feuvre, M., & Wieczorek, M. A. (2011). Nonuniform cratering of the Moon and a revised crater chronology of the inner Solar System. Icarus, 214, 1–20. (Pole-to-equator ratio about 0.80, D.1.)

[R11] Gallant, J., Gladman, B., & Ćuk, M. (2009). Current bombardment of the Earth–Moon system: emphasis on cratering asymmetries. Icarus, 202, 371–382. (Pole-to-equator ratio 0.912, D.1.)

[R12] Wang, Y., Wu, B., Xue, H., Li, X., & Ma, J. (2021). An improved global catalog of lunar impact craters (≥1 km) with 3D morphometric information and updates on global crater analysis. Journal of Geophysical Research: Planets, 126, e2020JE006728. (Catalogue LU1319373; depth-to-diameter ratios, D.5.)

Manuscrits de la série : Duran, Z., Papers 2 (LEF Mechanism in Principle), 3 (LEF Mechanism in Practice), 6 (The Craters), 4 (The Orbit) and 11 (The Requirements) of the LEF series, deposited at doi:10.5281/zenodo.22743615.

Déclarations

Auteur : Zohar Duran, chercheur indépendant, Ramla, Israël. ORCID 0009-0005-3194-9387. Aucun financement public ou commercial n’a été reçu. Une demande provisoire de brevet aux États-Unis (n° 64/116,844, déposée le 22 juillet 2026) existe pour un dispositif d’ingénierie lié au modèle.