LA OBSERVACIÓN DEL RADIANTE DE LOS NEOS PARA PREVENIR IMPACTOS El software NEORIES7P.EXE selecciona los neos que se aproximan a la órbita de la Tierra y genera un listado de radiantes de asteroides para ser observados entre unas fechas determinadas en previsión de su aproximación y poco probable impacto. La primera línea del listado que genera muestra: el nombre y clave del asteroide, su magnitud intrínseca H, su diámetro en metros y masa en toneladas, así como la energía mínima transformada si llegara a impactar (a velocidad infinita 0 y sólo por la atracción de la Tierra) expresada en megatones. En la segunda línea se obtiene: la fecha teórica de máxima aproximación del asteroide a la órbita terrestre, el punto de la esfera celeste (coordenadas de ascensión recta y declinación (radiante) desde donde incidiría un cuerpo que impactase sobre nuestro planeta si su órbita estuviera establecida, longitud eclíptica del punto de máxima aproximación, velocidad de incidencia o infinita en Km/s. , velocidad de impacto tras ser acelerado por la masa de la Tierra, el mes de paso de la Tierra y la diferencia en meses con el paso del asteroide. Cuando el asteroide supera unos parámetros establecidos de riesgo, obtenemos una tercera línea que indica: una escala de interés de 0 a 5 en función del riesgo, la fecha de aproximación de la Tierra al “punto de encuentro” y el periodo para la observación del radiante. Si la escala de interés es elevada (>1), es interesante estudiar más en profundidad la aproximación del asteroide a la órbita terrestre mediante la actualización de los elementos y su corrección consecuente de aproximación, con mayor precisión, además de consultar el riesgo obtenido de otros centros sobre el asteroide. Los estudios de riesgo de asteroides mayores de 50 metros (no se desintegran en el aire y alcanzan la superficie terrestre) son bastante precisos, pero el problema principal por su elevado número e imprecisión radica en los asteroides mayores de 15 metros (R>1) que inciden sobre la Tierra con un periodo medio próximo al año y pueden causar daños graves si caen sobre el cielo de las ciudades. La gran mayoría de asteroides de entre 15 y 50 metros no se ha detectado ya que son cientos de miles. No obstante, hay más de 5.000 que tienen órbitas establecidas pero con periodos de traslación muy poco precisos. Una media algo superior de un asteroide controlado con probabilidad aceptable de un impacto al año puede ser estudiada si disponemos del listado adjunto. El listado de 2013 nos muestra una decena de NEOS con escala de riesgo superior a 0, Después de analizados, se descartan los mayores ya que, si bien ofrecían un riesgo estadístico previo, son astros de órbitas precisas que los descartan de riesgo de choque. Hay varios asteroides menores de 15 metros (R<2), que podrían chocar con la Tierra y de hecho tanto Neodys, como JPL los confirman. No obstante, sus pequeñas dimensiones los excluyen como causantes de daños potenciales aunque fueran a parar sobre el cielo de una gran urbe. De cualquier forma, la lista de dichos asteroides se doblará cada pocos años y serán susceptible de análisis previo de impacto varios asteroides cada año, Si bien, un NEO de menos de 20 metros presenta una magnitud muy elevada>18 incluso en un radio inferior a medio millón de kilómetros, la observación del radiante facilita mucho su localización, tanto por su poco desplazamiento al venirnos de cara, lo que reduce hasta tres magnitudes su capacidad de detección, como al podernos centrar en un área de localización inferior a un grado al cuadrado, si los elementos orbitales no son muy dispares y el paso por el perihelio no diverge en más de unos pocos días. EJEMPLO DEL SOFTWARE SOBRE UN CASO REAL Mediante el software NEORIE7P.EXE, introducimos los elementos orbitales en la base de datos del asteroide 2008 TC3 que chocó con la atmósfera el 7 de octubre de 2008. La función 1 nos da un listado del radiante en sucesivos días con pocas variaciones que lo sitúan en: 23h16m de A.R. y +7º30´de declinación (eq2000). De haber tenido los elementos por un paso anterior previo al choque, se hubiera previsto el impacto observando las coordenadasdesde en días anteriores. Como se ve, el campo a observar es inferior a un grado, lo que facilita mucho la localización. Se obtiene además, la velocidad de impacto de 13 Km/s y los otros datos arriba indicados. NEORIE7P.EXE NEOS CON POSIBLE RIESGO. t.NEO A.R.(Radiante)Dec. Dis.orb. r. Long T. Vinf.Vimp. MesT. D.mes 2008 TC3 H, DIAM, MAS, MEGAT V0= 30.4 2 15 0 2008.1006 23.1637 +7.300 0.004 1.003 12.89 6.7 13.0 10.1 0.0 RIESGO (0-5), FECHA, 1 2008 10 5 <-0 7 > 2008 TC3 H, DIAM, MAS, MEGAT V0= 30.4 2 15 0 2008.1007 23.1566 +7.298 0.000 1.000 13.95 6.6 13.0 10.1 0.0 RIESGO (0-5), FECHA, 1 2008 10 6 < 1 8 > 2008 TC3 H, DIAM, MAS, MEGAT V0= 30.4 2 15 0 2008.1008 23.1491 +7.292 0.002 0.996 15.02 6.6 13.0 10.2 0.0 RIESGO (0-5), FECHA, 1 2008 10 7 < 2 9 > Como comprobación vemos las coordenadas de la circular MPEC 2008-T50 con el descubrimiento del asteroide cerca del radiante. M.P.E.C. 2008-T50 Issued 2008 Oct. 6, 14:59 UT The Minor Planet Electronic Circulars contain information on unusual minor planets and routine data on comets. They are published on behalf of Commission 20 of the International Astronomical Union by the Minor Planet Center, Smithsonian Astrophysical Observatory, Cambridge, MA 02138, U.S.A. Prepared using the Tamkin Foundation Computer Network MPC@CFA.HARVARD.EDU URL http://www.cfa.harvard.edu/iau/mpc.html ISSN 1523-6714 2008 TC3 Observations: K08T03C* C2008 10 06.27767 23 17 00.78 +07 49 22.7 18.9 V ET050G96 K08T03C C2008 10 06.28762 23 16 54.58 +07 49 25.8 18.8 V ET050G96 K08T03C C2008 10 06.29770 23 16 48.36 +07 49 27.6 18.8 V ET050G96 K08T03C C2008 10 06.30770 23 16 42.34 +07 49 28.7 19.1 V ET050G96 K08T03C C2008 10 06.32199 23 16 33.95 +07 49 28.0 18.8 V ET050G96 K08T03C C2008 10 06.33431 23 16 27.03 +07 49 25.5 18.8 V ET050G96 K08T03C C2008 10 06.34546 23 16 20.97 +07 49 29.8 18.1 R ET050854 K08T03C C2008 10 06.34667 23 16 20.48 +07 49 21.2 18.8 V ET050G96 K08T03C C2008 10 06.34885 23 16 19.30 +07 49 27.2 18.5 R ET050854 K08T03C C2008 10 06.35223 23 16 17.59 +07 49 26.9 19.0 R ET050854 K08T03C C2008 10 06.35904 23 16 14.38 +07 49 15.3 18.8 V ET050G96 K08T03C C2008 10 06.37547 23 16 07.17 +07 49 04.6 19.2 V ET050G96 K08T03C C2008 10 06.37584 23 16 07.03 +07 49 04.2 18.5 V ET050G96 K08T03C C2008 10 06.37620 23 16 06.88 +07 49 04.0 19.1 V ET050G96 K08T03C C2008 10 06.37656 23 16 06.74 +07 49 03.7 19.1 V ET050G96 K08T03C .C2008 10 06.61056 23 19 06.42 +09 15 20.0 17.8 R ET050D90 K08T03C .C2008 10 06.61326 23 19 04.70 +09 15 46.0 17.9 R ET050D90 K08T03C .C2008 10 06.61598 23 19 02.93 +09 16 12.4 17.9 R ET050D90 K08T03C C2008 10 06.63767 23 18 14.04 +09 15 39.0 18.2 V ET050E12 K08T03C C2008 10 06.63803 23 18 13.86 +09 15 42.4 17.4 V ET050E12 K08T03C C2008 10 06.63839 23 18 13.68 +09 15 45.7 18.0 V ET050E12 K08T03C C2008 10 06.63874 23 18 13.53 +09 15 49.5 18.0 V ET050E12 K08T03C C2008 10 06.63910 23 18 13.39 +09 15 52.2 17.6 V ET050E12 K08T03C C2008 10 06.64112 23 18 12.31 +09 16 09.5 17.9 V ET050E12 K08T03C C2008 10 06.64148 23 18 12.09 +09 16 13.0 17.7 V ET050E12 K08T03C C2008 10 06.64182 23 18 11.98 +09 16 15.9 17.7 V ET050E12 Observer details: 854 Sabino Canyon Observatory, Tucson. Observer J. E. McGaha. 0.36-m f/10.0 Schmidt-Cassegrain + CCD. D90 RAS Observatory, Moorook. Observers C. Jacques, E. Pimentel. Measurer C. Jacques. 0.25-m f/6.0 reflector + CCD. E12 Siding Spring Survey. Observer G. J. Garradd. Measurers E. C. Beshore, A. Boattini, G. J. Garradd, A. R. Gibbs, A. D. Grauer, R. E. Hill, R. A. Kowalski, S. M. Larson, R. H. McNaught. 0.5-m Uppsala Schmidt + CCD. G96 Mt. Lemmon Survey. Observer R. A. Kowalski. Measurers E. C. Beshore, A. Boattini, G. J. Garradd, A. R. Gibbs, A. D. Grauer, R. E. Hill, R. A. Kowalski, S. M. Larson, R. H. McNaught. 1.5-m reflector + CCD. Orbital elements: 2008 TC3 Earth MOID = 0.0000 AU Epoch 2008 Oct. 6.0 TT = JDT 2454745.5 MPC M 328.44590 (2000.0) P Q n 0.68871009 Peri. 233.94885 +0.37343159 -0.92760486 a 1.2699260 Node 194.13377 +0.86362066 +0.35153118 e 0.2847655 Incl. 2.32440 +0.33868601 +0.12639247 P 1.43 H 30.4 G 0.15 U 8 Residuals in seconds of arc 081006 G96 0.3+ 0.2+ 081006 854 0.2- 1.0+ 081006 E12 0.1- 0.3- 081006 G96 0.1- 0.1+ 081006 G96 0.1- 0.1- 081006 E12 0.1- 0.1- 081006 G96 0.4- 0.2- 081006 G96 0.3+ 0.2- 081006 E12 0.0- 0.0+ 081006 G96 0.0+ 0.1- 081006 G96 0.4+ 0.4- 081006 E12 0.4+ 0.7+ 081006 G96 0.1+ 0.1- 081006 G96 0.3+ 0.3- 081006 E12 1.0+ 0.2+ 081006 G96 0.0- 0.1- 081006 G96 0.4+ 0.3- 081006 E12 0.1- 0.5- 081006 854 1.0- 0.9+ 081006 D90 0.2+ 0.3- 081006 E12 0.8- 0.2- 081006 G96 0.0+ 0.2- 081006 D90 0.2+ 0.1+ 081006 E12 0.1+ 0.3- 081006 854 0.0+ 0.3- 081006 D90 0.7- 0.6+ Ephemeris: 2008 TC3 a,e,i = 1.27, 0.28, 2 q = 0.9083 Date TT R. A. (2000) Decl. Delta r Elong. Phase V 2008 10 06.50 23 19.43 +08 18.5 0.00246 1.002 157.4 22.5 18.4 2008 10 06.51 23 19.54 +08 19.1 0.00242 1.002 157.4 22.5 18.4 2008 10 06.52 23 19.65 +08 19.8 0.00239 1.002 157.5 22.5 18.4 2008 10 06.53 23 19.76 +08 20.5 0.00235 1.002 157.5 22.5 18.3 2008 10 06.54 23 19.88 +08 21.2 0.00231 1.002 157.5 22.5 18.3 2008 10 06.55 23 20.00 +08 21.9 0.00227 1.002 157.5 22.4 18.3 2008 10 06.56 23 20.13 +08 22.7 0.00224 1.002 157.5 22.4 18.2 2008 10 06.57 23 20.26 +08 23.5 0.00220 1.002 157.6 22.4 18.2 2008 10 06.58 23 20.40 +08 24.3 0.00216 1.002 157.6 22.4 18.2 2008 10 06.59 23 20.54 +08 25.1 0.00212 1.002 157.6 22.4 18.1 2008 10 06.60 23 20.68 +08 26.0 0.00209 1.001 157.6 22.3 18.1 2008 10 06.61 23 20.83 +08 26.9 0.00205 1.001 157.7 22.3 18.0 2008 10 06.62 23 20.99 +08 27.8 0.00201 1.001 157.7 22.3 18.0 2008 10 06.63 23 21.15 +08 28.8 0.00197 1.001 157.7 22.2 17.9 2008 10 06.64 23 21.32 +08 29.8 0.00193 1.001 157.7 22.2 17.9 2008 10 06.65 23 21.50 +08 30.9 0.00190 1.001 157.8 22.2 17.9 2008 10 06.66 23 21.68 +08 32.0 0.00186 1.001 157.8 22.1 17.8 2008 10 06.67 23 21.87 +08 33.1 0.00182 1.001 157.8 22.1 17.8 2008 10 06.68 23 22.07 +08 34.3 0.00178 1.001 157.9 22.1 17.7 2008 10 06.69 23 22.28 +08 35.5 0.00174 1.001 157.9 22.0 17.7 2008 10 06.70 23 22.49 +08 36.8 0.00171 1.001 158.0 22.0 17.6 2008 10 06.71 23 22.72 +08 38.2 0.00167 1.001 158.0 22.0 17.6 2008 10 06.72 23 22.96 +08 39.6 0.00163 1.001 158.1 21.9 17.5 2008 10 06.73 23 23.21 +08 41.1 0.00159 1.001 158.1 21.9 17.5 2008 10 06.74 23 23.47 +08 42.6 0.00155 1.001 158.2 21.8 17.4 2008 10 06.75 23 23.74 +08 44.2 0.00152 1.001 158.2 21.8 17.4 2008 10 06.76 23 24.03 +08 46.0 0.00148 1.001 158.3 21.7 17.3 2008 10 06.77 23 24.33 +08 47.8 0.00144 1.001 158.3 21.6 17.2 2008 10 06.78 23 24.65 +08 49.7 0.00140 1.001 158.4 21.6 17.2 2008 10 06.79 23 24.99 +08 51.7 0.00136 1.001 158.5 21.5 17.1 2008 10 06.80 23 25.35 +08 53.8 0.00133 1.001 158.5 21.4 17.1 2008 10 06.81 23 25.73 +08 56.0 0.00129 1.001 158.6 21.3 17.0 2008 10 06.82 23 25.87 +08 56.9 0.00125 1.001 158.7 21.3 16.9 2008 10 06.83 23 26.31 +08 59.5 0.00121 1.001 158.7 21.2 16.9 2008 10 06.84 23 26.79 +09 02.3 0.00117 1.001 158.8 21.1 16.8 2008 10 06.85 23 27.33 +09 05.4 0.00113 1.001 159.0 21.0 16.7 2008 10 06.86 23 27.91 +09 08.9 0.00109 1.000 159.1 20.9 16.6 2008 10 06.87 23 28.54 +09 12.6 0.00106 1.000 159.2 20.8 16.5 2008 10 06.88 23 29.23 +09 16.7 0.00102 1.000 159.4 20.6 16.5 2008 10 06.89 23 29.99 +09 21.1 0.00098 1.000 159.5 20.4 16.4 2008 10 06.90 23 30.80 +09 25.9 0.00094 1.000 159.7 20.3 16.3 2008 10 06.91 23 31.68 +09 31.0 0.00090 1.000 159.9 20.1 16.2 2008 10 06.92 23 32.64 +09 36.5 0.00086 1.000 160.1 19.9 16.1 2008 10 06.93 23 33.66 +09 42.5 0.00082 1.000 160.3 19.7 16.0 2008 10 06.94 23 32.61 +09 36.4 0.00080 1.000 160.1 19.9 15.9 2008 10 06.95 23 33.76 +09 43.0 0.00076 1.000 160.3 19.6 15.8 2008 10 06.96 23 35.11 +09 50.9 0.00072 1.000 160.6 19.4 15.7 2008 10 06.97 23 36.71 +10 00.1 0.00068 1.000 161.0 19.0 15.5 2008 10 06.98 23 38.57 +10 10.8 0.00063 1.000 161.4 18.6 15.4 2008 10 06.99 23 40.74 +10 23.2 0.00059 1.000 161.8 18.1 15.2 2008 10 07.00 23 43.27 +10 37.6 0.00054 1.000 162.4 17.6 15.0 2008 10 07.01 23 46.20 +10 54.1 0.00050 1.000 163.0 17.0 14.8 2008 10 07.02 23 49.60 +11 13.1 0.00045 1.000 163.7 16.3 14.5 2008 10 07.03 23 53.58 +11 35.1 0.00040 1.000 164.4 15.6 14.3 2008 10 07.04 23 58.26 +12 00.7 0.00035 1.000 165.3 14.7 14.0 2008 10 07.05 00 03.86 +12 30.7 0.00031 1.000 166.3 13.7 13.6 2008 10 07.06 00 10.69 +13 06.7 0.00026 1.000 167.4 12.6 13.2 2008 10 07.07 01 07.35 +17 29.9 0.00019 1.000 167.4 12.6 12.5 2008 10 07.08 01 31.80 +19 02.5 0.00015 1.000 163.3 16.7 12.1 2008 10 07.09 02 06.87 +20 51.7 0.00010 1.000 156.2 23.8 11.6 2008 10 07.10 03 08.55 +22 53.5 0.00007 0.999 142.8 37.2 11.0 2008 10 07.11 05 36.95 +21 35.3 0.00004 0.999 109.4 70.6 10.8 The nominal orbit given above has 2008 TC3 coming to within one earth radius around Oct. 7.1. The absolute magnitude indicates that the object will not survive passage through the atmosphere. Steve Chesley (JPL) reports that atmospheric entry will occur on 2008 Oct 07 0246 UTC over northern Sudan. Gareth V. Williams (C) Copyright 2008 MPC M.P.E.C. 2008-T50