Análisis orbital y espectral de meteoros

1. Elementos orbitales y lluvias de meteoros

1.1. Principales elementos orbitales de un meteoro

Para saber de dónde procede una roca espacial y cómo se mueve, es necesario definir su “documento de identidad” matemático: los elementos orbitales keplerianos, que describen la forma, el tamaño y la orientación de su trayectoria alrededor del Sol.

Elemento	Símbolo	Qué es
Semieje mayor	a	Define el tamaño de la órbita y su periodo (P = a3/2).
Excentricidad	e	Describe la forma de la órbita, desde circular (e = 0) hasta elíptica (0 < e < 1).
Distancia al perihelio	q	Distancia mínima al Sol, q = a(1 − e).
Inclinación	i	Ángulo de la órbita respecto al plano de la Tierra (eclíptica).
Longitud del nodo ascendente	Ω	Punto en el que la órbita cruza el plano de la Tierra.
Argumento del perihelio	ω	Orientación de la elipse dentro de su plano.

Tabla 1. Elementos orbitales keplerianos. Adaptada de Ceplecha (1987) y Espartero y Madiedo (2016).

Para ilustrar cómo estos elementos definen una órbita real, en el Anexo A se incluye la Figura A.1, con la trayectoria kepleriana del bólido SPMN051213.

1.2. Criterios de procedencia y fundamento matemático

La procedencia de un meteoroide se establece mediante dos criterios fundamentales, cuyo desarrollo matemático se recoge en los Anexos A.1.1 y A.1.2:

1. Parámetro de Tisserand (T_J). Es clave para distinguir entre cometas y asteroides. Si T_J > 3, el origen es asteroidal (como el bólido profundo de 2013, con T_J = 3,24). Si T_J < 2, el origen es cometario (cometas de tipo Júpiter).
2. Criterio de Southworth y Hawkins (D_SH). Mide la similitud entre la órbita de un meteoro y la órbita media de una lluvia conocida. Una órbita se asigna a una lluvia si D_SH < 0,15.

Aplicación al bólido SPMN051213. Para aplicar estos conceptos se utilizó la metodología de reducción de datos de Ceplecha (1987). Los parámetros orbitales y dinámicos del bólido profundo SPMN051213 (Espartero et al., 2014) dan un parámetro de Tisserand T_J = 3,24, lo que confirma su origen asteroidal. Los detalles de esta clasificación se encuentran en el Anexo A.2, y la geometría de sus elementos puede verse en la Figura A.1, que representa su órbita heliocéntrica y la posición de sus principales parámetros keplerianos.

1.3. Definición y origen de las lluvias de meteoros

Una lluvia de meteoros (o stream) se produce cuando la Tierra cruza una corriente de escombros en el espacio. Estos escombros son partículas que comparten órbitas similares, por lo que, desde nuestra perspectiva, parecen proceder de un mismo punto del cielo: el radiante.

La causa fundamental es el paso de la Tierra a través de esta corriente de polvo y pequeños fragmentos, que orbitan el Sol en una trayectoria similar a la de un cuerpo progenitor. La Unión Astronómica Internacional (IAU) cataloga actualmente cientos de corrientes, aunque las lluvias principales observadas cada año son entre 10 y 12.

Tipos de cuerpos progenitores (Ceplecha, 1985):

• Origen cometario. Es el más común. Los cometas subliman al acercarse al Sol y liberan polvo y gases que dejan un rastro de partículas a lo largo de su órbita. Por ejemplo, las ζ-Perseidas diurnas son un rastro ligado al complejo del cometa Encke (Madiedo et al., 2016).
• Origen asteroidal. El cuerpo progenitor es una roca sólida o un asteroide que se ha desintegrado. Estas partículas son más densas y menos frágiles. Por ejemplo, las ω-Escórpidas del Norte (NSC) se confirmaron como de origen asteroidal gracias a su T_J elevado y a su comportamiento dinámico en la atmósfera (Espartero y Madiedo, 2016).

Lluvia de meteoros	Objeto progenitor	Máximo aproximado	Tasa horaria zenital (THZ) máxima
Cuadrántidas	Asteroide 2003 EH1	Principios de enero	~120 meteoros/h
Líridas	Cometa Thatcher	Finales de abril	~18 meteoros/h
Eta Acuáridas	Cometa Halley	Principios de mayo	~60 meteoros/h
Perseidas	Cometa Swift-Tuttle	Mediados de agosto	~100 meteoros/h
Oriónidas	Cometa Halley	Finales de octubre	~23 meteoros/h
Leónidas	Cometa Tempel-Tuttle	Mediados de noviembre	~10–15 meteoros/h
Gemínidas	Asteroide 3200 Phaethon	Mediados de diciembre	~150 meteoros/h

