Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.puldb.ru/laza_gao/PUB/work_2.pdf Дата изменения: Sun Apr 12 12:30:34 2009 Дата индексирования: Mon Oct 1 20:37:36 2012 Кодировка: Поисковые слова: sojourner

Initial orbit determinations of NEO-2005 with Pulkovo AMP-method
O.P. Bykov
Central Astronomical Observatory of RAS, Russia, 196140, Saint-Petersburg, Pulkovskoye shosse 65/1

Abstract Pulkovo Apparent Motion Parameters Method (the AMP-method) was proposed by prof. A.A.Kiselev and developed by his colleagues for initial orbit determinations of Artificial Earth Satellites, Space Debris, Asteroids, and Comets on the base of very short observational arc. Observational data for the AMP-method are the celestial body position in a middle moment of its observations, its angular velocity and acceleration for this moment, positional angle of its motion and a curvature of its apparent tra jectory. These values are calculated from mathematical processing of crowded accurate coordinat sets of an observed celestial body. The AMP-method is a further development of the Classic Laplace's orbital method. Here the AMP-method was applied for the fast orbit determinations of the NEOs which were discovered in 2005. The first nights positions (1 or 3 dates) from Minor Planet Electronic Circulars were used for this job. Our preliminary osculating orbits were compared with the same orbits published in MPEC. Detailed results are given in (http://www.nm.neopage.ru and http://www.accuracy.puldb.ru)

Key words: Astrometric CCD observations. Initial orbit determination. Asteroids.

1. Intro duction In the beginnings of 1970th the method of the fast determination of Artificial Earth Satellites orbits was created and developed by Prof. A.A.kiselev and his colleagues at Pulkovo Astronomical Observatory (A.A. Kiselev, O.P.Bykov., 1973), (A.A. Kiselev, O.P.Bykov., 1976). Method was named Apparent Motion Parameter Method (the AMP-method) and successfully was applied to calculations of osculating AES orbits on the base of their photographic positional observations distributed along a short topocentric arc. But obtaining photographic observations with necessary density and accuracy were very difficult for using AMP-method in practice. Further the AMP-method was applied to determination of asteroid' and comet' orbits (O.P.Bykov., 1989). Now modern CCD positional asteroid and AES observations are resuscitating of the AMP and classical Laplace's methods for the operative orbit determination of observed ob ject in a real time, sometimes directly near observing telescope. The AMP-method requires a known six observational parameters for a given moment of U T , namely and - usual spherical coordinates of celestial body, µ and µ its angular topocentric velocity and acceleration, - poEmail address: oleg@OB3876.spb.edu (O.P. Bykov). Preprint submitted to Elsevier

