7. Los polos

El entorno computacional (a): datos del cuerpo y dinámica orbital

A continuación se definen los parámetros físicos, las distancias y los estados dinámicos que sirven de base para el cálculo del modelo de la Luna en formación (modelo LEF) en la fase acoplada:

El entorno computacional (b): datos del registro geométrico y del recuento empírico

A continuación se definen los indicadores, las bases de datos y los supuestos geométricos básicos utilizados para calcular las áreas superficiales y convertir los recuentos de cráteres en densidades normalizadas:

Dinámica de las fases de separación: evolución morfológica y balance de esfuerzos

A continuación se expone, en orden cronológico y mecánico, la secuencia de fases de desarrollo de la Luna joven —desde la consolidación de la masa inicial en el extremo de la estructura de la pera hasta su estado actual—, definiendo el balance de fuerzas, los cambios de densidad y la transformación geométrica:

Cronología de los acontecimientos y ventana de formación de los cráteres

El modelo geodinámico define una escala temporal total de unos 600 millones de años desde el comienzo de la formación del abultamiento original hasta la fase de liberación final (datada en 3.873 mil millones de años antes del presente), dividida en tres subintervalos que constituyen la ventana morfogenética exclusiva de los cráteres de este artículo:

• Fase inicial de la pera (unos 400 millones de años): elevación de la masa siálica sobre la prolongación del material de la Tierra y crecimiento del abultamiento.
• Fase de achatamiento dinámico (unos 100 millones de años): fase de deformación triaxial extrema y alargamiento de los ejes cerca de la línea de corrotación.
• Fase de redondeo acoplado (unos 100 millones de años): fase de liberación mecánica del estiramiento, espesamiento del cuerpo, rigidización de la corteza y separación final.

Definiciones básicas e indicadores mecánicos para la lectura de la tabla

• Eje de marea $a$: el eje espacial que recorre la línea de simetría que une directamente la Luna en formación con el centro de la Tierra.
• Eje de movimiento $b$: el eje perpendicular al eje de marea, orientado a lo largo del vector de velocidad orbital del cuerpo.
• Eje polar $c$: el eje geográfico perpendicular al plano del sistema, que une los polos (el eje más corto durante el achatamiento).
• “Se sostiene por sí mismo”: la capacidad mecánica del cuerpo para mantener su integridad y una forma geométrica cerrada dentro de las ecuaciones de un cuerpo fluido, bajo el efecto de la autogravedad y de la resistencia estructural (no significa la neutralización de las fuerzas de marea externas).
• “Separación”: separación física y apertura mecánica completa de la capa de contacto y de la interfaz entre la Luna y el material de la pera terrestre.
• Confirmación física de la densidad en el máximo ($\approx 1,700{\text{ kg/m}}^3$): este valor está respaldado por dos comprobaciones dinámicas independientes:
  1. La parte de la pita del cuerpo ($b/c\approx 1.4$) coincide matemáticamente con la forma de equilibrio exacta del esferoide de Maclaurin para un fluido homogéneo que gira con un periodo de 5.5 horas a una densidad de 1,600–1,750 ${\text{kg/m}}^3$.
  2. El límite mecánico superior para la contracción de los extremos del eje de marea durante el redondeo exige físicamente una densidad en el máximo inferior a 1,800 ${\text{kg/m}}^3$.
• Límite de Roche rígido ($10,314$ y $12,994$ km): calculado según la ecuación estructural $1.26R_{\oplus }({\rho }_{\oplus }/\rho {)}^{1/3}$. Un cuerpo más ligero ($1,300{\text{ kg/m}}^3$) posee un agarre gravitatorio propio más débil, por lo que su límite de Roche se sitúa más lejos del origen ($12,994$ km), donde logra por primera vez estabilizar una forma propia mientras sigue acoplado.

Criterio de estabilidad del fluido libre (índice $k$)

La evaluación del estiramiento mecánico se realiza mediante el índice del sistema:

• Umbral de estabilidad fluida: se fija en $k<0.09$. Cuando el valor del índice es inferior a este umbral, un cuerpo fluido libre es capaz de mantener una forma cerrada y estable sin desintegrarse bajo las fuerzas de marea de la Tierra ($M_E$).
• Sensibilidad paramétrica: dado que la distancia orbital ($r$) aparece elevada a la tercera potencia, el valor del índice $k$ cae rápidamente a medida que el cuerpo se aleja. Por otra parte, como la densidad ($\rho$) está en el denominador, un cuerpo de mayor densidad hará descender el valor de $k$ por debajo del umbral crítico en una fase más temprana y más cerca del origen. Un cuerpo fluido libre con la densidad de liberación habría necesitado una distancia de unos 20,000 km para alcanzar la estabilidad, distancia que la Luna no alcanzó en la fase acoplada; de ahí que la estructura se sostuviera únicamente gracias a su capa rígida.

Introducción — la pregunta y cómo se pone a prueba

1. Metodología de la investigación: definición de las regiones de recuento y de los indicadores empíricos

Para realizar la comparación morfométrica y evitar sesgos derivados de la incompletitud en la identificación de cráteres pequeños, se fijó un umbral mínimo estricto para el diámetro de los cráteres contados. Los datos se basan en el catálogo Robbins (2019) y se analizan según las siguientes unidades geométricas:

• Índice de densidad normalizada: se define como el número total de cráteres en una región dada dividido por su área estática. El índice se expresa en cráteres por millón de km², lo que permite una comparación neutral entre celdas de área de distintas dimensiones geométricas.
• Los casquetes polares: dos casquetes geométricos de igual área, definidos por un radio angular de $30^{\circ }$ desde el centro de la Luna alrededor de los polos de rotación (perpendiculares al plano ecuatorial).
  • El casquete norte: se extiende desde la latitud $60^{\circ }\text{ N}$ hasta $90^{\circ }\text{ N}$.
  • El casquete sur: se extiende desde la latitud $60^{\circ }\text{ S}$ hasta $90^{\circ }\text{ S}$.
• El cinturón ecuatorial: una banda global definida en el intervalo de latitudes entre $15^{\circ }\text{ S}$ y $15^{\circ }\text{ N}$.
• Asignación espacial: cada cráter se cataloga y se asigna a la región de recuento correspondiente según las coordenadas de su punto central geométrico.

2. Resultados del recuento bruto frente al procedimiento de recuento limpio

El recuento empírico inicial muestra un exceso de densidad significativo en las regiones polares, que se acentúa a medida que aumenta el diámetro de los cráteres examinados. Sin embargo, el análisis geológico muestra que el recuento bruto sufre sesgos estructurales agudos que requieren una normalización matemática:

• Sesgo del cinturón ecuatorial (maria y llanuras basálticas): las áreas basálticas y los maria (el sistema de alimentación del LROC) constituyen suelos jóvenes en los que los cráteres antiguos fueron fundidos o cubiertos debido a la actividad térmica y a una corteza fundida por el calor de la Tierra. Estas áreas representan el 22% del cinturón ecuatorial, pero solo el 1% de los casquetes polares. Su inclusión en el recuento bruto rebaja artificialmente la densidad ecuatorial e infla la razón polos/ecuador.
• Sesgo de los casquetes polares (cuenca South Pole–Aitken): la cuenca SPA borró físicamente los cráteres de su interior. Dado que el 24% del área de los casquetes polares (principalmente en el sur) se encuentra dentro de la cuenca, su inclusión hace que los polos parezcan más vacíos de lo que realmente están y reduce la razón.
• Procedimiento de recuento limpio: eliminación uniforme de las áreas de los maria y de la cuenca SPA en ambos lados de la comparación. Es el único procedimiento que permite comparar superficies antiguas (Highlands) con superficies antiguas, y de él se deriva el valor estricto de la razón del sistema en el que se apoya la investigación: $1.20\pm 0.05$.

3. Supuesto dinámico de partida y mecanismo de contracción mecánica

• Supuesto de distribución uniforme del material: en el marco del modelo LEF (Artículo 6), la población antigua de cráteres se formó por la caída de material terrestre retornado sobre la corteza mientras esta se extendía como un abultamiento sobre el cuerpo progenitor (la Tierra). Para los fines del análisis se adopta como hipótesis de trabajo una distribución completamente uniforme del material por unidad de área. Este supuesto neutraliza la cuestión de la asimetría de la fuente: toda falta de uniformidad en la densidad que se observa hoy en la superficie debe proceder de procesos mecánicos posteriores que actuaron sobre la corteza que transportaba los cráteres.
• Descarte de hipótesis alternativas: durante la investigación se examinaron seis mecanismos físicos distintos como posible explicación de la diferencia de densidad; las seis explicaciones fueron descartadas y dejadas de lado.
• La séptima explicación (redondeo y disminución de volumen): la Luna joven se comportó físicamente como un cuerpo achatado y aplastado (semejante a un globo sometido a la presión de unas manos que pierde volumen gradualmente). En el máximo de achatamiento, debido a los esfuerzos de rotación simétricos alrededor del eje de rotación, el cuerpo era ancho y plano en sus extremos y estrecho en su ecuador. Los cráteres quedaron impresos sobre esta morfología de manera uniforme.
• Apiñamiento de los polos: en el momento en que el centro de masa cruzó la línea de corrotación y el estiramiento se liberó mientras el cuerpo seguía acoplado, las “manos” mecánicas “se apartaron” y el cuerpo tiende a la forma esférica con disminución de volumen. En este proceso geométrico, las partes planas y anchas convergieron hacia dentro, hacia dos pequeños casquetes polares. Puesto que el eje de rotación pasa por los polos y el achatamiento era simétrico a su alrededor, cada meridiano retrocedió hacia los dos puntos extremos, que la rotación nunca desplazó dinámicamente (C.3).
• Conservación morfométrica de los cráteres: la capa exterior de la Luna era lo bastante rígida (C.7) para transportar los cráteres como elementos íntegros y aislados. Los propios cráteres conservaron su diámetro y tamaño originales; fue la superficie elástica entre ellos la que se redujo y contrajo (la superficie se cerró en una cuarta parte del área media y en alrededor de la mitad del área de los polos), lo que produjo un aumento directo del índice de densidad en las regiones polares en comparación con la zona ecuatorial (donde el borde convergió menos).

1. Predicciones del modelo geométrico y mecanismo de apiñamiento diferencial

El grado de apiñamiento mecánico de los cráteres no es función directa únicamente de la edad cronológica, sino un registro geométrico de la forma del cuerpo y de su achatamiento dinámico en el momento de la formación del cráter. Dado que se eyectó y retornó material a lo largo de todas las fases de la transición, el nivel de deformación estructural quedó conservado en la morfología de la corteza según los siguientes principios:

• Formación en el máximo de achatamiento: los cráteres formados cuando el cuerpo se hallaba en el máximo de la deformación triaxial ($1.8:1.4:1$) fueron transportados a lo largo de toda la secuencia de redondeo. En consecuencia, experimentaron el máximo mecánico de contracción de la superficie y presentan hoy una razón de densidades predicha de 1.49 a favor de los casquetes polares.
• Formación en fases tardías del redondeo: los cráteres formados cuando la Luna ya había empezado a colapsar hacia dentro y se había redondeado hasta una forma de 1.20 (cerca de la fase de rigidización de la capa) experimentaron menos contracción superficial residual, por lo que la razón de densidades predicha para ellos cae a 1.17.
• Formación tras la liberación: los cráteres formados sobre un cuerpo que ya había completado la transición a la forma esférica no experimentan contracción superficial alguna, por lo que la razón predicha para ellos es 1.00 (distribución completamente uniforme).
• Correlación entre el diámetro del cráter y el momento de su formación: la mayoría de los cráteres más grandes se formaron cerca del máximo de la fase de eyección y retorno del material terrestre. Por ello, el modelo predice de forma coherente que la razón de densidades aumentará cuanto mayor sea el diámetro del cráter.

2. Correlación empírica: predicción del modelo frente a los indicadores del recuento limpio

A continuación se presenta una comparación estricta y normalizada entre las razones de densidades teóricas predichas (derivadas de la forma de la Luna en el momento de la formación) y los indicadores empíricos medidos efectivamente en el recuento limpio (eliminación de maria y de la cuenca SPA), por clases de diámetro:

• Análisis de tendencias e incertidumbre: los valores empíricos medidos abarcan con precisión el intervalo teórico y conservan su orden jerárquico (disminución del exceso a medida que los cráteres son más pequeños). El grupo total ($\ge 20$ km) muestra una concordancia casi completa entre la predicción mecánica de los seis ejes (1.21) y el resultado en la superficie ($1.20\pm 0.05$). La clase de cráteres más grande (100–300 km) tiene una incertidumbre estadística amplia ($\pm 0.38$) debido a una muestra reducida.
• Productividad del método: la correspondencia fila a fila sirve como guía física para leer la tendencia morfológica, y no constituye una afirmación de que cada cráter haya sido asignado cronológicamente en exclusiva a la fase concreta.

3. Conservación del diámetro de los cráteres sobre una corteza en contracción

Una de las preguntas centrales en la elaboración del modelo geofísico es: si la superficie polar se contrajo mecánicamente durante el redondeo, ¿por qué los propios cráteres no se contrajeron con ella (lo que los habría llevado por debajo del umbral de recuento de 20 km y habría hecho desaparecer el efecto)?

• Mecanismo de la corteza rígida aislada (C.7): en el marco del modelo físico, los cráteres conservaron su integridad morfológica y su diámetro original. Fueron transportados como elementos rígidos sobre una corteza superior sólida, mientras que la contracción de la superficie se produjo mediante el cierre y la reducción de los espacios elásticos entre los cráteres. Por tanto, el diámetro de los cráteres medido hoy refleja con precisión el diámetro original en el momento del impacto.
• Prueba matemática a partir de la distribución de tamaños (D.2): si los cráteres se hubieran contraído físicamente junto con la superficie, el efecto de apiñamiento se habría vuelto estadísticamente invisible para cualquier umbral de recuento fijo. El aumento del número de cráteres por unidad de área habría quedado completamente neutralizado por la disminución del diámetro de los cráteres limítrofes por debajo del umbral de 20 km (los dos efectos opuestos en la distribución de tamaños de los cráteres lunares se anulan mutuamente en una contracción homogénea). La mera existencia del exceso significativo medido ($1.20$) constituye una prueba empírica directa de que los cráteres conservaron su tamaño original.

4. Estructura y orden de los capítulos del artículo

La continuación del análisis científico del artículo se organiza según la siguiente estructura:

• Parte A (la observación empírica): presentación de los datos completos del recuento por latitudes, por longitudes, dentro de los casquetes polares y alrededor de los seis ejes geométricos del cuerpo (tanto en la versión bruta como en la limpia).
• Parte B (descarte de hipótesis alternativas): análisis detallado de las seis explicaciones científicas alternativas examinadas y de las razones mecánicas de su descarte.
• Parte C (el mecanismo físico): el cálculo matemático de los cambios diferenciales de área que acompañan la transición de un cuerpo triaxial achatado a una esfera perfecta.
• Parte D (las cuatro pruebas estructurales): presentación de los resultados de las cuatro pruebas estadísticas para evaluar la validez del modelo (razones de densidades, diámetros medianos, poblaciones jóvenes y simetría de los polos).
• Parte E (síntesis y cuestiones restantes): definición de los temas y preguntas abiertas que quedan para la investigación futura.

Parte A — La observación empírica

La Parte A presenta los resultados empíricos brutos tal como se cuentan en la superficie, antes de aplicar ninguna interpretación física ni modelos de formación. La observación se basa en un recuento normalizado: el número de cráteres en cada celda de área se divide por el área geométrica de esa celda.

En esta parte se realizan cuatro recuentos independientes sobre la misma base de datos y con la misma metodología:

• A.1 (corte latitudinal): división de la Luna en bandas de latitud paralelas al ecuador, para examinar la diferencia de densidad entre los polos y el ecuador.
• A.2 (corte longitudinal): división de la Luna en sectores longitudinales de polo a polo, para examinar la diferencia de densidad entre las distintas caras de la Luna (un patrón separado y distinto).
• A.3 y A.4 (procedimiento limpio): repetición exacta de las secciones A.1 y A.2 tras eliminar los suelos jóvenes (maria, llanuras basálticas) y las grandes cuencas (SPA) y recalcular las áreas residuales.
• A.5 (síntesis de las observaciones): resumen y conclusiones de las cuatro pruebas espaciales.

A.1 Densidad de cráteres por latitud (datos brutos)

El análisis espacial de esta sección tiene por objeto determinar si existe falta de uniformidad en la densidad de cráteres a lo largo de los ejes de latitud geográfica.

Criterios metodológicos y calibración del recuento:

1. Completitud de la base (catálogo Robbins): el recuento se limita a cráteres con un diámetro ajustado de 20 km o más. Esta limitación es estricta y garantiza la completitud del recuento en todas las latitudes: neutraliza los sesgos ópticos (condiciones de baja iluminación y sombras largas en los polos que ocultan cráteres pequeños) y filtra los cráteres secundarios (cúmulos y cadenas de fragmentos de menos de 20 km de diámetro, que producirían una densidad local engañosa).

2. Grupo de control y confirmación de la tendencia: se realiza una segunda serie para cráteres de 40 km de diámetro o más, a fin de verificar que el patrón espacial no depende de la población de cráteres limítrofes más pequeños.

3. Ausencia inicial de máscaras: el primer recuento presenta la superficie tal como es, sin omitir artificialmente los maria, partiendo de la idea de que forman parte de la predicción estructural (el efecto de su eliminación se analiza por separado en la Parte A.3).

4. Cálculo de áreas geométricas diferenciales: la Luna se divide en 12 bandas de latitud de anchura constante de $15^{\circ }$. El área de cada banda se calcula según la ecuación espacial:

donde el radio medio adoptado es $R=1,737.4\text{ km}$ (estándar IAU) y ${\varphi }_1,{\varphi }_2$ representan los límites del sector de la banda.

Según esta ecuación, el área de una banda sobre una esfera no viene determinada por su anchura angular, sino por su extensión a lo largo del eje polar (como la altura de la etiqueta de un cilindro). Cerca del ecuador la banda es alta a lo largo del eje, mientras que cerca de los polos la superficie se curva bruscamente y la circunferencia del anillo es pequeña. En consecuencia, las áreas de las bandas varían considerablemente: la banda ecuatorial ($0-15^{\circ }$) abarca 4.909 millones de km², mientras que la banda polar ($75-90^{\circ }$) abarca solo 0.646 millones de km² (una razón de alrededor de 7:1). Esta diferencia geométrica obliga a utilizar el índice de densidad (cráteres por 10⁶ km²) y no números absolutos.

Resumen de los datos del recuento de cráteres por bandas de latitud (A.1)

A continuación se presentan los datos numéricos completos y relativos medidos en las 12 bandas de latitud globales, sin aplicar máscaras geológicas:

• Determinación del punto de referencia (normalización): el valor de base 1.00 se define como la densidad combinada de la banda ecuatorial completa ($15^{\circ }\text{ S}-15^{\circ }\text{ N}$). En el grupo $\ge 20\text{ km}$ la densidad de referencia es exactamente 166.436 cráteres por millón de km² (un total de 1,634 cráteres sobre un área total de 9.818 millones de km²). En el grupo $\ge 40\text{ km}$ la densidad de referencia es 53.170. Las dos semibandas ecuatoriales muestran una asimetría interna (0.93 en el norte frente a 1.07 en el sur), y solo su combinación da el valor de referencia normalizado.
• Intervalos de las escalas de densidad: para el grupo $\ge 20\text{ km}$ el intervalo de valores va de 151 a 300; para el grupo $\ge 40\text{ km}$ el intervalo va de 49 a 101 cráteres por millón de km². Las razones se calculan a partir de los valores brutos de recuento y de área, sin redondear.

Análisis y conclusiones del perfil espacial

De la lectura de las columnas de densidad y del perfil gráfico se desprenden los siguientes resultados principales:

• Curvatura bipolar marcada: los dos extremos de la tabla presentan los valores de densidad máximos del sistema. La banda junto al polo norte ($75-90^{\circ }\text{ N}$) alcanza una densidad máxima de unos 300 cráteres por millón de km², valor que equivale a 1.80 veces el del ecuador. La banda junto al polo sur ($75-90^{\circ }\text{ S}$) tiene una densidad de 215, equivalente a 1.29 veces el ecuador.
• Zonas de mínimo ecuatorial: las filas más pobres en cráteres se sitúan al norte del cinturón del ecuador geográfico, en los intervalos $0-15^{\circ }\text{ N}$ y $15-30^{\circ }\text{ N}$, con valores de densidad mínimos de 155 y 151 respectivamente.
• Acentuación de la tendencia en tamaños mayores: en el grupo de control de cráteres de diámetro $\ge 40\text{ km}$, la diferencia de densidad bipolar se agudiza y muestra razones extremas de 1.89 en el norte y 1.54 en el sur respecto a la referencia ecuatorial.
• Asimetría del hemisferio medio: el perfil no asciende de forma completamente lineal y suave del ecuador a los polos. El hemisferio sur es más rico en densidad de cráteres que el hemisferio norte en las latitudes medias (valor de 230 en el sector $45-60^{\circ }\text{ S}$ frente a 174 en el sector $45-60^{\circ }\text{ N}$). La gradación sube y baja en distintos lugares, pero la conclusión estructural es firme: las bandas extremas junto a ambos polos presentan un exceso de densidad anómalo que ninguna banda intermedia alcanza. Parte del valor bajo en el centro se debe a la concentración de maria en el ecuador, factor que se trata y neutraliza en la sección A.3.

Escala de color: a la izquierda, «low density», densidad baja (amarillo suave); a la derecha, «high density», densidad alta (rojo suave).

Figura 1. Densidad de cráteres por latitud, cráteres de 20 km o más, en doce bandas de 15° (A.1); el catálogo completo, sin máscara. Los dos polos son los más altos. (En el gráfico: eje vertical — cráteres ≥ 20 km por 10⁶ km², 100–300; eje horizontal — banda de latitud, el polo sur a la izquierda y el polo norte a la derecha, de 75–90S a 75–90N; la línea azul oscuro — la densidad en cada banda.)

A.2 Densidad de cráteres por longitud (datos brutos)

Mientras que la sección A.1 halló un exceso de densidad significativo en las regiones polares, esta sección examina la uniformidad de la distribución espacial en la dirección perpendicular, a lo largo de la circunferencia global de la Luna. El objetivo del análisis es aislar y evitar la mezcla estadística entre el efecto de latitud (el tema central de la investigación) y posibles asimetrías entre distintos hemisferios a lo largo de los ejes de longitud.

Criterios metodológicos y geometría de los sectores:

1. Estructura de las celdas de recuento: la superficie lunar se dividió en 12 bandas de longitud (sectores) de anchura constante de $30^{\circ }$. Cada banda se extiende de forma continua de polo a polo e incluye todas las latitudes (como un gajo de naranja).

2. Igualdad de áreas intrínseca: a diferencia de las bandas de latitud, todos los sectores de longitud sobre una esfera tienen exactamente la misma área. Esta forma geométrica hace innecesario aplicar factores de corrección o de normalización por el tamaño de la celda, y permite comparar directamente el número absoluto de cráteres (los recuentos brutos) con las razones de densidades.

3. Dimensiones de la superficie: el área de cada banda se calcula con precisión según:

Para un radio $R=1,737.4\text{ km}$, el área de cada sector de longitud es constante e igual a 3.161 millones de km² (1/12 de toda la superficie lunar).

4. El sistema de coordenadas dinámico: la banda que representa el sector $345^{\circ }-15^{\circ }$ este contiene el punto subterrestre ($0^{\circ }\text{ N},0^{\circ }\text{ E}$), centro de la cara orientada permanentemente hacia la Tierra debido al acoplamiento de marea.

5. Definición de los componentes orbitales (“delantero” frente a “de cola”): los términos se aplican según las definiciones mecánicas habituales:

• El punto delantero (Leading): el ápice del vector de movimiento orbital de la Luna, situado en ($0^{\circ },90^{\circ }\text{ W}$). Mira hacia delante en la dirección de la órbita (el frente del cuerpo).
• El punto de cola (Trailing): el punto que mira hacia atrás a lo largo de la dirección del vuelo orbital, situado en ($0^{\circ },90^{\circ }\text{ E}$).

6. Base de normalización fija: el punto de referencia (1.00) es exactamente el mismo que se fijó en la sección A.1: la densidad bruta de la banda ecuatorial completa ($15^{\circ }\text{ S}-15^{\circ }\text{ N}$). Su valor es 166.436 cráteres por millón de km² para el grupo $\ge 20\text{ km}$ y 53.170 para el grupo $\ge 40\text{ km}$. Mantener esta base fija permite una comparación directa entre las cuatro tablas de análisis de la Parte A (de A.1 a A.4).

Resumen de los datos del recuento de cráteres por bandas de longitud (A.2)

A continuación se presentan los datos numéricos completos medidos en los 12 sectores de longitud globales, sin aplicar máscaras geológicas:

Análisis y conclusiones del perfil longitudinal global

Del examen de los indicadores y del perfil gráfico se desprenden los siguientes resultados principales:

• Marcada falta de uniformidad en la circunferencia del cuerpo: la superficie lunar presenta una asimetría extrema a lo largo de las longitudes. La cara oculta y los sectores que rodean la longitud de cola ($90^{\circ }\text{ E}$) presentan densidades sistemáticamente altas, entre 1.22 y 1.51 veces el valor de referencia ecuatorial en el grupo $\ge 20\text{ km}$.
• Zonas de máximo y mínimo globales:
  • Densidad máxima: registrada en el sector oriental de la cara oculta, en el intervalo $105^{\circ }-135^{\circ }$ este, con 793 cráteres y una densidad de 251 por millón de km² (razón de 1.51).
  • Densidad mínima: registrada en la cara cercana occidental, en el intervalo $315^{\circ }-345^{\circ }$ este, con solo 267 cráteres y una densidad de 84 por millón de km² (razón de 0.51).