Tabla 2. Características de las principales lluvias anuales de meteoros. Elaboración propia a partir de la International Meteor Organization (IMO), el IAU Meteor Data Center y la información de NASA/AMS sobre lluvias de meteoros.

2. Ablación, bólidos y meteoritos

2.1. El proceso de ablación y la luminosidad

La ablación es la pérdida de masa por fusión y vaporización de las capas externas del meteoroide al entrar en la atmósfera. La entrada a velocidades hipersónicas (entre 11 y 72 km/s) genera una compresión extrema del aire y una onda de choque que convierte la energía cinética en calor. Este calentamiento funde y vaporiza la superficie del meteoroide y da lugar a fenómenos físicos clave: la compresión del aire, el calentamiento extremo y la emisión de luz.

La masa total ablatada (m_a) puede estimarse a partir de la intensidad luminosa I(t) del bólido, ya que la luz emitida es proporcional a la energía que el meteoroide disipa durante la ablación. La expresión formal, basada en la integración de la energía radiada (Madiedo et al., 2016), y la definición de sus parámetros se presentan en el Anexo B.1.1.

La ablación no es solo un proceso microscópico de vaporización: afecta directamente a la evolución dinámica del meteoro durante la entrada. Para ilustrar la relación entre la pérdida de masa y el descenso atmosférico se construyó un modelo simplificado (Figura 2).

2.2. Presión dinámica, fragmentación y clasificación de Ceplecha

No todos los meteoros son iguales. La diferencia entre ver un bólido o encontrar un meteorito depende de una batalla entre la resistencia de la roca y la presión del aire. El bólido se define por un criterio visual de brillo: cualquier meteoro que alcanza una magnitud absoluta de −4 o más brillante (Espartero et al., 2014). El meteorito se produce solo si el objeto logra desacelerar hasta unos 3 km/s antes de consumirse por completo y entra en la fase de vuelo oscuro (dark flight) hasta impactar contra el suelo (Ceplecha, 1987).

El criterio de presión dinámica. Que un meteoro llegue a ser meteorito depende del equilibrio entre la presión dinámica (p) que ejerce la atmósfera y la resistencia mecánica (tensile strength) de la roca. La presión máxima que sufre el meteoroide en el punto de mayor esfuerzo es:

p = ρ_aire · v²

donde ρ_aire es la densidad del aire y v la velocidad (Madiedo et al., 2016).

Clasificación física de Ceplecha (1985). Ceplecha estableció el estándar mundial para clasificar estos objetos según su capacidad de sobrevivir a la atmósfera. A partir de datos fotográficos, dividió los bólidos en cuatro grupos físicos:

Grupo	Composición típica	Densidad (ρ)	¿Cae meteorito?
Tipo I	Condrita ordinaria (roca sólida)	~3,7 g/cm³	Sí. Son rocas duras y suelen llegar al suelo.
Tipo II	Condrita carbonácea (roca porosa)	~2,1 g/cm³	A veces. Son más frágiles, pero las más grandes sobreviven.
Tipo IIIA	Cometario (bola de polvo y hielo)	~0,75 g/cm³	No. Se desintegran a gran altura.
Tipo IIIB	Cometario blando (Dracónidas)	< 0,3 g/cm³	No. Son extremadamente frágiles, como terrones de azúcar cósmicos.

Tabla 3. Clasificación de los bólidos según su resistencia estructural y su capacidad de supervivencia (Ceplecha, 1985).

La caída de un meteorito está determinada principalmente por la composición (tipos I o II) y la velocidad inicial: las velocidades bajas, por debajo de 25 km/s, favorecen la supervivencia. Los bólidos de origen cometario (tipos III), aunque sean brillantes, se consumen por completo en la atmósfera (Ceplecha, 1985).

3. El espectro de un meteoro

3.1. Método de determinación espectral y tipo de espectro

El espectro de un meteoro se determina mediante espectroscopía de alta velocidad, generalmente sin rendija (slitless), acoplando una red de difracción a cámaras de alta sensibilidad que registran el fenómeno en fracciones de segundo. Al atravesar la red, la luz del meteoro se dispersa en sus distintas longitudes de onda y produce una imagen espectral que puede analizarse fotograma a fotograma.

Debido a la extrema rapidez del fenómeno (decenas de km/s), la instrumentación requiere cámaras de alta sensibilidad y velocidad de lectura, una calibración precisa en longitud de onda y espectrógrafos o redes de difracción montados directamente sobre el sistema óptico.

