¿Cómo se calculan rigurosamente los mapas de probabilidad de impacto?

ingrese la descripción de la imagen aquí

Este es un ejemplo de un mapa de probabilidad de impacto. Se publicó en esta pregunta y es de Aerospace Corporation. Pero Mi pregunta es sobre el proceso, no sobre este mapa en particular, aunque es un buen ejemplo de trabajo para plantear esta pregunta.

Con una nave espacial dada, su TLE más reciente o los últimos, y una ventana de tiempo arbitraria definida por t 1 y t 2 , uno puede usar SGP4 para propagar la órbita entre esos dos tiempos para hacer un gráfico de seguimiento terrestre, que sería una serie de líneas delgadas onduladas.

Ese es más o menos el proceso descrito en este comentario , pero no es suficiente. Probablemente t 1 y t 2 debe elegirse para poner entre paréntesis de forma segura cualquier conjetura razonable sobre el momento del impacto. entonces en t 1 es casi seguro que todavía está en órbita, y en t 2 es casi seguro que está en el suelo (o vaporizado, o algo de ambos). Las trayectorias terrestres en esos puntos deberían contribuir mínimamente al histograma de probabilidad de área 2D final.

Un corolario de eso podría abordar lo que sucede con el histograma debajo de la trayectoria en el suelo en el momento más probable del impacto . Imagen a continuación: es de la publicación del blog ESA RocketScience Tiangong-1 Reentry Updates . Haga clic para tamaño completo.

Probabilidad de caída de la órbita de Tiangong-1

Mirando más de cerca el histograma de probabilidad 2D (ver más abajo), se destacan algunas cosas.

  1. ¡ Los bordes interiores entre la región verde ancha, interior (probabilidad inferior) y las bandas amarillas estrechas (probabilidad superior) cerca de los extremos de la latitud son irregulares! Dado que la atmósfera en sí misma no tiene esta estructura, esto debe ser algún tipo de efecto residual debido al número finito de órbitas consideradas. Ya no son líneas finas, pero se sigue viendo su efecto.
  2. Fuera de las estrechas bandas de alta probabilidad cerca de los extremos hay un segundo par de bandas aún más estrechas entre la alta probabilidad y la probabilidad cero. Esto no es lo que uno obtendría al observar solo la probabilidad del tiempo de permanencia de la pista terrestre.

Por lo tanto, parece haber mucho más que "lo que necesita es el tiempo que pasa sobre una unidad de área determinada de superficie en el transcurso de muchas órbitas".

Pregunta: ¿Cómo se calculan rigurosamente los mapas de probabilidad de impacto?

ingrese la descripción de la imagen aquí

arriba: Recortado de una versión más grande del mapa, que se encuentra aquí .

@Litho Hay más espacio para ampliar su comentario como una nueva respuesta aquí.
Mucho más para leer en el reingreso de Tiangong-1 de Spaceflight101 y el reingreso de Tiangong-1 de Aerospace Corporation .
Tres palabras, demasiado cortas para ser una respuesta: Simulación Monte Carlo. Dos palabras más, también demasiado cortas para ser una respuesta: Análisis Lincov.
@Litho También he agregado una recompensa.
@RussellBorogove respondí a tu comentario antes de que terminaras de editarlo. Si puede ampliar eso, en lo que respecta a este problema en particular, sería genial. Esa delgada franja verde en el borde exterior es de particular interés, ya que sería algo como esto del artículo de la ESA vinculado. i.stack.imgur.com/dqgnu.png
La franja verde representa un extremo exterior de las "barras de error" de incertidumbre en la posición/velocidad actual de la estación más un extremo exterior del posible rango de desviación aerodinámica acumulada a medida que cae en la atmósfera.
@RussellBorogove Los señalo porque resaltan algunos de los aspectos destacados que entrarían en el cálculo riguroso y serían importantes como parte de una respuesta. No estoy preguntando hey, ¿qué son esas cosas?

Respuestas (2)

No estoy totalmente de acuerdo con la explicación dada hasta ahora. Aquí están mis conjeturas adicionales:

2) Los bordes exteriores: este no es un artefacto de agrupamiento como lo indica @PearsonArtPhoto. Los escombros no solo caen a lo largo del camino del satélite, hay una región con un cierto ancho alrededor del camino donde tenemos que esperar que caigan escombros. Es razonable suponer que la distribución dentro de esta región tiene la forma de una curva de campana: cuanto más lejos del camino, más improbable es encontrar escombros. Eso significa que no hay un corte repentino sino una distribución suave entre alta probabilidad y cero.