• Amplitud del patrón longitudinal: la razón extrema entre el sector máximo y el sector mínimo es de casi tres veces ($251/84=2.97$ para el grupo $\ge 20\text{ km}$ y $90/31=2.88$ para el grupo $\ge 40\text{ km}$). Esta diferencia es considerablemente mayor que las diferencias brutas de densidad medidas entre los polos y el ecuador en el capítulo A.1. Este resultado confirma que se trata de un patrón geográfico independiente y completamente separado, que requiere un aislamiento metodológico.
• Contexto geofísico y morfológico (modelo LEF): el sector de mínimo extremo ($315^{\circ }-345^{\circ }$) alberga la cuenca del Mar de las Lluvias (Mare Imbrium). En el marco del modelo LEF, esta estructura actuó como la abertura principal de alimentación y conducción de calor desde la Tierra al interior de la Luna, y no como un lugar de expulsión y eyección de material. Este sector, junto con las extensas áreas de maria occidentales concentradas en él, se caracteriza por una corteza fundida y reiniciada térmicamente, por lo que hoy presenta pocos cráteres. En cambio, el notable exceso de cráteres en la cara de cola se interpreta en el modelo a través de la dispersión de los materiales eyectados que retornaron, quedaron rezagados respecto a la velocidad de rotación y cayeron detrás (Artículo 6).

Escala de color: a la izquierda, «low density», densidad baja (amarillo suave); a la derecha, «high density», densidad alta (rojo suave).

Figura 2. Densidad de cráteres por longitud, cráteres de 20 km o más, en doce bandas de 30° (A.2); el catálogo completo, sin máscara. El máximo se encuentra algo al este de la longitud de cola, y el mínimo algo al oeste del punto subterrestre. (En el gráfico: eje vertical — cráteres ≥ 20 km por 10⁶ km², 50–300; eje horizontal — banda de longitud (° este), desde el punto subterrestre hacia el este pasando por la longitud de cola, el punto antiterrestre y la longitud delantera, de 345–15 a 315–345; la línea azul oscuro — la densidad en cada banda.)

A.3 Densidad de cráteres por latitud tras excluir las cuencas y los maria

Mientras que las secciones A.1 y A.2 contaron la superficie global completa, esta sección aplica una máscara geológica estricta para eliminar las unidades de superficie cuyo registro de cráteres se vio afectado por factores endógenos (actividad térmica volcánica y borrado mecánico) no relacionados con la latitud geográfica.

Fundamento de la eliminación de las unidades geológicas:

1. El sistema de maria y llanuras basálticas (máscara LROC): los suelos de los maria representan aberturas de alimentación jóvenes sobre domos magmáticos (Artículo 5), en las que los cráteres antiguos fueron fundidos. Las llanuras basálticas (como Procellarum, Tranquillitatis, Fecunditatis, Insularum, Frigoris) representan corteza fundida por el calor global irradiado por la Tierra, sin un domo de alimentación local. En ambos casos, la escasez de cráteres se debe a un reinicio térmico del registro. Dado que estas unidades se concentran en la cara cercana y en latitudes bajas, rebajan artificialmente la densidad ecuatorial.

2. Las grandes cuencas ($\ge 500\text{ km}$): estas estructuras de impacto se extienden sobre amplias porciones de las bandas de latitud (principalmente en las regiones polares). Su suelo y su borde tienen una historia morfológica propia que afecta de forma independiente a la estadística de la banda.

3. Objetivo del procedimiento: aislar las porciones de corteza históricamente más estables (Highlands), en las que el registro se conservó uniforme, para obtener un indicador fiable de la distribución original de los cráteres.

Metodología de aplicación de la máscara geológica

• Definición de la máscara: combinación de los polígonos globales de maria del LROC (LROC_GLOBAL_MARE_180) y las 20 cuencas más grandes con diámetro de borde exterior $\ge 500\text{ km}$ del catálogo Neumann et al. (2015) (Tabla 1).
• Geometría de las cuencas: cada cuenca se representa como un casquete esférico centrado en las coordenadas de origen, con un radio angular de $D/(2R_{\text{Moon}})$.
• Regla de filtrado y área residual: un cráter se elimina del recuento si su punto central geométrico cae dentro de una cuenca o de un polígono de mare. Los solapamientos espaciales se eliminan una sola vez. El cálculo del área residual se realiza sobre una malla global de alta resolución de $0.05^{\circ }\times 0.05^{\circ }$, utilizando áreas exactas de celdas esféricas basadas en la distancia de círculo máximo.
• Mantenimiento del punto de referencia para la comparación: para que el lector pueda seguir el movimiento de las bandas, se mantiene la base de referencia estática de la sección A.1 (valor 1.00 = densidad bruta del cinturón ecuatorial completo sin máscara). Su valor es 166.436 (para $\ge 20\text{ km}$) y 53.170 (para $\ge 40\text{ km}$).
• Definición del índice de cambio de densidad ($\Delta \rho$): índice que mide cada banda respecto a sí misma antes y después de la máscara:

Un valor positivo indica que el área eliminada era pobre en cráteres en relación con el resto de la banda (y por tanto su eliminación aumentó la densidad residual).

Resumen de los datos de densidad de cráteres enmascarada por latitud (A.3)

A continuación se presentan los valores numéricos exactos obtenidos en las 12 bandas de latitud tras aplicar la máscara de filtrado:

• Escalas de colores numéricas de la tabla: en el grupo $\ge 20\text{ km}$ el intervalo va de 202 a 305; en el grupo $\ge 40\text{ km}$ el intervalo va de 65 a 124 cráteres por millón de km².

Análisis de los efectos espaciales y del sesgo intrínseco de la base estática

De la comparación del perfil enmascarado con el perfil bruto se desprenden las siguientes conclusiones científicas:

1. Comportamiento de los casos extremos de la columna de cambio:

• Cambio polar nulo en el norte ($75-90^{\circ }\text{ N}$): esta banda no toca físicamente la máscara (no contiene maria ni grandes cuencas). Su área ($0.646\times 10^6{\text{ km}}^2$) y su número de cráteres (194) permanecieron constantes, por lo que $\Delta \rho =0.0\%$.
• Cambio ecuatorial máximo en el norte ($15-30^{\circ }\text{ N}$): esta banda pierde más de la mitad de su área geométrica (caída de 4.574 a 2.153 millones de km²) debido a la concentración de maria de la cara cercana. La densidad residual de lo que queda se dispara un +48.4% para el grupo de 20 km y un +55.9% para el grupo de 40 km (el aumento más pronunciado del sistema).

2. Persistencia del exceso polar endógeno: la aplicación de la máscara eleva considerablemente las latitudes bajas y medias (especialmente en el norte), con lo que aplana parte de la diferencia espacial. Sin embargo, las dos bandas junto a los polos siguen siendo las regiones más densas de todo el cuerpo lunar. En el grupo $\ge 20\text{ km}$, el polo sur alcanza 305 y el norte 300 por millón de km², mientras que las 10 bandas intermedias son más bajas (valores de 202 a 246). En el grupo $\ge 40\text{ km}$, el polo sur encabeza todo el sistema con un valor máximo de 124 (razón de 2.32). Este resultado demuestra que el exceso polar no es un artefacto de la presencia de los maria en el ecuador; está inscrito en el propio terreno sólido.

3. Aclaración metodológica sobre la comparación de igual a igual: la columna de razones de esta tabla divide una banda enmascarada (un numerador que ha perdido sus maria) por un ecuador bruto (un denominador que los conservó), con lo que produce una exageración matemática del exceso predicho (razones de 1.40–1.48). Si se comparan los dos casquetes polares enmascarados (882 cráteres en 3.706 millones de km², lo que da una densidad de 238.0) con el cinturón ecuatorial enmascarado con exactamente la misma máscara (1,216 cráteres en 5.856 millones de km², lo que da una densidad de 207.7), la razón real es 1.15. Por ello, para la prueba científica final del artículo se aplica el recuento limpio de los capítulos A.11 y D.1 (eliminación simultánea de los maria y de la cuenca SPA en ambos lados), que da el valor de razón normalizado y estricto: $1.20\pm 0.05$.

Escala de color: a la izquierda, «low density», densidad baja (amarillo suave); a la derecha, «high density», densidad alta (rojo suave).

Figura 3. Las mismas doce bandas, tras excluir las veinte cuencas de 500 km o más y los maria y descontar su área (A.3). La máscara eleva el ecuador y deja los polos. (En el gráfico: eje vertical — cráteres ≥ 20 km por 10⁶ km², 100–300; eje horizontal — banda de latitud, el polo sur a la izquierda y el polo norte a la derecha; la línea roja — el perfil enmascarado.)

A.4 Densidad de cráteres por longitud tras excluir las cuencas y los maria

Esta sección aplica la máscara geológica estricta (las 20 grandes cuencas de Neumann y los polígonos de maria del LROC) a las 12 bandas de longitud de la sección A.2. Su objetivo principal es determinar qué parte de la asimetría extrema entre este y oeste procede directamente de la presencia de los maria jóvenes y de las cuencas, y qué parte está inscrita en la propia corteza montañosa antigua.

Metodología del cálculo de áreas y calibración del sistema:

1. Pérdida de la igualdad de áreas en las bandas de longitud: a diferencia del recuento bruto, en el que todos los sectores de longitud tienen el mismo tamaño, la aplicación de la máscara modifica las áreas de las bandas de forma diferencial. Una banda que contenía una porción amplia de maria o cuencas pierde más área. Las áreas residuales efectivas oscilan ahora en un amplio intervalo de entre 1.073 y 2.729 millones de km², lo que obliga a basarse exclusivamente en el índice de densidad y no en recuentos absolutos.

2. Coherencia y cierre matemático de la máscara: al sumar las 12 bandas se obtiene una concordancia completa con los datos de la sección A.3: los dos métodos de corte geométrico (latitudinal y longitudinal) dividen exactamente la misma área de corteza residual, que constituye algo más de la mitad de la Luna (21.570 millones de km²), y contienen exactamente el mismo número de cráteres (4,848 en el grupo $\ge 20$ km y 1,599 en el grupo $\ge 40$ km).

3. Mantenimiento de la base de referencia fija: el punto de referencia (valor 1.00) permanece completamente fijo: la densidad bruta del cinturón ecuatorial completo sin máscara (166.436 para $\ge 20$ km y 53.170 para $\ge 40$ km).

Resumen de los datos de densidad de cráteres enmascarada por longitud (A.4)

A continuación se presentan los valores numéricos exactos obtenidos en los 12 sectores de longitud tras aplicar la máscara de filtrado:

Análisis dinámico del efecto de la máscara sobre la asimetría longitudinal

Del examen de los indicadores residuales y de la columna de cambio geológico se desprenden las siguientes conclusiones:

1. Relleno de la depresión occidental extrema: la banda de longitud $345^{\circ }-315^{\circ }$ este —que alberga la cuenca Mare Imbrium y la mayor parte de los maria occidentales de la cara cercana— muestra el aumento más pronunciado del sistema: un salto de densidad del +108.5% (la densidad se duplica con creces, de 84 a 176 por millón de km²). Este resultado demuestra que el mínimo extremo medido en la sección A.2 se debió directamente al reinicio del registro por los suelos basálticos jóvenes.

2. Comportamiento anómalo en la zona de cola: el sector que rodea la longitud de cola ($105^{\circ }-75^{\circ }$ este) es la única banda del sistema en la que la densidad residual disminuye tras aplicar la máscara (caída del -7.8% en el grupo de 20 km y del -9.0% en el grupo de 40 km). Esto significa que las partes eliminadas en este sector (cuencas locales) eran más ricas y densas en cráteres que la corteza que quedó a su alrededor.

3. Reducción de la amplitud del patrón longitudinal: la máscara de filtrado reduce drásticamente la asimetría entre este y oeste y comprime el perfil en un intervalo más estrecho. La razón extrema entre el sector máximo y el sector mínimo cae de 2.97 a 1.68 para cráteres $\ge 20$ km (máximo de 266 en el sector $135^{\circ }-105^{\circ }$ frente a un mínimo de 158 en el sector $315^{\circ }-285^{\circ }$), y de 2.88 a 1.79 para cráteres $\ge 40$ km (92 frente a 51).

4. Persistencia del patrón residual de la circunferencia: a pesar de la reducción matemática, el patrón espacial básico no desaparece. La cara oculta más allá de la longitud de cola mantiene firmemente su condición de región más densa en ambos grupos de control, mientras que la cara cercana más allá de la longitud delantera sigue siendo la más pobre. Este patrón residual se analiza más adelante en el artículo mediante el modelo de los seis ejes (A.9). Para resolver el enigma polar, la propia medición aislada y la separación de los efectos longitudinales proporcionan la base metodológica necesaria.

Escala de color: a la izquierda, «low density», densidad baja (amarillo suave); a la derecha, «high density», densidad alta (rojo suave).

Figura 4. Las mismas doce bandas de longitud, tras excluir las veinte cuencas de 500 km o más y los maria y descontar su área (A.4). El mínimo occidental sube; el máximo de la cara oculta se mantiene. (En el gráfico: eje vertical — cráteres ≥ 20 km por 10⁶ km², 50–300; eje horizontal — banda de longitud (° este), desde el punto subterrestre hacia el este pasando por la longitud de cola, el punto antiterrestre y la longitud delantera; la línea roja — el perfil enmascarado.)

A.5 Síntesis y conclusiones de las cuatro pruebas de bandas

Esta sección reúne y analiza los resultados combinados de las cuatro pruebas de bandas espaciales realizadas en los capítulos A.1–A.4. El análisis pretende responder a una pregunta estricta y estructural: ¿es el exceso de densidad de cráteres en las regiones polares una propiedad geofísica real de la corteza lunar, o se trata de un sesgo estadístico artificial derivado del modo de dividir la superficie o de la presencia de grandes formaciones geológicas locales?

La metodología se basa en el cruce de cuatro pruebas complementarias:

1. Distribución por latitud sin máscara (A.1): examen de la tendencia ecuador frente a polos sobre la superficie bruta.

2. Distribución por longitud sin máscara (A.2): examen de la tendencia de la circunferencia global (este frente a oeste) sobre la superficie bruta.

3. Distribución por latitud enmascarada (A.3): limpieza de los efectos endógenos (maria, cuencas $\ge 500$ km) a lo largo de los ejes de latitud.

4. Distribución por longitud enmascarada (A.4): limpieza de los efectos endógenos a lo largo de los ejes de longitud.

La combinación de las cuatro pruebas permite aislar entre dos direcciones de corte perpendiculares (latitud frente a longitud) y, al mismo tiempo, controlar el grado de dependencia del registro respecto a las superficies basálticas jóvenes y a las cuencas gigantes.

A.5.1 Separación de los patrones: dos configuraciones espaciales independientes

La primera conclusión científica que se desprende de la síntesis es que no existe un perfil espacial único capaz de explicar la distribución global de los cráteres. Los ejes de latitud y los de longitud revelan dos patrones geográficos separados y completamente distintos:

1. Características del perfil de latitud (tendencia bipolar):

• La diferencia de densidad extrema: en el recuento bruto (sin máscara), en el grupo de diámetro $\ge 20$ km, la banda junto al polo norte ($75-90^{\circ }\text{ N}$ ◐) muestra una densidad máxima de unos 300 cráteres por millón de km², mientras que la banda mínima ($15-30^{\circ }\text{ N}$ ◐) se sitúa en unos 151 por millón de km²: una diferencia de casi el doble (razón 1.80) [A.1]. En el grupo de control $\ge 40$ km se registró una razón similar: unos 101 por millón de km² en el polo norte frente a unos 51 por millón de km² en las latitudes bajas ($15-30^{\circ }\text{ N}$) [A.1].
• Comportamiento no monótono: el perfil espacial no asciende en una gradación suave, continua y coherente de banda en banda en el camino del ecuador a los polos; incluye oscilaciones, subidas y bajadas locales [A.1].
• La conclusión estructural acotada: el argumento científico no habla de una gradación lineal suave, sino de una propiedad de los extremos firme: las bandas extremas junto a ambos polos alcanzan valores de densidad anómalos que ninguna banda intermedia del recorrido es capaz de alcanzar [A.1]. Un observador que se desplace del ecuador hacia el polo verá un recuento oscilante, pero la banda final, en la que se produce el corte polar, mostrará una densidad aproximadamente el doble de la de la banda más pobre medida en su recorrido [A.1].

2. Características del perfil de longitud (asimetría de la circunferencia):

• La diferencia de densidad extrema: en el recuento bruto, en el grupo $\ge 20$ km, la densidad va desde un mínimo absoluto de unos 84 por millón de km² en el sector occidental de la cara cercana ($315-345^{\circ }\text{ E}$ ◐) hasta un máximo global de unos 251 por millón de km² en el sector oriental de la cara oculta ($105-135^{\circ }\text{ E}$ ◐) [A.2]. En el grupo $\ge 40$ km el intervalo va de unos 31 a unos 90 cráteres por millón de km² [A.2].
• Índice de falta de uniformidad (razón máximo/mínimo): el factor de variación (densidad del sector más rico dividida por la densidad del sector más pobre) es de 2.97 para el grupo de 20 km y de 2.88 para el grupo de 40 km [A.2].
• Comparación de la magnitud de los efectos: el grado de falta de uniformidad a lo largo de la circunferencia (una diferencia de casi el triple entre este y oeste) es más agudo y amplio que el medido a lo largo de los ejes de latitud (una diferencia del doble entre los polos y el ecuador) [A.1, A.2]. Esto significa que un desplazamiento circunferencial alrededor del ecuador produce cambios de densidad más extremos que un desplazamiento vertical del ecuador a los polos [A.2].

Análisis metodológico: descarte de la mezcla de patrones

El hecho de que la asimetría a lo largo de las longitudes alcance un intervalo del triple podría llevar a afirmar que el exceso polar (que es menor, del doble) no es más que un subproducto o una distorsión estadística causada por el patrón mayor de longitud. La investigación descarta esta afirmación mediante la siguiente explicación geométrica:

• Perpendicularidad estructural absoluta: los dos perfiles se calculan y se miden en direcciones completamente perpendiculares entre sí.
• Mecanismo de compensación geométrica: cada banda de latitud es un anillo completo que rodea la Luna en $360^{\circ }$. Esto significa que cada banda de latitud absorbe exactamente todas las longitudes globales: las longitudes más ricas (la cara oculta y la de cola) y las más pobres (la cara cercana occidental) entran en el denominador estadístico de cada banda de latitud en la misma proporción geométrica fija.
• Aislamiento del efecto: la existencia simultánea de dos perfiles distintos demuestra que la presencia de una región rica o pobre en cráteres a lo largo del eje de la circunferencia no constituye por sí misma una explicación ni una fuente del exceso de densidad inscrito en las latitudes polares. La asimetría longitudinal es un patrón independiente y separado que se trata en el marco del modelo de los seis ejes (A.9) [A.2]. El presente análisis aísla con éxito el patrón polar para examinar con precisión qué queda de él en condiciones de control uniformes.

A.5.2 La máscara como prueba de control

Las secciones A.3 y A.4 no son un nuevo recuento de una población de cráteres nueva, sino una repetición metodológica exacta de las dos pruebas de bandas (latitudinal y longitudinal) en condiciones de filtrado con una máscara geológica uniforme. Esta máscara resta del cálculo las 20 cuencas de diámetro $\ge 500$ km y los polígonos de maria del LROC [A.3, A.4]. Un cráter cuyo centro geométrico cae dentro de los límites de la máscara se omite del recuento, y el área del sector se descuenta en consecuencia del denominador de la densidad [A.3, A.4].

Fundamento científico de la elección de las superficies filtradas:

• Los maria y las llanuras basálticas (reinicio térmico): estas unidades son suelos geológicamente jóvenes. El mare actúa como una abertura sobre un domo de alimentación magmático que fundió los cráteres tempranos, mientras que la llanura basáltica representa una corteza fundida bajo la carga térmica general irradiada por la Tierra. Su escasez de cráteres se debe a la renovación de la superficie y no a una baja tasa de impactos. Puesto que se concentran estadísticamente en la zona ecuatorial, su inclusión distorsiona el índice de latitud [A.3].
• Las cuencas gigantes (borrado mecánico): la formación de estas estructuras de impacto borró físicamente la población de cráteres que existía en su interior antes que ellas [A.5.2].
• La regla de comparación de igual a igual: para determinar si la diferencia de densidad procede de un efecto geodinámico dependiente de la latitud, la física exige comparar unidades de corteza con una historia uniforme (terreno antiguo frente a terreno antiguo). Incluir superficies parcialmente renovadas mediría el ritmo y la distribución de la actividad de renovación térmica, y no el efecto de la geometría del cuerpo [A.5.2].

La dinámica matemática del índice de densidad:

Una disminución del número absoluto de cráteres ($N$) tras aplicar la máscara no implica necesariamente una disminución de la densidad residual [A.3]. La ecuación de la densidad se rige por la relación variable entre el numerador y el denominador sobre el mismo terreno [A.3]:

Si la máscara resta una porción relativa de área mayor que la porción de cráteres que omite (situación típica de las áreas de maria vacías de cráteres), la densidad global de la superficie residual mostrará una tendencia creciente [A.3]. El mecanismo funciona de manera similar al cálculo de un índice de árboles por hectárea en un bosque que contiene un lago: restar el área del lago (vacía de árboles) del denominador de área total elevará matemáticamente el valor de densidad de la tierra restante, aunque no se haya añadido ni un solo árbol físico al sistema [A.5.2].

Aplicar la misma operación analítica a los dos ejes de corte perpendiculares produce una poderosa herramienta de control: si los dos efectos espaciales (latitudinal y longitudinal) fueran consecuencia exclusiva de la presencia de los maria y las cuencas, la aplicación de la máscara debería modificar ambos perfiles de forma coherente y similar. Los resultados empíricos lo descartan y muestran un comportamiento diferencial significativo [A.5.2].

A.5.3 Respuesta del perfil de longitud a la máscara de filtrado

A diferencia del perfil de latitud, el efecto de la máscara sobre el perfil de longitud es drástico y modifica de forma fundamental la amplitud (Mesa) bruta del sistema [A.5.3]:

1. Convergencia de los valores extremos y reducción de la amplitud a lo largo de la circunferencia:

• El sector mínimo occidental ($315^{\circ }-345^{\circ }$ este): en el recuento bruto (sin máscara), este sector presentaba una densidad mínima de unos 84 cráteres por millón de km² [A.2]. Al aplicar la máscara y limpiar las áreas de Mare Imbrium y del conjunto de maria occidentales, la densidad calculada sobre las superficies montañosas restantes se disparó a unos 176 por millón de km² (un aumento del +108.5% respecto a su estado bruto) [A.4].
• Índice global de falta de uniformidad (razón máximo/mínimo): como consecuencia del aumento de los valores mínimos, el factor de variación global a lo largo de los ejes de longitud se redujo considerablemente: para el grupo $\ge 20$ km, la razón cayó de 2.97 en el recuento bruto a 1.68 en el recuento enmascarado [A.2, A.4]. En el grupo de control $\ge 40$ km se registró un comportamiento paralelo, con una caída de la razón de 2.88 a 1.79 [A.2, A.4]. Esta coherencia en ambos grupos de tamaño confirma que la amplia amplitud original de la longitud se debía principalmente a un sesgo geológico de la distribución de los maria [A.5.3].

2. Colapso de la asimetría entre la longitud de cola y la delantera:

El análisis estadístico de las bandas de longitud completas (que se extienden de polo a polo) revela el colapso de la significación de la asimetría orbital a nivel de banda [A.5.3]:

• En el recuento bruto: la banda de cola ($75^{\circ }-105^{\circ }$ este) contaba 750 cráteres ($\rho =237$) frente a 487 cráteres ($\rho =154$) en la banda delantera ($255^{\circ }-285^{\circ }$ este) [A.2]. Esta diferencia daba una razón significativa de 1.54, valor que representa 7 desviaciones típicas del recuento ($7\sigma$) y descarta por completo la probabilidad de una fluctuación aleatoria (la variación natural de un recuento de cráteres es del orden de la raíz del recuento $\sqrt{N}$, y en una muestra de cientos de cráteres la incertidumbre es de unos pocos por ciento) [A.5.3].
• Tras aplicar la máscara: la banda delantera conservó un área residual mayor ($2.044\times 10^6{\text{ km}}^2$, 417 cráteres) que la banda de cola ($1.692\times 10^6{\text{ km}}^2$, 370 cráteres) debido a una concentración más extensa de maria en la cara de cola [A.4]. En consecuencia, los valores de densidad enmascarados convergieron a 219 frente a 204 cráteres por millón de km², lo que da una razón residual de solo 1.07 [A.4]. Como la incertidumbre intrínseca del recuento en esta razón es de alrededor del 7%, la diferencia restante equivale a una sola desviación típica ($1\sigma$) y no tiene significación estadística a nivel de banda (para el grupo $\ge 40$ km se registró un valor residual de 77 frente a 74) [A.4, A.5.3].

3. El efecto geométrico: promediado en bandas frente a concentración en casquetes:

A pesar del colapso de la significación a nivel de la banda larga, el patrón longitudinal no desaparece por completo de la corteza sólida [A.4]. La razón reside en la limitación geométrica de las bandas de longitud completas: una banda que se extiende de polo a polo tiende a diluir y enmascarar los excesos locales al promediarlos con áreas de terreno ordinarias a lo largo de toda la longitud [A.5.3].

Cuando se cambia el método geométrico y se pasa a examinar casquetes centrados con un radio de $30^{\circ }$ alrededor de los dos puntos extremos orbitales (sección A.9), la población enmascarada mantiene valores de densidad de 252 por millón de km² en la cara de cola frente a 174 en la cara delantera [A.5.3]. Esta diferencia da una razón de 1.45, que representa 5 desviaciones típicas ($5\sigma$) y tiene una significación fuerte. Este resultado demuestra que el exceso de cráteres que sobrevive a la máscara se concentra localmente alrededor de las coordenadas centrales de la longitud $90^{\circ }$ este, y no se distribuye de forma homogénea a lo largo de toda la banda [A.5.3].

Síntesis metodológica:

La aplicación de la máscara de control produce un doble resultado: modifica y reduce una parte considerable de la amplitud longitudinal bruta (y confirma el papel de los maria jóvenes), pero no elimina por completo la falta de uniformidad residual a lo largo de la circunferencia [A.5.3]. No obstante, para resolver el enigma polar y confirmar el modelo de redondeo y contracción mecánica del artículo, basta con que el patrón longitudinal haya sido medido, definido y separado estadísticamente, por lo que la explicación física de la asimetría entre este y oeste queda fuera del argumento central de la presente investigación [A.5.3].

A.5.4 Respuesta del perfil de latitud a la máscara de filtrado: estabilidad del exceso polar

A diferencia del perfil de longitud, en el que la máscara borra una parte considerable de la amplitud asimétrica [A.5.3], el perfil de latitud responde al procedimiento de filtrado de un modo que confirma la persistencia del efecto polar en la corteza rocosa antigua [A.5.4]. El comportamiento de las bandas extremas refleja la mecánica de la ecuación de la densidad residual [A.3, A.5.4]:

• La banda junto al polo sur (75°-90° S ◐): en el recuento bruto (sin máscara) se contaron en esta banda 139 cráteres (ρ ≈ 215 por millón de km²) [A.1]. La aplicación de la máscara restó al sector una porción de área enorme y dejó solo 84 cráteres [A.3]. Sin embargo, dado que la cuenca Polo Sur–Aitken (South Pole–Aitken - SPA) —la estructura gigante que cubre una parte considerable de esta banda— se caracteriza por un suelo pobre en cráteres debido al borrado mecánico (efecto de la “piscina seca”), la máscara restó un área relativa mucho mayor que la cantidad de cráteres omitidos [A.5.4]. En consecuencia, el área disponible se redujo a unos 0.275 × 10⁶ km², y la densidad calculada sobre ella se disparó a unos 305 cráteres por millón de km² [A.3].
• La banda junto al polo norte (75°-90° N ◐): en el recuento bruto se contaron en este sector 194 cráteres, lo que da una densidad de unos 300 por millón de km² [A.1]. Dado que esta banda septentrional apenas se ve afectada geográficamente por los polígonos de la máscara (ausencia de grandes cuencas o de maria extensos cerca del polo norte), los datos de área y numerador permanecieron sin cambios, y la densidad residual queda fijada en unos 300 cráteres por millón de km² [A.3].
• Convergencia de los valores extremos bipolares: la limpieza de los registros geológicos conduce a una corrección matemática fundamental. En el recuento bruto, el hemisferio sur presentaba una apariencia engañosa de escasez respecto al norte (215 frente a 300) [A.1]. Tras aplicar la máscara, esta asimetría artificial colapsa, y las dos bandas junto a los polos convergen casi al mismo valor firme y elevado: 305 en el sur y 300 en el norte [A.3].

A.5.5 Conclusiones científicas de la comparación de los polos y prueba de igualdad de superficie

La convergencia de los valores entre norte y sur constituye uno de los resultados metodológicos principales de la prueba de control, y proporciona una lección firme en el tratamiento de datos planetarios [A.5.5]:

• Colapso de la asimetría norte-sur bruta: la comparación bruta original incurría en el error de enfrentar una banda polar norte compuesta por corteza virgen y antigua apenas alterada con una banda polar sur compuesta en su mayor parte por el suelo de una cuenca de impacto gigante que borró el registro temprano [A.5.5]. Esta comparación no representaba una propiedad real de las latitudes, sino una propiedad del tipo de terreno [A.5.5].
• Establecimiento de una comparación de igual a igual: en cuanto se aplica una máscara de control uniforme que omite las mismas unidades geológicas en ambos lados, las dos bandas se miden sobre el mismo tipo de superficie (Highlands vírgenes), y la asimetría interna entre ellas se anula en el grupo de ≥ 20 km [A.5.5].
• La conclusión morfológica: la neutralización de la diferencia bruta no produce un polo denso y un polo opuesto pobre, sino que revela dos casquetes polares de densidad similar, estable y considerablemente más alta que el resto de la superficie [A.5.5].

A.5.6 Síntesis de la primera fase de las observaciones espaciales y limitaciones de la geometría de bandas

Las cuatro pruebas de bandas (A.1–A.4) cierran por completo la cuestión de la validez del hallazgo morfológico y permiten reducir el problema científico a cuatro conclusiones firmes y una limitación estructural [A.5.6]:

Conclusiones establecidas:

1. El efecto es global y bipolar: la alta densidad polar no es una característica local exclusiva del hemisferio norte; es común a las dos bandas polares extremas una vez eliminados los sesgos [A.5.6].

2. El exceso no depende de la máscara: el efecto no desaparece ni se disipa tras eliminar las grandes cuencas y los maria. Al contrario, en el polo sur la limpieza de la superficie puso de relieve y elevó la densidad real de la corteza restante [A.5.6].

3. Aislamiento total respecto a la asimetría longitudinal: la densidad polar no es una expresión indirecta de la asimetría entre este y oeste. El perfil de longitud muestra una amplitud distinta, un comportamiento geométrico separado y una respuesta completamente opuesta a la máscara (colapso de la mayor parte de la amplitud longitudinal a nivel de banda frente a la conservación y el realce de la densidad en los extremos polares) [A.5.3, A.5.6].

4. No se requiere un perfil monótono suave: los datos empíricos de la superficie no respaldan la afirmación de una gradación suave que asciende de forma continua del ecuador a los polos. El hallazgo científico preciso —y físicamente más sólido— establece que las bandas junto a los polos se sitúan de forma significativa en el extremo superior absoluto de la distribución global de latitudes [A.5.6].

La limitación geométrica intrínseca de las bandas (The Ring Averaging Limitation):

A pesar de cerrar la cuestión inicial y confirmar la persistencia del hallazgo, el uso de bandas de latitud completas (75°-90°) choca con una barrera analítica clara:

• Una banda de latitud rodea la Luna en 360° alrededor del eje de rotación, por lo que comprime y promedia en un único número estadístico todas las longitudes globales que contiene [A.5.6].
• La geometría de la banda no es capaz de indicar cómo se organiza y distribuye la densidad dentro de la propia celda polar [A.5.6]. Un valor medio de 300 por millón de km² puede proceder físicamente de varios escenarios morfológicos distintos: una acumulación puntual muy densa junto al vértice del polo, una estructura de anillo periférico a cierta distancia angular, una concentración asimétrica en un flanco del casquete, o una distribución uniforme y homogénea [A.5.6].
• El método analítico es comparable a determinar el número medio de personas por m² en una sala redonda: el dato global proporciona el grado de densidad general, pero no permite distinguir si el público está concentrado en el centro de la sala, repartido a lo largo de las paredes o distribuido de manera uniforme [A.5.6]. Resolver esta limitación y descifrar la estructura interna de los casquetes requiere pasar a la siguiente metodología de la investigación: el análisis de los casquetes centrados y el recuento alrededor de los seis ejes del cuerpo (A.6–A.9) [A.5.3].

A.5.7 — Fin de la fase de bandas y paso a los casquetes

En este punto, el trabajo analítico que las bandas de latitud y longitud pueden realizar en el sistema ha llegado a su término completo [A.5.7].

Síntesis de los resultados en una línea:

Se ha establecido la existencia de una densidad de cráteres alta y anómala en los dos extremos polares del perfil de latitud, completamente aislada del patrón espacial entre este y oeste, y que se mantiene estable y firme incluso tras eliminar los efectos de las cuencas y los maria [A.5.7].

Con un detalle metodológico más amplio, el análisis de bandas cumplió cuatro objetivos de control principales [A.5.7]:

1. Localización espacial: definición de la ubicación exacta de la anomalía de densidad respecto a los ejes de latitud.

2. Separación de patrones: aislamiento total del efecto polar respecto a la asimetría orbital a lo largo de las longitudes.

3. Respuesta geológica: examen cuantitativo de la variación de ambos patrones bajo la restricción de una máscara geológica fija y uniforme.

4. Equilibrio bipolar: revelación de que, tras eliminar los sesgos geológicos, se obtienen dos bandas polares (norte y sur) con valores de densidad similares, altos y equilibrados en el sistema.

Cambio de método analítico: el paso a la pregunta geométrica

A pesar de la persistencia de los resultados, la mecánica de las bandas choca con una barrera estructural y no puede responder a la siguiente pregunta científica [A.5.7]. Hasta este punto, las preguntas examinadas eran: ¿se caracterizan las latitudes polares por una mayor densidad de cráteres?, y ¿sobrevive este resultado a los procedimientos de control espacial [A.5.7]?

De aquí en adelante, el argumento científico cambia de naturaleza y se convierte en un problema geométrico y morfológico estricto [A.5.7]:

¿Cómo se dispone y organiza la densidad de cráteres dentro del propio espacio polar?

Para descifrar la estructura interna de la anomalía, la investigación cambia de forma fundamental el modo de medición y el corte geométrico de la corteza [A.5.7]:

• Supresión del promediado circunferencial: abandono de la división en bandas globales cerradas que rodean el cuerpo en 360° y promedian las longitudes [A.5.6, A.5.7].
• Método de casquetes centrados: definición de cada polo geográfico como el centro estructural exclusivo de una celda esférica (casquete polar) [A.5.7].
• Análisis del vector de distancia: examen y cartografía de la distribución de la densidad de cráteres en función directa de la distancia angular y radial desde el vértice del polo [A.5.7].

Por tanto, el capítulo siguiente —la sección A.6 (los casquetes polares)— no constituye una repetición estadística ni un recuento adicional de la misma población [A.5.7]. Representa la siguiente fase de desarrollo del argumento, que se deriva directamente del límite mecánico de la prueba de bandas: un paso nítido de la demostración de la mera existencia de una alta densidad polar al desciframiento de la geometría espacial y de la disposición interna de los cráteres dentro del casquete polar [A.5.7].

A.6 Los casquetes polares, sin máscara

En esta sección, la investigación cambia formalmente la geometría de la medición y pasa del uso de bandas circunferenciales al análisis de seis casquetes esféricos concéntricos (anidados) alrededor de cada polo, con radios angulares de: $10^{\circ },15^{\circ },20^{\circ },30^{\circ },45^{\circ },$ y $60^{\circ }$ [A.6]. El objetivo de la prueba es descifrar la estructura interna de la anomalía y examinar cómo varía la densidad acumulada a medida que el espacio analítico converge hacia dentro, hacia el punto extremo del polo de rotación [A.6].

Principios de la estructura geométrica de los casquetes:

1. Definición del radio angular ($\theta$): un casquete de radio $\theta$ incluye todos los cráteres cuya distancia de arco desde el vértice del polo cumple $\le \theta$ (por ejemplo: un casquete de $30^{\circ }$ se extiende desde la latitud $60^{\circ }$ hasta el $90^{\circ }$ del polo) [A.6].

2. Ecuación diferencial del área: el área de cada celda esférica se calcula con precisión según:

Para un radio $R=1,737.4\text{ km}$, el factor de área crece de forma no lineal: el casquete central reducido ($10^{\circ }$) abarca 0.288 millones de km², mientras que el casquete amplio ($60^{\circ }$) abarca 9.483 millones de km² (una extensión de área más de 30 veces mayor) [A.6].

3. Estructura acumulativa (anidada): las filas de esta tabla no representan zonas anulares aisladas e independientes, sino un perfil de suma acumulada [A.6]. El mecanismo es similar a los anillos de una diana de dardos: el casquete de $30^{\circ }$ no es solo el anillo exterior, sino toda la diana desde el centro ($0^{\circ }$) hasta $30^{\circ }$ [A.6].

• Si los cráteres están físicamente concentrados cerca del polo: los casquetes pequeños mostrarán valores de densidad altos, y la densidad disminuirá de forma constante a medida que el casquete se amplía e incorpora terreno ordinario [A.6].
• Si la distribución es homogénea: todos los casquetes anidados mostrarán valores de densidad similares, con independencia del radio [A.6].

4. Mantenimiento de la base de referencia estática: el punto de referencia (valor 1.00) permanece completamente fijo: la densidad bruta del cinturón ecuatorial completo sin máscara ($15^{\circ }\text{ S}-15^{\circ }\text{ N}$). Su valor es 166.436 (para $\ge 20\text{ km}$) y 53.170 (para $\ge 40\text{ km}$), lo que permite una comparación directa con los datos de las partes A.1–A.4 [A.6].

5. Sensibilidad de la muestra y factor de tamaño: los casquetes interiores ($10^{\circ }$) abarcan celdas de área pequeñas y contienen una cantidad reducida de cráteres, por lo que sus intervalos de confianza son amplios y conllevan la mayor incertidumbre estadística (el análisis completo se presenta en D.4) [A.6]. El grupo de control ($\ge 40\text{ km}$) representa una población de cráteres más antigua; una razón más alta en comparación con el grupo de 20 km indica una acentuación del efecto en los cráteres antiguos (se discute en A.8.5) [A.6].

Resumen de los datos de densidad de cráteres en seis casquetes anidados (A.6)

A continuación se presentan los valores numéricos y relativos medidos en los casquetes concéntricos alrededor de ambos polos, sin aplicar máscaras geológicas:

• Escalas de colores: en el grupo $\ge 20\text{ km}$ el intervalo va de 189 a 300; en el grupo $\ge 40\text{ km}$ el intervalo va de 63 a 108 cráteres por millón de km² [A.6].

Análisis morfológico del perfil radial bruto

De la lectura de los indicadores y del comportamiento gráfico de los dos hemisferios se desprenden los siguientes resultados:

1. Predominio del exceso polar global: ningún casquete polar, de ningún diámetro y en ningún hemisferio, muestra escasez respecto al cinturón ecuatorial [A.6]. La razón mínima del sistema es 1.14 (el casquete norte más amplio, de $60^{\circ }$), lo que demuestra que el interior de todos los casquetes presenta un exceso de densidad intrínseco respecto al ecuador bruto [A.6].

2. Comportamiento diferencial del polo norte: el perfil norte (la línea azul oscuro) muestra una convergencia y concentración mecánica clásica hacia el interior [A.6]. A medida que se reduce el radio de medición, la densidad asciende de forma continua: de 189 en el casquete amplio ($60^{\circ }$), pasando por 227 en el casquete de prueba medio ($30^{\circ }$), hasta un máximo de 300 cráteres por millón de km² en el casquete de $15^{\circ }$ [A.6]. En el casquete central extremo ($10^{\circ }$) el perfil se estabiliza en un valor casi idéntico (298), que da una razón de 1.79 [A.6]. En el grupo de control de los cráteres grandes ($\ge 40\text{ km}$), la tendencia de concentración septentrional se acentúa y alcanza una razón máxima de 2.02 en el casquete de $10^{\circ }$ [A.6].

3. Comportamiento plano del polo sur bruto: el perfil sur (la línea roja) muestra a primera vista una tendencia completamente distinta [A.6]. A lo largo de la mayor parte del intervalo de medición ($15^{\circ }$ a $60^{\circ }$), la curva de densidad permanece plana, horizontal y casi independiente de la distancia al polo, con valores constantes en el estrecho intervalo de 207–221 por millón de km² [A.6]. Solo en el casquete central más interior ($10^{\circ }$) se registra un aumento brusco hasta un valor de 250 (razón de 1.50) [A.6].

Límite de la interpretación geofísica en esta fase:

Los resultados empíricos brutos muestran aparentemente una profunda diferencia estructural entre los dos polos: el norte se caracteriza por una concentración interior marcada, mientras que el sur muestra una distribución plana y homogénea [A.6]. Sin embargo, la física del modelo reconoce que el perfil sur está cortado y promediado sobre el suelo de la cuenca gigante South Pole–Aitken, que borró una parte amplia del registro [A.6, A.11]. Para determinar qué parte de la variación procede de la deformación de la forma del cuerpo y qué parte es una distorsión debida a una única cuenca de impacto, debe aplicarse la máscara de filtrado geológica [A.6].

Leyenda de colores: transición continua del amarillo suave («densidad baja») al rojo suave («densidad alta»).

Figura 5. Densidad de cráteres de 20 km o más en casquetes de radio creciente alrededor de ambos polos (A.6), el catálogo completo. El norte sube hacia su polo; el sur es plano: el suelo de South Pole–Aitken, del que el casquete sur de 30° está compuesto en su mayor parte (A.3, D.4). (En el gráfico: eje horizontal — “radio del casquete alrededor del polo”, 10°, 15°, 20°, 30°, 45°, 60°; eje vertical — “cráteres ≥ 20 km por 10⁶ km²”, 150–300; la línea azul oscuro — “polo norte”; la línea roja — “polo sur”.)

A.7 Los casquetes polares tras excluir las cuencas y los maria

Esta sección examina si el comportamiento diferencial entre los dos hemisferios observado en la sección A.6 (perfil norte concentrado frente a perfil sur plano) es una propiedad de la población original de cráteres, o si se trata de una distorsión causada por el tipo de terreno y por los suelos geológicos renovados [A.7]. Para ello, se aplica la máscara de filtrado estricta (20 cuencas de Neumann de diámetro $\ge 500$ km y los polígonos de maria del LROC) a los seis casquetes concéntricos [A.7]. Este procedimiento garantiza una comparación precisa y de igual a igual sobre una superficie antigua y uniforme, fuera de los lugares de borrado térmico y mecánico [A.7].

Indicadores y calibración metodológica:

1. Cálculo del área y del numerador residuales: los cráteres cuyo punto central cae dentro de los dominios de la máscara se omiten, y el área de la máscara se descuenta con precisión del denominador [A.7]. La malla de medición se basa en celdas esféricas exactas con una resolución de $0.05^{\circ }\times 0.05^{\circ }$, eliminando los solapamientos una sola vez [A.7]. Allí donde las geometrías de corte coinciden (los casquetes de $30^{\circ }$ frente a las bandas de latitud $60^{\circ }-90^{\circ }$ en A.3), se obtiene una coherencia y un cierre matemático completos de los datos de área y de los numeradores [A.7].

2. Fijación de la regla de referencia constante: el punto de referencia (valor 1.00) permanece sin cambios: la densidad bruta del cinturón ecuatorial completo sin máscara ($166.436$ para 20 km y $53.170$ para 40 km) [A.7]. Esta elección es deliberada y sirve como regla estática constante para seguir cuantitativamente el grado de variación geológica de cada casquete, y no como indicador del exceso polar final [A.7].

3. Índice de cambio de densidad ($\Delta \rho$): compara cada casquete concéntrico consigo mismo en su estado bruto de A.6:

Un valor positivo indica un aumento de la densidad del área restante tras eliminar unidades pobres en cráteres [A.7].

Resumen de los datos de densidad de cráteres enmascarada en seis casquetes anidados (A.7)

A continuación se presentan los indicadores numéricos y los valores relativos obtenidos en los casquetes concéntricos enmascarados:

• Escalas de colores: en el grupo $\ge 20$ km el intervalo va de 232 a 348; en el grupo $\ge 40$ km el intervalo va de 74 a 135 cráteres por millón de km² [A.7].

Análisis geodinámico del perfil radial enmascarado

Del examen de los indicadores residuales y del comportamiento gráfico conjunto (Figura 6) se desprenden las siguientes conclusiones:

Figura 6. Los mismos seis casquetes alrededor de ambos polos tras excluir las cuencas y los maria (A.7), cráteres de 20 km o más. La línea norte no cambia en el intervalo de 20°; la línea sur sube en sus casquetes interiores, y ambas se aproximan. (En el gráfico: eje horizontal — “radio del casquete alrededor del polo”, 10°, 15°, 20°, 30°, 45°, 60°; eje vertical — “cráteres ≥ 20 km por 10⁶ km²”, 150–350; la línea azul oscuro — “polo norte”; la línea roja — “polo sur”.)

1. Efecto geológico nulo en el núcleo norte: los tres casquetes interiores del norte ($10^{\circ },15^{\circ },20^{\circ }$) muestran valores de cambio constantes de $0.0\%$ [A.7]. Este dato demuestra matemáticamente que no hay cuencas de Neumann ni maria basálticos en un radio de $20^{\circ }$ desde el vértice del polo norte; por tanto, el alto perfil central medido allí (300 por millón de km²) refleja con precisión el estado virgen de la corteza [A.7].

2. El efecto extremo en el núcleo sur y el borrado de la cuenca SPA: los casquetes más pequeños del sur muestran el efecto del borrado mecánico extremo de la cuenca SPA [A.7]. En el casquete sur de $10^{\circ }$, el área de medición se reduce en más de la mitad (caída de 0.288 a 0.104 millones de km²), lo que indica que la mayor parte del núcleo sur interior estaba cubierta por el suelo de la cuenca [A.7]. La limpieza de este suelo eleva la densidad real de la corteza montañosa residual en $+39.1\%$ en el grupo de 20 km (valor de 348, razón de 2.09) y en $+49.8\%$ en el grupo de 40 km (valor de 135, razón de 2.54) [A.7].

3. Convergencia morfológica hacia una curva de concentración uniforme: la comparación de la Figura 6 con la Figura 5 revela un vuelco morfológico [A.7]. Mientras que en el estado bruto la línea sur era completamente plana, la aplicación de la máscara de control igualó el tipo de corteza y elevó los casquetes interiores del sur [A.7]. Ahora la curva sur (roja) cambia por completo de forma y asciende de manera pronunciada hacia el interior, hacia el polo, de modo similar a la curva norte (azul oscuro) [A.7]. En cuanto ambos polos se cuentan sobre el mismo tipo de superficie, los dos muestran una única correlación firme: cráteres cada vez más densos a medida que nos acercamos al vértice del polo, hallazgo que confirma el modelo geodinámico de contracción de la superficie [A.7].

Aclaración metodológica sobre los valores del exceso limpio:

La columna de razones de esta tabla compara un numerador enmascarado con un denominador bruto (que conserva los maria ecuatoriales vacíos), por lo que sirve únicamente como regla de seguimiento [A.7]. Para derivar el exceso polar limpio real, el artículo aplica una comparación de igual a igual (eliminación simultánea de los maria y de la cuenca SPA en ambos lados) sobre los casquetes estándar de $30^{\circ }$ frente al cinturón ecuatorial ($|\varphi |<15^{\circ }$) [A.7]:

• Catálogo bruto completo (incluidos los maria): razón de 1.32 [A.7].
• Eliminación solo de los maria en ambos lados: razón de $1.11\pm 0.04$ [A.7].
• Eliminación de los maria y de la cuenca SPA en ambos lados (el recuento limpio): razón de $1.20\pm 0.05$ [A.7].

La desviación típica ($\pm 0.05$) refleja la fluctuación aleatoria del recuento; el intervalo de confianza de $95\%$ es el doble de amplio y establece con plena confianza estadística que el valor real de la razón se sitúa en el intervalo entre 1.10 y 1.30, valor situado claramente por encima de la base de distribución uniforme (1.00), y sobre él se apoya la conclusión científica de la investigación [A.7].

A.8 Síntesis y conclusiones de la prueba de casquetes

Las secciones A.6 y A.7 examinaron las regiones polares con una geometría concéntrica (anidada) esencialmente distinta del promediado espacial de las pruebas de bandas [A.8]. En lugar de medir la densidad media en anillos de latitud globales que rodean el cuerpo en 360°, cada polo de rotación se definió como el centro estructural de seis casquetes esféricos concéntricos de radios angulares $10^{\circ },15^{\circ },20^{\circ },30^{\circ },45^{\circ }$ y $60^{\circ }$ [A.5.6, A.8]. La prueba se realizó en una primera fase sobre el catálogo Robbins bruto sin máscaras, y en la fase siguiente bajo la restricción estricta de una máscara geológica (eliminación de las 20 grandes cuencas de Neumann y de los polígonos de maria del LROC) [A.8].

El cruce de los dos procedimientos separa metodológicamente dos preguntas centrales [A.8]:

1. ¿Cuál es la geometría espacial observada de la densidad de cráteres a medida que el vector de distancia converge hacia el punto extremo del polo de rotación [A.8]?

2. ¿Qué parte de esta geometría radial está inscrita en la corteza virgen antigua, y qué parte es una distorsión residual causada por una actividad endógena tardía (cuencas y maria) [A.8]?

El resultado principal establece que la aplicación de la máscara de control cambia de raíz la morfología estadística y el modo en que debe definirse físicamente la diferencia estructural entre los dos polos [A.8].

A.8.1 El estado bruto (sin máscara): asimetría de la morfología interna

En el recuento empírico inicial sin máscara (A.6), los dos hemisferios muestran dos perfiles radiales distintos y claramente diferenciados [A.8.1]:

• La concentración interior del polo norte: para el grupo de cráteres $\ge 20\text{ km}$, la densidad de cráteres asciende de forma constante a medida que se reduce el radio de medición: de 189 cráteres por millón de km² en el casquete amplio ($60^{\circ }$), pasando por 198 ($45^{\circ }$), 227 ($30^{\circ }$), 270 ($20^{\circ }$), hasta estabilizarse en un máximo de unos 300 y 298 en los casquetes centrales interiores ($15^{\circ }$ y $10^{\circ }$ respectivamente) [A.8.1]. El perfil no exige una definición monótona suave hasta el propio vértice (ya que los valores de 10° y 15° son estadísticamente idénticos), pero la tendencia espacial es firme: el núcleo interior es considerablemente más denso que su entorno amplio [A.8.1].
• Confirmación del grupo de control en el norte: el grupo de cráteres grandes ($\ge 40\text{ km}$) muestra un comportamiento coherente: la densidad asciende de 63 en el casquete amplio ($60^{\circ }$), pasando por 71 ($45^{\circ }$) y 87 ($30^{\circ }$), hasta un intervalo máximo de 100–108 por millón de km² en los casquetes centrales ($20^{\circ }-10^{\circ }$) [A.8.1]. Este resultado demuestra que la convergencia de los cráteres en el norte no es una distorsión estadística inherente únicamente al umbral de 20 km [A.8.1].
• El perfil plano del polo sur: en su estado bruto, el hemisferio sur muestra un comportamiento opuesto [A.8.1]. En el grupo de $\ge 20\text{ km}$, la curva de densidad permanece plana, amplia y casi independiente de la distancia al polo, manteniendo valores constantes de 207–221 por millón de km² a lo largo de todo el intervalo entre $15^{\circ }$ y $60^{\circ }$ [A.8.1]. Solo dentro del casquete central más reducido ($10^{\circ }$) se registra un salto aislado a 250 por millón de km² [A.8.1]. El grupo de $\ge 40\text{ km}$ refleja un perfil amplio similar (72–82 por millón de km² en el intervalo de $15^{\circ }-60^{\circ }$ frente a 90 en el casquete de $10^{\circ }$) [A.8.1].

Si el análisis científico se detuviera en esta fase (A.6), la única conclusión física sería que los dos polos difieren esencialmente tanto en la intensidad de su densidad absoluta como en su geometría estructural interna (norte concentrado frente a sur homogéneo y plano) [A.8.1]. Los datos de la sección A.7 demuestran que esta descripción es incompleta y procede de un sesgo geológico agudo [A.8.1].

A.8.2 Comportamiento del perfil norte bajo la restricción de la máscara

Uno de los resultados más claros e importantes de la aplicación de la máscara de control geológica (A.7) es la persistencia absoluta del núcleo estructural del polo norte [A.8.2]:

• Cambio nulo en el núcleo virgen septentrional: en los tres casquetes centrales interiores del hemisferio norte ($10^{\circ },15^{\circ },$ y $20^{\circ }$), la máscara de filtrado no modifica el número de cráteres ($N$) ni el área geométrica a nivel de cálculo [A.8.2]. Los valores de densidad permanecen completamente constantes en 298 (a 10°), 300 (a 15°) y 270 (a 20°) por millón de km² para el grupo de 20 km [A.8.2]. En el grupo de 40 km se conservan exactamente los mismos valores brutos de origen: 108, 101 y 104 por millón de km² respectivamente [A.8.2].
• Significado geofísico: este resultado demuestra que la convergencia de la densidad interior en el norte no es un artefacto de la presencia de cuencas o maria [A.8.2]. En un radio de $20^{\circ }$ desde el vértice del polo norte, la máscara no resta nada, pues ningún polígono de mare ni cuenca de Neumann se acerca tanto a estas coordenadas [A.8.2]. El núcleo norte es una corteza antigua virgen (Highlands) cuyo registro estaba limpio y normalizado desde el principio [A.8.2].
• Efecto de la máscara sobre el fondo periférico septentrional: el efecto del filtrado solo empieza a medirse cuando el radio de los casquetes se amplía más allá del núcleo [A.8.2]. En el casquete norte de $30^{\circ }$, la densidad de 20 km sube ligeramente de 227 a 233 (cambio de $\Delta \rho =+2.8\%$) [A.8.2]. En el casquete de $45^{\circ }$ asciende de 198 a 232 ($+16.7\%$), y en el casquete amplio de $60^{\circ }$ de 189 a 235 ($+24.2\%$) [A.8.2]. Esta tendencia se reproduce exactamente en el grupo de $\ge 40\text{ km}$ [A.8.2].

Así pues, la máscara modifica el fondo espacial que rodea la región polar norte, pero deja intacto el núcleo [A.8.2]. Antes de la máscara, el contraste estructural entre el núcleo ($15^{\circ }$) y el entorno ($60^{\circ }$) era extremo (300 frente a 189) [A.8.2]. Tras la máscara, mientras el núcleo se mantiene estable en 300, el entorno de $60^{\circ }$ subió matemáticamente a 235 por la limpieza de los maria situados en él [A.8.2]. Por tanto, la concentración hacia el interior sobrevive firmemente, pero su amplitud (Contrast) se reduce [A.8.2]. La máscara demuestra que parte del contraste bruto original se debía al sesgo de distribución de las superficies jóvenes fuera del núcleo, y no a una asimetría de impactos [A.8.2].

A.8.3 El vuelco geodinámico del perfil sur

A diferencia del norte, donde el núcleo interior era inmune a la máscara, en el hemisferio sur la máscara de filtrado modifica de forma fundamental la imagen morfológica, ya dentro de los casquetes más reducidos [A.8.3]:

• El núcleo sur de $10^{\circ }$: la aplicación de la máscara retira una porción enorme del área del casquete y de la población bruta de cráteres [A.8.3]. El número de cráteres de 20 km cae exactamente a la mitad (de 72 a 36) [A.8.3]. Sin embargo, dado que el área geométrica disponible se reduce relativamente mucho más (caída brusca de 0.288 a 0.104 millones de km²), la densidad residual de lo que queda se dispara un $+39\%$ y alcanza un valor máximo de 348 cráteres por millón de km² [A.7, A.8.3].
• El casquete sur de $15^{\circ }$: en el recuento bruto se contaron en este sector 139 cráteres ($\rho \approx 215$) [A.6]. Tras la limpieza de la máscara quedan 84 cráteres sobre un área residual reducida de 0.275 millones de km², y la densidad real de la corteza se dispara un $+42\%$ y fija un valor de 305 cráteres por millón de km² [A.7, A.8.3].
• Acentuación de la tendencia en los cráteres antiguos ($\ge 40\text{ km}$): en el casquete de $15^{\circ }$, el grupo de cráteres grandes muestra el cambio más agudo del sistema: la densidad residual salta de 82 a 124 por millón de km², un salto estructural de $\Delta \rho =+51\%$ [A.7, A.8.3].
• Atenuación del efecto de filtrado al ampliar el radio: en el casquete sur de $20^{\circ }$ se registra un aumento considerable (el de 20 km sube de 207 a 280, y el de 40 km de 73 a 98) [A.7, A.8.3]. A medida que el radio se amplía e incorpora terreno montañoso ordinario, el efecto de la máscara se atenúa: un aumento de alrededor del $17\%$ en el casquete de $30^{\circ }$, y una estabilización en un aumento moderado de alrededor del $11\%-12\%$ en los casquetes amplios de $45^{\circ }$ y $60^{\circ }$ [A.7, A.8.3].

El mecanismo geofísico:

En un radio de $20^{\circ }$ desde el vértice del polo sur, el único componente estructural de la máscara es la cuenca gigante South Pole–Aitken (SPA) [A.8.3]. El suelo de esta cuenca se caracteriza por una escasez extrema de cráteres debido al borrado mecánico completo del registro temprano [A.8.3]. La exclusión de las áreas de la cuenca resta matemáticamente una porción de área enorme (el denominador) y muy pocos cráteres (el numerador), por lo que revela la alta densidad real de la corteza montañosa virgen residual que la rodea [A.8.3]. Este resultado firme descarta por completo la posibilidad de leer el perfil sur plano original como el registro de una población de cráteres virgen; se trata de una distorsión local debida a una única cuenca de impacto, que se neutraliza con éxito [A.8.3].

A.8.4 Convergencia bipolar e igualdad estadística tras la máscara

En cuanto ambos polos se miden y calculan sobre el mismo tipo de superficie cortical (Highlands vírgenes fuera de los lugares de reinicio térmico y mecánico), la gradación morfológica cambia y muestra una imagen uniforme y coherente [A.8.4]:

• Comparación con el estado bruto sin máscara: la limpieza geológica hace desaparecer la asimetría aparente en la que los casquetes del sur mostraban valores bajos (215 frente a 300 en $15^{\circ }$, y 207 frente a 270 en $20^{\circ }$) [A.8.4].
• Cierre coherente con la prueba de bandas: estos resultados detallan y confirman geométricamente los resultados de la sección A.3 (en la que las bandas extremas $75^{\circ }-90^{\circ }$ se igualaron por completo tras la máscara) [A.8.4]. Los casquetes concéntricos demuestran que la convergencia bipolar no es producto del promediado de una banda circunferencial completa; penetra hacia el interior y fija una igualdad de densidad estrecha también dentro de los casquetes centrales interiores de $15^{\circ }$ y $20^{\circ }$ [A.8.4].
• Falta de coincidencia en los detalles del grupo de control ($\ge 40\text{ km}$): a pesar de la convergencia general, los polos no se vuelven completamente idénticos en cada indicador local [A.8.4]. En el grupo de $\ge 40\text{ km}$, en el casquete de $15^{\circ }$ se midieron 101 en el norte frente a 124 en el sur; a 10° se midieron 108 en el norte frente a 135 en el sur; y a 20° se midieron 104 en el norte frente a 98 en el sur [A.8.4]. Los valores mantienen firmemente el mismo orden de magnitud físico, pero su estructura geométrica local muestra una variación interna moderada [A.8.4].

Conclusión metodológica final:

El análisis estadístico completo de los intervalos de confianza ($95\%$) detallado en el capítulo D.4 establece que, dado que el valor de simetría perfecta (1.00) está contenido en los intervalos de incertidumbre de las muestras en todos los casquetes y para ambos umbrales de diámetro, los dos polos se definen como iguales y completamente equilibrados dentro de los límites del error aleatorio del recuento [A.8.4]. Por tanto, la formulación científica precisa no afirma que “la máscara hace que los polos sean completamente idénticos”, sino que: la máscara de filtrado geológica neutraliza y elimina la mayor parte del contraste bruto aparente que existía entre los perfiles interiores del norte y del sur, y los lleva a converger en intervalos de densidad cortical similar, alta y completamente equilibrada [A.8.4].

A.8.5 Comparación de los umbrales de diámetro lineales y factor de tamaño del cráter (reloj cronológico aproximado)

El uso diferencial de dos umbrales de diámetro estrictos ($\ge 20\text{ km}$ y $\ge 40\text{ km}$) constituye una herramienta de control metodológico esencial, que verifica que las tendencias geodinámicas observadas no son producto de una anomalía estadística inherente a una población de tamaño específica [A.8.5]:

• Sincronización morfológica en el norte: en ambos umbrales de diámetro se registró un perfil bruto idéntico de concentración hacia el interior (densidad mayor en el núcleo que en los casquetes amplios) [A.8.1, A.8.5]. Bajo la restricción de la máscara, el núcleo central ($10^{\circ }-20^{\circ }$) permanece inmune y sin cambio numérico en ambas muestras, mientras que el fondo periférico experimenta un aumento en ambas [A.8.2, A.8.5].
• Sincronización morfológica en el sur: ambos umbrales responden al colapso de la cuenca SPA con un salto brusco de densidad en los casquetes interiores [A.8.3, A.8.5]. En el casquete sur de $15^{\circ }$, se registra un aumento residual de $+42.0\%$ en el grupo de 20 km y un aumento agudo de $+50.7\%$ en el grupo de 40 km [A.7, A.8.5]. En el casquete de $20^{\circ }$ el aumento se estabiliza estrechamente en $+35.2\%$ y $+34.7\%$ respectivamente [A.7, A.8.5].

Esta coherencia en ambos hemisferios demuestra que la gradación morfológica general no depende de la elección del umbral de diámetro [A.8.5]. No obstante, se registra una diferencia constante de profundo significado físico: las razones respecto al cinturón ecuatorial son sistemáticamente más altas en el grupo de $\ge 40\text{ km}$ que en el grupo de $\ge 20\text{ km}$ [A.8.5]. En el casquete norte bruto de $30^{\circ }$, la razón es de 1.64 (a 40 km) frente a 1.36 (a 20 km); en el casquete sur enmascarado de $15^{\circ }$, se acentúa hasta 2.32 frente a 1.83 [A.6, A.7, A.8.5].

El mecanismo del reloj cronológico del diámetro del cráter:

En el marco del modelo físico, el diámetro de un cráter actúa como un reloj cronológico aproximado: los cráteres gigantes ($\ge 40\text{ km}$) se asocian estadísticamente a las fases más tempranas del registro de impactos, mientras que los cráteres más pequeños se acumularon a lo largo de la secuencia posterior [A.8.5].

• Si el exceso de densidad se debiera a un factor estático independiente del tiempo, las razones de densidades serían idénticas en todas las clases de tamaño [A.8.5].
• El hecho de que el exceso polar aumente en función directa del diámetro (como se confirma con mayor finura en el recuento limpio de D.2, en el que la razón de los casquetes de $30^{\circ }$ respecto al cinturón asciende de 1.09 en la clase de 20–30 km hasta 1.53 en la clase de 60–100 km) demuestra que el proceso de apiñamiento mecánico se produjo de forma simultánea y paralela a la continuación de la caída de la población de cráteres [A.7, A.8.5].

Cuanto antes se formó un cráter (cerca del máximo del achatamiento geométrico de la Luna joven), mayor fue la porción de la secuencia posterior de redondeo y disminución de volumen que experimentó, por lo que sufrió el máximo de la contracción mecánica de la superficie que convergió hacia los polos [A.8.5].

A.8.6 Integración estructural: cruce de los datos de bandas y casquetes

La combinación de los dos niveles de corte geométrico de la Parte A permite formar una imagen física coherente y aislada de la corteza del cuerpo [A.8.6]:

• La prueba de bandas (A.1–A.4): respondió a la pregunta sobre la mera existencia del fenómeno y estableció firmemente que las latitudes polares se caracterizan por un exceso de densidad real que permanece estable y equilibrado entre ambos hemisferios (norte y sur) tras limpiar los suelos basálticos jóvenes y las cuencas [A.8.6].
• La prueba de casquetes (A.6–A.7): respondió a la pregunta morfológica precisa y descifró cómo se organiza la densidad hacia el interior, hacia los puntos extremos de los polos [A.8.6]. Mientras que en el estado bruto se obtuvo una ilusión de asimetría (norte concentrado frente a sur plano), la aplicación de la máscara eliminó el sesgo de la cuenca SPA y elevó los valores del núcleo sur, llevando así los casquetes interiores de ambos polos a una convergencia morfológica similar (especialmente en el grupo de $\ge 20\text{ km}$) [A.8.1, A.8.3, A.8.6].

El cruce revela una distinción crítica entre dos tipos de información: el exceso polar en las latitudes altas es una propiedad física estable que sobrevive a la máscara, mientras que la variación geométrica inicial entre norte y sur es una distorsión local pasajera que no sobrevive al control [A.8.6]. Esta distinción evita un grave error metodológico y permite separar el estado actual de la corteza conservada de las propiedades primarias del fenómeno geodinámico [A.8.6].

A.8.7 Síntesis de las conclusiones empíricas de la primera fase

Tras el tratamiento de las secciones A.1–A.7, se establecen cinco conclusiones observacionales firmes que delimitan el fenómeno que el mecanismo físico debe explicar [A.8.7]:

1. Falta de uniformidad latitudinal significativa: la población global de cráteres no está distribuida de forma homogénea según la latitud; las bandas polares se sitúan de forma absoluta en el extremo superior del perfil de densidad [A.8.7].

2. Persistencia endógena robusta: el exceso de densidad polar no es función de la presencia de las 20 grandes cuencas ni de las superficies de maria, y sobrevive por completo a la aplicación de la máscara geológica [A.8.7].

3. Simetría bipolar residual: la limpieza de los sesgos geológicos revela valores de densidad similares y equilibrados entre ambos polos, tanto a nivel de la media global de las bandas como dentro de la estructura interna de los casquetes interiores ($\ge 20\text{ km}$) [A.8.7].

4. Inestabilidad de la asimetría norte-sur: la diferencia morfológica entre la curva norte concentrada y la curva sur plana colapsa bajo el control, y se demuestra como producto de un impacto de cuenca local en el sur y no como una propiedad primaria [A.8.3, A.8.7].

5. Coherencia entre umbrales: la tendencia general se reproduce y se observa simultáneamente tanto en el umbral de $\ge 20\text{ km}$ como en el de $\ge 40\text{ km}$, y no depende de la elección de un único diámetro [A.8.7].

A.8.8 Límites de la prueba de control y fronteras del conocimiento en la fase actual

Para mantener el rigor científico, la investigación define tres límites analíticos estrictos para la interpretación de los datos en esta fase [A.8.8]:

• Limitación de la corrección estática de la máscara: la aplicación de la máscara geológica es únicamente una herramienta de control y filtrado espacial, y no pretende realizar una reconstrucción paleogeográfica completa de la superficie lunar primitiva [A.8.8]. No corrige mecánicamente los procesos graduales de renovación de la corteza, no reconstruye cráteres completamente borrados y no constituye una prueba absoluta de que la corteza montañosa residual haya conservado íntegramente la población original [A.8.8].
• Limitación de la medición acumulada de los casquetes: dado que los casquetes geométricos están anidados unos dentro de otros (el casquete de $30^{\circ }$ incluye matemáticamente los casquetes estrechos de $20^{\circ },15^{\circ },$ y $10^{\circ }$), un cambio en el valor de densidad acumulado entre dos filas no refleja directamente la densidad local dentro del anillo aislado entre ellas [A.8.8]. Determinar la ubicación exacta del cambio radial requerirá en el futuro un análisis de anillos independientes, pero no es necesario para establecer los resultados acumulados presentados [A.8.8].
• Límite de resolución de las muestras reducidas (incertidumbre aleatoria): a medida que el radio del casquete converge hacia el interior, la celda de área se reduce y el número de cráteres cae, especialmente en el umbral de $\ge 40\text{ km}$ [A.8.8]. En el núcleo de $10^{\circ }$ se contaron solo 86 cráteres en el norte y 36 en el sur [A.8.8]. Debido a la muestra reducida, los intervalos de incertidumbre (como se detalla en la matriz de $95\%$ del capítulo D.4) se amplían considerablemente [A.8.4, A.8.8]. Por tanto, está científicamente vedado utilizar pequeñas fluctuaciones numéricas entre 10° y 15° como prueba de una estructura interna fina sin el respaldo de pruebas de significación estadística adecuadas [A.8.8].
• Aislamiento de los componentes de la máscara en el núcleo sur: en el estrecho radio de $20^{\circ }$ desde el vértice del polo sur, la cuenca SPA es el único componente de la máscara, y en él existe una coincidencia completa entre los datos de la máscara total y los de la cuenca aislada [A.8.3, A.8.8]. La estructura interna particular de la población de cráteres dentro de la propia cuenca se separa y analiza aparte en el capítulo A.10 [A.8.8].

A.8.9 Fin de la prueba de casquetes y definición del fenómeno científico

Esta sección resume y cierra formalmente la pregunta morfológica que condujo del análisis de bandas a la fase de casquetes anidados, y consolida la base de datos virgen de la corteza residual [A.8.9].

Síntesis de la dinámica analítica:

• Limitación de la primera fase: las pruebas de bandas (A.1–A.4) demostraron que las latitudes polares son más densas y que el efecto sobrevive a la máscara de filtrado, pero debido al promediado circunferencial en $360^{\circ }$ no podían indicar cómo se organiza la densidad dentro del propio espacio polar [A.5.6, A.8.9].
• Desciframiento de la segunda fase (sin máscara): la sección A.6 cubrió esta laguna y reveló una asimetría geométrica aguda [A.8.9]. El norte mostró una concentración clásica hacia el interior (caída de 300 en el núcleo de $15^{\circ }$ a 189 en el casquete de $60^{\circ }$ para el grupo de 20 km) [A.8.1, A.8.9]. El sur, en cambio, mostró un perfil plano, horizontal y homogéneo (valores constantes de 207–221 por millón de km² en la mayor parte del intervalo) [A.8.1, A.8.9]. Detenerse en esta fase habría llevado a la conclusión errónea de que los dos polos difieren geométricamente de raíz [A.8.1, A.8.9].
• El vuelco de la tercera fase (aplicación de la máscara): la sección A.7 reveló que la asimetría morfológica era consecuencia de sesgos geológicos [A.8.9]. En el núcleo norte ($10^{\circ }-20^{\circ }$), la máscara no restó nada y demostró que la corteza allí era virgen y limpia desde el principio (conservación de los valores de 300 y 270 por millón de km²) [A.8.2, A.8.9]. En el sur, la máscara penetró hasta los casquetes más pequeños; la exclusión del suelo pobre de la cuenca SPA borró una porción relativa de área (el denominador) mayor que la cantidad de cráteres omitidos (el numerador) [A.8.3, A.8.9]. En consecuencia, la densidad residual en el sur se disparó: de 215 a 305 en el casquete de $15^{\circ }$, y de 207 a 280 en el casquete de $20^{\circ }$ [A.7, A.8.9].

Colapso de la asimetría y convergencia de los datos residuales

La limpieza de los registros geológicos condujo a una reducción drástica de las diferencias bipolares y a una convergencia estructural de los valores de densidad dentro de los núcleos interiores [A.8.9]:

El resultado acotado y preciso establece que la mayor parte de la asimetría geométrica observada inicialmente entre los casquetes interiores no es estable y se disipa en condiciones de control uniformes [A.8.9]. Este resultado complementa y confirma los resultados de las bandas: a nivel de bandas, la máscara igualó los sectores de $75^{\circ }-90^{\circ }$ en valores de 300 y 305; a nivel de casquetes, la misma convergencia firme se reproduce y no desaparece aunque el espacio analítico se reduzca gradualmente hacia el vértice polar [A.3, A.7, A.8.9]. Las dos geometrías de corte distintas apuntan exactamente al mismo vector físico: el alto exceso de densidad no es una característica local de un solo polo [A.8.9].

Definición de los límites del control mecánico de la cuenca SPA

La investigación mantiene una separación científica estricta entre la descripción empírica y el mecanismo mecánico [A.8.9]:

1. Identificación de los componentes de la máscara: la máscara de filtrado se define como la unión de los polígonos de maria y las 20 grandes cuencas de Neumann [A.7, A.8.9].

2. Papel geodinámico de la cuenca SPA: la cuenca Polo Sur–Aitken (South Pole–Aitken) es el componente gigante principal que penetra profundamente en el espacio polar sur, y su geometría morfológica proporciona el fundamento físico de que la máscara actúe y eleve los valores en el sur de forma anómala respecto al norte virgen [A.8.3, A.8.9].

3. Aislamiento de componentes paso a paso: el cálculo enmascarado total no constituye una prueba de la cuenca SPA por sí sola [A.8.9]. Por ello, está metodológicamente vedado atribuir todo el vuelco morfológico del sur exclusivamente a la SPA sin aislar su contribución residual de los demás componentes de la máscara [A.8.9]. Esta separación se realiza y se confirma en dos planos estrictos:

• Dentro del radio de $20^{\circ }$ del núcleo sur: la SPA es el único componente enmascarado del sistema, por lo que allí existe una coincidencia numérica absoluta entre ella y la máscara general [A.8.8, A.8.9].
• En la fase de la prueba física final (A.11 y D.1): se aplica el procedimiento explícito de “recuento limpio”, que omite de forma activa y equilibrada las áreas de los maria y la cuenca SPA en ambos lados de la comparación simultáneamente (tanto de los casquetes como del cinturón ecuatorial) [A.7, A.8.9].

La base de datos consolidada: definición de la anomalía cortical

Al término de la prueba de casquetes, se define el fenómeno físico virgen que las demás partes del artículo deben explicar mediante las ecuaciones de la geodinámica y el modelo de contracción mecánica [A.8.9]:

• ${\mathbf{\Delta }\rho }_{\text{polar}}>1.00$ — exceso de densidad bipolar significativo, presente en ambos polos simultáneamente
• $\rho \propto \frac{1}{\theta }$ — apiñamiento radial continuo que se acentúa hacia el vértice de cada polo de rotación
• ${\rho }_{\text{North}}{\equiv \rho }_{\text{South}}$ — igualdad mecánica equilibrada entre ambos polos dentro de los límites del error aleatorio del recuento
• ${\mathbf{\Delta }\rho }_{\ge 40}{>\mathbf{\Delta }\rho }_{\ge 20}$ — coherencia entre umbrales, con un exceso estructural que se acentúa en las clases de diámetro grandes

Sobre la base de estos cuatro indicadores centrales, la investigación pasa ahora a la fase de integración geométrica y de prueba física de los ejes del cuerpo [A.8.9].

A.9 Los seis ejes y la definición de la geometría triaxial de la Luna joven

Esta sección amplía los límites del análisis y plantea una prueba geométrica combinada. Mientras que las secciones A.1–A.8 examinaron la superficie a través de una sola dirección de recuento cada vez (latitud o longitud por separado), esta sección combina los datos sobre la base de los tres ejes naturales de un cuerpo con acoplamiento de marea [A.9]:

1. El eje de rotación (polar): pasa por ambos polos. Este eje gobierna todos los fenómenos afectados por los esfuerzos centrífugos, por el achatamiento rotacional y por su relajación geodinámica (el mecanismo central propuesto para la reducción y el apiñamiento de la superficie en los polos) [A.9].

2. La línea de la Tierra (radial): va del punto subterrestre (la cara cercana, orientada permanentemente hacia la Tierra) al punto antiterrestre (centro de la cara oculta). Este eje gobierna los esfuerzos de marea externos y la interacción mecánica de dos cuerpos (como la zona de contacto primitiva con la pera del cuerpo progenitor) [A.9].

3. La dirección de la órbita (tangencial): va del punto delantero (el ápice del movimiento en el espacio, $90^{\circ }\text{ W}$) al punto de cola (la retaguardia orbital, $90^{\circ }\text{ E}$). Este eje refleja asimetrías dinámicas derivadas del movimiento del cuerpo a través del vuelo de los fragmentos o materiales eyectados [A.9].

La metodología se basa en comparar la densidad dentro de seis casquetes de igual área con un radio angular de $30^{\circ }$ (un radio de unos 910 km sobre la superficie) alrededor de los seis extremos de estos ejes, tras aplicar la máscara de filtrado geológica estricta (eliminación de las 20 cuencas gigantes de Neumann y de los polígonos de maria del LROC) [A.9]. La igualdad de áreas intrínseca permite una comparación cuantitativa directa y libre de sesgos de dimensión [A.9].

Resumen de los datos de densidad de cráteres en los seis casquetes principales (A.9)

A continuación se presentan los indicadores numéricos exactos, ordenados estrictamente según la gradación de la densidad residual, del casquete menos denso al más denso del sistema [A.9]:

• Base de normalización fija: el índice de razón (1.00) se calcula de forma coherente respecto a la densidad bruta del cinturón ecuatorial completo sin máscara del capítulo A.1 (valor de 166.436 para 20 km y 53.170 para 40 km) [A.9].
• Valores brutos en el catálogo completo (antes de la máscara): los valores numéricos originales utilizados para construir la Figura 7 son 94 (subterrestre), 129 (delantero), 211 (polo sur), 216 (antiterrestre), 227 (polo norte) y 253 (de cola) cráteres por millón de km² [A.9].

Figura 7. Los seis ejes (A.9): densidad de cráteres de 20 km o más en casquetes de 30° alrededor de los puntos subterrestre, delantero, polo sur, antiterrestre, polo norte y de cola. Los polos se sitúan con la cara oculta; el punto subterrestre es el menos denso. Las columnas corresponden al catálogo completo, antes de la máscara; leídas de izquierda a derecha, ascienden desde el punto que miraba hacia la Tierra hasta la cara de cola. (En el gráfico: eje vertical — “cráteres ≥ 20 km por millón de km²”; columnas de izquierda a derecha — subterrestre 94, delantero 129, polo sur 211, antiterrestre 216, polo norte 227, de cola 253; leyenda inferior — “casquete de 30° alrededor de cada uno de los seis puntos principales, catálogo completo”.)

Análisis dinámico del efecto de la máscara y confirmación del mínimo subterrestre

De la lectura de la columna de cambio de densidad ($\Delta \rho$) y de la gradación estadística se desprende una respuesta clara a una objeción metodológica previsible [A.9]:

• La objeción científica inicial: dado que la cara cercana alberga la inmensa mayoría de las superficies de maria jóvenes vacías de cráteres, podría argumentarse que el valor mínimo extremo del punto subterrestre no es más que un reflejo indirecto de la presencia de los maria, y nada más [A.9].
• Refutación del sesgo geológico: la aplicación de la máscara descarta esta hipótesis [A.9]. El casquete del punto subterrestre obtiene efectivamente el aumento más pronunciado del sistema tras la limpieza de los suelos basálticos ($\Delta \rho =+57.8\%$), pero, a pesar del salto teórico, permanece claramente en la base de la gradación como el casquete menos denso de la superficie lunar ($148$) [A.9]. No son los maria los que producen el mínimo subterrestre; está inscrito físicamente en la propia corteza [A.9].
• Comportamiento de otros casos extremos: el casquete de cola apenas se ve afectado por el filtrado geológico ($-0.3\%$), lo que indica que su registro estaba limpio desde el principio [A.9]. En cambio, el punto delantero muestra un aumento considerable del +34.7% (de 129 a 174) [A.9]. El examen a nivel de los componentes de la máscara muestra que una parte decisiva de la escasez del casquete delantero en su estado bruto se debía al borrado térmico de la cuenca aislada Orientale, cuya limpieza eleva el valor residual a 162 [A.9].

Lectura de la gradación empírica como registro de un cuerpo triaxial achatado

La gradación numérica conservada tanto antes como después de la máscara (subterrestre $\to$ delantero $\to$ antiterrestre / polos $\to$ de cola) demuestra que la distribución de los cráteres no está gobernada por un único eje, sino que es reflejo de una geometría triaxial de un cuerpo deformado [A.9]. El modelo geofísico interpreta el patrón espacial del siguiente modo [A.9]:

1. Equilibrio del eje de rotación (polos frente a ecuador): los polos presentan un exceso de densidad marcado sobre los puntos ecuatoriales (subterrestre y delantero), como dicta el modelo de redondeo y contracción del achatamiento rotacional [A.9]. Tras la máscara, la media de los dos polos es 1.12 veces superior a la media de las dos caras ecuatoriales (delantera y de cola) [A.9].

2. Asimetría del eje de la línea de la Tierra (cabeza libre frente a disco de contacto): los dos extremos del eje de la línea de la Tierra muestran un comportamiento completamente distinto [A.9]. El extremo antiterrestre (la cara oculta) presenta un exceso de densidad alto ($221$) y se iguala por completo con el nivel de los polos [A.9]. En cambio, el extremo subterrestre (la cara cercana) es la región menos densa del sistema ($148$) [A.9].

• El mecanismo físico: esta diferencia concuerda exactamente con la estructura morfológica de la Luna joven mientras estaba acoplada a la Tierra [A.9]. El extremo antiterrestre actuó como una “cabeza libre” (la cima de la cúpula libre) y no sufrió contracción ni borrado mecánico durante el redondeo, por lo que se sitúa sistemáticamente al nivel de los polos [A.9]. En cambio, el extremo subterrestre coincidía con el centro del disco y del sistema de contacto físico de la estructura de la pera con la Tierra, por donde fluía una presión de origen elevada que provocó el empobrecimiento del registro original por razones mecánicas no relacionadas con la forma global del cuerpo (analizado en A.9.2) [A.9].

La combinación de estos esfuerzos describe un cuerpo sometido a un achatamiento rotacional simétrico alrededor del eje polar ($c$) y, al mismo tiempo, extremadamente estirado a lo largo del eje radial ($a$) entre la gravedad terrestre hacia dentro y la fuerza centrífuga hacia fuera [A.9]. Esta estructura define una forma triaxial precisa que cumple la gradación física: un eje largo en dirección a la Tierra ($a$), un eje intermedio a lo largo del vector orbital ($b$) y un eje corto y comprimido a través de los polos ($c$) [A.9]. La forma teórica de un cuerpo de marea libre (C.2) reproduce exactamente una gradación de razones idéntica (razones de 1.66, 1.56 y 1.06 para polos frente a subterrestre, caras frente a subterrestre y polos frente a caras, respectivamente) [A.9].

Determinación de los indicadores de la forma geométrica a partir del registro de cráteres

En lugar de suponer una forma hipotética, las ecuaciones de redondeo y disminución de volumen (C.3) permiten leer las dimensiones originales del cuerpo directamente a partir del registro empírico del exceso de cráteres [A.9]:

• Reconstrucción de la forma de la población general ($\ge 20\text{ km}$): el ajuste de los valores de exceso de densidad en los dos extremos de la línea de la Tierra (que combinan el efecto de cabeza libre con valores normalizados, y muestran razones de exceso de 1.12 y 1.30 frente a las caras y al subterrestre) da una forma triaxial original de $1.37:1.15:1$ en las fases tardías del redondeo [A.9].
• La prueba de integración independiente: el redondeo teórico de la forma $1.37:1.15:1$ hacia una esfera perfecta predice de forma completamente independiente una razón de densidades limpia de 1.21 entre los casquetes polares de $30^{\circ }$ y el cinturón ecuatorial [A.9]. Este valor predicho concuerda de forma absoluta y sin ajustes con las razones de densidades empíricas medidas en la superficie en el recuento limpio del capítulo D.1 ($1.19-1.20$) [A.9]. Dos regiones geográficas completamente distintas de la corteza lunar se leen y se concilian con éxito mediante una única ecuación geométrica [A.9].
• Reconstrucción de la forma en el máximo cronológico a partir de los cráteres grandes: dado que el diámetro del cráter es un reloj cronológico aproximado, los cráteres más grandes reflejan las fases de deformación más tempranas y extremas [A.8.5, A.9]. El tratamiento de los datos de los seis ejes para las clases de tamaño de 60–100 km y 60–300 km da formas de alrededor de $2.1:1.2:1$ y alrededor de $2.0:1.3:1$ [A.9]. A pesar de la incertidumbre local derivada de una muestra reducida en cada casquete (solo 13–68 cráteres), la combinación analítica conjunta de los seis ejes y del cinturón limpio en la matriz de dispersión estadística (las tablas detalladas en C.3) converge exactamente alrededor del valor de forma del máximo de achatamiento del sistema: $1.8:1.4:1$ (la forma del esferoide de Maclaurin con una rotación de 5.5 horas) [A.9].

A.9.1 La asimetría orbital: la cara de cola frente a la cara delantera

Los dos extremos del eje de la dirección orbital muestran una diferencia estadística coherente y aguda en todos los indicadores de análisis espacial aplicados en la Parte A; la dirección de la anomalía se mantiene firme sin inversión de signo, mientras que la amplitud cuantitativa varía según la geometría de corte [A.9.1]:

• El grupo de control (≥ 40 km): en los casquetes de 30° se midieron valores de 84 frente a 47 por millón de km² en el estado bruto, y de 92 frente a 61 en el estado enmascarado, dato que confirma la estabilidad del patrón en los cráteres antiguos [A.9.1].
• Descarte del mecanismo de barrido frontal (The Shielding/Sweep Effect): los datos empíricos invierten la expectativa física simple: un cuerpo planetario que se mueve en el espacio a través de una corriente de partículas o polvo debería recibir una densidad de impactos mayor en el frente de su vuelo (la cara delantera), como el parabrisas de un coche que circula a gran velocidad y recoge insectos en comparación con la ventanilla trasera [A.9.1]. La Luna muestra una anomalía completamente opuesta, ya que la cara delantera es la más dispersa y vacía del sistema, hallazgo que descarta una captura espacial directa y exige un mecanismo geofísico interno y dinámico [A.9.1].
• Explicación de la diferencia entre casquetes y bandas (efecto de dilución espacial): el colapso de la razón enmascarada en las bandas (de 1.54 a 1.07) frente a su supervivencia en los casquetes (1.45) se debe a que una banda de longitud de 30° se extiende de forma continua de polo a polo, con lo que promedia y diluye la anomalía local con áreas de corteza ordinaria en latitudes altas [A.9.1]. En cambio, el casquete delimita con precisión la celda de área relevante con un radio de unos 910 km alrededor del punto central [A.9.1]. Por tanto, se establece que el exceso no es una propiedad global de toda la longitud, sino una concentración geográfica localizada en torno a las coordenadas de 90° este [A.9.1].

A.9.2 La correlación geológica: distribución de las superficies de maria y basalto

El aislamiento de la anomalía local en torno a 90° este plantea la cuestión de la variación morfológica: ¿qué existe en esta celda de área que no existe en la cara opuesta [A.9.2]? Un corte preciso de los polígonos de maria del LROC con la geometría de los casquetes de 30° da los siguientes indicadores cuantitativos [A.9.2]:

• Inventario de las estructuras geológicas en cada flanco:
  • En el flanco de cola (90° E): se concentran las estructuras Mare Marginis, Mare Smythii, Mare Undarum, Mare Spumans, junto con pequeñas aberturas en los cráteres Artamonov, Joliot, Hubble, Lomonosov [A.9.2].
  • En el flanco delantero (90° W): se concentran las estructuras del sistema de la cuenca Orientale (Mare Orientale, Lacus Veris) y los márgenes occidentales extremos de Oceanus Procellarum [A.9.2].
• Neutralización del sesgo de la cuenca de contacto (Oceanus Procellarum): la estructura Procellarum representa una llanura basáltica gigante (corteza fundida por el calor irradiado por la Tierra, sin un domo de alimentación magmático propio) que pertenece geográficamente a la zona de contacto de la Luna joven con el cuerpo progenitor, y no a la cara delantera orbital [A.9.2]. El centro de esta zona de contacto (con una anchura de base de 3,000 km) se sitúa en 24.1° N, 335.2° E, y contiene la cuenca Mare Imbrium, la abertura de alimentación principal desde la Tierra [A.9.2]. Dado que el área de Procellarum penetra en el casquete delantero únicamente por un estiramiento geométrico, su omisión revela la razón real de las aberturas de las caras: la cara de cola presenta un exceso de 2.36 veces en el área de maria propios y de 2.10 veces en el número de aberturas respecto a la cara delantera [A.9.2].
• Identidad de los indicadores vírgenes: si se limpian del casquete de cola los márgenes de los maria de la cara cercana que se han extendido hasta él (Mare Crisium y la llanura de Fecunditatis), el área residual es de 170,839 km² frente a 87,336 km² en el delantero [A.9.2]. Este dato da una razón de áreas de aberturas de 1.96, valor absolutamente idéntico a la razón de densidades de cráteres bruta original en los mismos casquetes (1.96) [A.9.1, A.9.2]. Esta proximidad numérica constituye una observación firme que dicta una relación física directa: la cara que alberga el exceso de cráteres es exactamente la cara que alberga el exceso de maria propios [A.9.2].

A.9.3 La ecuación del balance morfológico: los maria como aberturas de eyección y de alimentación

Hay que resolver la siguiente paradoja científica: ¿por qué el sector más rico en cráteres es el sector que alberga el exceso de maria, cuando los propios maria son suelos jóvenes vacíos de cráteres [A.9.3]? La respuesta se basa en la mecánica del modelo LEF [A.9.3]:

1. El mecanismo de eyección y retorno (Artículo 6): en el marco del modelo, cada mare es una abertura cortical sobre un domo de alimentación magmático activo [A.9.3]. Estas aberturas expulsaron y lanzaron hacia arriba cantidades ingentes de material [A.9.3]. El material eyectado conservó su energía cinética y su momento angular, y describió trayectorias balísticas que regresaron y cayeron sobre la corteza a velocidades de 2 a 5 km/s (exactamente la velocidad con la que había salido de la abertura) [A.9.3]. Estos impactos se produjeron en un radio de 50 a 1,000 km alrededor de la abertura de origen en las fases tempranas, y se redujeron más tarde con el redondeo y la compactación del cuerpo [A.9.3].

2. El efecto bidireccional de una abertura de eyección: el mare actúa sobre el registro de la corteza con dos vectores opuestos y simultáneos: produce cráteres de impacto a su alrededor mediante el material retornado, mientras que su propio suelo permanece caliente, funde y reinicia los cráteres impresos en él (el suelo se congeló y se fijó solo tras la separación) [A.9.3]. El registro geométrico de la superficie se expresa con precisión mediante la ecuación del balance morfológico:

Fuera de los límites de los maria, el término de borrado térmico se anula mientras que el término del material retornado es máximo; dentro de ellos, el término de borrado domina y anula el numerador [A.9.3].

Asimetría de la acción de la máscara:

La ecuación explica por qué la aplicación de la máscara (A.3) modifica los valores de forma no simétrica [A.9.3]:

• En el flanco de cola (90° E): la densidad residual permanece completamente constante (un movimiento menor de 253 a 252) [A.9.3]. Este resultado demuestra que la superficie fuera de los maria en este sector ya estaba saturada y densamente poblada de cráteres exactamente como el casquete entero, ya que constituye la zona principal de caída del material retornado [A.9.3]. La máscara no reconstruye lo que se borró en los maria; solo aísla el área sólida de caída [A.9.3].
• En el flanco delantero (90° W): la densidad residual salta de 129 a 174 [A.9.3]. Este aumento se debe a que la máscara limpió de él las áreas vacías de Procellarum y del suelo de Orientale, restando así un denominador de área grande en relación con el reducido numerador de cráteres del sector [A.9.3].

A.9.4 Integración de los mecanismos de apiñamiento de la cara de cola

El exceso de cráteres conservado en la cara de cola se explica mediante la combinación de dos mecanismos físicos independientes que actúan en la misma dirección dinámica [A.9.4]:

1. El mecanismo de dispersión local del material retornado:

La escala de las trayectorias balísticas del material (hasta unos 1,000 km) concuerda exactamente con las dimensiones geométricas de los casquetes de 30° (≈ 910 km), y explica por qué el efecto se conserva en el casquete y colapsa en las bandas largas (que promedian el área con terreno virgen al que no llegaron los fragmentos) [A.9.1, A.9.4]. La relación mecánica se define en la siguiente cadena lógica estricta [A.9.4]:

Mayor concentración de aberturas (maria) en la estructura → mayor volumen de material eyectado y retornado → densidad de cráteres anómala en la corteza circundante

El cálculo cuantitativo completo de la historia de las eyecciones, la producción de cráteres y las geometrías de vuelo se desarrolla y detalla en la Parte E.3 del artículo y en la Parte K del Artículo 6 [A.9.4].

2. El mecanismo de rezago rotacional dinámico (Lag Effect):

Durante la fase acoplada, la Luna mantuvo el acoplamiento rotacional y giró simultáneamente con la Tierra [A.9.4]. El material lanzado con fuerza hacia arriba desde aberturas de eyección situadas en otras zonas de la superficie del cuerpo experimentó un rezago físico respecto a la rápida velocidad de rotación del sistema (periodo de 5.5 horas) [A.9.4]. La dinámica funciona como lanzar una pelota directamente hacia arriba desde el suelo de un tiovivo en movimiento: el cuerpo avanza y el material que cae aterriza, en sentido figurado, detrás de su punto de partida [A.9.4]. Por tanto, la retaguardia orbital (la cara de cola) actuó como una trampa dinámica de momento, que captó y atrajo hacia sí cantidades residuales de material expulsado desde todas las partes del cuerpo que aterrizaron detrás, además del material retornado desde sus propias aberturas [A.9.4].

A.10 La cuenca Polo Sur–Aitken (SPA) como objeto espacial independiente

Durante el tratamiento de los datos bipolares, se estableció que la cuenca South Pole–Aitken no es solo un componente de filtrado pasivo; posee una estructura interna propia y una falta de uniformidad considerable en la distribución de la población en su interior [A.10]. La cuenca se extiende a lo ancho de unos 2,400 km sobre la cara oculta y penetra profundamente en el casquete polar sur (cubre el 24% del área de todos los casquetes polares) [A.10].

La investigación descarta el uso de una única media global para la cuenca, ya que aplicar un único número monolítico produciría una grave distorsión metodológica: una media combinada de un centro completamente vacío y unos márgenes montañosos ordinarios presentaría la ilusión cortical de una cuenca moderadamente pobre y homogénea en cada punto [A.10]. Por ello, para mantener un rigor académico estricto, la sección A.10 aísla la estructura de la cuenca y realiza una serie de mediciones radiales continuas que se extienden de forma escalonada desde su punto central geométrico hacia fuera, hacia los márgenes, con el fin de cerrar las cuestiones abiertas del registro de cráteres del Artículo 7 [A.10].

A.10.1 Metodología del análisis radial de la cuenca South Pole–Aitken (SPA)

Para descifrar la falta de uniformidad estructural interna de la cuenca South Pole–Aitken, las coordenadas del centro geométrico de la cuenca se fijaron según los mismos parámetros utilizados para construir la máscara de filtrado: $53^{\circ }$ sur, $191^{\circ }$ este, con una definición estricta del diámetro del borde exterior de 2,400 km [A.10.1].

• Radio angular aproximado: un diámetro de 2,400 km da un radio físico de 1,200 km, cuya división por el radio lunar ($1,737.4\text{ km}$) da un ángulo de arco de $39.57^{\circ }$ desde el centro de la cuenca hasta sus márgenes [A.10.1].
• Estructura de las celdas de recuento (división en anillos): en lugar de promediar la cuenca en un único valor que distorsionaría la realidad, el espacio se dividió en cuatro anillos esféricos interiores y dos anillos de control exteriores: $0-10^{\circ },10-20^{\circ },20-30^{\circ },$ y $30^{\circ }-39.57^{\circ }$ (hasta el borde); y más allá del borde: $39.57^{\circ }-50^{\circ }$ y $50^{\circ }-60^{\circ }$ [A.10.1].
• Zonas aisladas y cálculo del área esférica: a diferencia de la serie de casquetes anidados, aquí los anillos no están anidados y cada fila representa una celda de área aislada e independiente [A.10.1]. El área se calcula según la ecuación espacial:

El factor $\cos {\theta }_1-\cos {\theta }_2$ compensa la curvatura de la esfera y hace que un anillo de anchura $10^{\circ }$ alejado del centro tenga mayor área que un anillo de igual anchura cercano a él [A.10.1].

• Regla de referencia local: el valor de referencia (1.00) de este capítulo se define como la densidad medida en el primer anillo de control situado inmediatamente fuera del borde de la cuenca ($39.57^{\circ }-50^{\circ }$) [A.10.1]. Esta elección es local y deliberada: el terreno residual fuera del borde se encuentra en las mismas latitudes y en la misma cara oculta del cuerpo, con lo que neutraliza y elimina los efectos del exceso polar y de la asimetría longitudinal de los capítulos A.1–A.9, y aísla exclusivamente la anomalía geológica de la propia cuenca [A.10.1].

Resumen de los datos de densidad de cráteres en anillos alrededor del centro de la SPA (A.10)

A continuación se presentan los indicadores numéricos exactos obtenidos en los seis anillos independientes, calculados para el grupo de $\ge 20\text{ km}$ y el grupo de $\ge 40\text{ km}$ [A.10.1]:

• Escalas de colores numéricas: para $\ge 20\text{ km}$ el intervalo va de 132 a 222; para $\ge 40\text{ km}$ el intervalo va de 52 a 75 cráteres por millón de km² [A.10.1]. Los valores relativos se calculan a partir de datos de recuento y área sin redondear [A.10.1].

A.10.2 La zona central como foco del mínimo (estructura en escalón)

Del análisis de las columnas de densidad y del perfil espacial (Figura 8) se desprenden los siguientes resultados científicos principales:

Leyenda de colores: del amarillo («bajo») al rojo («alto»).

Figura 8. Densidad de cráteres de 20 km o más en seis anillos separados alrededor del centro de South Pole–Aitken (A.10). La zona central es el mínimo del perfil; la densidad sube ya en el anillo de 10–20° y se aproxima a la del entorno hacia el borde. La parte sombreada es el interior de la cuenca; la línea discontinua es su borde. (En el gráfico: eje vertical — “cráteres ≥ 20 km por millón de km²”; eje horizontal — “distancia al centro de la SPA (puntos medios de los anillos)”, con marcas en 0°, 10°, 20°, 30°, 50°, 60° y “borde”; leyendas — “dentro de la SPA”, “fuera de la SPA”; sobre la curva roja — “132 — el mínimo” y 222.)

• El disco del mínimo del núcleo: el mínimo absoluto del registro se sitúa dentro del anillo central geográfico ($0-10^{\circ }$), que se extiende como un recinto amplio sobre un área de unos 288,000 km² (y no como un punto singular) [A.10.2]. Dentro de esta celda de área se contaron solo 38 cráteres de diámetro $\ge 20\text{ km}$, lo que da una densidad de 132 cráteres por millón de km² [A.10.2].
• El vector del salto periférico: fuera del disco central, la densidad no asciende con una pendiente lineal moderada (rampa), sino que muestra un comportamiento mecánico de escalón marcado [A.10.2]. Ya en el paso al segundo anillo ($10-20^{\circ }$), la densidad salta bruscamente a 193 y se estabiliza en valores casi iguales y altos de 200 y 202 en el resto del interior de la cuenca ($20^{\circ }$ hasta el borde) [A.10.2]. Al cruzar el borde de 2,400 km, la densidad residual sube a 222 por millón de km² [A.10.2].
• El factor de déficit cuantitativo: la densidad del núcleo central constituye una razón de solo 0.59 respecto al terreno de control fuera del borde, dato que refleja un déficit relativo extremo de alrededor del 41% en el registro de cráteres del centro de la cuenca [A.10.2].

A.10.3 Calibración del grupo de control ($\ge 40\text{ km}$) e intervalos de incertidumbre

El grupo de cráteres grandes ($\ge 40\text{ km}$) reproduce exactamente la tendencia morfológica y confirma la coherencia de la dirección del resultado [A.10.3]:

• La densidad en el núcleo es de 52 y asciende a un intervalo constante de 68–71 en las partes exteriores de la cuenca, hasta alcanzar un valor de 75 fuera del borde [A.10.3]. La razón del centro respecto al entorno residual se fija en 0.69 [A.10.3].
• Límite de resolución de la muestra reducida: en el núcleo central se contaron solo 15 cráteres de diámetro $\ge 40\text{ km}$ [A.10.3]. Dado que el error de dispersión aleatoria y la fluctuación natural de un recuento crecen como la raíz cuadrada del recuento ($\sqrt{N}$), el peso estadístico relativo del error aumenta a medida que la muestra se reduce [A.10.3]. Una muestra de 15 cráteres conlleva una incertidumbre mucho más amplia que una muestra de 38 cráteres [A.10.3]. Por tanto, el grupo de 40 km sirve únicamente como herramienta de control cualitativa para confirmar el vector de la curva, mientras que la determinación cuantitativa se apoya en el grupo de 20 km [A.10.3].

A.10.4 Implicaciones dinámicas para el procedimiento de recuento limpio (A.11)

Los resultados empíricos demuestran que definir la cuenca SPA como un bloque único bajo el término monolítico “déficit” falta al rigor científico, ya que la escasez de cráteres está limitada y concentrada físicamente solo dentro del recinto central [A.10.4]. Los tres anillos exteriores del interior de la cuenca ($10^{\circ }$ hasta el borde) muestran valores de densidad considerablemente más altos que el centro (193–202) y próximos al nivel del fondo montañoso ordinario del entorno [A.10.2, A.10.4].

Esta conclusión estructural es esencial para comprender el procedimiento de recuento limpio de los capítulos A.11 y D.1 [A.10.4]:

• En el recuento limpio equilibrado, la cuenca SPA se elimina por completo del cálculo porque su centro es pobre y distorsiona los datos [A.10.4].
• La aplicación de este filtro monolítico tiene un coste estadístico: excluye del cálculo el centro pobre (el término de borrado necesario) y, al mismo tiempo, omite con él los tres anillos exteriores, que no son pobres en cráteres y conservan un registro cortical normal [A.10.4].
• La demostración de que el déficit se limita al corazón de la cuenca delimita el alcance de la conclusión observacional: dentro de la cuenca SPA existe un amplio recinto central en el que la densidad de cráteres es baja y fue claramente reiniciada, mientras que las partes periféricas y el terreno de los márgenes sobrevivieron sin un borrado agudo [A.10.4]. Las cuestiones de la geometría no circular detallada del recinto, la definición de sus límites físicos exactos y el desciframiento de su origen geodinámico exceden el objetivo del Artículo 7, y la explicación completa de su formación se traslada a los capítulos del Artículo 8 [A.10, A.10.4].

A.10.5 Interfaz de la investigación: el punto de contacto con el Artículo 8

El hallazgo morfológico de la existencia de un recinto central pobre en cráteres dentro de la cuenca South Pole–Aitken surgió de forma independiente de la prueba de control de los cráteres polares en el presente artículo (Artículo 7) [A.10.5]. Por tanto, la identificación de la anomalía y su cuantificación numérica forman parte inseparable del registro de la observación empírica de esta investigación [A.10.5].

No obstante, el establecimiento de este mínimo central plantea al sistema una nueva cuestión física y geodinámica: ¿por qué el núcleo geográfico interior de la cuenca SPA presenta una densidad de cráteres tan anómalamente baja en comparación con los anillos periféricos y las áreas exteriores de los márgenes de la cuenca [A.10.5]?

La división estructural de funciones entre las partes de la serie científica

Para mantener una estructura académica sobria y contenida, se ha establecido una división de funciones definida y estricta entre las partes de la serie, que delimita el alcance del Artículo 7 [A.10.5]:

• La función del Artículo 7 (el presente artículo): la función analítica de este artículo se apoya en este punto y concluye en él por completo [A.10.5]. El artículo identifica el fenómeno morfológico, lo cuantifica con precisión en anillos separados y establece una afirmación observacional demostrada: la amplia zona central de la cuenca SPA constituye el punto mínimo absoluto del perfil radial de densidad de cráteres medido [A.10.5]. El Artículo 7 no pretende ni intenta proponer un mecanismo físico para descifrar este enigma [A.10.5].
• La función del Artículo 8 (investigación de la cuenca SPA): la resolución de la cuestión mecánica y dinámica pasa exclusivamente al Artículo 8, en el que la estructura South Pole–Aitken se convierte en el objeto central de investigación y análisis [A.10.5]. El Artículo 8 desarrollará y examinará los siguientes datos complementarios:
  1. La geometría detallada y no circular de la zona pobre (indicadores de elipticidad y asimetría estructural de la cuenca) [A.10.5].
  2. Las coordenadas exactas del punto de mínima densidad absoluta [A.10.5].
  3. Posibles relaciones antipodales (Antipodal Relations) sobre la superficie del cuerpo lunar [A.10.5].
  4. Los mecanismos geofísicos y dinámicos capaces de explicar la formación de la estructura anómala de la cuenca [A.10.5].

Aclaración metodológica estricta

Todos los valores y datos numéricos presentados en este capítulo se derivaron y procesaron de forma coherente a partir del catálogo Robbins (2019), que sirve como base de datos fija a lo largo de toda la Parte A [A.10.5]. Se subraya estrictamente que la definición del límite de $10^{\circ }$ del núcleo central es únicamente un elemento de división analítica para definir el anillo interior, y no se presenta ni se interpreta como un límite físico o geológico demostrado de la zona de déficit de cráteres [A.10.5].