El resultado es un espectro de emisión: líneas brillantes sobre un fondo oscuro. Estas líneas proceden del plasma generado por la ablación, en el que los átomos y moléculas del meteoroide (y del aire) se excitan y emiten luz en longitudes de onda discretas. Constituyen la huella dactilar química del meteoroide.

3.2. Información espectral y relación físico-química

Las líneas de emisión permiten identificar los elementos presentes en el meteoroide. Entre los multipletes más relevantes destacan Na I-1 (sodio), Mg I-2 (magnesio) y Fe I (hierro).

La intensidad relativa de estos multipletes es clave para distinguir entre condritas, meteoroides asteroidales y fragmentos cometarios. Por ejemplo, una marcada depleción de sodio (Na/Mg bajo) caracteriza a los meteoroides de la corriente de las ω-Escórpidas del Norte (NSC), mientras que las condritas ordinarias muestran relaciones Na/Mg y Fe/Mg cercanas a la unidad.

Para visualizar estas diferencias composicionales de forma cuantitativa se generó un espectro sintético comparativo (Figura 3). Los espectros simulados comparan un meteoroide condrítico con relaciones Na/Mg ≈ 1 y Fe/Mg ≈ 1 frente a un meteoroide de la corriente NSC con Na/Mg = 0,47 y Fe/Mg = 0,80, según los valores de Espartero y Madiedo (2016), y muestran con claridad la menor contribución del multiplete Na I-1 en este último.

3.3. Procesos físicos observables en el espectro

El espectro no solo permite conocer la composición, sino también inferir procesos físicos que ocurren durante la ablación:

• La ionización y la recombinación del plasma.
• La fragmentación del meteoroide (líneas ensanchadas o con múltiples componentes).
• La interacción con la atmósfera (líneas de oxígeno y nitrógeno ionizados).
• La temperatura del plasma, a partir de la intensidad relativa de las líneas.

Estos procesos ayudan a caracterizar el tipo físico del meteoroide y su comportamiento estructural durante la entrada atmosférica.

Relación con la dinámica y la masa ablatada. La espectroscopía también es fundamental para estimar la masa perdida por ablación. La intensidad luminosa total (I) registrada está relacionada físicamente con la tasa de pérdida de masa (dm/dt), la velocidad (v) y la eficiencia luminosa (τ):

I = −½ · τ · v² · (dm / dt)

donde I es la intensidad luminosa total, dm/dt la tasa de pérdida de masa por ablación, v la velocidad del meteoro y τ la eficiencia luminosa (Ceplecha, 1985). Esta relación permite, a partir del espectro y de la fotometría del meteoro, estimar cuánta masa se pierde por vaporización y evaluar su comportamiento físico en la atmósfera.

Conclusión

El análisis orbital, dinámico y espectral permite caracterizar de manera integral el origen y la naturaleza física de los meteoroides. Los elementos keplerianos y los criterios de Tisserand y de Southworth–Hawkins identifican su procedencia; la ablación y la presión dinámica describen su evolución atmosférica y su posible supervivencia; y el espectro de emisión revela su composición química y los procesos físicos que tienen lugar durante la entrada.

En conjunto, estos métodos ofrecen una reconstrucción completa del comportamiento de un meteoro, desde su órbita heliocéntrica hasta su interacción final con la atmósfera terrestre.

Anexo A. Fundamentos físico-matemáticos y datos del bólido SPMN051213

A.1.1. Cálculo del parámetro de Tisserand (T_J)

El parámetro de Tisserand respecto a Júpiter (a_J ≈ 5,204 UA) es:

T_J = a_J / a + 2 · cos i · √[(a / a_J) · (1 − e²)]

Sustitución y resultado para SPMN051213 (datos de Espartero et al., 2014). Con a = 2,36 UA, e = 0,619 e i = 8,99°:

T_J = 5,204 / 2,36 + 2 · cos(8,99°) · √[(2,36 / 5,204) · (1 − 0,619²)] ≈ 3,25

La pequeña diferencia respecto al valor publicado de 3,24 se debe al redondeo de los datos orbitales.

Los valores de Tisserand de otros eventos analizados confirman el régimen asteroidal (T_J > 3):

• ZPE 2012 (bólido rozador de las ζ-Perseidas): T_J = 3,31.
• NSC (promedio de las ω-Escórpidas del Norte): T_J = 3,52.
• Bólido de 2013 (SPMN051213): T_J = 3,24.

Todos estos valores caen dentro del régimen asteroidal (T_J > 3) y se diferencian del origen cometario.

