The following information shows the result of the orbit fit based on Gary Bernstein's method. Most of the information should be self-explanatory. Take special note that while the original Bernstein software works with barycentric coordinates, we convert these results into a heliocentric coordinate system.
# Object: 98WS31
# Created Sun Nov 22 02:37:10 2009
# Orbit generated from Bernstein formalism
# Fitting 24 observations of 24
# Arc: 10.01y
# First observation: 1998/11/18
# Last observation: 2008/11/21
Preliminary a, adot, b, bdot, g, gdot:
-0.003234 0.042357 -0.000323 0.004661 0.029623 0.000000
# Chi-squared of fit: 69.54 DOF: 42 RMS: 0.30
# Min/Max residuals: -1.42 0.73
# Exact a, adot, b, bdot, g, gdot:
1.945178E-05 3.990220E-02 3.216591E-06 4.381018E-03 3.272323E-02 -1.088043E-03
# Covariance matrix:
3.1613E-13 -2.1845E-13 2.0613E-14 -2.2072E-14 2.2501E-13 -4.4831E-13
-2.1845E-13 2.2117E-13 -2.0484E-14 2.3505E-14 -1.3870E-13 4.7988E-13
2.0613E-14 -2.0484E-14 9.6921E-14 -1.4745E-14 1.7295E-14 -4.5024E-14
-2.2072E-14 2.3505E-14 -1.4745E-14 7.0788E-15 -1.4063E-14 5.2148E-14
2.2501E-13 -1.3870E-13 1.7295E-14 -1.4063E-14 4.2959E-13 -2.7868E-13
-4.4831E-13 4.7988E-13 -4.5024E-14 5.2148E-14 -2.7868E-13 1.0654E-12
# lat0 lon0 xBary yBary zBary JD0
2.641138 37.660085 -0.309438 0.042821 -0.933604 2451135.598851
# Heliocentric elements and errors
Epoch: 2451130.5000 = 1998/11/13
Mean Anomaly: 356.53551 +/- 0.006
Argument of Peri: 27.80828 +/- 0.010
Long of Asc Node: 15.90765 +/- 0.000
Inclination: 6.72435 +/- 0.000
Eccentricity: 0.20752251 +/- 0.0001
Semi-Major Axis: 39.71965852 +/- 0.0051
Time of Perihelion: 2452010.4205 +/- 1.6
Perihelion: 31.47693536 +/- 0.0053
Aphelion: 47.96238169 +/- 0.0070
Period (y) 250.3321 +/- 0.05
# Ecliptic coordinates at JD0 (AU and AU/d)
Ecliptic X 24.72539828 +/- 0.0005
Ecliptic Y 19.46691984 +/- 0.0004
Ecliptic Z 1.40830379 +/- 0.0000
Ecliptic XDOT -0.00212569 +/- 0.0000
Ecliptic YDOT 0.00258614 +/- 0.0000
Ecliptic ZDOT 0.00036193 +/- 0.0000
# Distances at JD0 (AU)
Heliocenter to KBO 31.50062869 +/- 0.0004
Geocenter to KBO 30.55932948 +/- 0.0006
# Hcoef: 8.26
The following table shows the complete astrometric record for 98WS31. The first three columns show the date of observation. The next six columns are RA and DEC. The next column (when provided) is the observed magnitude and filter. The next column is the object name (98WS31) followed by the observatory code and reference code for the source of the astrometry.
1998 11 18.09812 02 17 37.85 +16 33 43.0 22.5R 98WS31 695 a3608 1998 11 18.16698 02 17 37.44 +16 33 41.4 98WS31 695 a3608 1998 11 18.24713 02 17 36.96 +16 33 39.1 98WS31 695 a3608 1998 12 17.84201 02 15 13.98 +16 21 24.9 23.1R 98WS31 950 a3608 1998 12 17.87639 02 15 13.85 +16 21 24.2 98WS31 950 a3608 1998 12 17.91186 02 15 13.72 +16 21 23.7 98WS31 950 a3608 1999 11 30.19125 02 25 14.98 +17 27 05.5 22.8R 98WS31 695 Ca8499 1999 11 30.19836 02 25 14.98 +17 27 05.1 22.7R 98WS31 695 Ca8499 1999 11 30.30921 02 25 14.37 +17 27 02.3 22.4R 98WS31 695 Ca8499 1999 11 30.31631 02 25 14.34 +17 27 02.0 22.6R 98WS31 695 Ca8499 2000 01 02.15602 02 23 07.83 +17 15 47.1 23.3R 98WS31 697 Ca8499 2000 01 02.20304 02 23 07.62 +17 15 46.0 23.2R 98WS31 697 Ca8499 2000 01 02.22104 02 23 07.72 +17 15 46.0 98WS31 568 Ca8499 2000 01 02.23592 02 23 07.69 +17 15 45.7 23.2R 98WS31 568 Ca8499 2000 01 02.24998 02 23 07.65 +17 15 45.5 98WS31 568 Ca8499 2000 01 03.23913 02 23 05.61 +17 15 33.3 23.1R 98WS31 697 Ca8499 2000 01 05.11914 02 23 02.19 +17 15 10.8 23.2R 98WS31 697 Ca8499 2000 01 05.18084 02 23 02.09 +17 15 10.0 23.1R 98WS31 697 Ca8499 2002 09 11.39596 02 59 14.29 +20 39 18.5 22.2R 98WS31 662 Cg2658 2002 09 11.49551 02 59 14.06 +20 39 17.8 98WS31 662 Cg2658 2008 11 21.65583 03 48 50.40 +25 01 56.8 24.1R 98WS31 250 S~06AK 2008 11 21.66208 03 48 50.35 +25 01 56.7 98WS31 250 S~06AK 2008 11 21.66902 03 48 50.31 +25 01 56.5 98WS31 250 S~06AK 2008 11 21.67527 03 48 50.26 +25 01 56.5 98WS31 250 S~06AK
The following table shows the residuals to the orbit fit. The first coumn is the point number. The second column is the time, in years, measured from the first observation. The third and fifth columns are the regularized positions used in the orbit fit. The fourth and sixth columns are the residuals, in arc seconds, for RA and Dec respectively.
1 0.0000 0.00 -0.04 0.00 -0.23 2 0.0002 -6.09 0.05 0.43 0.10 3 0.0004 -13.37 0.00 0.52 0.09 4 0.0814 -2197.33 0.08 -13.94 0.09 5 0.0815 -2199.33 0.00 -13.98 0.09 6 0.0816 -2201.26 0.05 -13.83 0.27 7 1.0324 7243.90 -0.38 905.83 0.33 8 1.0324 7243.77 0.08 905.45 -0.05 9 1.0327 7234.59 0.02 905.58 0.10 10 1.0328 7234.09 0.10 905.43 -0.05 11 1.1227 5300.26 0.19 846.77 -0.04 12 1.1228 5297.06 -1.42 846.70 0.01 13 1.1229 5298.41 0.35 846.23 -0.42 14 1.1229 5297.91 0.35 846.09 -0.53 15 1.1229 5297.30 0.22 846.08 -0.50 16 1.1256 5265.69 -0.13 843.94 0.15 17 1.1308 5212.03 0.14 838.41 0.04 18 1.1310 5210.41 0.28 838.11 -0.10 19 3.8146 38716.55 -0.23 3369.57 0.73 20 3.8149 38713.14 0.25 3369.78 0.09 21 10.0111 86579.76 0.06 9588.83 -0.06 22 10.0111 86578.95 -0.05 9588.86 -0.02 23 10.0111 86578.28 0.05 9588.75 -0.12 24 10.0112 86577.49 -0.04 9588.89 0.02
The following table comes from a 10My integration of the orbit of the object. Three columns are shown. The first column is the result of integrating the nominal orbit. The other two columns are based on clones of the nominal orbit that are +/- 3 sigma from the nominal orbit. If all three types agree then the classificiation is deemed secure. The basis for these calculations is described in more detail in AJ, 129, 1117 (2005). Any use made of these calculations should refer to and credit this publication and the Deep Ecliptic Survey Team.
98WS31 quality flag:3 Type: 3:2E 3:2E 3:2E axisobj 39.786 39.801 39.771 ecceobj 0.209 0.209 0.209 incobj 6.727 6.727 6.727 qmin 31.264 31.283 31.246 qmax 47.733 47.694 47.846 amean 39.455 39.455 39.455 amin 39.382 39.395 39.368 amax 39.535 39.522 39.549 emean 0.195 0.195 0.195 emin 0.179 0.179 0.179 emax 0.207 0.208 0.210 imean 6.871 6.880 6.859 imin 6.257 6.272 6.241 imax 7.726 7.732 7.713 excite_mean 0.229 0.229 0.229 fracstop 1.000 1.000 1.000 cjmean 2.992 2.992 2.992 libcent 0 180.0 179.8 179.9 libamp 0 30.4 23.3 38.3 libcent 1 -180.0 -180.0 -180.0 libamp 1 -180.0 -180.0 -180.0 libcent 2 -180.0 -180.0 -180.0 libamp 2 -180.0 -180.0 -180.0 libcent 3 -180.0 -180.0 -180.0 libamp 3 -180.0 -180.0 -180.0 libcent 4 -180.0 -180.0 -180.0 libamp 4 -180.0 -180.0 -180.0 kozaimean 189.6 189.1 190.2 kozaiamp 180.0 180.0 180.0