Con una cuenta desbloqueas el cuestionario, el foro de dudas y el tutor con IA de esta lección, además del registro de progreso y el certificado verificable al terminar «Mecánica Orbital: de Kepler al TLE».
Las cuatro lecciones anteriores han construido un edificio impecable sobre una hipótesis falsa: que la Tierra es una masa puntual y que no hay más fuerzas que su gravedad. En el mundo real cinco de los seis elementos orbitales no son constantes. Cambian lentamente, con ritmos que van del grado por día a la fracción de milímetro por revolución, y esos cambios lentos son responsables de la mitad del diseño orbital moderno.
El enfoque estándar es tratar las desviaciones como perturbaciones sobre la órbita kepleriana: se conserva la parametrización de la lección 3 y se escribe cómo derivan sus elementos. Ordenadas por magnitud relativa a la aceleración gravitatoria principal:
Perturbación
Magnitud relativa
Dónde domina
Efecto principal
Achatamiento J2
10−3
Todas las órbitas terrestres
Precesión de Ω y de ω
Arrastre atmosférico
10−5 a 10−9
LEO por debajo de 800 km
Decaimiento de a, reentrada
Tercer cuerpo (Luna y Sol)
10−6
MEO, GEO, HEO
Deriva de i, unos 0.85°/año en GEO
Armónicos superiores (J3, J22)
10−6
GEO y órbitas de precisión
Deriva en longitud de los GEO
Presión de radiación solar
10−7
Alta relación área/masa
Excitación de la excentricidad
Todas son pequeñas frente al término dominante. Ninguna es despreciable si el horizonte de predicción pasa de horas a semanas, y por eso ningún propagador serio —el SGP4 de la lección 7 incluido— es puramente kepleriano.
La Tierra no es una esfera: el término J2
La rotación achata el planeta: el radio ecuatorial supera al polar en 21.4 km. Ese exceso de masa acumulado en el ecuador deforma el campo gravitatorio, que deja de ser central. El potencial se desarrolla en armónicos esféricos, y el primer término no trivial domina abrumadoramente sobre todos los demás:
con ϕ la latitud geocéntrica. El coeficiente J2 es unas mil veces mayor que cualquier otro armónico, hasta el punto de que en la mayoría de los análisis se habla simplemente de "la perturbación J2" como si fuera la única.
Su efecto es un par de fuerzas: el abultamiento ecuatorial tira del satélite hacia el plano del ecuador cuando este se encuentra fuera de él. Igual que en una peonza, ese par no endereza la órbita: la hace precesar, girando el plano orbital alrededor del eje polar sin cambiar su inclinación. Adiós al plano congelado de la lección 1.
Regresión nodal
Promediando el efecto sobre una revolución, el nodo ascendente deriva a razón de:
Ω˙=−23a2(1−e2)2nJ2R⊕2cosi
y en la forma práctica que conviene tener a mano, con a en kilómetros y el resultado en grados por día:
Ω˙≈−2.06474×1014a−7/2(1−e2)−2cosi
Tres lecturas inmediatas. El signo negativo con cosi>0 significa que las órbitas prógradas regresan hacia el oeste. Las órbitas polares (i=90°) no precesan, porque cos90°=0. Y las retrógradas (i>90°) precesan hacia el este, porque el coseno cambia de signo: esa es la puerta que abre las heliosíncronas.
Casi cinco grados diarios: el plano de la estación completa una vuelta entera en unos 73 días. Esto no es un tecnicismo, es planificación de misión pura. Explica por qué las ventanas de lanzamiento hacia la ISS son de minutos —hay que insertarse en un plano que se mueve—, por qué los periodos de iluminación de sus paneles varían a lo largo del año, y por qué las oportunidades de acoplamiento se repiten con una cadencia de dos meses largos.
Órbitas heliosíncronas: convertir el defecto en función
La Tierra recorre su órbita alrededor del Sol a 360°/365.2422=0.9856 grados por día. Si conseguimos que el plano orbital de un satélite precese exactamente esa misma cantidad y en el mismo sentido, el ángulo entre el plano y la dirección al Sol permanecerá constante todo el año. Es la definición de órbita heliosíncrona, y el resultado observable es que el satélite cruza el ecuador siempre a la misma hora solar local.
Como necesitamos Ω˙=+0.9856 grados/día, y el signo de la fórmula es negativo para las prógradas, hace falta cosi<0: la órbita tiene que ser retrógrada. Despejando para 550 km de altitud (a=6928 km):
Ese es el origen del número que aparece en casi todos los TLE de satélites de observación. Y como Ω˙ depende también de a, la inclinación heliosíncrona cambia con la altitud:
Altitud
a (km)
i heliosíncrona
Ejemplo real
400 km
6 778
97.0°
Misiones cortas de observación
550 km
6 928
97.6°
Capa SSO de Starlink, cubesats
600 km
6 978
97.8°
Sentinel-3 (aprox.)
800 km
7 178
98.6°
NOAA, Metop
1 000 km
7 378
99.5°
Órbitas de vigilancia
La elección del ángulo respecto al Sol define además dos familias operativas. La mediodía-medianoche (LTAN 12:00) maximiza la iluminación para instrumentos ópticos; la amanecer-atardecer (LTAN 6:00 o 18:00) mantiene el satélite sobre la línea del terminador, con los paneles iluminados de forma casi permanente y un ambiente térmico muy estable: es la preferida por los radares de apertura sintética, que no necesitan luz solar pero sí mucha potencia.
Avance del perigeo y la inclinación crítica
El mismo J2 hace girar la elipse dentro de su plano:
ω˙≈+1.03237×1014a−7/2(1−e2)−2(5cos2i−1)
El factor (5cos2i−1) es la clave y cambia de signo. Se anula cuando cos2i=1/5, es decir:
icrıˊtica=63.4349°y116.5651°
Ahí está el origen exacto del número que la lección 4 dejó pendiente. Una Molniya con i=63.4° mantiene su apogeo clavado sobre el hemisferio norte durante toda la vida de la misión, sin gastar una gota de propelente. Con i=45°, el mismo satélite vería su argumento del perigeo derivar unos +3.4 grados diarios, y en poco más de un mes el apogeo habría emigrado 100 grados: el sistema de comunicaciones dejaría de funcionar para su público objetivo. Diseñar en la inclinación crítica es el ejemplo canónico de explotar una perturbación en lugar de combatirla.
Arrastre atmosférico: el enemigo de la órbita baja
Por debajo de unos 800 km la atmósfera no ha terminado: hay una exosfera tenuísima cuyo rozamiento, aunque minúsculo, actúa sin descanso y siempre en el mismo sentido. La aceleración de frenado es:
aD=−21ρmCDAv2v^
donde CDA/m es el coeficiente balístico inverso del objeto. La densidad ρ cae de forma casi exponencial con la altura, con una altura de escala de 40 a 80 km:
Altitud
ρ mínimo solar (kg/m³)
ρ máximo solar (kg/m³)
Vida orbital típica
200 km
2.5×10−10
4×10−10
Días
300 km
2×10−11
6×10−11
Semanas
400 km
3×10−12
2×10−11
Meses a 1 año
550 km
1×10−13
2×10−12
3 a 8 años
800 km
2×10−15
5×10−14
Décadas a siglos
1 000 km
3×10−16
8×10−15
Más de mil años
Las dos columnas de densidad no son un exceso de celo. La actividad solar calienta y expande la termosfera, y su indicador operativo es el flujo de radio a 10.7 cm, el índice F10.7, medido en unidades de flujo solar. Oscila entre unos 65 sfu en el mínimo y más de 250 sfu en el máximo del ciclo de once años, y la densidad a 400 km puede variar un factor diez entre ambos extremos. Predecir el decaimiento de un satélite exige, por tanto, predecir el Sol, que es un problema abierto: es la mayor fuente de incertidumbre en toda la predicción orbital operativa.
Ejemplo resuelto 2: cuánto pierde la ISS cada día
El decaimiento del semieje mayor por revolución en órbita casi circular vale aproximadamente:
Δarev≈−2πmCDAρa2
Para la ISS, con CDA/m≈0.011 m²/kg, ρ≈2×10−12 kg/m³ y a=6.778×106 m:
Δarev≈−2π×0.011×2×10−12×4.594×1013=−6.4m
Con 15.5 revoluciones diarias son unos 100 metros al día, es decir del orden de 2 a 3 km al mes. Por eso la estación necesita reboosts periódicos; por eso su altitud dibuja el característico diente de sierra en las gráficas históricas; y por eso la cadencia de esos reboosts se acelera notablemente en los máximos solares.
La paradoja del arrastre
Hay un detalle que desconcierta a todo el mundo la primera vez. El rozamiento retira energía, luego a disminuye, luego —por la tercera ley de Kepler— el periodo se acorta y el satélite acelera. Frenar a un satélite lo hace ir más rápido. No hay contradicción: la energía total disminuye, pero al caer a una órbita más baja convierte energía potencial en cinética con creces. Es el mismo mecanismo del fasing de la lección 6, aquí impuesto por la atmósfera en lugar de por un motor.
En órbitas excéntricas el arrastre actúa casi solo cerca del perigeo, donde la densidad es alta y la velocidad máxima. El resultado es que baja el apogeo dejando el perigeo casi intacto: la órbita se circulariza progresivamente. Solo cuando la elipse ya es casi una circunferencia baja empieza el descenso final, que es rápido y de duración muy difícil de predecir. Esa es la firma característica de toda reentrada no controlada.
Casos reales
En febrero de 2022 SpaceX lanzó 49 satélites Starlink a una órbita de despliegue de 210 km, deliberadamente baja para poder desorbitar rápido cualquier unidad defectuosa. Una tormenta geomagnética aumentó la densidad atmosférica en torno a un 50 % respecto a lo previsto, el arrastre superó la capacidad de maniobra de los satélites en modo de seguridad y 38 de los 49 reentraron en días. Es la demostración más cara de que el clima espacial es un parámetro de diseño y no un detalle.
A 550 km, la altitud operativa de la constelación, un Starlink sin propulsión activa reentra de forma natural en unos cinco años: por diseño, y no por accidente, cumple así el requisito de mitigación de basura espacial sin depender de un sistema de desorbitado. En el extremo opuesto, la constelación Iridium a 780 km tardaría más de un siglo, lo que convierte cada satélite muerto en un problema heredado por varias generaciones. El tema se retoma a fondo en la lección 9.
B*: el arrastre empaquetado en el TLE
Todo lo anterior tiene que caber en los 69 caracteres de un TLE. La forma en que se hace es el parámetro B∗ (B-star), un coeficiente balístico modificado que el propagador SGP4 usa como única entrada de arrastre:
B∗=2ρ0mCDA
expresado en inversos de radio terrestre. Un B∗ grande identifica un objeto ligero y de gran superficie —un panel suelto, un cubesat con velas, un fragmento de basura— que decaerá rápido; un B∗ pequeño, un objeto denso y compacto. En la práctica B∗no es una propiedad física medida, sino un parámetro ajustado por mínimos cuadrados contra las observaciones de rastreo, y absorbe todo lo que el modelo no explica. De ahí que un satélite maniobrando pueda aparecer con B∗ negativo, lo que físicamente sería un arrastre que empuja: no es un error del catálogo, es el ajuste haciendo su trabajo.
Cuando en la lección 7 abras un TLE real y leas la secuencia 30777-3 en las columnas 54 a 61, estarás leyendo exactamente esto: el resumen de toda esta sección en ocho caracteres. Y cuando en Sistemas APT un pase predicho se desvíe treinta segundos respecto al horario calculado con un TLE de hace una semana, la causa más probable será que el B∗ ajustado ya no describe la atmósfera que el satélite está atravesando hoy.
Las que faltan: SRP y tercer cuerpo
Dos perturbaciones más completan el cuadro, y aunque su tratamiento cuantitativo excede esta lección, conviene saber cuándo dominan.
La presión de radiación solar transfiere momento de los fotones al satélite, con una presión de 4.56×10−6 N/m² a la distancia terrestre. La aceleración resultante en un GEO típico es de unos 10−7 m/s², despreciable en un día pero capaz de excitar la excentricidad de forma apreciable en meses. Escala con la relación área/masa, de modo que castiga sobre todo a los satélites de grandes paneles y a los objetos ligeros; y como desaparece durante los eclipses, produce un efecto periódico que complica el modelado.
La atracción de Luna y Sol es la perturbación dominante por encima de MEO. En GEO inclina el plano orbital entre 0.75 y 0.95 grados por año, y combatirla es el 90 % del presupuesto de station-keeping que cuantificaremos en la lección 6. Curiosamente, si se deja actuar libremente durante décadas, la inclinación de un GEO abandonado oscila entre 0° y unos 15° con un periodo de unos 53 años: la naturaleza no destruye la órbita, solo la saca de servicio.
Laboratorio
El simulador de esta lección separa las dos perturbaciones dominantes en dos escenas conectadas. La primera muestra el plano orbital girando por efecto de J2: eliges altitud, inclinación y excentricidad, y la escena dibuja la precesión de la RAAN acelerada miles de veces, con la dirección al Sol superpuesta para que puedas ver si el plano la sigue, la adelanta o se queda atrás. Un panel muestra Ω˙ y ω˙ calculados con las fórmulas prácticas de esta lección.
La segunda escena integra el decaimiento por arrastre con un modelo exponencial de densidad y un control de actividad solar F10.7. Verás simultáneamente la evolución del semieje mayor, del apogeo y del perigeo, y podrás comprobar de primera mano la circularización previa a la reentrada: el apogeo desciende mientras el perigeo se mantiene, hasta que ambos se juntan y la caída se dispara.
Presta atención al acoplamiento entre ambas: al bajar la altitud, la precesión nodal se acelera (depende de a−7/2) al mismo tiempo que el decaimiento se dispara. Un satélite que pierde altura pierde también su sincronía heliosíncrona, y esa es la razón operativa por la que las misiones de observación gastan propelente en mantener la altitud aunque su instrumento funcionaría igual de bien un poco más abajo.
Tres retos concretos:
Encuentra la heliosíncrona a tres altitudes. Para 400, 700 y 1 000 km, ajusta la inclinación hasta que Ω˙ marque +0.9856 grados/día y compara con la tabla de la lección. Explica por qué la inclinación necesaria aumenta con la altitud.
Verifica la inclinación crítica. Con e=0.7 y periodo de 12 h, mide ω˙ para i=45°, 63.4° y 80°. Estima cuántos días tardaría el apogeo de la Molniya de i=45° en desplazarse 90 grados de latitud útil.
Compara ciclos solares. Sitúa un cubesat con CDA/m=0.02 m²/kg a 450 km y simula su vida orbital con F10.7 = 70 y con F10.7 = 200. Cuantifica el factor entre ambas vidas y comenta qué implica para el cumplimiento de la regla de desorbitado en 25 años.
1 de 8
Comprueba la idea que hace posibles las heliosíncronas:
¿Por qué una órbita heliosíncrona debe ser retrógrada, con i>90°?
Con esto cerramos el bloque de física orbital: sabemos qué forma tiene una órbita, cómo se describe con seis números, dónde se sitúa sobre el mapa y cómo la deforma el mundo real. En la lección 6 el satélite deja de ser espectador y enciende motores.
Laboratorio: Perturbaciones · Precesión J2 y decaimiento por arrastre
8. Del TLE al pase: la geometría de la visibilidad