recently , @xcite reported on the use of minor planets to assess gravity in the outer solar system . that paper was mostly devoted to the use of astrometry of asteroids to investigate gravity in the that region and the potential use of such observations to ascertain the reality of the pioneer effect , an unexplained acceleration towards the sun that perturbs the motion of pioneer 10 and 11 beyond a distance of 20 au @xcite
however , some statements were made that applied to comets and the outer planets .
in particular , with respect to the outer planets , @xcite commented that the ephemerides of the outer planets are almost entirely based on optical observations and are much less accurate than those for the inner planets .
this , coupled with the pioneer effect not being observed in the inner solar system argues that the major planets are not good candidates for investigating gravity in the outer solar system .
@xcite , carried out detailed calculations on the impact that the pioneer effect would have on orbital motions of the outer planets .
they found a pioneer - like acceleration would produce significant secular and periodic effects and conclude that the absence of such evidence indicates that the pioneer effect is ruled out as a phenomenon effecting the outer planets .
the approach of these authors involves analytically and numerically investigating the gaussian rate equations for rates of change in orbital elements as a function of time and a perturbing force .
this approach will be discussed more fully later in section [ results ] .
additionally , @xcite observes that @xcite offered no calculations to substantiate its conclusions about the outer planets and indicates that other authors reach different conclusions about the use of the outer planets in this role .
@xcite also sought to illuminate the issues involved by using a simplified four parameter model of orbital motion .
referencing several papers @xcite , @xcite calls the dispute about the observability of pioneer effect - like perturbations in the motion of the outer planets indicative of a lack of consensus in discussions of the planet s motion .    both @xcite and @xcite
refer to @xcite to buttress their contention that a pioneer - like acceleration should be observable in the motion of the outer planets . in discussing options for a non - dedicated spacecraft investigation of the pioneer effect
, @xcite briefly considers whether the effect would be detectable in the motion of the outer planets and conclude that the pioneer effect would have an impact on outer planet motion that would be at variance to that observed .
this conclusion is largely based on the change in the @xmath0 product that would be required to result in the pioneer effect at the distance of neptune and the comparison of that value with the otherwise known uncertainty of the @xmath0 value .
the basis of their conclusion seems to be similar to that expressed in @xcite , where modifications to kepler s law are used to assess the possibility of several modifications to gravity .
@xcite restricts its attention to the case of the inner planets , however , where ranging data is available from both earth - based radar and spacecraft .
however , @xcite s conclusion that the pioneer effect would have been observed in the motion of the outer planets is too all - encompassing .
for example , some alternative theories of gravity provide for a variation in the force that a spacecraft would experience , as would distributed mass densities concentrated in the outer solar system .
the issue here is whether such accelerations are really observationally detectable independent of their source .
this paper s purpose is twofold .
first , we wish to deal with the problem of astrometric determination of the orbital elements of pluto , in particular by assessing the ability of astrometric observations to detect perturbations of the sort expected from the pioneer effect .
secondly , we wish to show that while using simplified models of the dynamics can lead to some insights , one must be careful to not over - simplify the issues involved lest one be misled by the analysis onto false paths .
specifically , this paper shows that , contrary to recent assertions in the literature , the current ephemeris for pluto does not preclude the existence of the pioneer effect .
we show that the orbit of pluto is currently not well enough characterized to make such an assertion .
thus we address the views alluded to above that have now propagated through a number of papers ( ) .
we believe that the conclusions presented in those papers are not established unambiguously and that caution needs to be used in drawing further inferences about the dynamics of the outer solar system .    notwithstanding our comments about the dangers of using simplified methods
, we note that our approach itself is a simplification .
the outer planets must be dynamically assessed as a system . in order to rigorously determine whether a small perturbation like the pioneer effect is detectable by astrometry
, we should simultaneously include the changes it produces in the orbits of uranus , neptune , and pluto . only in this way
can all the second order perturbations to the system be taken into account .
however , in order to illustrate the ideas concerned and the weaknesses of the approaches outlined above , we restrict our attention here to manipulating the orbit of pluto .
this approach is further discussed in section [ disc ] .
the remainder of this paper is divided into four further sections .
section [ m&m ] describes several `` back of the envelope '' approaches to understanding the nature of the orbital changes that a pioneer effect would cause , and a detailed discussion of the analysis methodology we have chosen to make our case .
section [ results ] addresses our results in terms of the relationships between orbital parameters and observation arc length on the observability of small perturbations .
finally , section [ disc ] provides a discussion of our results and section [ concl ] presents conclusions .
to begin , we must describe the working definition of the pioneer effect used in this paper .
following @xcite , we take the pioneer effect to be manifested by a radial acceleration , directed sunward , of magnitude @xmath1 m s@xmath2 .
since the primary purpose of this paper is to investigate the dynamical consequences of the pioneer effect , which apparently begin at about 20 au from the sun , and because there are no data indicating a more gradual onset of the pioneer effect , we will assume that the anomalous acceleration begins abruptly at a heliocentric distance of 20 au .
we recognize that this is a simplistic model of the pioneer effect .
alternative mechanisms exist that cause the acceleration to vary with object mass , orbital eccentricity , radial distance , and other parameters of the motion .
as further observations are made of the pioneer effect , they can be used to investigate different force models to explore various alternatives until the pioneer effect is either ruled out or its origin is found . however , the current status of knowledge of the pioneer effect argues that this simple model should be investigated first .
additionally , a perturbation beginning somewhat closer to the sun at a more gradual pace would be more easily seen in the motion of the planets . as shown recently by @xcite ,
if the pioneer effect occurred at shorter distances from the sun , its effect on planetary ephemerides would have already been detected .
thus , our assumptions about the pioneer effect represent the minimum plausible perturbation given the available data .
notwithstanding our concentration on the pioneer effect , the analysis presented here should be valid for any small constant perturbation to gravity in the outer solar system .
there are several ways to approach assessing whether small perturbations like the pioneer effect can be observed in the motions of the outer planets . for clarity , however , let us make clear that we are talking of perturbations to the motion of the outer planets that are spherically symmetric and directed towards the sun .
perturbations due to localized mass concentrations ( i.e. , planet x ) are specifically excluded .
however , the latter can be instructive .
some years ago , there were allegations that the motion of the outer planets contained anomalies that indicated the presence of a large mass concentration in the outer solar system .
@xcite showed how these presumed anomalies in uranus motion vanish when the orbital elements are adjusted while using correct values of neptune s mass as determined by spacecraft .
the important point here is that it is not enough to merely compare projected positions ; rather , one must adjust the orbital elements , and even other parameters defining the problem , to best fit observational data .
there are a number of approaches to addressing this problem .
one approach was taken by @xcite in which a variation in the value of the astronomical unit is investigated for the inner planets .
their analysis assumes a small eccentricity and makes use of ranging data available from radar and spacecraft observations .
this shows that quite stringent constraints are placed on the nature of gravity at scales approximating that of the inner solar system . in particular , any pioneer - like acceleration at those scales would long since have been detected .
another approach is to consider the relative magnitudes of the acceleration due to the sun and that due to the pioneer effect .
this is the approach taken in @xcite , although it was not documented in that paper .
certainly the pioneer effect should not be expected to have a dominating impact on the motion of the outer planets .
the ratio of the pioneer acceleration to that produced by the sun at a distance equal to the semimajor axis of the outer planets is 0.005 , 0.013 , and 0.023 percent for uranus , neptune , and pluto , respectively .
if we integrate an orbit without any perturbations to newtonian gravity other than the pioneer effect , and compare with the keplerian case , we find the orbital periods of these objects are systematically shorter .
uranus period shortens by 5.8 days and neptune s by 24.1 days , while pluto s period drops by 79.7 days .
these intervals correspond to 0.02 , 0.04 , and 0.09 percent of the periods of uranus , neptune , and pluto , respectively .
differentiating kepler s third law implicitly , we obtain a relationship between a small change in an orbiting object s period and a corresponding change in the orbit s semimajor axis .
in particular , @xmath3 where @xmath4 is the semimajor axis ( in au ) and @xmath5 is the orbital period ( in years ) . for uranus ,
the fractional period shortening due to the pioneer effect is equivalent to a fractional change in semimajor axis of approximately 0.013 percent .
similarly , the orbital period shortening is equivalent to a reduction in semimajor axis of 0.027 and 0.059 percent for neptune and pluto , respectively .
these simple calculations imply an equivalent change in aphelion distance of @xmath6 , @xmath7 , and @xmath8 km for uranus , neptune , and pluto . in the cases of uranus and neptune , this is less than the approximate maximum errors in range of @xmath9 km @xcite ) . in the case of pluto ,
the change in semimajor axis that would correspond to the shortening of its orbital period with the pioneer effect is about twice the radial distance uncertainty .
however , pluto has completed less than one - third of an orbit since its discovery and its orbital elements are even less well - determined than the other outer planets .
thus , small changes in other orbital elements could easily obscure any orbital changes due to the pioneer effect .      the values of orbital elements
are not directly observable and to reduce observational data to orbital elements it is necessary to proceed by using numerical approximations . the normal manner in which one proceeds is to first determine a `` preliminary orbit '' from a small number of observations and then to refine it by successive approximations into a `` definitive orbit '' as many more observations are added .
the method used to improve the orbit and obtain the definitive orbit through additional observations is called `` differential correction '' and is well described in standard celestial mechanics texts ( e.g. , @xcite , @xcite , or @xcite ) .
differential correction uses a least squares approach to iteratively refine the estimates of the elements as more observations become available .
additionally , statistical information on the errors of the elements naturally results from the differential correction process .
the process results in a set of orbital elements , along with error estimates for the elements and covariance parameters showing the degree to which the elements are correlated .
it is worth remembering that the equations governing celestial mechanics are nonlinear and normally a linearized version of the problem is used to determine the covariances and elements .
however , the nonlinearity manifests itself in occasional difficulties in determining elements and some of these difficulties will be discussed later in section [ sect3-b ] .
the whole chain of analysis outlined above for orbital element determination is associated with a number of errors that must be understood and characterized .
some are associated with the description of the solar system dynamics and calculation of the orbits , while some originate in the observations themselves .
however , in determining elements from observations , there is no prior warning about the pathological situations that can arise . for example
, if the inclination is small , the error in the longitude of the ascending node will be large .
similarly , if the eccentricity is small , the error in the argument of perihelion will be large and the time of perihelion passage or mean anomaly will be poorly defined .
such problems can be avoided through the use of alternative orbital elements ( e.g. , equinoctial elements ) .
thus , the orbital characterization process outlined above must be carried out to enable predictions of the future positions of objects on the sky .
the errors associated with the predictions must be compared statistically to determine whether a truly observable difference can be asserted between the alternative models of gravity .
importantly , since orbital elements are not directly observable , it is not enough to simply integrate an object s equations of motion forward in time from a set of initial conditions , with and without the perturbation ( as was done in @xcite ) , and then compare the differences .
one must first determine element values by adjusting them to a set of observations and a gravity model to obtain a new set of elements with associated error metrics .
this can be understood physically by noting that the orbital elements describe the conserved mechanical energy of the body and if the potential field changes due to the existence or absence of a perturbation the elements need to be redetermined .    in order to address these issues
, we use the orbfit software package @xcite .
this program uses observational data and data on the dynamics of the solar system to determine orbital elements and predict ephemerides for minor bodies .
since orbfit is available in source code , we added a simple option to include forces arising from the pioneer effect . by comparing orbits resulting from synthetic observations with and without the pioneer effect ,
we explore the expected effects of such a perturbation on asteroid orbits and examine when this effect can be detected astrometrically .
orbfit is a complex piece of software that performs a number of tasks that relate to this analysis .
the first task of interest here is that of orbit determination .
this process begins with astrometric observations of the body of interest .
an initial orbit is determined from a few of the available observations and a process of differential correction is applied to correct the initial elements and take into account any additional observations .    the second task of interest performed by orbfit is orbital position prediction .
to predict the position of bodies of interest , orbfit takes the position and velocity of the body of interest , and integrates the equations of motion of that body with many additional perturbations due to the other bodies in the solar system , solar oblateness , etc .
additionally , the process used by orbfit naturally provides measures of the errors in elements that result from the unavoidable errors in position measurement and the associated element covariances .
these data can be mapped on to positions in the sky and the error associated with the prediction , which permit our analysis to take place .    by using orbfit for our calculations
, we are able to include planetary and other perturbations and factors that affect the determination and prediction of orbits . these factors are dealt with in orbfit to create orbital predictions and also conduct orbital analyses in a manner that very closely mimics the process by which data from observations are reduced to orbital predictions and elements in the real world .
in particular , orbfit is able to deal with various errors and noise that are both unavoidable and necessary for calculations at a level of accuracy sufficient for distinguishing between orbits affected by the pioneer effect and those subject to gravity without the additional perturbation .
the calculations outlined above represent standard techniques of celestial mechanics and can be performed by any number of systems .
orbfit was chosen because of its availability in source code .
this permitted addition of a simple option to include the force arising from the pioneer effect , which were simply added to the forces exerted on the bodies of interest by the major planets .
rather than integrating the orbits of the main bodies in the solar system , their dynamics are introduced via the jpl de405 ephemeris . in this context
the planetary ephemeris is a sort of lookup table containing positions of the planets so that their effect on the motion of a body of interest might be determined .      in this paper
, we will parametrically vary the eccentricity of hypothetical test bodies , generate synthetic observations for them with and without a pioneer - like acceleration , and determine the conditions under which the pioneer effect can be observed . with the exception of eccentricity , the orbital elements of our test bodies are those of pluto @xcite .
the eccentricity of the bodies are 0.001 , 0.005 , 0.01 , 0.05 , 0.1 , 0.2 , and 0.3 .
these values are chosen to bracket the actual eccentricity of pluto while extending to nearly circular orbits similar to those found in the other outer planets .
the nominal elements at the epoch j2000 ( jd 2,451,545.0 ) , are referred to the mean ecliptic and equinox of jd 2000.0 and are given in table [ elements ] .
ll semimajor axis & 39.48168677 au + eccentricity & see text + inclination & 17.14175 degrees + longitude of ascending node & 110.30347 degrees + argument of pericenter & 113.76329 degrees + mean anomaly & 14.86205 degrees +    the observation cadence we have chosen to use roughly reflects that which exists for pluto .
if one reviews the observations for pluto that were used to construct the de405 ephemeris , we find approximate observation frequencies as indicated in table [ obsfreqs ] . with the exception of the 1914 - 1929 period
, we have adopted these cadences and generate synthetic observations accordingly .
the pre - discovery observations occuring from 1914 - 1929 are aggregated and an overall rate of one observation per year is used for this time period .
lrrl 1914 - 1919 & 11 & 6 & 2 per year + 1920 - 1929 & 6 & 10 & 1 per 2 years + 1930 - 1939 & 431 & 10 & 1 per week + 1940 - 1949 & 233 & 10 & 2 per month + 1950 - 1959 & 113 & 10 & 1 per month + 1960 - 1969 & 113 & 10 & 1 per month + 1970 - 1979 & 364 & 10 & 3 per month + 1980 - 1989 & 361 & 10 & 3 per month + 1990 - 1999 & 1,125 & 10 & 1 per 3 days + 2000 - 2006 & 962 & 7 & 1 per 3 days + 2007 - 2011 & est .
& 5 & 1 per 3 days + 2012 - 2163 & est .
& 152 & 1 per day +    in the analysis , we use observation arc lengths running up to 250 years that are evaluated at 50 year intervals . since the orbital period of these objects is approximately 250 years
, our synthetic observations span a complete orbit .
the total number of synthetic observations occurring over the observation arc lengths are shown in table [ arclength ] .
lrr 1914 - 1963 & 952 & 952 + 1964 - 2013 & 4,212 & 5,164 + 2014 - 2063 & 18,262 & 23,426 + 2064 - 2113 & 18,262 & 41,688 + 2114 - 2163 & 18,262 & 59,950 +    the general approach is to use orbfit with a set of elements defined as described in table [ elements ] to generate ephemerides separately in both the perturbed ( that is , with the pioneer perturbation present ) , and unperturbed ( without the pioneer effect ) cases .
the predicted positions of the test body on the sky represent `` perfect '' observations with no uncertainties in either observations or elements . in each case
, we can take these predicted ephemeris positions and add gaussian observational error .
these randomly altered positions then represent the results of synthetic astrometric observations .
the random observational error applied is equal to 0.3 seconds of arc , representing the results of good quality ccd astrometry reduced against modern star catalogs .
this value of error is assumed to be isotropic on the sky ; thus , we generate and apply synthetic errors in right ascension and declination equal to this value divided by the square root of two with appropriate adjustment to the right ascension for the cosine of the declination .
the resulting different sets of synthetic observations and alternative gravity models can then be used with orbfit to determine elements and errors associated with the elements as a function of the eccentricity and the observation arc length .
it is noted that this is the procedure used to investigate the motion of minor planets in @xcite and @xcite .
we recognize that this assumed error represents a level of accuracy much better than that found in observations taken until relatively recently . indeed
, there are many potential sources of position error that might be modeled , some of a known nature and some unknown .
the data were taken at different observatories and different techniques may have been used to reduce the observations , leading to systematic errors in the derived positions .
the positions have been reduced with various catalogs and transformed , sometimes repeatedly , to different reference frames . even for modern observations
, there may be errors in the adopted precession values for the j2000 epoch .
additional error sources include errors due to uncertainties of the initial conditions of the other planets ( e.g. , ephemeris errors ) , errors due to corrections to earth rotation , and unmodeled instrumental corrections .
since the primary impact of increasing the a priori positional error is to de - emphasize those observations with larger errors , we chose to equally weight all our synthetic observations . thus , we assume much more optimistic error estimates for early observations than are representative in the actual observation archive and results in equal weighting for all the synthetic observations .
this should result in our calculations providing an optimistic estimate of the detectability of the pioneer perturbation .
an additional reason for electing to model a uniform observational accuracy across the span of the synthetic observations in order to minimize variation in our results due to another source of observational noise . for parts of the synthetic observational arc
prior to the advent of modern catalogs and ccds , this error will be considerably less than that found in the actual observational archive .
interestingly , the orbfit program is capable of varying a priori observational error across time or across different observatories , but we did not make use of this capability in this paper .    given these assumptions , we have four cases , which we illustrate in table [ cases ] .
the first pair of cases ( the `` gravity with pe '' column in table [ cases ] ) , and the focus of our analysis , is for observations that are generated with a pioneer - effect perturbation present meaning that the pioneer effect exists in nature and determines an object s motion .
these same observations are analyzed with two gravity models .
the first model is that of newtonian gravity with the addition of a constant radial acceleration with a value commensurate with that associated with the pioneer effect ; the second is standard newtonian gravity .
the comparison of the orbital solutions for these two cases correspond to a situation in which the universe is one where the pioneer effect actually exists and we analyze it to determine which of our two gravity models is correct .
this is the situation portrayed in the first column of table [ cases ] .
the second pair of cases are comparable except that they are associated with a universe where the pioneer effect does not exist .
again , we compare two theories of gravity to see if we can observationally distinguish between the two . this situation is portrayed in the second column of table [ cases ] .
this paper investigates the first case only ; as shown in @xcite , the second case is nearly symmetric with the first , and produces similar conclusions .
lcc gravity with pe & `` matched '' & `` mismatched '' + only gravity & `` mismatched '' & `` matched '' +
as outlined above , our plan is to take ephemerides created with various known orbital elements and predict positions on the sky when the motion is governed by newtonian gravity with and without a pioneer - like perturbing acceleration .
this process involves integrating the equations of motion of a body subjected to two different force laws .
the results of this calculation provide sky positions as a function of time .
since the angular differences between the predicted positions are small , we separately consider the difference in right ascension and declination resulting from the two situations .
[ fig1 ] shows the resulting position difference .
since the epoch of the elements is jd 2,451,545.0 , the angular differences are zero at that date and diverge as one moves forward and backwards in time from the element epoch .
it is worth emphasizing that the two ephemerides being compared are generated with two separate gravity models with identical and exactly specified elements .
figure [ fig1 ] can be compared with figs . 9 and 10 in @xcite .
they drew the conclusion that the divergence in position , since not observed over the past 90 years , indicates the pioneer effect does not exist .
however , such a conclusion is based on the idea that the elements ( i.e. , the initial conditions ) are known precisely .
the reality is that these position predictions are those related to the elements derived from the assumption of newtonian gravity coupled with observed positions .
if there is a `` mismatch '' between the gravity model used to derive orbital elements and the reality that determines the actual motion of the orbital bodies , these will not be the correct elements . even if the elements are derived from a `` matching '' reality and gravity model
, it is remarkable how rapidly orbital predictions can degrade outside the observation arc , especially when that arc is short relative to the orbital period .
these aspects will be discussed in detail in the next subsection .
if we take the ephemeris positions associated with one `` reality '' ( e.g. , the one containing the pioneer perturbation  see table [ cases ] ) consonant with the observation cadence shown in table [ obsfreqs ] , and add an isotropic , normally distributed random position error of 0.3 arc seconds to each of the positions , the result is a set of synthetic positions that are illustrative and representative of those that might have been obtained as observations of the bodies .
this set of observations can be divided into five different arcs covering 250 years ( approximately one orbital period ) ; thus , each arc covers an additional 50 years from the start of the observations and represents an incremental one - fifth of an orbital revolution .
each synthetic observation arc can be analyzed when a pioneer - like perturbation is included and when it is not , corresponding to the two rows in the first column of table [ cases ] .
the observation fitting process serves to determine the elements of the orbit and these are the elements that can then be used to predict sky position . how do the observed positions differ when this procedure is followed ?
[ fig2 ] through [ fig6 ] shows the difference in sky position as a function of time for the orbital fits associated with the two gravity models . in all cases ,
the left hand graph shows the right ascension difference while the right hand graph shows the declination difference . comparing these figures with fig .
1 above and figs .
9 and 10 in @xcite shows that for shorter arcs ( those ranging up to 100 years long ) the difference in predicted motion for the two gravity models is quite small for the duration of the observation arc and only begins to substantially diverge when the position is being extrapolated to times beyond those for which observations exist . for a 150 year arc , small irregularities in the position difference begin to be noticeable during the observational arc .
as the arcs grow still longer ( to 200 years and beyond ) , more substantial irregularities are seen . also , it is generally true that positional differences are greater the greater the eccentricity , all other things being equal .
recall what the differences plotted in figs .
[ fig2 ] through [ fig6 ] represent .
synthetic observations for a universe with a pioneer effect were generated and orbital elements were determined by fitting the same observations to two alternative gravity models .
one of these models is correctly `` matched '' to the observations ; that is , the gravity model used for orbit fitting is the same one used to generate the observations .
the other model does not correspond to the one used to generate the observations .
thus , we would a priori expect that the former would fit the data better than the latter .    what we see
, however , is surprising . for observation arcs
as long as a century , which at best is the situation obtaining for pluto at the current time , there are only relatively small differences in the positions predicted by the correct gravity model and those predicted by the mismatched model . in all cases , however , the predicted positions outside the range of observations rapidly diverge from one another .
this behavior is a manifestation of two inter - related factors .
first , the orbital fitting problem is inherently nonlinear and is normally solved in the linear approximation .
even if not mathematically chaotic , the system of equations is certainly sensitively dependent upon initial conditions ; small changes in elements can result in large changes in predicted position outside the range of observations .
secondly , this sensitivity is exacerbated by the problem of a short observation arc .
the length of the entire observational archive for pluto is no more than about one - third of its orbital period . even in as simple a case as linear least squares curve fitting ,
a limited amount of independent variable data ( corresponding to a short observation arc ) will lead to relatively large errors in the fitting parameters ( corresponding to the orbital elements ) .
together , these factors conspire to potentially generate large position errors outside the observation arc , while increasing the length of the observation arc can markedly reduce error over the whole of the arc and even beyond it .
some of these issues can sometimes be alleviated by a change in variables .
use of other than keplerian orbital elements ( e.g. , equinoctial elements ) can allow one to avoid some problems with determining an orbital solution , in particular in cases with small eccentricities or small inclinations .
similarly , a linear combination of keplerian orbital elements can be used to improve the accuracy of some elements ( but see the footnotes on pages and ) . however , those types of variable change do nothing for the issues associated with the ill - conditioned nature of the linear matrix and a short observational arc .    in any case
we see that we must fit observations to particular gravity models and adjust orbital elements before predicted positions on the sky can be compared . noting that deriving orbital elements and predicting an orbit beyond a short observation arc provides little information about the actual motion that will be seen once new observations are made , we conclude that drawing conclusions from such extrapolations is totally unwarranted .
thus , we can conclude nothing regarding the pioneer effect and the motion of pluto as it is currently known , at least subject to the limitations of our methodology .
now let us examine more closely the relatively small position differences at the beginning of the observation arcs in figs .
[ fig2 ] through [ fig6 ] .
[ fig12 ] is similar to figs .
[ fig2 ] through [ fig6 ] , but shows an expanded view of the difference in position of a body with an eccentricity of 0.3 ( e.g. , similar to pluto ) and a 100 year observation arc . at the beginning of the arc , the difference in right ascension for the two gravity models is nearly six seconds of arc and the difference in declination is about 1.5 arcsec .
one would certainly think that this level of difference would be observable ; however , we argue below that it is not a distinguishable difference given our current knowledge of pluto s orbit .
figure [ fig13 ] shows a subset of the residuals in both right ascension and declination of actual observations with respect to a recent jpl ephemeris , de414 .
the full set of residuals show a number of clear outliers and we have removed from the set any residuals of magnitude
greater than ten seconds of arc .
prior to about 1960 there is clearly a greater spread in the residuals than at later times .
additionally , prediscovery images taken prior to 1930 and dating to as early as 1914 have been found in the archive .
these observations seem to contain serious irregularities manifested by large and systematically biased residuals . supporting this assertion , @xcite find that systematic differences between the observed and calculated orbit are a continuing issue . additionally , they find that only about half of pluto s motion since its discovery has been observed in a systematic and organized manner .
they also comment that it is impossible to extrapolate pluto s position more than a few years into the future .
this is indicative of a poorly characterized orbit .
the residuals of the synthetic observations with respect to a `` mismatched '' gravity model are similar in character and magnitude to the differences between predictions for different gravity models that are shown in fig .
7 . as observed previously , there seems to be a significant difference between the two sets of predictions .
a standard comparison of these models @xcite would calculate and compare @xmath10 values for the two models .
if the resulting value is about unity , the fit is good .
the overall issues associated with model goodness - of - fit are complex and will be discussed later in section [ quality ] ; however , anticipating this complexity let us look more deeply at the data .
if we compare our `` mismatched '' orbit fit with actual observational data and their associated fit , we will see that there are comparable systematic errors in both the synthetic and the real datasets .
systematic trends in residuals can be visually discerned more readily by examination of normal points rather than individual data points . in a normal point plot , means and standard deviations of the residuals are collected and plotted for each planetary opposition . fig .
[ fig14 ] shows the normal points for the de414 residuals in the `` high - low - close '' portion of the plot , where the `` closing '' value indicated by the small horizontal mark is the average residual for each opposition and the length of the vertical line indicates a one standard deviation variation above and below the mean residuals .
the comments above about the systematic bias of the prediscovery observations and the larger dispersion of residuals before 1960 are borne out in the normal plot .
now , let us compare this data with normal points for the synthetic observations of pluto . in developing the observation cadence for the synthetic observations , we did not consider the relative positions of the earth , the sun , and pluto and thus did not concentrate synthetic observations around particular yearly oppositions . since our primary purpose was to illustrate broader issues with the data , this
was done to avoid introducing additional noise through a varying observational cadence .
thus , for the post - discovery residuals associated with our synthetic observations we found the average and standard deviations for each calendar year . for prediscovery images , our assumed observation cadence of one observation per year clearly made the dispersion calculation impossible .
in these cases , we accumulated data by decade ; thus , we calculated the mean and standard deviation for the 1914 - 1919 period and the 1920 - 1929 period and plotted the resulting values accordingly .
the normal points associated with the synthetic observations are also shown in fig .
[ fig14 ] .
the average is shown by the solid line and the plus and minus one standard deviation limits are shown as dotted lines .
several comments are in order about these graphs .
first , the width of the standard deviation spread on the synthetic observations is different from that of the de414 residuals .
this is due to our assumption of a constant 0.3 arc second astrometric accuracy of the synthetic observations .
the actual observations that contribute to the dispersion of the de414 residuals have varying accuracies , ranging from large values for the early observations to small values for current ccd observations . as pointed out earlier , the level of astrometric accuracy used in the synthetic observations is much better than that found in the actual observation archive until relatively recently .
we chose to model a uniform observational accuracy across the span of our synthetic observation arc in order to make an optimistic estimate of the detectability of the pioneer perturbation as well as to minimize variation in our results due to another source of noise ( e.g. , the observation error ) .
thus , rather than the uniform difference between the upper and lower limits , the normal points of the actual de414 residuals vary in width . if we had modeled a varying error , such a variation in width would be expected for the synthetic observations as well .
another point to be made about fig .
[ fig14 ] is that while the prediscovery residuals for the synthetic observations and those with respect to de414 are of opposite signs , their magnitudes are of the same order . in the former case ,
the large residuals are due to fitting a `` mismatched '' gravity model consisting of only nominal gravity to observations that result from gravity plus an anomalous pioneer acceleration ; in the latter case , unknown observational errors have resulted in large and biased residuals .
what are we to make of the difference in sign ?
we would argue that since there are errors of an unknown nature present in these early observations , we have absolutely no basis is even assuming the sign of the residual is indicative of reality .
although the prediscovery residuals have characteristics that indicate that they might represent a constant offset in right ascension , an effect that can be caused by many errors in the analysis of the observations , the limits on the residual s variation can also accomodate a linear trend as is seen in the corresponding synthetic observations .
we can only observe that the approximate magnitude of the residuals at early times are roughly the same for both real observations relative to the de414 ephemeris and the synthetic observations with respect to their ephemeris .
thus , we can not draw any conclusions about the existence or nonexistence of the pioneer effect from our current knowledge of the orbit of pluto ; its observation arc is too short for such a determination .    a third point with respect to fig .
[ fig14 ] is that post - discovery residuals are similar for both the synthetic observations and those of de414 and very close for post-1960 observations .
overall , the only place where the two sets of normal points diverge is in the prediscovery era .
thus , neither set seems to indicate the existence or nonexistence of an additional gravitational anomaly .
both sets of residuals possess patterned irregularities , but it is impossible to say that one is superior to the other .    as a final comment on the comparison of the two sets of residuals , fig .
[ fig15 ] shows the total rms residual by epoch for the two cases .
this quantity is the square root of the sum of the squares of the residuals in right ascension and declination with an appropriate adjustment for the cosine of the declination .
as above , de414 s residuals are accumulated by opposition , while the residuals of the synthetic observations relative to their ephemeris are accumulated by decade for prediscovery observations and by year thereafter .
the quality of the residuals in the two cases is quite high over most of the observation arc . even in the prediscovery era ,
the trends of the residuals are similar .
as before , there is a linear trend in the synthetic residuals , but the entire prediscovery regime only includes two points . the trend of those two points is not dissimilar to that of the de414 residuals . even if the prediscovery observations possessed a systematic error of as large as two seconds of arc , the trend in rms residual is not substantially different from those of the synthetic observations .    from the discovery of pluto to about 1960 ,
the de414 total rms residuals are greater than those found for the synthetic observations .
this difference is expected due to the conservative assumptions of astrometric accuracy that were made in generating the synthetic observations . however , the trend of the residuals in both cases is similar , indicating a comparable fit to the observations . from about 1960 onwards ,
both sets of residuals fit quite well .    given the likely , but uncharacterized , errors in pre - discovery observations , we are led to the observation that real observations are explained about as well by de414 as the synthetic observations are explained by their ephemeris .    the context of this comment must be kept in mind .
the de414 residuals relative to its ephemeris represent a fit to the actual observation archive by a normal gravity model without a pioneer anomaly .
the synthetic observations are created with a pioneer anomaly present and thus their underlying motion is not in keeping with normal gravity .
when the synthetic observations are fit to a gravity model without a pioneer anomaly included in the dynamics , the fit is almost indistinguishable from that of the de414 ephemeris . especially given the uncharacterized errors apparently present in the prediscovery observation archive , we can not distinguish between the two gravity models on the basis of the current observational record .
however , we note that this is at least as much due to the shortness of the observation arc for pluto as due to any difference in gravity models
.    what adds special interest to the pluto case , however , is that due to its relatively large eccentricity pluto shows potentially detectable differences between the two gravity models if the entire one hundred year observation arc were uniformly reliable .
thus , we are led again to the desirability of continuing to gather astrometric data on pluto and expect that an astrometric difference in position might soon be able to falsify the existence of the pioneer anomaly . nevertheless , for the length of the currently available observation arc and
because of the uneven quality of the observational archive , we can not currently distinguish between the two gravity models on the basis of orbital position
.      there are six orbital elements corresponding to the six degrees of freedom of the classical dynamical problem in three dimensions . given that elements are determined from observations by a model fitting process
, we must ask whether this large number of degrees of freedom is somehow redundant and whether a simpler model of the dynamics can be used to illuminate the problem .
this is the approach taken by @xcite , where synthetic observations are generated for a simplified dynamical model with four degrees of freedom .
conclusions are then drawn about the degree to which a pioneer effect perturbation can be fit to this simplified model and thus whether the pioneer effect exists .
the problem with this approach is that the elements that are ignored or suppressed in @xcite are precisely those with the greatest uncertainty resulting from the orbital fitting process .
further , as will be discussed below , care must be taken when using rms residuals as a primary measure of goodness of fit as far as model selection is concerned .
the underlying cause of the problems outlined above was discussed earlier .
briefly , this behavior is due to two inter - related factors : first , the orbital fitting problem is inherently nonlinear and is normally solved in the linear approximation ; second , the sensitivity of the orbital solution to initial conditions is exacerbated , especially in the case of pluto , by the problem of a short observation arc .
these issues can sometimes be alleviated by a change in variables .
use of other than keplerian orbital elements can result in an improved determination of orbital elements @xcite . however , this type of variable change does nothing for the issues associated with a short observational arc .
the uncertainties arising from these sources make the approach in @xcite problematic . in simplifying the problem by assuming away two or more degrees of freedom , the full parameter space , which has ample room to conceal the differences in predicted position , is not available for that purpose .
this provides the false impression that perturbation effects are observable when , in fact , they are not
.    one might think that the pioneer acceleration , assumed to be radial , would only effect the motion of an orbiting body in its orbital plane . in a theoretical sense , this is true .
the pioneer effect exerts no torque on the orbiting body , and the plane of the orbit remains unchanged as time progresses .
an analytical investigation of such an orbit can be considered as two dimensional and the number of degrees of freedom of the motion can be reduced accordingly .
this is the essence of the approach of @xcite .
however , in an observational context the motion is not restricted to a plane .
there are at least two factors that force orbital motion to be considered to occur in three dimensional space with the associated degrees of freedom .
the first is that unless the observer s location ( e.g. , earth ) remains in the plane of the orbit of the body , there will be parallax introduced which forces the motion into three dimensions .
similarly , observational errors will unavoidably make the object s position vary from its theoretical planar motion . since both of these are a priori unknown and are only determined in the context of a least squares orbital solution for the observed object , the object s motion must necessarily be considered in its complete dimensionality .
of course , this requires all six degrees of freedom of the classical dynamical problem
. there can be large errors associated with orbital elements and the extra degrees of freedom can provide ample room to `` hide '' orbital variations .
this will be discussed in detail below .
as shown in any numerical analysis text ( for example , @xcite ) , an upper limit on the relative error of a fitted model parameter is related to the relative error of the independent variables by a condition number .
the condition number can be estimated by the square root of the ratio of the largest to the smallest eigenvalues of the associated normal system of equations @xcitethe differential correction problem can be stated as a linear matrix equation @xmath11 , where @xmath12 is an appropriate jacobian matrix , @xmath13 is a vector of elements , and @xmath14 is a vector of observations .
the normal equations are formulated by multiplying this equation by the transpose of @xmath12 ( denoted by @xmath15 ) , giving @xmath16 .
this system of equations is solved formally by inverting the @xmath17 product and multiplying again , giving @xmath18 . in general , the @xmath12 matrix is not square .
it has a column for each orbital element ( e.g. , six columns ) and a row for each measured sky coordinate ( e.g. , twice the number of observations ) .
however , the @xmath17 matrix is square and possesses a corresponding set of eigenvalues .
the ratio of the largest to the smallest eigenvalues of the @xmath17 matrix provides the condition number of the normal equations . however , the condition number of the normal equations is the square of the condition number of the @xmath12 matrix .
thus , the original equation @xmath11 has a condition number equal to the square root of the condition number of the normal equations . ] .
it might be thought that the ill - posed nature of the orbital fitting problem could be obviated through a coordinate transformation .
indeed , there are alternative orbital elements that alleviate some difficulties with orbits of too small eccentricity or inclination , for example , equinoctial elements @xcite .
alternatively , new variables that are linear combinations of the elements considered here can be used to alter the mix of errors among them .
thus , rather than specifying the time of perihelion passage , we can specify some instant of time and then specify either the mean anomaly , the mean longitude ( the sum of the longitude of the ascending node , the argument of perihelion , and the mean anomaly ) , or the true longitude ( the sum of the longitude of the ascending node , the argument of perihelion , and the true anomaly ) at that moment .
however , it can be shown that the condition number of a problem like this is invariant under linear coordinate transformations @xcitematrices @xmath12 and @xmath19 are said to be similar when @xmath20 , where @xmath21 is nonsingular .
such transformations arise from changes of variables . now , suppose we have an eigenvalue problem , @xmath22 and let us perform this similarity transformation
. then @xmath22 implies @xmath23 , which gives @xmath24 which can be stated as @xmath25 so that @xmath12 and @xmath19 have the same eigenvalues and the eigenvectors of @xmath12 and @xmath19 are related by @xmath26 .
thus , similarity transformations preserve eigenvalues but not eigenvectors .
since the condition number of a problem is the ratio of the largest to smallest eigenvalues , we see that condition number is preserved under linear coordinate transformations . ]
although nonlinear coordinate transformations might alleviate this problem , there are none known to the authors that are used in practical astrometric work using orbfit or other standard software packages .    fig .
[ fig7 ] shows the condition number of our orbital fitting problem as a function of orbital eccentricity and observation arc length . for low eccentricities , the condition number can become quite large . even for eccentricities
similar to that of pluto , the condition number for a short arc problem implies large errors in orbital elements .
ordinarily , the condition number is important in the context of numerical roundoff with finite precision arithmetic in a computer .
this is not the situation here ; however , our precision constraint is the accuracy of the position observations that drive the orbital fitting process .
thus an angular precision of 0.3 seconds of arc in right ascension and declination represents a relative error on the order of @xmath27 percent or less . with a condition number of 100 ,
this translates into an error of about 30 seconds of arc in orbital parameters like mean anomaly and argument of pericenter .
the corresponding results when we fit our synthetic observations bear out these theoretical considerations .
[ fig8 ] shows the relative error ( one sigma estimated error in the linear approximation ) in orbital elements obtained with different observation arc lengths .
the four panels show a representative set of orbital eccentricities .
what is generally true for these figures is that for short arcs , the relative error in the elements is much greater than it is for the longer arcs .
similarly , the uncertainty in the location of the perihelion ( as measured by the argument of the pericenter and the mean anomaly ) is quite high and remains so for relatively long observation arcs .
additionally , since these errors are derived in the linear approximation , the real errors are likely to be greater , perhaps much greater .
+    this implies a large error in mean anomaly or argument of pericenter , both associated with the spatial and temporal location of perihelion .
this , coupled with the observational nature of the astrometric problem necessitating a full three dimensional treatment , forces us to conclude that a substantial amount of error can be absorbed into a multidimensional parameter space and , since it is unclear how the error can be allocated across the parameter space dimensions , the full parameter space must be included to properly reflect the uncertainties of the motion .
thus in order to drive the error in orbital elements down to levels allowing us to determine whether the pioneer effect exists , we need a full dynamical model in three - dimensional space and either a relatively high eccentricity object or a long observation arc or both .    as observed before , our approach is itself an approximation . a `` full dynamical model '' should include not only adjustment of the orbital elements of pluto , but also simultaneous adjustment to those of uranus and neptune as well .
it is only in this way that all the second order perturbations to the system can be taken into account .
however , as observed above , in order to illustrate the ideas concerned and the weaknesses of the approaches outlined above , we restrict our attention here to manipulating the orbit of pluto . this approach will be justified further in section [ disc ] .
the problem of assessing model validity and comparing alternative models is huge and any kind of complete treatment is well beyond the scope of this paper .
however , given the comments above , it is appropriate to offer some observations on this large and interesting topic .
the issue of goodness - of - fit for models is not a simple one .
often , a @xmath28 is calculated and compared to determine which of competing models is the preferred . however , there are complexities involved .
for example , @xcite discusses goodness - of - fit and points out that while use of the @xmath10 statistic as a goodness - of - fit measure for binned data is justifiable and is often done , it has flaws .
he further observes that for unbinned data ( such as we have here ) , there is no equivalent popular method for measuring goodness - of - fit .
indeed , @xcite gives several examples of problematic goodness - of - fit cases .
similarly , @xcite advocates investigating the structure of residuals to find patterns , biases , and sytematic differences between the model and the data
.    we can divide metrics for model quality into two broad categories that , although useful , are certainly not mutually exclusive .
first , we can consider what might be termed point estimates of fit quality .
often , a measure of merit is used that is related to the rms residuals between the model and the data .
[ fig9 ] compares our rms residuals for various eccentricities .
one panel of this figure shows the results of fitting the `` mismatched '' case of the pioneer - perturbed observations being fit to a strictly newtonian gravity model ; the other panel shows the corresponding `` matched '' case where the gravity model includes a pioneer perturbation . in both cases
the underlying synthetic observations contain a pioneer perturbation .
the most striking feature of this figure is how slowly the residual degrades as the arc length increases , especially for low eccentricity objects . indeed , until the observation arc is over 100 years long , there are only small differences in the quality of the fit as measured by rms residuals for any of the eccentricities evaluated . even in the `` mismatched '' case ,
the rms residual is still less than half an arc second when the arc is 150 years long . in keeping with the discussion above about the ill - conditioned nature of the orbital fitting problem
, we see that strictly relying on rms residuals as a measure of merit for model selection can be problematic .
another measure of merit that is sometimes invoked is the uncertainties in the orbital elements themselves .
the orbital fitting process automatically provides mathematically well - justified estimates of the elements , although these are in a linear approximation that has its own difficulties since the orbital fitting problem is intrinsically nonlinear .
this approximation provides error estimates similar to those shown in fig .
[ fig8 ] . across the eccentricities shown ,
there is substantial error spread throughout the parameter space . in particular
, it might be instructive to compare panels a and d in fig .
[ fig8 ] . in panel a , a very low eccentricity case , the errors in semimajor axis , inclination , and longitude of the ascending node drop to a level allowing seven or eight significant digits to be present in the corresponding element values . at the same time
, however , the uncertainties in eccentricity , argument of the pericenter , and mean anomaly remain relatively large . in the case with
the largest eccentricity evaluated , panel d shows the uncertainty to be more uniformly spread across elements .
trading off the uncertainties of one element against those of another to compare two models is therefore at best arbitrary and could be misleading .
position error in the sky plane is yet another point metric for fit quality . in this case ,
uncertainties in the orbital elements are mapped directly onto the sky . here ,
all the comments above about the uncertainty of the elements are valid , plus the observation that the mapping from the six - dimensional space of orbital elements to the two - dimensional sky plane is highly degenerate ; many sets of elements can map to the same region of sky .
thus , the sky plane position , by itself , can have problems as a point measure of model fit quality .
in particular , comparing the nominal solution for two models can result in a large error volume about the calculated positions of objects .    all three approaches to providing point estimates of model fit quality have one common weakness .
they are all based on a linearized form of the orbital fitting problem and make the assumptions associated with least squares fitting like normally distributed errors , independence , and no bias or systematic errors .
possibly the biggest manifestation of this is that we should not expect to be able to accurately extrapolate very far beyond the available observation arc unless there is at least a full orbital revolution s volume of data and even then , extrapolation is dangerous .
a possibly better approach to comparatively assessing model alternatives is to use some type of global comparison .
of broadest applicability in this context is the idea of testing for the normality of the residuals between the model and the data .
the overall basis for these tests is to look at the structure of the residuals and to perform statistical tests for lack of fit ( see , for example , ( * ? ? ?
* section 4.4.4.6 ) ) .
the basis for these tests is to search for patterns , biases , and systematic differences between the model and the data .
another approach that could be used to address the global quality of fit is to use bootstrap techniques .
this approach could involve , as done here , synthesizing observations whose error characteristics reflect those found actual observations and then using the synthetic observations to assess actual variations in the fit quantities .
another approach of this sort would be to remove some of the data from the fit and see how the fitted solution extrapolates to the times that were removed .
a similar approach would be simply to remove data randomly from the observation set and analyze the resulting variations in fitting parameters .
one form that this approach could take would be to conduct a full bayesian analysis of competing models .
the distributions of orbital elements resulting from synthetic observations that arise from a monte carlo process could be used to compare the probability of obtaining the observations in light of the competing models ( see , for example , @xcite ) .
although beyond the scope of this paper , the findings discussed here show that there are not currently enough data to warrant this type of analysis ; however , the advent of pan - starrs and lsst will change that situation in the near future .
in the analysis described above we first showed that one can not simply take orbital elements resulting from the fit of observations to a particular force model and use them to predict positions resulting from motion under the influence of another force . rather , we must refit the orbits to the observations under the new force model . in the case of pluto , this produces a well - fitting orbit that is indistinguishable in a practical sense , at least as long as the observation arc is short enough , from the unperturbed motion .
thus , without redetermining the elements we can not make sweeping generalizations about whether or not the outer planets orbits show that the pioneer effect does or does not exist .
similarly , although making simplifying assumptions about a physical situation in order to draw conclusions is a time - honored theoretical mode of attack , if the physical model is oversimplified we can be misled into erroneous conclusions . as seen above , in an observational context ,
both observer position and observational errors lead to the necessity of introducing the third spatial dimension with its associated degrees of freedom .
if we do not keep an appropriate number of degrees of freedom , the problem can be oversimplified too much and mislead us into unwarranted conclusions . in particular , we must take care in using such a simplified model to conclude that the pioneer effect does not exist .    once again , it should be emphasized that our approach is itself an approximation . to properly conduct an analysis of the sort outlined here , the orbital parameters of the entire system of outer planets should be included
. this would bring into the calculation all second order perturbations .
however , we have taken a simplified approach and believe it accurately and fairly addresses the uncertainties in the orbit of pluto .
one reason we believe our simplified approach is valid is illustrated in fig . [ fig16 ] .
[ fig16 ] shows the force per unit mass ( in units of the pioneer acceleration ) exerted by uranus and neptune at pluto s position over the period of time we have considered .
perhaps surprisingly , due to the relative positions of the planets in their orbits , neptune exerts less force on pluto than uranus over most of this period .
the magnitude of the force is of the order of the pioneer acceleration for both planets .
if the pioneer effect exists , we would expect the orbital elements of uranus and neptune to change , but their positions would change by very little . the magnitude of the forces they exert on pluto would therefore change by an amount much less than the current magnitude of those forces .
thus , we argue that the approach used in this paper , while not as accurate as a full second order calculation involving all the outer planets , is accurate enough for our purposes .    given the comments above , how then are we to compare alternative models ? there are a number of methods that we can use to compare and assess models and their results . in the orbit fitting context , these have been discussed above in terms of point estimates and global estimates of goodness of fit .
the simplest point estimate is to simply calculate the rms residual of the fit compared with the input observations .
however , as we ve seen , the ill - conditioned nature of the orbital fitting problem can sometimes make the residual a poor candidate for this role .
errors , in the linear approximation , to the orbital elements can be assessed to determine the quality of the fit .
further physical insight can be gathered from inspecting the sky position errors that result from element errors .
global assessments of model fit primarily revolve around the normality of residuals .
there are lack of fit tests that can be used to test the residuals and determine if there are any indications of a deficient model .
generally , a lack of fit is manifested as patterns , biases , and systematic variations in the residuals , which would indicate a poorly fitting model .    which of these metrics is best ?
the reality is that determining the quality of a model and comparing model effectiveness is a complex problem .
we need all these measures of model fit .
statistical tests should be performed comparing predicted positions , taking proper account of the associated errors , to test the hypothesis that the positions predicted under the two dynamical models are different . only in this way
can the existence of a perturbation like the pioneer effect be falsified through astrometric methods .    as discussed earlier
, pluto is very interesting in that orbital fits to the two gravity models show deviations in predicted position that are probably detectable if the entire observation arc consisted of reliable measurements .
however , prior to about 1960 , there is a significantly greater dispersion in residuals that at later times .
only about half of pluto s motion since its discovery has been the result of systematic and organized observing campaigns @xcite .
what adds special interest , however , is that because of its relatively large eccentricity pluto is likely to show differences in position predictions for the two gravity models in the relatively near future as more observations are accumulated .
the analysis described above shows two major things .
first , we must fit observations to a particular dynamical model and adjust orbital elements before predicted positions on the sky can be compared .
orbital parameters are derived from observations which have associated an unavoidable error .
the determination of orbits is a model fitting process which has its own associated error sources .
extrapolating sky positions very far past the end of an observation arc can result in predicted observations becoming inaccurate so rapidly as to be worthless .
the implication of these findings is that `` matched '' and `` mismatched '' gravity models can not be distinguished on the basis of observable sky positions for observation arc lengths similar to those currently obtaining for pluto .
similarly , in order to draw conclusions about differences in position in the sky , we must be careful not to oversimplify the dynamical model used to draw the conclusions .
suppressing degrees of freedom in the dynamics simplifies the orbital determination problem to just such a degree .
the orbital determination problem is nonlinear and the customary solution methods are approximations .
thus , any missing or ignored degrees of freedom can , if present , conceal dynamical effects associated with differing gravity models ; we are forced to make use of the full dimensionality of the dynamical problem . in particular
, a substantial amount of variation can be absorbed into a multidimensional parameter space and the full parameter space must be considered to properly reflect differences in motion of the outer planets due to the pioneer effect .
the problem with the simplified approach is a two - headed one .
first , the orbital fitting problem is inherently nonlinear and is normally solved in the linear approximation .
even if not mathematically chaotic , the system of equations is sensitively dependent upon initial conditions .
thus , small changes in elements can result in large changes in predicted position outside the range of observations .
this sensitivity is exacerbated by the problem of a short observation arc .
the length of the entire observational archive for pluto is less than about one - third of a complete revolution .
together , these factors conspire to potentially generate large errors outside the observation arc , while increasing the length of the observation arc can markedly reduce error over the whole of the arc and even beyond it .
we must conclude that we do not know the orbit of pluto as well as we might have thought
. we must continue to perform astrometry on it in order to be able to comment on the accuracy with which we know its orbit . using current data
, we can not assert that the motion of pluto demonstrates that the pioneer effect does not exist .
that jury is still out . of course
, this does not mean that the pioneer effect exists .
it does mean that we can not deny the existence of the pioneer effect on the basis of motions of the pluto as currently known .
further observations are required before such an assertion can be made with confidence .
it should again be emphasized that our approach is itself an approximation .
the dynamical system that should be analyzed to provide a comprehensive answer to the question of the detectabilty of the pioneer effect should include not only adjustment of the orbital elements of pluto , but also simultaneous adjustment to those of uranus and neptune as well .
it is only in this way that all the second order perturbations to the system can be taken into account .
our approach here , however serves to illustrate the ideas concerned and the weaknesses of the approaches outlined above .
finally , it should be pointed out that , in addition to the observations of individual objects discussed in this paper , there are other related approaches to assessing gravity in the outer solar system . for example , recently @xcite have investigated the use of ensembles of trans - neptunian objects ( tnos ) to ascertain whether their motion reflects unknown additional perturbations and showed that the pioneer effect was not consistent with the motion of tnos . on the other hand , in the area of the observation of individual objects as discussed in this paper , the advent of pan - starrs and lsst in the next several years should provide sufficient data to determine whether the motion of outer solar system bodies reflect the action of unknown forces . this determination should occur over time frames discussed in this paper .
however , a combination of the techniques of @xcite and the considerations presented here should provide definitive answers more quickly .
the authors wish to thank dr .
myles standish for providing the residuals of jpl s de414 ephemeris .
we also wish to thank the anonymous referee for a close reading of the manuscript and comments that resulted in significant improvements to the paper .            ,
p.  r. , & robinson , d.  k. 2003 , data reduction and error analysis for the physical sciences ( data reduction and error analysis for the physical sciences , 3rd ed .
, by philip r.  bevington , and keith d.  robinson .
boston , ma : mcgraw - hill , isbn 0 - 07 - 247227 - 8 , 2003 . )                                          , p.  k. 1992 , explanatory supplement to the astronomical almanac ( explanatory supplement to the astronomical almanac , by p.  kenneth seidelmann .  published by university science books , isbn 0 - 935702 - 68 - 7 , 752pp , 1992 . )