A.11 Metodología del recuento limpio y límites de la medición empírica

Esta sección reúne las reglas mecánicas y los criterios estructurales en los que se basan todos los recuentos de la Parte A, y define los términos de la ecuación de la razón principal y oficial del artículo [A.11]. El análisis pretende normalizar la superficie, separar los sesgos geológicos y enmarcar las limitaciones metodológicas de la base de datos [A.11].

1. Definición de las regiones de comparación (la geometría espacial):

• Los casquetes polares (las celdas objetivo): se definen de forma continua y estricta desde la latitud $60^{\circ }$ hasta $90^{\circ }$ (tanto en el norte como en el sur) [A.11].
• El cinturón ecuatorial (la celda de referencia): una banda global estrecha definida en el intervalo de latitudes de $|\varphi |<15^{\circ }$ alrededor del ecuador [A.11]. Este cinturón es el sector más estrecho que contiene una muestra de cráteres suficientemente amplia; produce el contraste (Contrast) más nítido y limpio del sistema, ya que allí el vector de redondeo y disminución de volumen de la Luna joven alcanza su máximo mecánico (los polos son los puntos de encuentro críticos de los meridianos, mientras que el ecuador es la región más alejada de ellos) [A.11].
• Calibración del control angular: la ampliación del cinturón ecuatorial al intervalo de $|\varphi |<30^{\circ }$ (azul claro en la Figura 9) da un resultado global casi idéntico ($1.12$ frente a $1.11$), pero diluye y enmascara el contraste estructural, ya que incorpora latitudes medias en las que el efecto de contracción es moderado [A.11]. La concordancia numérica entre los dos cinturones demuestra que la mera existencia de la anomalía es robusta y no depende de una elección arbitraria de la anchura del sector [A.11].

Figura 9. Dónde se realizan los recuentos: los casquetes polares (naranja) frente al cinturón ecuatorial de ±15° (azul oscuro; el cinturón de ±30° en azul claro, a efectos de comparación). Los maria (gris) se eliminan en ambos lados; South Pole–Aitken (contorno rojo) se elimina del recuento limpio. Puntos: cráteres de 60 km o más.

El procedimiento de recuento limpio equilibrado (eliminación de sesgos)

Para realizar una comparación justa y neutral, la física exige comparar unidades de corteza con una historia evolutiva similar (terreno antiguo frente a terreno antiguo) [A.11]. Para ello se aplica una limpieza simultánea y equilibrada de dos superficies geológicas que sesgan los registros en direcciones opuestas [A.11]:

1. El sistema de maria y llanuras basálticas (gris): suelos jóvenes en los que el registro original fue borrado y reiniciado térmicamente [A.11]. Estas áreas se concentran en latitudes bajas y representan el 22% del cinturón ecuatorial, pero solo el 1% de los casquetes polares [A.11]. Su inclusión en el recuento bruto rebaja artificialmente la densidad del denominador ecuatorial e infla erróneamente la razón de los polos (razón bruta de 1.32) [A.7, A.11].

2. La cuenca Polo Sur–Aitken (contorno rojo): una estructura de impacto gigante que borró mecánicamente la población temprana de cráteres en el corazón de su núcleo (A.10) [A.11]. La cuenca se encuentra en la cara oculta y cubre el 24% del área de todos los casquetes polares (principalmente en el sur) [A.10, A.11]. Su inclusión en el recuento bruto rebaja artificialmente el numerador polar y sesga la razón a favor del ecuador [A.11].

3. El balance normalizado: restar estos dos elementos de ambos lados de la ecuación a la vez neutraliza los dos sesgos opuestos y deja una corteza montañosa virgen y de historia estable (Highlands) tanto en los polos como en el ecuador [A.11]. En este recuento limpio y equilibrado se apoya firmemente la conclusión del artículo (D.1–D.2), que fija el valor de la razón del sistema en $1.20\pm 0.05$ [A.7, A.11].

Los cinco límites del conocimiento y las limitaciones analíticas del recuento

Para mantener un estándar académico sobrio y contenido, la investigación define explícitamente cinco reglas básicas que limitan y delimitan lo que un recuento de cráteres es capaz de demostrar [A.11]:

• La regla de acumulación histórica (Accumulative Record): el registro de cráteres visible hoy en la superficie es una función acumulada [A.11]. Todo cráter impreso desde que la capa comenzó a congelarse y a recoger su registro sobre la Tierra se ha conservado en ella, a menos que una superficie posterior (un mare o una cuenca) lo haya borrado físicamente [A.11]. El recuento moderno refleja la suma de las poblaciones de todas las épocas juntas, y no una única fase cronológica; el exceso de densidad polar se analiza como un patrón mecánico añadido sobre este fondo total [A.11].
• Límite del umbral de diámetro y riesgo de saturación (Saturation Limit): la completitud de la identificación y la calidad de los datos varían drásticamente en función del tamaño del cráter [A.11]. Por debajo del umbral de 40 km, la completitud del catálogo se ve afectada por el tipo de terreno y por la intrusión de cráteres secundarios (cúmulos de fragmentos posteriores al impacto) [A.11]. Por ello, la clase de $\ge 20\text{ km}$ se utiliza para derivar los valores numéricos, mientras que la clase de $\ge 40\text{ km}$ sirve como herramienta de control para confirmar la dirección del signo [A.11]. Por debajo del umbral estricto de 20 km no se realiza ninguna comparación, ya que allí dominan los cráteres secundarios y la superficie se encuentra en estado de saturación (un equilibrio en el que los cráteres nuevos borran a los antiguos al mismo ritmo y el recuento deja de crecer) [A.11].
• Neutralización de los sesgos ópticos y de la iluminación solar: en los polos, los rayos del Sol inciden con un ángulo bajo sobre el horizonte y proyectan sombras largas y enormes que realzan cada surco y cada cráter pequeño, lo que produce la ilusión visual de una superficie mucho más craterizada que el cinturón ecuatorial (efecto de atardecer) [A.11]. Para bloquear este sesgo óptico, la investigación descarta el uso de imágenes o fotografías y se basa exclusivamente en datos topográficos y mapas de elevación digitales que miden el terreno real y no sus sombras [A.11].
• Presencia de fondo de la cuenca de contacto: la anomalía del Artículo 11 (el mínimo de la cara cercana en la zona de contacto) está presente como fondo básico bajo todo el conjunto presentado [A.11]. El exceso polar se mide respecto a ella y no la sustituye; el modelo de los seis ejes (A.9) presenta y equilibra ambos patrones simultáneamente [A.11].
• Obligación de uniformidad de la base de datos: todos los datos de densidad proceden del catálogo geométrico global de Robbins (2019), caracterizado por una resolución uniforme sobre toda la Luna [A.11]. La investigación descarta el uso de mapas geológicos clásicos, compilados a partir de fuentes históricas con un detalle desigual (varias veces más densos en la cara cercana que en la oculta), para evitar una distorsión comparativa [A.11, D.3]. Una comparación científica válida solo existe cuando ambas mitades del cuerpo se cuentan con el mismo método y el mismo instrumento de medición [A.11].

Cierre de la fase descriptiva y paso a las pruebas explicativas

Con la conclusión de las secciones A.1–A.11, queda cerrada por completo la parte descriptiva del artículo [A.11]. Las cinco etapas espaciales establecieron la siguiente base de datos:

• Identificación de un exceso de densidad marcado en las latitudes altas [A.11].
• Separación y aislamiento completos del patrón polar respecto a la asimetría longitudinal (este/oeste) [A.11].
• Confirmación de la persistencia del efecto bipolar incluso bajo la restricción de una máscara de filtrado geológica estricta [A.11].
• Desciframiento de la estructura interna de los casquetes y demostración de la convergencia de los núcleos hacia la igualdad estadística [A.11].
• Integración de las tendencias en una gradación geométrica coherente que confirma la existencia de un cuerpo triaxial primitivo [A.11].
• Aislamiento y cuantificación del mínimo central de la cuenca SPA y su traslado al análisis del Artículo 8 [A.11].

De aquí en adelante, el argumento científico cambia de rostro y pasa de la cuestión de la descripción empírica (qué hay en la superficie y dónde) a la cuestión del mecanismo físico (qué es capaz de producir un patrón así) [A.11]. Todo modelo o hipótesis alternativa que se examine en los capítulos siguientes del artículo tendrá que cumplir seis condiciones límite estrictas establecidas por la observación, sin posibilidad de excepción [A.11]:

• 1. Densidad de cráteres alta y anómala en ambas regiones polares simultáneamente
• 2. Reproducción coherente de la tendencia en los dos umbrales de diámetro aislados ( 20 km y 40 km)
• 3. Supervivencia y persistencia completas del exceso de los extremos tras aplicar la máscara de maria y cuencas
• 4. Explicación de la variación diferencial en la respuesta del norte frente al sur a la máscara (colapso de la cuenca SPA en el sur)
• 5. Convergencia morfológica hacia valores de densidad similares en los núcleos interiores residuales de los casquetes
• 6. Conservación de la gradación de los seis ejes, con el punto subterrestre fijado como el casquete menos denso

Parte B — Lo que intentamos antes de la respuesta

La Parte A definió cuantitativa y espacialmente la anomalía de densidad: los cráteres lunares son significativamente más densos en las regiones de ambos polos que en el cinturón ecuatorial, y el efecto se organiza de forma global en torno a un anillo completo de la superficie y no como un único punto singular [B]. La Parte C desarrollará y presentará las ecuaciones geodinámicas del mecanismo de redondeo y contracción mecánica propuesto como solución científica [B].

Este capítulo (Parte B) constituye la fase de prueba mecánica preliminar, y presenta las seis hipótesis e interpretaciones alternativas examinadas y rechazadas durante la investigación [B]. Su descarte mecánico es lo que condujo la investigación a la solución del redondeo, ya que el análisis conjunto de los seis fracasos demuestra que la física del sistema deja lugar a un único tipo de causa generadora [B].

La matriz de las cuatro preguntas estructurales:

Cada hipótesis alternativa se examina de forma coherente y sobria frente a cuatro preguntas estrictas, cuyos valores se derivan directamente de los datos de observación de la Parte A [B]:

1. Prueba del signo (Sign Test): ¿crea realmente el mecanismo físico propuesto un exceso de densidad en los polos respecto al ecuador, o actúa en sentido contrario y los empobrece [B]? (Un mecanismo que empobrece los polos queda descartado de inmediato) [B].

2. Prueba de la magnitud (Magnitude Test): si el signo dinámico es correcto, ¿es capaz el mecanismo de producir las razones empíricas medidas —$1.32$ en el recuento bruto, $1.20\pm 0.05$ en el recuento limpio—, junto con la tendencia creciente significativa a medida que aumenta el diámetro del cráter (D.1–D.2) [B]?

3. Prueba de la forma (Morphology Test): ¿sitúa y organiza el mecanismo el exceso de densidad en el espacio geográfico exacto en que se observa —presente en ambos polos en medida similar y ordenado según la gradación triaxial de los seis ejes (A.9)— [B]?

4. Prueba de la máscara (Mask Test): ¿es capaz el exceso de densidad producido por el mecanismo de sobrevivir a la aplicación de la máscara geológica (eliminación de los maria, las grandes cuencas y la cuenca SPA en ambos lados de la ecuación), como sobrevive de hecho en los recuentos residuales de la Parte A [B]?

Las pruebas se ordenan en una gradación metodológica desde la más gruesa (signo) hasta la más fina (máscara), de modo que la mayoría de las hipótesis se deciden y descartan en fases tempranas del análisis [B].

B.1 Conservación en tierras altas antiguas (Highlands Preservation Hypothesis)

Definición de la interpretación:

Esta hipótesis sostiene que el exceso de densidad en los polos no es más que el reflejo del grado de conservación de la superficie y de una edad geológica variable [B.1]. Las regiones polares son una corteza montañosa antigua y virgen (Highlands) que no experimentó la actividad de inundación volcánica tardía de las llanuras basálticas y los maria oscuros, que borraron y cubrieron el registro temprano en las latitudes bajas de la cara cercana [B.1]. Por tanto, los polos son más densos porque son más antiguos y no se vieron afectados por alteraciones endógenas, y no porque algún mecanismo dinámico concentrara o reuniera activamente cráteres en ellos [B.1]. Es la explicación inicial aceptada en la comunidad planetaria, ya que la densidad de cráteres sirve como indicador clásico para determinar la edad relativa de las superficies [B.1].

La predicción morfológica de la hipótesis (qué predice):

El modelo exige que la curva de densidad siga linealmente la edad cronológica y el grado de conservación de la celda de área [B.1]. La densidad debe ser completamente uniforme y homogénea en todas las unidades de corteza montañosa no inundadas, y mostrar valores bajos (depresiones) únicamente donde se extiende un mare o una cuenca joven [B.1]. La hipótesis descarta la existencia de un centro o vértice anómalo, ya que una meseta antigua de la misma edad es uniforme en el espacio [B.1]. Si la conservación térmica fuera el único factor, el polo mostraría valores idénticos a los de cualquier otro espacio montañoso de la misma edad en la superficie del cuerpo, y el mapa global de densidad estaría compuesto por mesetas planas y depresiones locales, sin una curva de máximo concentrado [B.1].

Análisis de los resultados empíricos de los recuentos (qué muestra la superficie):

El efecto de conservación geológica y el impacto térmico de los maria son completamente reales, y de hecho moldean la superficie lunar a la mayor escala (como se refleja en las zonas de exposición del Artículo 11, donde la cara oculta presenta una densidad bruta de 219 cráteres por 10⁶ km², el cinturón marginal 164 y el disco de contacto subterrestre solo 99) [B.1]. Sin embargo, tres observaciones firmes de la Parte A descartan la interpretación y conducen a su rechazo [B.1]:

1. Persistencia del exceso bajo la prueba decisiva de la máscara: la máscara de filtrado de la sección A.3 omite del cálculo precisamente los elementos en los que se apoya la hipótesis: los polígonos de maria y las cuencas gigantes [B.3, B.1]. Este procedimiento normaliza el denominador y garantiza una comparación limpia de meseta montañosa frente a meseta montañosa en ambos lados de la ecuación [B.1]. Si las diferencias de conservación fueran la única causa del fenómeno, la razón de densidades residual caería exactamente al valor de distribución uniforme ($1.00$) al eliminar la corteza sustituida [B.1]. En la práctica, la razón bruta cae de $1.32$ a $1.11$ (limpieza de maria) y vuelve a subir a un valor firme de $1.20\pm 0.05$ en cuanto se neutraliza también la cuenca SPA [B.1]. La razón no se reduce a 1 y sigue mostrando una tendencia creciente en función del diámetro del cráter hasta un valor de $1.53$ en las clases grandes [B.1]. Ninguna diferencia de densidad que sobreviva a este procedimiento de igual a igual puede atribuirse a la diferencia entre terreno inundado y terreno virgen [B.1].

2. El máximo bipolar por encima de la meseta más antigua del sistema: la corona de la cara oculta (alrededor de las coordenadas del punto antiterrestre) representa en el modelo dinámico la unidad de corteza más antigua y segura, que las superficies de maria apenas tocaron geográficamente [B.1]. Tras aplicar la máscara, este sector presenta una densidad de 221 cráteres por 10⁶ km² [B.1]. En cambio, los dos casquetes polares residuales muestran valores considerablemente más altos: 233 en el norte y 246 en el sur [B.9, B.1]. Los polos se sitúan morfológicamente por encima de la meseta más antigua, y no a su mismo nivel [B.1]. Si la edad cronológica de la corteza fuera el factor determinante, la corona de la cara oculta debería encabezar la gradación global; los resultados de la superficie lo descartan [B.1].

3. Existencia de una concentración radial (centro) en el núcleo norte virgen: dentro del espacio del polo norte, donde ningún polígono de mare ni cuenca de Neumann penetra en un radio de $20^{\circ }$ desde el vértice, la densidad no es homogénea ni plana, sino que muestra una curva de concentración marcada [B.1]. El índice de densidad asciende de forma continua desde 189 en el casquete amplio ($60^{\circ }$) hasta un valor máximo de 300 cráteres por millón de km² en el casquete interior de $15^{\circ }$ [B.6, B.1]. Una corteza montañosa uniforme y estática no es capaz de producir una convergencia radial de la densidad hacia un punto geográfico definido [B.1]. El núcleo norte constituye el lugar de observación más limpio para refutar la hipótesis, ya que no contiene ningún elemento geológico ajeno que requiera filtrado con la máscara [B.1].

Veredicto:

La hipótesis de conservación en tierras altas antiguas supera únicamente la prueba del signo, pero fracasa por completo en la prueba de la máscara y en la prueba de la forma [B.1]. La amplia variación espacial a escala global (como la gradación cara cercana frente a cara oculta) se explica mediante la geometría de las zonas de contacto del extremo de la pera (Artículo 3) y la distribución de los vectores de las caídas balísticas más duras, concentradas en la cara oculta (Artículo 6), y no en función de una conservación estática [B.1].

La investigación mantiene el mecanismo de conservación únicamente como la causa física de que las áreas de maria y las cuencas de impacto jóvenes muestren valores de densidad bajos (reinicio del registro), y no como fuente o explicación del exceso de densidad polar [B.1]. La presencia de las tierras altas antiguas de la cara cercana al sur del cráter Tycho, que muestran una densidad anómala de 302 cráteres en el intervalo de $15^{\circ }$, constituye un recordatorio mecánico adicional de que la densidad de una superficie antigua no es función únicamente de la edad cronológica [B.1]. Lo que la hipótesis de conservación deja sin respuesta física es precisamente el núcleo de la investigación del artículo: un exceso de densidad significativo conservado en la comparación de meseta frente a meseta, y que posee un centro geométrico anclado en el eje de rotación del cuerpo [B.1].

B.2 Achatamiento rotacional durante el acoplamiento (Rotational Flattening Hypothesis)

Definición de la interpretación:

Esta hipótesis se basa en el principio mecánico de que un cuerpo viscoso/blando que gira rápidamente alrededor de su eje experimenta una deformación achatada: se expande a lo largo del cinturón ecuatorial por los esfuerzos centrífugos y se acorta a lo largo del eje polar, como la masa de una pizza que se estira y se aplana cuando se la hace girar en el aire [B.2]. Dado que la Luna joven en la fase acoplada giraba en rápida sincronización con la Tierra (periodo de unas 5.5 horas), experimentó un achatamiento global marcado (C.1) [B.2]. La idea inicial sostiene que estos esfuerzos de compresión empujaron físicamente las regiones polares hacia dentro, acercando así mecánicamente entre sí los cráteres que transportaban y aumentando su densidad [B.2].

La predicción morfológica (qué predice):

La geometría física de una deformación achatada a volumen constante (situación en la que no se añade ni se resta material al sistema, sino que solo cambia su configuración) muestra un comportamiento completamente opuesto a la predicción de apiñamiento [B.2]:

• En el espacio polar: las superficies superior e inferior se aplastan y obligan al material a extenderse lateralmente, lo que provoca un estiramiento bidireccional de la superficie en el plano del polo [B.2].
• En el espacio ecuatorial: la superficie se estira a lo largo de la circunferencia del ecuador, pero se acorta a lo largo de los meridianos [B.2].

En consecuencia, el número de cráteres por unidad de área (${\rho }_{\text{crater}}$) disminuye y se empobrece en la región polar, y la ecuación diferencial de la extensión de la superficie se define con precisión según la razón de ejes del esferoide ($a/c$, donde $a$ es el radio ecuatorial y $c$ es el radio polar) [B.2]:

Esta ecuación dicta que cuanto más achatado está un cuerpo, más dispersos y vacíos de cráteres por unidad de área quedan sus polos (exactamente lo inverso de la ecuación de redondeo desarrollada en C.3) [B.2]. Para las razones de ejes del esferoide de Maclaurin en el máximo de achatamiento ($a/c=1.69$, densidad 1,300 ${\text{kg/m}}^3$, rotación de 5.5 horas) y en las fases tardías del redondeo ($a/c=1.20$, densidad 2,600 ${\text{kg/m}}^3$, periodo de 5.77 horas), el mecanismo de achatamiento habría dejado los polos en una posición de inferioridad teórica de solo 0.59 y 0.83 respecto al ecuador (déficit estructural de $41\%$ y $17\%$ respectivamente) [B.2].

Lo que dicen los recuentos empíricos:

Los datos físicos de la superficie muestran un exceso de densidad marcado en los polos (razón de $1.32$ en el recuento bruto y $1.20\pm 0.05$ en el recuento limpio) [B.2]. Por tanto, el signo mecánico de la hipótesis es erróneo de raíz; el achatamiento rotacional empobrece los cráteres en el vértice, mientras que las observaciones de la superficie los muestran con densidad máxima [B.2].

Veredicto:

La hipótesis del achatamiento rotacional fracasa por completo en la prueba del signo [B.2]. Si se comprime una pelota elástica desde arriba y desde abajo, la superficie del extremo superior no se contrae: se extiende y se expande, y los puntos impresos en ella se alejan unos de otros [B.2].

No obstante, la hipótesis capta correctamente la geometría general del cuerpo cuando se lee en sentido inverso: una superficie cortical que se estiró a su pesar durante la fase de achatamiento, y que más tarde queda liberada de la restricción del esfuerzo, se contraerá bruscamente en las regiones polares [B.2]. Esta comprensión inversa constituye el núcleo mecánico del capítulo C.3 [B.2]. Este hallazgo provocó un cambio metodológico fundamental durante la investigación: en lugar de buscar un mecanismo dinámico que produzca achatamiento, comenzó la búsqueda del proceso geodinámico que anula y libera el achatamiento (el proceso de redondeo de vuelta a la esfera) [B.2].

B.3 El eje polar acortado por la marea — efecto Roche (Roche Tidal Stretching Hypothesis)

Definición de la interpretación:

Durante la fase acoplada, la Luna joven estaba sometida al efecto Roche: esfuerzos de marea extremos que estiraron su volumen a lo largo de la línea de la Tierra (el eje radial $a$) debido a la diferencia dinámica entre la gravedad terrestre, que tiraba de la cara cercana hacia dentro, y la fuerza centrífuga de la corrotación, que tiraba de la masa hacia fuera, mientras el interior del cuerpo permanecía sujeto al extremo de la estructura de la pera [B.3]. Dado que el volumen del cuerpo permanecía casi constante, el alargamiento del eje radial se produjo a costa de la reducción de los dos ejes perpendiculares [B.3]. Dentro del sistema en rotación, el término centrífugo actuó para estabilizar y conservar el eje de la dirección del movimiento orbital ($b$), por lo que la mayor parte del esfuerzo de acortamiento y colapso mecánico recayó sobre el eje polar ($c$) [B.3]. Los polos fueron comprimidos y acercados entre sí, y la idea alternativa sostiene que los cráteres situados sobre ellos fueron arrastrados y se apiñaron con ellos [B.3].

La predicción morfológica (qué predice):

La geometría del cuerpo fluido libre de 2,600 ${\text{kg/m}}^3$ de densidad bajo estos esfuerzos de marea se calcula con precisión con una razón de ejes triaxial de $1:0.637:0.601$ (para los ejes $a:b:c$) [B.3]. Esta estructura produce un cuerpo elipsoidal muy largo en dirección a la Tierra y muy corto en los polos (un huevo apoyado sobre el eje de marea) [B.3]. Cuando se aplica esta deformación de estiramiento de marea a una corteza que transporta cráteres con distribución uniforme, el índice de densidad en cada celda de área varía en proporción inversa a los factores de contracción (exactamente lo inverso de los valores de C.3) [B.3]:

La lectura transversal de la matriz demuestra que un eje polar acortado mecánicamente bajo esfuerzos de estiramiento de marea extiende y expande la superficie situada sobre él, exactamente como en el mecanismo de achatamiento rotacional [B.3]. El polo es físicamente el sector más estirado de la forma, por lo que los cráteres en él se alejan unos de otros y se diluyen ($0.83$) [B.3]. En cambio, los extremos de la línea de la Tierra (el punto subterrestre) son las regiones en las que la superficie converge y se apiña hasta el máximo del sistema ($1.38$) [B.3].

Lo que dicen los recuentos empíricos:

Las observaciones reales de la superficie muestran un panorama completamente opuesto: los polos son las regiones más densas del sistema, mientras que el punto subterrestre es el más disperso y vacío de los seis casquetes principales (valor bruto de 94 y valor enmascarado de solo 148) [A.9, B.3]. El mecanismo de estiramiento de marea predice exactamente las inversiones opuestas en estos dos sectores [B.3].

Veredicto:

La hipótesis del estiramiento de Roche fracasa en la prueba del signo como mecanismo generador de apiñamiento [B.3].

No obstante, el efecto físico en sí es real y su magnitud se conserva firmemente en las ecuaciones del modelo: define la forma triaxial de partida de la sección C.2, y el redondeo geodinámico desde ella hacia una esfera perfecta es la respuesta científica de la investigación [B.3]. El efecto Roche es responsable de fijar el máximo de deformación y estiramiento del cuerpo acoplado; el redondeo posterior (la liberación del estiramiento y el retorno a la esfera) es el que impulsa y acerca físicamente los cráteres [B.3]. La lectura de la columna derecha de la tabla (el vector de redondeo) reproduce con plena concordancia la gradación de los recuentos: polos máximos, caras intermedias y punto subterrestre en la base del mínimo [B.3].

B.4 Gravedad más intensa en el polo (Polar Gravity Excess Hypothesis)

Definición de la interpretación:

En un cuerpo planetario de forma achatada, el vértice del polo se sitúa geométricamente más cerca del centro de masa del cuerpo que el cinturón ecuatorial [B.4]. Además, en el polo geográfico (situado directamente sobre el eje de rotación) la velocidad lineal se anula, por lo que no existe el componente de empuje centrífugo hacia fuera [B.4]. La combinación de estos dos factores produce en los polos un valor de gravedad efectiva (la diferencia entre la atracción propia y la contribución centrífuga) mayor que en el ecuador [B.4]. La interpretación alternativa propone que este campo gravitatorio reforzado en los polos actuó como una bomba mecánica que atrajo y desplazó material superficial elástico y suelto a lo largo de la corteza hacia las regiones polares, apiñando así el registro a su alrededor [B.4].

La predicción morfológica (qué predice):