A.1.2. Criterio de similitud orbital de Southworth y Hawkins (D_SH)

Este criterio permite asociar la órbita de un meteoroide individual (q, e, i, ω, Ω) a una corriente meteórica conocida, comparándola con la órbita media de la lluvia (q′, e′, i′, ω′, Ω′). La disimilaridad se define como:

D_SH² = (e − e′)² + (q − q′)² + [2 sen(I/2)]² + [(e + e′)/2]² · [2 sen(Π/2)]²

donde I es el ángulo entre los planos de ambas órbitas y Π la diferencia entre las longitudes de sus perihelios, medidas desde la intersección de dichos planos.

Criterio de asociación. El meteoroide se asocia a la lluvia si la disimilaridad es baja (D_SH < 0,15). Algunas aplicaciones:

• Bólido rozador ZPE 2012: D_SH = 0,13, asociado a la lluvia de las ζ-Perseidas.
• ω-Escórpidas (núcleo de la corriente): D_SH ≤ 0,07, valor que confirma una asociación muy estrecha.

A.2. Parámetros del bólido SPMN051213

La tabla recoge no solo parámetros orbitales, sino también parámetros dinámicos o cinemáticos esenciales para el análisis de la entrada en la atmósfera.

Parámetro	Símbolo	Valor	Unidad
Parámetro de Tisserand	TJ	3,24	—
Semieje mayor	a	2,36	UA
Excentricidad	e	0,619	—
Distancia al perihelio	q	0,898	UA
Velocidad geocéntrica	vg	29,5 ± 0,5	km/s

Tabla A.2. Datos adaptados de Espartero et al. (2014), procesados con Python (Matplotlib y Astropy).

Anexo B. Fundamentos físico-matemáticos, dinámica y clasificación estructural

B.1.1. Cálculo de la masa ablatada (m_a)

La masa que pierde el meteoro puede calcularse a partir de su brillo. Según Madiedo et al. (2016), la masa ablatada se obtiene integrando la energía radiada:

m_a = 2 ∫ I(t) / (τ · v²) dt

donde:

• I(t) es la intensidad luminosa (el brillo) en cada instante, que se obtiene de la magnitud absoluta M mediante I = I₀ · 10^−0,4·M, siendo I₀ la intensidad de un meteoro de magnitud cero.
• v es la velocidad del meteoro.
• τ es la eficiencia luminosa, que depende de la velocidad (Ceplecha y McCrosky, 1976).

Ejemplo del artículo sobre ZPE 2012. La aplicación de esta integración a eventos reales permite cuantificar la pérdida de masa. En el caso del bólido rozador de las ζ-Perseidas de 2012, Madiedo et al. (2016) determinaron que el meteoroide perdió una masa ablatada de m_a ≈ 260 ± 30 g.

Referencias

• American Meteor Society (s. f.). Meteor Showers. AMS.
• Ceplecha, Z. (1985). Fireball information on meteoroids and meteorites. Bulletin of the Astronomical Institutes of Czechoslovakia, 36(4), 237–242.
• Ceplecha, Z. (1987). Geometric, dynamic, orbital and photometric data on meteoroids from photographic fireball networks. Bulletin of the Astronomical Institutes of Czechoslovakia, 38(4), 222–234.
• Ceplecha, Z., & McCrosky, R. E. (1976). Fireball end heights: A diagnostic for the structure of meteoric material. Journal of Geophysical Research, 81(35), 6257–6275.
• Espartero, F., Madiedo Gil, J. M., & Castro Tirado, A. J. (2014). A deep-penetrating fireball produced by an asteroidal meteoroid. 45th Lunar and Planetary Science Conference.
• Espartero, F. A., & Madiedo, J. M. (2016). The Northern ω-Scorpiid meteoroid stream: orbits and emission spectra. Earth, Moon, and Planets, 118(2), 81–89.
• International Meteor Organization (2024). Meteor Shower Calendar 2025 (IMO INFO(3-24)). IMO.
• International Astronomical Union, Meteor Data Center (s. f.). IAU Meteor Shower Database. IAU.
• Jet Propulsion Laboratory (s. f.). Diagrama de la órbita del cometa 1P/Halley ilustrando los elementos cometarios.
• Madiedo, J. M., Espartero, F., Castro-Tirado, A. J., Pastor, S., & de los Reyes, J. A. (2016). An Earth-grazing fireball from the Daytime ζ-Perseid shower observed over Spain on 2012 June 10. Monthly Notices of the Royal Astronomical Society, 460(1), 917–922.