1) Estoy de acuerdo en que esto proviene de la cantidad de órbitas que se tuvieron que considerar dada la longitud de la ventana de reingreso. Pero: ¿Por qué las ondas no son simétricas alrededor de la órbita dibujada del satélite? A medida que la órbita decae, el período orbital se acorta. Por otro lado, como no conocemos la tasa exacta de descomposición, no podemos determinar con precisión la trayectoria orbital relativa a la superficie de la Tierra.

Sospecho que las curvas dibujadas son de una órbita constante, pero el cálculo de la probabilidad de impacto incluyó el cambio (desconocido y aleatorio) en el período orbital. Si la órbita se vuelve más corta, el satélite pasará sobre áreas más al este con cada revolución, pero nunca más al oeste que el camino "nominal". Y esto es lo que vemos: En el lado derecho del camino las probabilidades son más altas que en el lado izquierdo.

En resumen, todo esto parece lo que debe esperar después de ejecutar una serie de predicciones orbitales que involucran fluctuaciones aleatorias en la tasa de desintegración y suponiendo una cierta distribución de desechos alrededor del punto de desintegración calculado.

Es lamentable que la otra respuesta comience con "Supongo que...", pero no estoy seguro de que eso signifique que las conjeturas ahora son respuestas. Si todos empezáramos a publicar conjeturas como respuestas...
@uhoh Si tiene fuentes definitivas, adelante. Siempre que no tenga acceso a las fuentes originales, lo mejor que puede obtener es una conjetura informada. Decir que es una conjetura y no ponerlo como un hecho es solo una buena práctica.
Eliminé "rigurosamente" del título y "Así que me gustaría saber cómo los expertos hacen esto como una cuestión de procedimiento estándar". del cuerpo Si realiza una edición simbólica de su respuesta, puedo revertir mi voto negativo a positivo.

Mi conjetura es que la forma en que hicieron esto fue tomar los caminos orbitales que podrían ocurrir y averiguar por dónde volvería a entrar si ese fuera el caso.

Lo más probable es que el artefacto en los bordes sea del tamaño de un binning. Solo la parte central (superior e inferior) es una posibilidad, pero el contenedor todavía está allí, por lo que la cantidad real es pequeña. Una forma más sencilla de pensar en esto es imaginar que tiene contenedores para un número cada grado verticalmente de 0 a 90. Digamos que el más alto fue 50,1. El pequeño contenedor de 0,1 lo mostraría como una pequeña posibilidad, pero depende de en qué parte del contenedor se encuentre, el 50,0 bien podría ser la zona más probable de impacto. También podría tener en cuenta los vientos y otras cosas similares que podrían desdibujar un poco el borde.

Las olas parecen tener un período de unos 12 grados, o unos 60 picos. Eso corresponde a aproximadamente 4 días de órbitas, más o menos un poco, lo que sospecho corresponde a la región de incertidumbre.

La conclusión es crear un modelo de la ruta que incluya las incertidumbres de los parámetros clave. Varíe ese modelo con las diferentes incertidumbres variadas. Tenga algún tipo de imagen que muestre los posibles puntos de reingreso, y cada vez que algo golpee ese punto, agréguele un punto. Continúe iterando hasta que tenga una buena idea de ello. Es posible que hayan incluido algunas cosas, como que los escombros caigan sobre un arco en lugar de un solo punto también.

La conclusión es que se realizó algún tipo de simulación con condiciones aleatorias modeladas, y el gráfico es lo que salió.

No puedo entender lo que esto está tratando de decir; " Lo más probable es que el artefacto en los bordes sea del tamaño de un contenedor. Solo la parte central (superior e inferior) es una posibilidad, pero el contenedor todavía está allí, por lo que la cantidad real es pequeña ". aunque tengo el presentimiento de que es una idea útil. Creo que está sugiriendo que si un límite de contenedor cae cerca del centro de una distribución de pico estrecho, ¿esto puede mejorar artificialmente el ancho aparente de esa distribución?
Sí, eso es lo que estoy tratando de decir. Las palabras pueden ser difíciles...
Eliminé "rigurosamente" del título y "Así que me gustaría saber cómo los expertos hacen esto como una cuestión de procedimiento estándar". del cuerpo Si realiza una edición simbólica de su respuesta, puedo revertir mi voto negativo a positivo.