Short communication
, Volume: 6( 3)The Potentially Dangerous Asteroid (29075) 1950 DA
- *Correspondence:
- Wlodarczyk Ireneusz , 41-500 Chorzow, Poland, Tel: +48627221279; E-Mail: astrobit@ka.onet.pl
Received: July 10 2018; Accepted: August 01 2018; Published: August 08 2018
Citation: Wlodarczyk Ireneusz . The Potentially Dangerous Asteroid (29075) 1950 DA. J Phys Astron. 2018; 6(3):161.
Abstract
I present behavior of the dynamically interesting potentially dangerous asteroid (29075) 1950 DA. Using the public available OrbFit software I studied orbit of asteroid 900 years forward in the future searching for close approaches with the Earth, which lead to possible impact in 2880.
Keywords
Asteroids; Dynamic; Celestial mechanic; Near-Earth Objects; Orbit determination
Introduction
The asteroid (29075) 1950 DA was discovered on February 22, 1950 by the (IAU 662) Lick-Observatory, Mount Hamilton (IAUC 1258). It was observed for 19 days in 1950 till March, 12 and then lost because of the short observational arc [1]. On December 31, 2000, it was observed by J.E. Rogers from (IAU 670) Camarillo Observatory, recovered as 2000 YK66 and 2 hours later was recognized as 1950 DA( MPEC 2001-A-26). Between these dates the asteroid (29075) 1950 DA was observed on September 30, 1981 by Siding Spring Observatory-DSS ( IAU 260). Asteroid (29075) 1950 DA belongs to the Apollo group, comprising 9194 members as of July 5, 2018), and is one of 18345 known Near-Earth Asteroids at this time http://www.minorplanetcenter.net/iau/lists/Unusual.html. According to the JPL NASA: https://ssd.jpl.nasa.gov/sbdb.cgi#top, its absolute magnitude is 17.1, suggesting a diameter of about 2.0 km (see also next Section). Asteroid (29075) 1950 DA is dynamically interesting. According to my computattions, it will pass close to the Earth on March 02.28506, 2032 at about 0.0758 au (11.3×106 km).
Computations of Giorgini from JPL NASA [1] are based on observations till 2002 (214 optical observations and 12 radar observations). Now, July 5, 2018 we have 576 and 12 observations, respectively, and the longer observational arc. Sitarski [2] determined the asteroid’s orbit based on 231 optical observations from the period February 22, 1950 to December 18, 2004 and applied his cracovian methods to find an impact orbit for 2880.
The Initial Orbital Elements of Asteroid (29075) 1950 DA
The orbit of the asteroid (29075) 1950 DA was computed using the freely distributed OrbFit software, version 5.0.5 (http://adams.dm.unipi.it/~orbmaint/orbfit/ ). This new version includes the new error model and the debiasing and weighting scheme described by [3] and the DE431 version of JPL's planetary ephemerides. The NEODyS site (http://newton.dm.unipi.it/neodys/) and the OrbFit software used above given star catalog debiasing and an error model with assumed astrometric RMS errors computed from the tests in the above paper. We can compute the Yarkovsky non-gravitational effect as a model uncertainty, that is solves for 7 parameters instead of 6, i.e. only orbital parameters. This applies both to the orbit determination (at the initial epoch) and to the propagations used in the impact monitoring procedure.
Nongravitational effects in orbital motion of asteroid can be taken into account in various ways. The simplest one is to assume the existence of a secular change of semimajor axis of the orbit, da/dt. We can also compute the non-gravitational parameter A2. Negative value of A2 for asteroid denotes that the mean semimajor axis drift da/dt<0 and hence asteroid can be retrograde rotator. Observational arc length is 24824 days (almost 68 years) so it is possible to compute non-gravitational parameter A2.
The orbital elements given in TABLE 1 were computed using the JPL Planetary and Lunar Ephemerides DE431 and additional perturbations of the 16 massive asteroids were used by [4-6]: (1)Ceres, (2)Pallas, (4)Vesta, (10)Hygiea, (3)Juno, (9)Metis, (15) Eunomia, (16) Psyche, (29) Amphitrite, (31) Euphrosyne, (52) Europa, (65) Cybele, (87) Sylvia, (511) Davida, (532) Herculina and (704) Interamnia. And additionally I used perturbation of Pluto.
semimajor axis | 1.69853729475696 ± 1.08961E-9 au |
eccentricity | 0.5079907792308 ± 3.20793E-8 |
inclination | 12.16964176754 ± 4.56135E-6 deg |
longitude node | 356.67753537197 ± 8.10688E-6 deg |
argument of perihelion | 224.65716201213 ± 9.28745E-6 deg |
mean anomaly | 234.67603804158 ± 2.35067E-6 deg | ±
non-gravitational parameter | A2 (-6.03004 ± 1.24958)E-15 au/d^2 |
Table 1: Keplerian orbital elements.
TABLE 1 lists the orbital parameters of asteroid (29075) 1950 DA, which were computed using the OrbFit software based on all 576 optical observations (of which 6 are rejected as outliers). from February 22.23014, 1950 to February 09.19461, 2018, and also on 12 radar observations from March 3, 2001to 2012 May 1, 2012.
The TABLE 1 also includes their 1-σ uncertainties, which are the square roots of the diagonal elements of the computed covariance matrix. RMS is the quadratic mean of the residuals between the calculated orbit and all 588 observations used in the solution.
I computed directly from the OrbFit software an absolute magnitude H of (17.064 ± 0.507) mag. Hence, I calculated the diameter D of the asteroid from its H and albedo pv using the formula of [7]:
D = 1329 * 10-H/5 Pv-1/2,
where D is in km.
For H = (17.064 ± 0.507) mag and albedo of 0.07 (EAR_A_COMPIL_5_NEOWISEDIAM_V1_0[neos] (http://adsabs.harvard.edu/abs/2011ApJ...743..156M), I find the diameter of (29075) 1950 DA to be between 1.5 km and 2.5 km.
Close Approaches Analysis of Asteroid (29075) 1950DA till 2300
FIG. 1 presents close approaches (CA) between the Earth and asteroid (29075) 1950 DA between 2000 and 2900. I showed only CA<0.02 au. They are computed using orbital elements and non-gravitational parameter A2 from TABLE 1 The first CA is in 2641 at 0.0155 au. The deepest CAs begin in 2860 at 0.0033 au with possible impact with the Earth in 2880.
Possible Impact Solutions
Next, using starting orbital elements from TABLE 1 and 601 Vitual Asteroids (clones) (VAs) with 3-σ uncertainties, and propagating them till 2900 I searched for impacts with the Earth. TABLE 2 presents possible impact solution for 2880.
Date of impact | March 16.992, 2880 |
Impact Probability | 2.00E-05 |
Mass | 3.05E12 kg |
Impact velocity | 17.99 km/s |
Impact energy | 1.18E5 MT |
1 MT | 4.2E15 J |
Table 2: Possible impact solution for 2880.
Conclusion
The main goal of my work is to present methods which allow us to follow the orbit of the potentially dangerous asteroid (29075) 1950 DA. I computed the minimum distances between the asteroid and the Earth in the time-span 900 years forward from now. I found impact solution for 2880. It is possible to recover asteroid (29075) 1950 DA in 2032 when the distance to the Earth will be of about 0.076 au. Optical and radar observations made during this close approach can refine the orbit of the asteroid and give more precise possible impact solutions.
References
- Giorgini J,. Asteroid 1950 DA?s Encounter with Earth in 2880. Physical Limits of Collision Probability Prediction. Science. 2002;296:132-136.
- Sitarski, G. Generating of ?clones? of an Impact Orbit for the Earth-Asteroid Collision. Acta Astronomica. 2006;56:283-292.
- Farnocchia, D., Chesley, S. R., Chamberlin, A.B., Tholen, D. J. Star catalog position and proper motion corrections in asteroid astrometry. Icarus. 2015;245:94-111
- Farnocchia, D., Chesley, S. R., Chodas, P.~W., Micheli, M., Tholen, D.~J., Milani, A., Elliott, G.~T., Bernardi, F. Yarkovsky-driven impact risk analysis for asteroid (99942) Apophis.Icarus 2013;224:192-200.
- Farnocchia, D., Chesley, S. R., Vokrouhlick'y, D., Milani, A., Spoto, F., Bottke, W.~F. Near Earth Asteroids with measurable Yarkovsky effect. Icarus. 2013;224:1-13.
- Wlodarczyk, I. The Potentially Hazardous Asteroid (410777) 2009 FD. Acta Astronomica. 2015;65:215-231.
- Fowler, J. W., Chilemi, J. R. IRAS asteroid data processing, in The IRAS Minor Planet Survey, ed. E. F. Tedesco, Phillips Laboratory, Hanscom AFB, Technical Report. 1992. PL-TR-92-2049