La ley física dicta que el vector de la gravedad actúa hacia abajo (radialmente hacia el centro) y no lateralmente (tangencialmente a lo largo de la superficie) [B.4]. Para que el material cortical experimente un movimiento lateral activo a lo largo de la superficie hacia el polo, el polo debe funcionar físicamente como el punto más bajo en términos del potencial efectivo del cuerpo (el fondo de un valle de potencial, hacia el que el material se desliza y fluye como el agua pendiente abajo) [B.4].

En el sistema dinámico del cuerpo acoplado que gira rápidamente con la Tierra, la situación física es completamente la opuesta [B.4]. En el sistema en rotación, el potencial efectivo en el cinturón ecuatorial es inferior al del polo en un factor estricto de $\frac{1}{2}{\Omega }^2{R}^2$ (donde $\Omega$ es la velocidad de rotación y $R$ es el radio), valor que representa la energía cinética de movimiento excedente que posee el cinturón ecuatorial y de la que carece por completo el eje polar [B.4].

En las fases tardías del redondeo (densidad 2,600 ${\text{kg/m}}^3$, periodo de 5.77 horas, radio medio $R\approx 1,890\text{ km}$), con un valor de gravedad superficial de $g\approx 1.37{\text{ m/s}}^2$ (C.7), la ecuación de la altura teórica se define según [B.4]:

Esta ecuación demuestra que, sobre un cuerpo achatado que aún no ha completado su colapso hacia su forma de equilibrio hidrostático, las regiones polares se alzan como colinas topográficas gigantes de unos 120 km de altura “por encima” del cinturón ecuatorial en sentido energético [B.4]. Son los puntos de mayor potencial (la cima de una colina), por lo que todo material suelto o corteza elástica experimentará un vector de empuje opuesto y se deslizará alejándose de ellos hacia el valle de potencial ecuatorial [B.4]. Este es el mecanismo que produce achatamiento (B.2), y no apiñamiento polar [B.4].

Lo que dicen los recuentos empíricos:

El exceso se sitúa en los polos [B.4]. Mientras el cuerpo giraba, la combinación de la gravedad y la aceleración centrífuga actuó para alejar y diluir material de los polos, y no para concentrarlo en ellos [B.4].

Veredicto:

La hipótesis de la gravedad polar reforzada fracasa por completo como mecanismo de desplazamiento lateral mientras el cuerpo se encuentra en la fase acoplada y en rotación [B.4].

El efecto ocupa su lugar físico correcto solo en el momento de la transición dinámica: la autogravedad de la Luna es la que la mantuvo como un cuerpo íntegro e impidió su desintegración más allá del límite de Roche rígido ($\approx 13,000\text{ km}$ para una densidad de 1,300) [B.4]. En el momento en que el centro de masa de la Luna cruzó la línea de corrotación y la compresión mecánica sobre la base se anuló, el cuerpo quedó liberado de la restricción de la rotación síncrona excedente [B.4]. En esta situación, el achatamiento que la rotación sostenía hasta entonces se convierte en una carga anómala que el cuerpo no puede soportar [B.4]. Solo en ese momento el polo del cuerpo achatado colapsa y pasa de ser una cima energética a ser el fondo de un valle de potencial profundo [B.4]. Ahora la autogravedad actúa como el motor principal que impulsa el flujo y la contracción de la superficie hacia dentro, hacia los polos: este es el mecanismo de redondeo de vuelta (C.4) [B.4]. La diferencia geométrica de 120 km representa la escala exacta de la energía que la rotación mantenía arriba, y su liberación mecánica en el redondeo es la que apiña los cráteres [B.4].

B.5 Sumidero de drenaje en el polo (Polar Drain/Vortex Hypothesis)

Definición de la interpretación:

Esta hipótesis se basa en la analogía mecánica intuitiva de la “succión hacia dentro”: la existencia de un mecanismo de drenaje, de hundimiento cortical o de un vórtice dinámico interno bajo cada polo geográfico [B.5]. Según este modelo, el sistema magmático o geodinámico actúa como el agua en el fondo de un fregadero que fluye y se reúne hacia el desagüe, atrayendo, concentrando y recogiendo activamente la superficie y los cráteres situados sobre ella hacia un único punto extremo [B.5].

La predicción morfológica (qué predice):

La dinámica de un sumidero de drenaje dicta una concentración dirigida exclusivamente hacia un punto singular (The Vertex Point) [B.5]. Por tanto, el mapa de densidad debe mostrar una tendencia creciente pronunciada e ilimitada (una singularidad estadística) a medida que nos acercamos físicamente al propio punto del polo, y la topografía geofísica debe mostrar una depresión o un hoyo profundo en el punto central [B.5]. El modelo requiere la acción de un vector de fuerza mecánica lateral dirigido directamente hacia las coordenadas del polo [B.5].

Lo que dicen los recuentos empíricos:

Los datos medidos en la superficie descartan la estructura de singularidad puntual [B.5]:

1. Detención y estabilización del núcleo norte virgen: dentro del espacio limpio del polo norte, la densidad de cráteres ($\ge 20\text{ km}$) asciende efectivamente de forma radial: de 189 en el casquete de $60^{\circ }$, a 227 ($30^{\circ }$) y a 270 ($20^{\circ }$) [A.6, B.5]. Sin embargo, al penetrar en el núcleo extremo, la curva de densidad se detiene por completo y se estabiliza: un valor de 300 en el casquete de $15^{\circ }$ frente a un valor de 298 en el casquete de $10^{\circ }$ [A.6, B.5]. Los datos revelan una concentración global hacia una región polar amplia (Highlands Plateau) y no una convergencia puntual hacia un sumidero de drenaje [B.5].

2. Sensibilidad de la muestra sur residual: en el sur, tras aplicar la máscara, los casquetes interiores muestran valores de 305 (a 15°) y 348 (a 10°) [A.7, B.5]. Sin embargo, el análisis del sistema muestra que, en un radio de $20^{\circ }$ desde el vértice del polo sur, la máscara se compone únicamente de la unión coincidente con la cuenca SPA, y el casquete de $10^{\circ }$ se apoya en una muestra reducida de solo 36 cráteres [B.5]. En una muestra tan pequeña, una fluctuación aleatoria de un puñado de cráteres desplaza el valor numérico en proporciones amplias, por lo que el aumento en el sur refleja la geometría local de la cuenca de impacto y los márgenes de error del recuento, y no evidencia de un sumidero de drenaje [B.5]. El peso científico de la prueba se apoya en la secuencia norte limpia [B.5].

3. Ausencia de evidencias topográficas: los datos topográficos planetarios reales no muestran ningún par de depresiones gravitatorias o estructurales coincidentes con los polos geográficos [B.5].

Veredicto:

La hipótesis del sumidero de drenaje supera la prueba del signo (positivo), pero fracasa por completo en la prueba de la forma [B.5].

Dentro de un cuerpo planetario en rotación, todo vector de fuerza física debe conservar la simetría axial alrededor del eje de rotación, y tal simetría produce una estructura geométrica de anillo (Ring) y no de punto único; no existe ninguna fuerza física capaz de converger en el propio punto del polo [B.5]. El exceso de densidad real se extiende a lo largo de un anillo completo —el anillo situado a una distancia de $90^{\circ }$ del punto subterrestre, que pasa de forma continua por ambos polos y ambas caras (A.9, C.3)— y los polos encabezan la gradación únicamente porque el eje de rotación determina el centro de su convergencia geométrica [B.5]. La interpretación del sumidero de drenaje tiene el signo correcto (exceso), pero un patrón espacial que produce una distorsión: exige un máximo puntual, mientras que las observaciones muestran un nivel constante sobre una celda de área amplia y una banda de apiñamiento global que rodea todo el cuerpo [B.5].

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

Definición de la interpretación:

Esta hipótesis propone explicar el exceso de densidad mediante procesos isotrópicos (uniformes) globales, como una expansión de volumen uniforme, una contracción térmica uniforme de todo el cuerpo durante sus fases de enfriamiento, o la acción de una autogravedad intensa que atrae por igual todos los materiales de la corteza hacia el centro de masa; es decir, un cambio mecánico de la dimensión y el tamaño del cuerpo ($R$) sin cambio en su forma geométrica (estructura) [B.6]. Cada uno de estos acontecimientos constituye un proceso físico real y documentado en la historia temprana de un cuerpo planetario joven y caliente, y cada uno de ellos desplaza y modifica la posición de toda la superficie [B.6].

La predicción morfológica (qué predice):

Un cambio uniforme en la dimensión del cuerpo modifica cada elemento y cada unidad de área cortical exactamente en el mismo factor y la misma proporción [B.6]. La separación geográfica entre los cráteres aumenta o disminuye de forma homogénea en cada coordenada de la superficie del cuerpo, por lo que el índice de densidad varía en un factor constante y uniforme tanto en los polos como en el cinturón ecuatorial [B.6]:

Para una superficie craterizada uniformemente, un cambio global de tamaño multiplica el numerador (el número de cráteres) y el denominador (el área superficial) por el mismo factor matemático exacto, por lo que el valor de la razón estadística entre ellos no puede cambiar [B.6]. El mecanismo funciona como el inflado uniforme de un globo: la operación hace que todos los puntos impresos sobre él estén más dispersos, pero mantiene la igualdad y no es capaz de producir una región más densa o más dispersa que otra [B.6].

Lo que dicen los recuentos empíricos:

Las observaciones de la superficie fijan una diferencia morfológica marcada: una razón de $1.32$ en el recuento bruto, $1.20\pm 0.05$ en el recuento limpio y una tendencia ascendente constante en función del diámetro del cráter hasta alcanzar un valor de $1.53$ en las clases grandes [B.1, B.6].

Veredicto:

La hipótesis del cambio uniforme fracasa de forma absoluta tanto en la prueba del signo como en la prueba de la magnitud, ya que no es capaz de producir ningún exceso de densidad ni variación espacial alguna [B.6].

No obstante, es la interpretación metodológica que trazó el camino analítico de la investigación [B.6]. Una vez descartado físicamente un proceso de cambio uniforme de tamaño, la física dictó una conclusión firme: un cambio de tamaño no puede producir un exceso polar; solo un cambio en la forma del cuerpo (Anisotropic Distortion) es capaz de producir tal anomalía [B.6]. Este cambio de forma debe ser no uniforme en el espacio de manera definida: debe contraer mecánicamente la superficie en las regiones polares en la máxima proporción, y debe conservar una simetría axial completa alrededor del eje de rotación [B.6].

Dado que las hipótesis B.2 y B.3 representan cambios de forma reales pero con el signo mecánico erróneo (dilución de los polos), la solución científica del enigma exige su inversión geodinámica exacta [B.6]. Esta fase redujo la pregunta científica: la búsqueda se orientó exclusivamente a encontrar el cambio de forma específico que contrae la superficie polar más que cualquier otra región, y en la historia del modelo existe un único proceso geodinámico que cumple esta condición [B.6].

B.7 La respuesta, y por qué es la que permanece (The Flattening Release / Spherical Rebound)

Un único cambio de forma geodinámico en el sistema da una respuesta precisa y completa a todas las condiciones límite establecidas por las observaciones y en las que fracasaron las interpretaciones anteriores: la transición mecánica de un cuerpo triaxial muy achatado que se libera de sus restricciones y se redondea de vuelta a la forma de una esfera perfecta con una disminución espontánea de volumen [B.7].

Cierre de las condiciones y ajuste al sistema físico:

1. El signo mecánico (positivo): el mecanismo constituye la inversión geométrica y mecánica exacta de las hipótesis B.2 y B.3 [B.7]. La liberación del estiramiento y la contracción posterior actúan con el signo correcto y apiñan la corteza [B.7].

2. El motor físico: lo accionan las fuerzas de autogravedad de la sección B.4 en su lugar físico exacto: en el momento en que el centro de masa de la Luna cruza la línea de corrotación y la fuerza de compresión neta del cuerpo sobre la base de la pera se anula [B.4, B.7]. El cuerpo se libera de la restricción centrífuga excedente mientras sigue acoplado y girando en sincronización con la Tierra [B.7].

3. La geometría espacial: produce exactamente la estructura de anillo simétrico alrededor del eje de rotación (B.5), en la que las regiones polares encabezan el máximo del apiñamiento y la contracción mecánica porque el eje de rotación pasa por ellas y todos los meridianos retroceden hacia ellas (los polos son los únicos puntos que la rotación nunca desplazó) [B.5, B.7].

4. Prueba de la magnitud y coherencia cronológica: las dimensiones físicas se derivan directamente de la forma que los propios cráteres registran alrededor de los seis ejes (A.9) y de las pruebas de ajuste a la forma de equilibrio hidrostático del esferoide de Maclaurin (C.1) [B.7]. El modelo de redondeo no predice un valor estático, sino una gradación cronológica dinámica que concuerda de forma absoluta con los resultados empíricos del recuento limpio [B.7]:

• Cráteres gigantes antiguos ($60-100\text{ km}$): se formaron en el máximo de achatamiento sobre el cuerpo más deformado ($1.8:1.4:1$, densidad $1,700{\text{ kg/m}}^3$), experimentaron la secuencia de redondeo completa y por ello muestran el máximo apiñamiento: ~1.5 predicho frente a $1.53\pm 0.19$ empírico [A.8.5, B.7].
• Cráteres intermedios/tardíos ($20-30\text{ km}$): se formaron en fases tardías del colapso hacia dentro sobre un cuerpo que ya había empezado a redondearse ($\approx 1.20$, densidad $2,600{\text{ kg/m}}^3$), experimentaron solo el final de la contracción de la superficie y por ello muestran un efecto residual moderado: ~1.2 predicho frente a $1.09\pm 0.07$ empírico [A.8.5, B.7].
• Toda la población ($\ge 20\text{ km}$): el modelo de los seis ejes reconstruye una forma combinada de $1.37:1.15:1$, que predice de forma independiente una razón de 1.21, en concordancia exacta con el resultado del recuento limpio residual: $1.20\pm 0.05$ [A.9, B.7].

5. Prueba de la máscara: el efecto es robusto y sobrevive por completo a la limpieza de los maria, las cuencas y el sistema SPA, ya que está inscrito en la mecánica de la propia corteza sólida [B.7].

6. Momento y ventana de formación definidos: el proceso se sitúa cronológicamente en la fase de acoplamiento temprana (antes de la separación mecánica final, en la que la Luna se liberó en un estado redondo y rígido con una razón de ejes de $\le 1.04$) [B.7]. El interior del cuerpo era lo bastante blando y viscoso para permitir el flujo y el redondeo, mientras que la capa exterior era lo bastante rígida para transportar los cráteres como elementos íntegros y aislados (la contracción cerró los espacios elásticos entre ellos y generó la presión de C.4) [B.7].

El significado del aumento con el diámetro del cráter:

La tendencia ascendente de la razón de densidades a medida que aumenta el diámetro del cráter constituye la prueba empírica más clara y directa de la existencia de un proceso de redondeo que se desarrolla a lo largo del eje temporal [B.7]. Un cráter cortical solo puede conservar y registrar la porción de contracción superficial que se produjo después de su formación [B.7]. Los cráteres más grandes (los más antiguos) experimentaron toda la transformación, mientras que los cráteres pequeños (los tardíos) captaron solo el acorde final geométrico [B.7]. Ninguna hipótesis alternativa (de B.1 a B.6) es capaz de producir una dependencia estructural entre el exceso de densidad y el diámetro del cráter; el modelo de redondeo la exige físicamente [B.7].

La matriz de síntesis de la prueba mecánica

A continuación se presenta un resumen cuantitativo y sobrio de las seis interpretaciones rechazadas frente a la hipótesis principal del redondeo, que define la contribución residual de cada mecanismo descartado a la formulación científica final [B.7]:

Lean la matriz de arriba abajo: cada interpretación rechazada dejó tras de sí un componente físico firme, y la fila inferior los reúne en una única estructura coherente [B.7]. El estiramiento de B.2 y B.3 proporciona los indicadores de deformación que deben anularse, la autogravedad de B.4 proporciona el vector motor, la simetría axial de B.5 define la geometría, B.6 establece la ley de un cambio de forma necesario, y B.1 aísla el efecto de las superficies jóvenes [B.7].

Disponemos ahora de una causa física exclusiva que supera con éxito las cuatro pruebas empíricas, y de seis hipótesis rigurosamente descartadas [B.7].

Parte C — El mecanismo geodinámico y sus números

La Parte B de la investigación descartó las seis explicaciones alternativas [B.7]. El descarte analítico dejó una única causa mecánica posible: la transición geométrica de un cuerpo triaxial muy achatado que se libera de sus restricciones de esfuerzo y se redondea de vuelta a la forma de una esfera perfecta, con disminución espontánea de volumen y aumento de densidad [B.7, C]. Esta parte convierte esa respuesta en un mecanismo físico cuantitativo: define las dimensiones de la forma de la Luna joven, los vectores de esfuerzo que la modificaron, el efecto espacial sobre la corteza y la energía que impulsó el proceso [C].

La secuencia de fuerzas deformantes y la evolución morfológica:

La geometría de la Luna joven quedó determinada por una ecuación de equilibrio dinámico entre tres vectores de fuerza principales [C]:

1. La atracción de la autogravedad: tendía a comprimir todas las partes del cuerpo uniformemente hacia el centro, y por sí sola habría producido una esfera perfecta [C].

2. La aceleración centrífuga de la rotación rápida: en la fase de acoplamiento síncrono con la Tierra (periodo de unas 5.5 horas), la rotación empujó el cinturón ecuatorial hacia fuera y rebajó las regiones polares [C].

3. Los esfuerzos de estiramiento de marea (efecto Roche): estiraron el cuerpo de forma extrema a lo largo del eje radial entre la Tierra y la Luna [B.3, C]. La gravedad de la Tierra tiraba de la cara cercana hacia dentro, mientras que la fuerza centrífuga de la corrotación tiraba de la masa hacia fuera, con la base del cuerpo sujeta físicamente al extremo de la estructura de la pera [C]. Este estiramiento acortó los dos ejes perpendiculares y obligó al eje polar a absorber el máximo acortamiento mecánico [C].

La cronología de la secuencia de forma se despliega a lo largo de las fases mecánicas de separación [C]:

• Fase de agregación inicial: una masa siálica ligera de 1,300 ${\text{kg/m}}^3$ de densidad se consolidó en el extremo de la estructura de la pera de la Tierra y comenzó a desplazarse hacia fuera sobre una base de 6,000 km de anchura, con un achatamiento inicial bajo la restricción de la Tierra [C].
• Fase de estabilización del cuerpo y achatamiento máximo: a medida que la masa se alejaba, la base de la pera se redujo a una anchura de 3,000 km [C]. Cuando la pera alcanzó la corrotación, las dos caras del cuerpo fueron arrastradas en direcciones opuestas y el cuerpo alcanzó el máximo de su deformación geométrica cerca de la base de 3,000 km [C]. Esta forma máxima es la que analizan los capítulos C.1–C.2 [C].
• Fase de redondeo acoplado y disminución de volumen: en el momento en que el centro de masa de la Luna cruzó la línea de corrotación, la atracción neta que la presionaba sobre la base de la pera se anuló [C, C.6]. La corrotación permitió al cuerpo separarse por debajo y comenzar a redondearse [C]. La liberación del estiramiento de marea y rotacional provocó una disminución espontánea y marcada del volumen; el cuerpo se volvió más pequeño, más compacto y más denso, mientras su capa exterior se iba rigidizando [C].
• Fase de liberación final (la separación): en el momento de la separación mecánica final, el cuerpo ya era completamente redondo y su capa era rígida [B.7, C]. Tras la liberación, la ralentización de la rotación acoplada solo eliminó el pequeño abultamiento elástico del cuerpo rígido, y no el achatamiento temprano, que había desaparecido antes [C, C.5]. El proceso de redondeo es lo que la propia gravedad del cuerpo produjo en el momento en que los esfuerzos que sostenían la deformación se relajaron [C].

C.1 La forma inicial en el máximo de achatamiento: rotación y densidad

La pregunta de la sección cuantifica matemáticamente cuán achatada y deformada estaba la Luna en el máximo de la fase de achatamiento, y cómo cambió su forma a lo largo de las fases de colapso hacia dentro [C.1]. Para calcular el componente rotacional se necesitan solo dos parámetros físicos: la velocidad angular de rotación ($\Omega$) y la densidad media del material ($\rho$) [C.1].

Para un fluido homogéneo autogravitante (un cuerpo de densidad uniforme mantenido unido solo por su gravedad, sin resistencia mecánica inicial), la forma exacta de equilibrio hidrostático se define como el esferoide de Maclaurin (Maclaurin Spheroid) [C.1]. El grado de desviación del esferoide respecto a una esfera perfecta se mide mediante su excentricidad, $e=\sqrt{1-c^2/a^2}$ (donde $a$ es el radio ecuatorial y $c$ es el radio polar; $e=0$ da una esfera perfecta y $e\to 1$ da un disco aplanado) [C.1]. El cuerpo se mantiene estable en reposo cuando la velocidad angular $\Omega$ (expresada como $2\pi /T$, donde $T$ es el periodo de rotación) cumple la ecuación estructural [C.1]:

• La mecánica de la ecuación: el miembro izquierdo compara la fuerza centrífuga hacia fuera (${\Omega }^2$) con la autogravedad hacia dentro de un cuerpo uniforme ($G\rho$, donde $G$ es la constante de gravitación universal) [C.1]. El miembro derecho es una función trascendente que depende únicamente de la forma geométrica [C.1]. El miembro derecho asciende desde cero (en el estado esférico) hasta un valor máximo absoluto de 0.4493, alcanzado con una excentricidad de $e=0.930$ (dato que da una razón de ejes marcada de $a/c=2.72$) [C.1]. Más allá de este límite mecánico, no existe físicamente ninguna forma de Maclaurin estable [C.1].
• Sensibilidad paramétrica: cuanto más rápida es la velocidad de rotación en relación con la densidad, más achatada es la forma de equilibrio [C.1]. Un cuerpo de baja densidad (ligero) es más sensible al achatamiento que un cuerpo denso a la misma velocidad de rotación, porque su campo gravitatorio débil ofrece menos resistencia al esfuerzo centrífugo [C.1].

Matriz de solución de la razón de ejes dinámica ($a/c$):

La solución de la ecuación de Maclaurin para los valores de los parámetros de las fases de separación da la siguiente matriz de razones entre el radio ecuatorial y el radio polar [C.1]:

• Lectura horizontal (el vector de densidad): la forma geométrica experimenta un redondeo continuo (el valor $a/c$ cae hacia 1.00) a medida que aumenta la densidad del material, ya que una velocidad de rotación constante no es capaz de sostener una deformación achatada frente a un campo de autogravedad que se intensifica y atrae hacia dentro [C.1].
• Lectura vertical (el vector de rotación): la forma se redondea ligeramente a medida que la velocidad de rotación se ralentiza de 5.5 a 5.77 horas [C.1].
• Calibración de los extremos del redondeo cortical (los valores destacados): la matriz define los dos extremos de la transformación de la parte de la pita del cuerpo: el cuerpo estirado e hinchado en el máximo de achatamiento ($1.41$ con una densidad de unos 1,700 ${\text{kg/m}}^3$) y el cuerpo denso en las fases tardías del redondeo ($1.20$ con una densidad de 2,600 ${\text{kg/m}}^3$) [C.1].
• El valor de referencia de la pera ($1.300{\text{ kg/m}}^3$): el valor de 1.69 refleja la forma que la masa ligera acumulada en el extremo de la pera habría adoptado como fluido libre bajo la sola rotación; este dato sirve únicamente como contexto físico, y no como valor máximo medido del cuerpo autónomo [C.1].

El límite inferior de estabilidad para la densidad del material:

Con un periodo de rotación de 5.5 horas y una densidad de 1,000 ${\text{kg/m}}^3$, la ecuación no da solución de equilibrio, porque el esfuerzo centrífugo supera el punto máximo admisible del miembro derecho [C.1]. En esta situación hipotética, el cuerpo se desintegraría y perdería masa por los bordes de su ecuador [C.1]. Este punto define un límite inferior estricto para la densidad del material necesaria para mantener la integridad del cuerpo a cualquier velocidad de rotación:

Para un periodo de rotación de 5.5 horas, el límite mínimo de densidad se fija matemáticamente en 1,070 ${\text{kg/m}}^3$; un cuerpo planetario más ligero que este valor no es capaz de sostenerse como fluido con tal velocidad de rotación [C.1].

Ubicación del máximo de deformación y la geometría triaxial

• El momento del achatamiento máximo: en las fases iniciales, la masa actuaba como parte integral del extremo de la pera (densidad 1,300, distancia $\approx 8,991\text{ km}$), y su forma estaba dictada directamente por la Tierra [C.1]. El achatamiento morfológico propio se desarrolló en el momento en que la masa se convirtió en un cuerpo autónomo capaz de estabilizarse a sí mismo, condición que se cumple al cruzar el límite de Roche rígido [C.1]. Más allá de este límite, la autogravedad y la resistencia del material superan las fuerzas de marea disgregadoras de la Tierra [C.1]. Para un material ligero ($1,300{\text{ kg/m}}^3$), el límite de Roche se sitúa más lejos, a una distancia de 12,994 km (el límite varía en relación funcional con ${\rho }^{-1/3}$; para la densidad de liberación de 2,600 ${\text{kg/m}}^3$ el límite es más cercano y se sitúa en 10,314 km) [C.1]. El cuerpo se desplazó hacia fuera más allá de 13,000 km mientras su base se reducía de 6,000 a 3,000 km, y alcanzó el máximo de deformación geométrica y de achatamiento al aproximarse a la base de 3,000 km, muy cerca del cruce de la línea de corrotación ($15,818\text{ km}$ con un periodo de 5.5 horas) [C.1].
• La estructura triaxial (un huevo que es también una pita): la forma en el máximo de deformación no era un esferoide simétrico biaxial [C.1]. Mientras la rotación achataba el cuerpo alrededor de su eje, los esfuerzos de marea de dos cuerpos de la Tierra lo estiraban activamente a lo largo del eje de la línea de la Tierra [C.1]. La combinación de los esfuerzos impuso una forma triaxial con tres ejes distintos: un eje largo orientado hacia la Tierra ($a$), un eje intermedio orientado a lo largo del vector de movimiento orbital ($b$) y un eje corto y comprimido que une los polos ($c$) [C.1].