sitional angle of its apparent motions and c - a curvature of its visible tra jectory. All these values may be calculated from mathematical processing of crowded accurate positions distributed on a short observational arc. System of equations of AMP-method allows us also to obtain these parameters as an exact function of celestial body oculating elements and topocentric observer's coordinates. Table 1 shows these apparent motion parameters for Numbered asteroid 3317 Paris for dates of the Workshop meeting. Osculating elements of this asteroid was taken from Bowell's orbital catalog (Lowell Astronomical Observatoty) via Internet. Parameters were calculated with our Pulkovo EPOS Software (L'vov V.N. et. al., 2001). EPOS (Ephemeris Program for Ob jects of the Solar system) is the ephemeris program system with wide possibilities for professional and amateur observers of small celestial bodies - from calculation of the ob jects' ephemerides of various type and for control the accuracy of positional observations of the Solar system's ob jects, including their identifications with the use of Apparent Motion Parameters. The ephemerides of minor bodies are based on numerical integration of equations of their perturbed motion. The EPOS software package works as a Windows application on the IBM PC compatibles. Table 1 demonstrates a possibility to get new ephemeris data. It is our modern approach to the Ephemeris Service of asteroid and comet observations.
15 January 2007


Table 1 Apparent motion parameters for NMP 3317 Paris Parameters 2006 November µ, per day µ c
h

Table 2 Potentially Hazard Asteroids discovered in the first half of 2005 No Nomination Jan.-Febr. 1 2 3 4 5 6 7 2005 AD13 2005 AN26 2005 AV27 2005 AY28 2005 BC 2005 BY2 2005 BG14 2005 CJ 2005 CL 2005 CZ36 2005 DD March 12 13 14 15 16 17 18 19 20 2005 EA 2005 EE 2005 EO33 2005 EK94 2005 EY95 2005 EJ225 2005 ED318 2005 FH 2005 FE3 2005 GL 2005 GU 39 40 41 42 43 32 33 34 35 36 37 38 23 24 25 26 27 28 29 30 31 No Nomination April 2005 GJ8 2005 GY8 2005 GO21 2005 GP21 2005 GO22 2005 GE59 2005 GD60 2005 GH81 2005 GC120 May 2005 JU1 2005 JQ5 2005 JF21 2005 JE46 2005 JU81 2005 JS108 2005 KJ10 June 2005 LW3 2005 LX36 2005 LW39 2005 LO40 2005 MO13

Date, DT 16.20000 12 23 54.06
m

17.20000 12 24 37.446 14 09 35.02 626.52867 -36.76324 93.58182 1.774134

18.20000 12 25 20.529 14 08 56.81 621.84039 -36.65173 93.20344 1.714423

5s

14 10 17.60" 631.13530" -36.86803" 93.95019


1.831310

A problem of an identification of observed asteroids may be easy solved with the use of Apparent Motion Parameters. With EPOS Software we also can calculate full list of ephemerides: for example, a phase angle (8.8 ), visual magnitude (16.4m ), elongation (56.2 W), topocentric distance and its the first and the second derivatives. Horizontal coordinates are also available for MPC codes observatories. We present possibilities of the AMP-method of osculating orbit calculations for the NEOs discovered in the first half of 2005. Forty fast moving asteroids were fixed in this period by MPC. Table 2 presents a list of these NEOs. Our investigation of these NEOs included also an estimation of accuracy of their CCD observations made by various MPC code' stations. Here we cannot give all our results on this theme (see http://www.mn.neopage.ru) but they are usual, i.e. approximately in two time worse then that of for numbered asteroids. 2. Problem of the sort arc asteroid' observations Now a wast expancion of CCD matrices in astronomy allows to elaborate new approach to traditional problem of useing modern observational information obtained from very short arc positions. One can easy get coordinate' sets of any moving celestial ob ject with any dense distribution along short topocentric arc. These data can give us the Apparent Motion Parameters in a middle moment of observations, at least coordinates , , µ, . Several minutes of CCD asteroid' accurate observations are enough for getting these values from mathematical processing of coordinates by means of their linear approximation. Consequently a circular AMP-orbit may be calculated immediately, and it will be suitable to continuation of further observations of fixed ob ject during several close nights by observer himself. Of course, an observer, especially amateur, must inform the Minor Planet Center about his asteroid discovery but he will not wait new ephemerides from MPC for new observations. But Minor Planet Center does not take into consideration one night asteroid positions, up to now only two nights observations are processed in MPC for preliminary orbit determination. If an accuracy of CCD observations and their distribution along topocentric arc allow to calculate additionally values angular accelerations µ and curvature c of visible tra 2

8 9 10 11

21 22

jectory we can obtained preliminary elliptic orbit. For confirmation these possibilities of AMP-orbits we had considerd results of NEOs CCD observations in the first nights of their discoveries for determinations of our orbits as it could be in real time of observations. Initial orbits for these asteroids were taken from Minor Planets Electronic Circulars (MPEC). Our AMP-orbits were compared with the MPCorbits obtained on the base the same CCD observations. As a rule our data practically coincidence with MPEC elements. In Tables below we give an example of our results. 3. Example of initial orbit determination Here we present one fragment from our orbital calculations. Original CCD observations were taken from (T.B.Spahr, 20051). Table 3 illustrates a used observational data, namely the Normal Places (NPLs) of asteroid 2005 EE. The NPLs were obtained by polynomial approximation of real CCD observations made by observatories with given MPC codes. Only MPC code 703 is a professional astronomical observatory, other are amateurs of Astronomy with small CCD telescopes. All theire data are given in this Table. Only 25 separate positions were


obtained by these telescopes for 3 close nights during 1-3 March. We got Normal Places (NPL) for each observatory by means of linear approximations of spheric coordinate' sets. This procedure are smooth for improvement of observations. Sometimes it is very important for useing the AMP method in a case of many participating observatories with different telescopes and CCD matrices. Here only 4 or 3 positions per night (symbol "j" in the Table) were used for NPLs getting. The NPLs are presented in the Table 3 together with (O-C)s and their errors. The (O-C) values were calculated by EPOS Software with Bowell's orbital elements which were derived from all available observations of 2005 EE (46 positions from March to May). The errors of (O-C)s were obtained as errors of the average meanings for considered night. Naturally, the (O-C) values are reflecting differencies between observed and calculated positions. As seen an observational set are good for solving orbital task. It would be noted that the seven circular orbits can be calculated for each MPC code observatory with the use of one night observations but these circular orbits (real eccentricity is equal 0.33) may be available for the next observational night only. Next Table 4 presents the Apparent Motion Parameters calculated from Table 3 data. These AMPs are compared with their meanings calculated by EPOS with the use of Bowell's orbital elements.
Table 4 Apparent Motion Parameters for 2005 EE. Ep o ch 2005 03 02.19464 Calculated = 11h 07
m

Table 5 Osculating orbit calculated with various metho ds Method Positions Epo ch AMP Improvement 25 2005 03 02.19464 a e i N w M 1.129298 0.326745 6.14621 110.41971 284.90132 90.23289


MPEC 21 2005 02 02.00000 1.143650 0.345471 6.49906 115.60292 282.75954 75.13349

46 2005 03 02.19464 1.130414 0.328254 6.17246 110.98637 284.68501 89.61010

see that the (O-C) values calculated with MPEC orbit are very large in positions and and almost normal in , . These special features are used in our asteroid identification process.
Table 6 Orbital presentations of used observations No AMP for 1 2 3 4 5 6 7 Impr. for AMP for Impr. for 703 649 651 H06 703 448 J95 co de

+0.08" -0.29" -0.20 -0.37 +0.53 +0.71 -0.61 -0.04 +0.21 -0.04 +0.82 +0.96 -0.16 +0.14

+0.35" -0.17" -0.17 +0.20 -0.40 -0.10 +0.46 -0.42 -0.42 -0.07 -0.68 -0.34 +0.28 -0.24

Exact 27.253 00.57 2425.924 -92.818 283.60982 1.841095

26.650s 02.21" 2425.756" -91.634" 283.6160


= 29 57 µ, per day µ c

1.838965

Table 8 Verification of AMP-orbits of NEOs-2005 Nomin. Obs. c a e 2005 arc 7.4 1.3 1.4 2.3 0.2 (a.e) 0.000 -0.001 -0.000

i

N

(deg.) (deg.) 0 .0 0 .0 1 .7

We can see that NEO 2005 EE had angular velocity near 0.7 per day that is not high for such celestial ob jects. Our AMPs and their "exact" values are coinciding. But small differencies dive the various elements of orbits which are presented in Table 5. MPEC elements presented in the last column of Table 5 were calculated by T.B.Sphar (T.B.Spahr, 20051). His orbit is worse than AMP-orbit but he did not use 3 positions obtained with MPC code J95. And Spahr's result is good for 2 nights observations and his method of orbit determination. I don't know what a method was used by Dr.T.Sphar for these orbital calculations. Table 6 shows a presentation of used NPLs by calculated orbits. Table 7 demonstrates an ephemeris service with elements from Table 5. Our experiense allows us to state that the errors of initial orbital elements have an influence on the AMPs much less than on spherical coordinates. Data presented in Table 7 illustrates this statement. We can 3

AD13 AN26 AV27 AY28 BG14 CL DD EA EE EK94 EY95 EJ225 ED318 FH

0.004 -0.082 -0.022 -0.3

0.001 -0.117 -0.022 -0.1 -0.4 0.000 -0.001 0.001 0.080 0.157 0.001 0.002 0.001 0.013 0.006 0.030 0.000 0.002 0.001 0 .3 - 0 .0 2 .9 4 .8

4.0 -0.000 1.3 -0.000 1.6 -0.000 1.2 -0.002 1.3 1.4

0 .4 - 0 .5 0.0 0.0

0.0 -0.0 0.0 -0.5

0.011 -0.079 -0.015 -0.4 -2.2 0.000 0.034 0.043 0.3 -1.7 0.2 -0.3

1.1 -0.011 -0.054 -0.005 1.1 1.3

0.089 -0.075 -0.023 -0.1 -0.4 0.017 -0.071 -0.009 -0.3 -0.3


Table 3 Normal places and (O-C)s for PHA 2005 EE and their errors No Date 2005 1 03 01.32761 (O-C) (O-C) MPC code Posit. 703 4 Telescop es Catalina Sky Survey, Arizona, USA D=0.4, FL=1.2, FOV=175x175', Scale=2.6" 649 4 Powell Observatory, Louisburg, USA D=0.3, FL=3.8, FOV=30x30', Scale=1.8" 651 4 Grasslands Observatory, Tucson, USA D=0.6, FL=3.1, No data H06 3 New Mexico Skies Observatory, Germany D=0.25, FL=0.8, FOV=56x37', Scale=2.2" 703 4 Catalina Sky Survey, Arizona, USA D=0.4, FL=1.2, FOV=175x175', Scale=2.6" 448 3 Desert Mo on Observatory, Las Cruces, USA D=0.3, FL=2.0, FOV=21x21', Scale=2.4" J95 3 Great Shefford, Berkshire, England D=0.3, FL=1.8, FOV=25x25', Scale=1.5"

11 10 05.970 -0.29" +29 48 17.06 -0.17" ±0.54 ±0.09 +29 56 50.53 -0.42 ±0.10 +29 56 52.71 -0.07 ±0.15 +29 57 10.10 -0.68 ±0.09 +29 57 53.20 -0.34 ±0.24 +29 57 59.05 +0.28 ±0.28 +30 04 47.19 -0.24 ±0.11

2

03 02.17446

11 07 30.305 +0.21 ±0.21

3

03 02.17765

11 07 29.707 -0.04 ±0.28

4

03 02.20918

11 07 24.054 +0.82 ±0.39

5

03 02.28467

11 07 10.390 +0.96 ±0.64

6

03 02.29397

11 07 08.622 -0.16 ±0.22

7

03 03.06168

11 04 51.454 +0.14 ±0.38

Table 7 Ephemerides with initial orbits from Epo ch 2005 03 02.19464 to Ep och 2005 04 15.01153 MPC co de Observations Calculations (O-C) Calculations Apr. 15.01153 J50 Date Apr. 15.01153 Apr. 15.01153
h m

10 23 18.9
s

+2 5 27 31"

29.7
s



Orbit

-774" Improv. -770 -4" -758 MPEC -16" AMP

10 22 57.7 21.2 10 25 40.7 - 141.8

+25 28 54 - 83" +25 29 59 - 148"

29.3 0 .4
s

31.0 - 1 .3
s

Table 8 summarizes results of our job. It consists columns of "absolute" errors of observational and calculational values of curvatures and orbital elements. These meanings were obtained with formulae: ""AMP elements" minus "Improve elements"". If an accuracy of CCD used observations and topocentric arc are optimal for applying of AMP method, these tabular errors are small or close to zero. As a hole our initial AMP-orbits give a good ephemeris prognosis during one or two week. It makes them practically important.

References A.A. Kiselev, O.P.Bykov. The determination of the satellite orbit by a single photograph with many satellites trails. 1973. Soviet Astron. J., vol. 53, 1298 - 1308. A.A. Kiselev, O.P.Bykov. The determination of a satellite elliptical orbit with the use of parameters of the satellite's apparent motion. 1976. Soviet Astron. J., vol. 56, 879 888. O.P.Bykov. A determination of celestial bodies orbits with Direct methods. 1989. In.: "Problems of a constraction of coordinates system in Atronomy". 328 - 355. L'vov V.N., Smekhatchiova R.I., Tseikmeister S.D. EPOS as the program package for Solar system ob jects investigations. 2001. In proceedings of the conference "Near Earth astronomy in XXI century". 235 - 240. M.P.E.C. 2005-E06

4. Conclusion We are sure that Pulkovo Apparent Motion Parameters Method may be recommended for initial orbit determinations to observers especially amateurs of Astronomy. The matter is only in expansion of according SOFT for example, Pulkovo EPOS. It seems to us the circular orbits may be used for identification of moving ob jects in Gaia observational campaign. But this question requares special investigation. 4