El registro empírico del exceso de cráteres reproduce y confirma con precisión estas dimensiones (secciones A.9, C.3): la forma máxima original tenía una razón de aproximadamente $1.8:1.4:1$, dentro de un intervalo de dispersión morfométrica de entre $1.64:1.45:1$ y $1.95:1.33:1$ [C.1]. La población total ($\ge 20\text{ km}$) converge en una media combinada de $1.37:1.15:1$ [C.1]. Se trata de una deformación moderada y físicamente clara: el eje largo es de 1.6 a 1.95 veces mayor que el eje polar, y el eje intermedio es de 1.33 a 1.45 veces mayor que este [C.1].

Confirmación de los indicadores de densidad en el máximo ($\approx 1,700{\text{ kg/m}}^3$):

La determinación del valor de densidad del cuerpo en el máximo de achatamiento se basa en el cruce de dos comprobaciones dinámicas independientes, que dan una concordancia completa [C.1]:

1. Confirmación de la parte de la pita según la relación de Maclaurin: la parte de la pita de la forma ($b/c\approx 1.4$) viene determinada exclusivamente por los esfuerzos de rotación [C.1]. En primera aproximación, la diferencia entre el eje intermedio y el corto de un cuerpo estirado por la marea depende de la velocidad de rotación propia, por lo que puede leerse directamente de la ecuación de Maclaurin [C.1]. Con un periodo de rotación de 5.5 horas, la relación de Maclaurin da un valor exacto de $b/c=1.41$ para un cuerpo de 1,700 ${\text{kg/m}}^3$ de densidad (el valor destacado en la tabla) [C.1]. El intervalo de dispersión empírico ($1.33-1.45$) acota la densidad entre 1,615 y 1,960 ${\text{kg/m}}^3$, con el núcleo concentrado en 1,600–1,750 [C.1].

2. Límite geométrico del cierre de la superficie: la transición de un cuerpo estirado a una esfera perfecta implica una liberación del estiramiento y una disminución global de volumen [C.1]. Las ecuaciones geométricas (C.3) demuestran que, si la densidad del cuerpo en el máximo hubiera sido superior a 1,800 ${\text{kg/m}}^3$, el volumen de estiramiento disponible habría sido demasiado pequeño; en tal caso, los extremos del eje de marea se habrían visto obligados a expandirse y crecer durante el redondeo en lugar de contraerse y cerrarse [C.1].

Ambas pruebas fijan el valor del cuerpo en el máximo como un cuerpo estirado, hinchado y con una densidad de unos 1,700 ${\text{kg/m}}^3$, lo que da un radio de esfera equivalente (el radio de una esfera de igual volumen) de unos 2,180 km [C.1]. Este valor contrasta con el radio de 1,890 km en la liberación y el radio actual de 1,737 km [C.1]. La masa siálica original en el extremo de la pera ($1,300{\text{ kg/m}}^3$, radio equivalente de unos 2,380 km) constituye la base de partida primitiva, antes de la activación de la fase de achatamiento dinámico [C.1].

Las fases del redondeo mecánico

Cuando el centro de masa cruzó la línea de corrotación, el cuerpo inició el proceso de redondeo de vuelta a la esfera, que se desarrolló en dos fases mecánicas separadas [C.1]:

1. Liberación de la parte del huevo (eje de marea $a$): se liberó primero y con rapidez, porque sobre la línea de corrotación el vector de atracción neta que presionaba el cuerpo contra la base de la pera cayó a cero [C.1, C.6]. Esta parte del estiramiento a lo largo de la línea de la Tierra colapsó hacia dentro [C.1].

2. Liberación de la parte de la pita (achatamiento rotacional $b/c$): avanzó en paralelo al espesamiento y la compactación del material bajo un periodo de rotación casi constante (C.4) [C.1]. Esta fase se define a lo largo del aumento de los indicadores de densidad ($2,000\to 2,300\to 2,600{\text{ kg/m}}^3$) y de un ligero alargamiento del día terrestre (hasta alcanzar $5.77$ horas en el momento de la liberación) [C.1]. Con esta dinámica, la razón de ejes cayó y se fijó en un valor de alrededor de 1.20 en las fases tardías del redondeo [C.1].

Un cuerpo fluido que completa una rotación propia cada pocas horas debe presentar físicamente una deformación achatada y no puede existir como esfera perfecta en estado viscoso [C.1]. La transición final y rígida desde una razón de alrededor de 1.20 hasta un cuerpo esférico perfecto en el momento de la liberación se estabiliza y refuerza mediante la rigidización de la capa exterior [C.1]. Este es un requisito estructural del modelo: una corteza superior con resistencia mecánica y rigidez de $\mu \ge 10\text{ GPa}$ es capaz de absorber los esfuerzos de rotación residuales y estabilizar un cuerpo redondo con una razón de ejes de 1.04 como máximo [B.7, C.1]. En el momento de la liberación la Luna es redonda y su capa es rígida [C.1]. El achatamiento residual menor que quedó tras la liberación es solo un pequeño abultamiento elástico, que decae y desaparece por completo a medida que el cuerpo se aleja y su velocidad de rotación en la órbita se ralentiza (C.5) [C.1].

La ecuación de Maclaurin es una ley física fundamental, y las soluciones de la matriz no dependen de nada salvo de los valores de la velocidad angular ($\Omega$) y de la densidad del material ($\rho$) [C.1].

C.2 El componente de marea: la forma libre y el efecto Roche

Mientras que la sección C.1 calculó el efecto de los esfuerzos de rotación únicamente, esta sección incorpora el segundo vector de fuerza que moldeó morfológicamente el cuerpo: las fuerzas de marea de la Tierra (efecto Roche) y el balance de potenciales compuesto del sistema [C.2]. El análisis sirve como herramienta geométrica para ilustrar la contribución del estiramiento radial, y no constituye la base de la prueba estadística de densidades [C.2].

La dinámica de dos cuerpos en el sistema en rotación:

Un cuerpo planetario situado muy cerca de la Tierra y con acoplamiento rotacional síncrono (periodo orbital igual al periodo de rotación propia $\Omega$) se estira de forma extrema a lo largo del eje radial entre él y la Tierra [C.2]. En la fase acoplada, la masa está sujeta al extremo de la estructura de la pera, bajo un equilibrio de esfuerzos opuestos: la gravedad terrestre tira del flanco cercano hacia dentro y la fuerza centrífuga tira de él hacia fuera [C.2].

El modelo físico clásico reproduce esta estructura mediante el cálculo de la forma de un fluido homogéneo libre (sin restricción de contacto), que adopta la configuración de un elipsoide de Roche (Roche Ellipsoid) con tres ejes distintos ($a:b:c$): el eje de marea hacia la Tierra ($a$), el eje de la dirección del movimiento orbital ($b$) y el eje corto que une los polos ($c$) [C.2]. A diferencia de la forma biaxial de Maclaurin (que se ve como un círculo perfecto desde encima del polo), el elipsoide de Roche está deformado y tiene forma de huevo desde cualquier ángulo de observación [C.2].

La ecuación del balance de potenciales hidrostático:

La superficie de un cuerpo fluido en reposo debe constituir físicamente una superficie equipotencial (Equipotential Surface), para impedir el flujo de material a lo largo de ella [C.2]. Dentro de un elipsoide uniforme, el potencial total se compone de la suma matemática de tres componentes de energía [C.2]:

1. La autogravedad del cuerpo: la energía crece con el cuadrado de la distancia al centro a lo largo de los ejes ($x,y,z$, paralelos a los ejes $a,b,c$), y los coeficientes estáticos $A_1,A_2,A_3$ (que cumplen la ecuación de índices de Chandrasekhar: $A_1+A_2+A_3=2$) determinan la pendiente del aumento [C.2]. Al eje más corto ($c$) le corresponde el coeficiente máximo ($A_3$), ya que la masa está concentrada más cerca de él [C.2].

2. El potencial de marea de la Tierra: se expresa como $-\mu (x^2-\frac{1}{2}{y}^2-\frac{1}{2}{z}^2)$, donde el factor de estiramiento es $\mu =\frac{GM}{r^3}$ ($M$ es la masa de la Tierra y $r$ la distancia orbital) [C.2]. Este componente rebaja el potencial a lo largo del eje de marea ($x$) y lo eleva transversalmente ($y,z$), acción que ejerce un estiramiento radial y una compresión tangencial [C.2].

3. El potencial centrífugo de la órbita: se expresa como $-\frac{1}{2}{\Omega }^2(x^2+y^2)$ (según la tercera ley de Kepler en el sistema: ${\Omega }^2=\frac{G(M+m)}{r^3}$) [C.2]. Este componente rebaja el potencial hacia fuera, lejos del eje de rotación, y deja sin cambio el componente polar ($z$) [C.2].

Para que la superficie se equilibre en una forma cerrada, la atracción neta hacia dentro (la gravedad menos los factores de rotación y de marea) multiplicada por el cuadrado de la longitud del eje debe ser exactamente igual en los tres extremos [C.2]:

Esta ecuación combinada da dos sistemas de condiciones para las dos razones de ejes independientes: $b/a$ y $c/a$ [C.2].

Matriz de solución de las ecuaciones del elipsoide de Roche

A continuación se presentan los valores numéricos exactos obtenidos al resolver el sistema de ecuaciones, calculados a una distancia de 20,595 km (algo más allá del límite de Roche fluido de un cuerpo de 2,600 ${\text{kg/m}}^3$ de densidad), que sirve como herramienta de ilustración morfológica [C.2]:

• Desciframiento de las dimensiones de la forma estirada ($1:0.637:0.601$): se establece una asimetría extrema, en la que el eje de marea orientado hacia la Tierra es dos tercios más largo ($\approx 66\%$) que el eje polar corto [C.2].
• Normalización al volumen de una esfera equivalente ($R_s$): el radio medio de una esfera de igual volumen se define según $R_s=(abc{)}^{1/3}$ [C.2]. Las dos últimas columnas cuantifican el grado de desplazamiento de la superficie durante la transición geométrica (el redondeo) hacia esta esfera equivalente [C.2]:
  • El componente del polo: el valor $R_s/c=1.208$ dicta que, durante el redondeo, el vértice del polo debe ascender y desplazarse hacia fuera en alrededor de una quinta parte de su altura ($\approx 21\%$) [C.2].
  • El componente del eje de marea: el valor $R_s/a=0.726$ dicta que el punto subterrestre debe colapsar hacia dentro y hundirse en alrededor de una cuarta parte de su longitud ($\approx 27\%$) [C.2].
• Confirmación matemática independiente: una nueva resolución de las condiciones de contorno desde un punto de partida libre y analítico reproduce exactamente las formas (valores de $1:0.634:0.598$ para una densidad de 2,600 y $1:0.762:0.725$ para una densidad de 3,344), con una desviación de coherencia mínima inferior al medio por ciento ($\le 0.5\%$) y una precisión de potenciales de dos partes por mil [C.2].

Implicaciones geodinámicas y límites del límite de Roche rígido

El cálculo matemático fija dos leyes físicas fundamentales que limitan las fases de separación [C.2]:

1. Descarte de la liberación fluida ligera: con la densidad de origen de la pera ($1,300{\text{ kg/m}}^3$), el límite de Roche fluido se sitúa muy lejos, a una distancia de unos 25,200 km [C.2]. Esto significa que el cuerpo no puede físicamente alcanzar su fase de liberación dinámica manteniendo esta baja densidad, ya que se desintegraría por la marea; este dato confirma el supuesto básico del Artículo 2, que estableció que la densidad debe subir a 2,600 ${\text{kg/m}}^3$ para recortar los límites de Roche hacia dentro [C.2].

2. Colapso del límite de Roche fluido exacto: con una densidad de 2,600, la solución de las ecuaciones muestra que a 19,660 km la forma se desintegra, mientras que a 20,595 es estable; el límite de Roche fluido exacto se sitúa entre ambos: en 19,971 km según el coeficiente clásico, y en 20,100 km ($k=0.09$) con la solución exacta de los potenciales [C.2].

El límite de la corteza rígida en la secuencia dinámica:

En el modelo LEF actualizado, la Luna se libera y se separa mecánicamente mucho antes y mucho más cerca de la Tierra: a una distancia de unos 16,400 km en el ecuador ($16,399\text{ km}$ con un periodo de 5.77 horas) y de 17,287 km en las latitudes de $22.5^{\circ }$ [C.2]. A estas distancias cercanas, la física dicta que no existe ni puede existir una forma fluida libre estable con una densidad de 2,600 (la región se encuentra muy dentro de los dominios del límite de Roche fluido) [C.2].

El significado físico es firme: la Luna joven resistió y conservó su integridad estructural en esta región gracias a su capa exterior, que se rigidizó y se congeló pronto, y no en virtud de un equilibrio hidrostático fluido [C.2]. Por tanto, el modelo triaxial de este capítulo sirve únicamente como indicador y herramienta de ilustración que muestra cómo los esfuerzos de marea externos impusieron un estiramiento radial marcado a lo largo del eje $a$, además del achatamiento rotacional del eje $b/c$ (C.1) [C.2]. El exceso de densidad y la anomalía estadística medidos en la superficie son consecuencia exclusiva de esta liberación combinada de esfuerzos durante la secuencia de redondeo y disminución de volumen (C.1, C.5, D.2) [C.2].

Mientras que las secciones C.1 y C.2 definieron las dimensiones iniciales de la deformación, esta sección calcula la transformación de la superficie deformada durante su colapso mecánico hacia una forma esférica perfecta bajo la condición de disminución de volumen (compactación y endurecimiento del núcleo) [C.3]. El objetivo es determinar matemáticamente qué celdas de área cortical se cierran y contraen en la máxima proporción, cuáles en la mínima, y dónde se apiñan los cráteres existentes [C.3]. El resultado del cálculo es un índice de razón de densidades diferencial entre los polos y el ecuador, que es robusto y no depende de la tasa de contracción isotrópica general del cuerpo [C.3].

• Distinción cronológica en la serie: el modelo geodinámico define tres fases de deformación sucesivas:
  1. El redondeo hacia fuera: la fase inicial en la que la masa siálica se separó del extremo de la pera y se convirtió en una esfera fluida aislada (Artículo 3).
  2. La fase de achatamiento: la fase de alargamiento y achatamiento extremo de dos cuerpos cerca de la línea de corrotación.
  3. El redondeo hacia dentro (el redondeo de este artículo): la fase de colapso de la estructura estirada de vuelta a la esfera, que se produce tarde en la secuencia dinámica, después de que el centro de masa cruzara la línea de corrotación [C.3].

La analogía del globo aplastado:

Para hacer accesible la mecánica, imaginemos un globo redondo cubierto de puntos (cráteres) distribuidos uniformemente por unidad de área, que es presionado y aplastado con fuerza entre dos palmas planas (los esfuerzos de dos cuerpos) [B.5, C.3]. Debido a la presión, sus partes superior e inferior se vuelven anchas, achatadas y planas, mientras que el borde periférico (el ecuador) se estira bruscamente [B.5, C.3].

En el momento en que los esfuerzos externos se relajan (cruce de la línea de corrotación), el globo libera la energía elástica residual y recupera de golpe su forma esférica [C.3]. En este proceso geométrico, las superficies anchas y planas de los extremos se recogen rápidamente hacia dentro y convergen en un recinto reducido: los puntos impresos en ellas se acercan físicamente entre sí en la máxima proporción, mientras que el borde medio experimenta el mínimo cambio de área [B.5, C.3].

La mecánica geofísica de la fijación de los puntos extremos en los polos

¿Por qué son precisamente las regiones polares geográficas las que funcionan como los centros principales de convergencia y apiñamiento de la corteza [C.3]?

• Refutación de las hipótesis de espesor o gravedad: las observaciones de los satélites GRAIL descartan la existencia de una corteza engrosada o de anomalías gravitatorias laterales en las regiones polares respecto al cinturón de latitud en que se encuentran (B.4) [C.3].
• La restricción del eje de rotación (The Axial Geometric Anchor): los polos son los puntos extremos estructurales por los que pasa físicamente el eje de rotación del cuerpo [C.3]. Dado que los esfuerzos de achatamiento rotacional actúan de forma completamente simétrica alrededor de este eje, todos los meridianos globales se deformaron y estiraron de manera uniforme, mientras que los dos puntos extremos del eje permanecieron estáticos (no tienen un componente de momento “hacia fuera” propio) [C.3].
• El mecanismo de los radios de la rueda: en el momento en que el achatamiento rotacional se libera, el material cortical que antes había sido empujado regresa hacia dentro exactamente a lo largo de los mismos meridianos [C.3]. Dado que todos los meridianos de un cuerpo esférico convergen y se encuentran en los polos, estas regiones extremas funcionan como un cubo estructural que reúne en sí la superficie de todos los meridianos [C.3]. Esta propiedad geométrica explica por qué en el sistema existen exactamente dos focos de apiñamiento polares, situados en posición diametralmente opuesta el uno frente al otro [C.3]. El modelo fija los polos del redondeo primitivo como los polos geográficos observados hoy en la superficie [C.3].

Distribución de los esfuerzos de tracción diferenciales sobre la superficie

Durante la fase de deformación, la corteza en las regiones polares era la más plana, estirada y extendida del sistema [C.3]. En cambio, el estiramiento de marea de la Tierra (efecto Roche) estiró el cuerpo de forma extrema a lo largo del eje $a$ (la línea radial), con su parte cercana sujeta al extremo de la estructura de la pera [C.3]. En consecuencia, la superficie en los dos extremos de la línea de la Tierra (los puntos subterrestre y antiterrestre) se apiñó y estaba menos estirada que cualquier otro sector [C.3].

Por tanto, en el momento de la liberación de los esfuerzos y la relajación del cuerpo:

• Máxima contracción de la superficie: se produce en la superficie polar (retorno del estiramiento rotacional máximo) y genera el mayor exceso de densidad de cráteres [C.3].
• Mínima contracción de la superficie: se produce en los extremos del eje de marea (el punto subterrestre) y deja el registro de cráteres en su estado más disperso ($94$ cráteres por 10⁶ km²) [A.9, C.3].

La ecuación del balance geométrico: el factor de densidad combinado

El apiñamiento mecánico no se limita a los extremos polares; afecta a un anillo global completo situado a una distancia angular de $90^{\circ }$ del eje de la línea de la Tierra (el anillo polar-tangencial que pasa por ambos polos y ambas caras orbitales) [C.3]. La razón es que los dos extremos del eje de marea, al colapsar hacia dentro durante el redondeo, empujan y desplazan su porción de superficie lateralmente, hacia el anillo medio [C.3].

A continuación se presenta una calibración matemática de los factores bidireccionales de cambio de área, utilizando el elipsoide de Roche (C.2) como modelo ilustrativo [C.2, C.3]:

• Desciframiento de la dinámica matemática: para cada pequeña celda de área cortical se definen dos direcciones de esfuerzo perpendiculares a lo largo del terreno [C.3]. El redondeo acorta o alarga las dimensiones de la celda a lo largo de estos dos ejes, mientras que el número total de cráteres sobre ella permanece sin cambio [C.3]. Por tanto, la densidad final de cráteres se calcula como el producto de los inversos de estos dos factores de cambio de área (${\rho }_{\text{combined}}={\text{Factor}}_1\times {\text{Factor}}_2$) [C.3].
• Normalización del sistema: los valores se expresan en relación con el cambio de área global medio [C.3]. Un factor $\ge 1.00$ indica que la celda de área se cerró y contrajo en esa dirección más que la media (los cráteres se acercaron entre sí); un factor $<1.00$ indica una contracción inferior a la media (los cráteres se alejaron relativamente) [C.3]. En el polo, la ecuación establece: $1.38\times 0.88=1.21$ [C.3].
• La gradación del anillo: los polos y las caras experimentan un factor de contracción idéntico desde el eje de marea ($1.38$), pero los polos reciben material de las caras en una proporción superior al material que les transfieren [C.3]. En consecuencia, los polos se sitúan en la cabeza de la gradación ($1.21$) y las caras después de ellos ($1.14$) [C.3]. El punto subterrestre se expande en sus dos ejes respecto a la media, por lo que muestra el valor de densidad mínimo del sistema ($0.73$) [C.3].

Confirmación de la gradación a partir de los datos limpios de los seis ejes

A continuación se presentan los datos empíricos reconstruidos a partir del catálogo Robbins en los casquetes de $30^{\circ }$, tras aplicar el procedimiento de “recuento limpio” (eliminación completa de los maria y de la cuenca SPA en ambos lados de la ecuación), que confirman con precisión las predicciones de la ecuación del balance geométrico [A.9, C.3]:

Los resultados empíricos son firmes: todo punto situado a lo largo del anillo polar muestra un exceso de densidad significativo en comparación con el valle de densidad del punto subterrestre, que contiene menos de la mitad de la cantidad de cráteres de los polos [C.3]. La asimetría interna entre la cara de cola y la delantera se explica mediante la distribución de las aberturas de eyección locales y el flujo de momento del material retornado que quedó rezagado respecto a la velocidad de rotación y cayó detrás (A.9) [C.3]. Las regiones polares son el lugar geométrico en el que se produce el máximo apiñamiento y en el que el fenómeno se mide en su forma más virgen y limpia, y de ahí toman su nombre los capítulos del artículo [C.3].

C.3 El cálculo matemático de la contracción diferencial de la superficie y las predicciones morfométricas

Esta sección constituye el núcleo analítico estricto del artículo [C.3]. Traduce los cambios de forma geométricos del cuerpo en ecuaciones de transformación de la superficie y calcula con precisión las razones de apiñamiento diferenciales entre las distintas regiones de la corteza [C.3]. El cálculo se basa en una separación física absoluta entre dos componentes de colapso independientes:

1. El cambio de forma (Anisotropic Deformation): se calcula bajo la restricción estricta de volumen constante hacia una esfera equivalente de radio $R_s=(abc{)}^{1/3}$ [C.3].

2. El cambio de tamaño isotrópico (Isotropic Volume Contraction): la Luna experimentó una contracción térmica global y una liberación del estiramiento que redujeron su volumen en un factor de alrededor de 1.5 durante el redondeo [C.3]. Dado que esta contracción común multiplica y reduce cada celda de área de la superficie del cuerpo por el mismo factor constante, se cancela matemáticamente y no modifica ninguna razón de densidades [B.6, C.3]. Por tanto, el índice de apiñamiento de los polos frente al ecuador es completamente robusto frente a la tasa de contracción general del cuerpo [C.3].

La ecuación de la transformación elipsoidal homogénea:

La transición geométrica se representa como una deformación elástica homogénea que transforma la forma del elipsoide en una esfera, como si la superficie estuviera dibujada sobre una lámina de goma uniforme [C.3]. Cada coordenada a lo largo del eje largo, el intermedio y el corto se multiplica por los factores de estiramiento correspondientes: $R_s/a,R_s/b,$ y $R_s/c$ respectivamente [C.3]. Una pequeña celda de área cuyo vector normal exterior (perpendicular a la superficie y de longitud unidad) se define como $\hat{n}=(n_x,n_y,n_z)$ cambia su área en el factor geométrico combinado [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}

Los cráteres impresos de manera uniforme sobre la superficie original cambian sus densidades en proporción inversa a este factor de área ($J$), ya que la misma cantidad de cráteres comprimida en una celda de área más pequeña mostrará una densidad mayor [C.3].

El significado de la fórmula en el espacio de la corteza:

La ecuación establece que el área de una mancha orientada a lo largo de un eje espacial está gobernada por el producto de los estiramientos de los dos ejes perpendiculares a él [C.3]. Dado que el producto de los tres estiramientos es exactamente igual a 1 (conservación del volumen en la fase de forma), el cambio de área es exactamente igual al inverso del estiramiento a lo largo del eje normal de la mancha [C.3]:

• En el vértice del polo ($n_z=1$): el factor se reduce a su valor mínimo, $c/R_s$ [C.3]. La superficie experimenta la máxima contracción geométrica, por igual en todas las direcciones tangentes, y la densidad residual de cráteres se dispara en la proporción máxima de $R_s/c$ [C.3].
• En el punto subterrestre ($n_x=1$): el factor adopta su valor máximo en el sistema, $a/R_s$ [C.3]. La superficie experimenta la mínima contracción, y la densidad de cráteres allí es la más baja, $R_s/a$ [C.3].
• La transición radial: entre el polo y el ecuador, el factor de contracción se debilita de forma continua a medida que nos alejamos del vértice del eje [C.3]. Esta propiedad mecánica explica geométricamente por qué el perfil de densidad del polo norte virgen (A.6) asciende y se acentúa hacia el interior, directamente hacia el punto del polo geográfico [C.3].

Calibración del modelo biaxial (el esferoide achatado)

Para un cuerpo achatado simétricamente en el que el eje largo y el intermedio son iguales ($a=4b$), la condición de conservación del volumen dicta $R_s^3=a^{2}c$ [C.3].

• La superficie polar cambia en un factor de $c/R_s=(c/a{)}^{2/3}$, que representa un acortamiento biaxial uniforme de $(c/a{)}^{1/3}$ en todas las direcciones sobre el terreno [C.3].
• La superficie ecuatorial cambia en un factor de $a/R_s=(a/c{)}^{1/3}$, que representa un acortamiento este-oeste de $(c/a{)}^{1/3}$ frente a una contracción mínima en el eje norte-sur: es la celda de área que menos se contrae de toda la superficie del cuerpo [C.3].

La razón estadística entre estas dos regiones define la ecuación del apiñamiento absoluto del polo frente al ecuador, en la que el radio de la esfera ($R_s$) se cancela por completo [C.3]:

Se establece una ley firme: el apiñamiento teórico máximo en el punto del polo frente al ecuador es exactamente igual a la razón de ejes ecuatorial/polar ($a/c$) de la forma achatada original [C.3]. Un cuerpo que era 1.69 veces más ancho que alto en el máximo de deformación dejará los cráteres impresos en él con una densidad exactamente 1.69 veces mayor en el vértice del polo que en el ecuador [C.1, C.3]. Sin embargo, cuando se calculan las medias sobre los casquetes empíricos finales, el valor de la razón es menor, porque la contracción se atenúa con la distancia angular al polo [C.3].

Calibración del modelo triaxial y ajuste a los seis ejes

Para un cuerpo triaxial achatado y estirado ($a:b:c$), la misma ecuación diferencial establece una gradación de densidades triple: la densidad asciende en un factor $R_s/c$ en los polos, en un factor $R_s/b$ en las caras ecuatoriales y en un factor $R_s/a$ en la base de la línea de la Tierra [C.3]. Esta gradación (polos $\to$ caras $\to$ línea de la Tierra) reproduce y concuerda de forma absoluta con la gradación empírica de los seis ejes medida en la sección A.9 [C.3].