对火星轨道变化问题的最后解释 (第5/5页)

cmap.Similarly，thelongitudeerrorofPlutocanbeestimatedas12°.ThisvalueforPlutoismuchbetterthantheresultinKinoshita&Nakai(1996)wherethedifferenceisestimatedas60°.

    3Numericalresults–I.Glanceattherawdata

    Inthissectionwebrieflyreviewthelong-termstabilityofplanetaryorbitalmotionthroughsomesnapshotsofrawnumericaldata.Theorbitalmotionofplanetsindicateslong-termstabilityinallofournumericalintegrations:noorbitalcrossingsnorcloseencountersbetweenanypairofplanetstookplace.

    3.1Generaldescriptionofthestabilityofplanetaryorbits

    First，webrieflylookatthegeneralcharacterofthelong-termstabilityofplanetaryorbits.Ourinterestherefocusesparticularlyontheinnerfourterrestrialplanetsforwhichtheorbitaltime-scalesaremuchshorterthanthoseoftheouterfiveplanets.AswecanseeclearlyfromtheplanarorbitalconfigurationsshowninFigs2and3，orbitalpositionsoftheterrestrialplanetsdifferlittlebetweentheinitialandfinalpartofeachnumericalintegration，whichspansseveralGyr.Thesolidlinesdenotingthepresentorbitsoftheplanetsliealmostwithintheswarmofdotseveninthefinalpartofintegrations(b)and(d).Thisindicatesthatthroughouttheentireintegrationperiodthealmostregularvariationsofplanetaryorbitalmotionremainnearlythesameastheyareatpresent.

    Verticalviewofthefourinnerplanetaryorbits(fromthez-axisdirection)attheinitialandfinalpartsoftheintegrationsN±1.Theaxesunitsareau.Thexy-planeissettotheinvariantplaneofSolarsystemtotalangularmomentum.(a)TheinitialpartofN 1(t=0to0.0547×109yr).(b)ThefinalpartofN 1(t=4.9339×108to4.9886×109yr).(c)TheinitialpartofN1(t=0to0.0547×109yr).(d)ThefinalpartofN1(t=3.9180×109to3.9727×109yr).Ineachpanel，atotalof23684pointsareplottedwithanintervalofabout2190yrover5.47×107yr.Solidlinesineachpaneldenotethepresentorbitsofthefourterrestrialplanets(takenfromDE245).

    ThevariationofeccentricitiesandorbitalinclinationsfortheinnerfourplanetsintheinitialandfinalpartoftheintegrationN 1isshowninFig.4.Asexpected，thecharacterofthevariationofplanetaryorbitalelementsdoesnotdiffersignificantlybetweentheinitialandfinalpartofeachintegration，atleastforVenus，EarthandMars.TheelementsofMercury，especiallyitseccentricity，seemtochangetoasignificantextent.Thisispartlybecausetheorbitaltime-scaleoftheplanetistheshortestofalltheplanets，whichleadstoamorerapidorbitalevolutionthanotherplatheinnermostplanetmaybenearesttoinstability.ThisresultappearstobeinsomeagreementwithLaskar‘s(1994，1996)expectationsthatlargeandirregularvariationsappearintheeccentricitiesandinclinationsofMercuryonatime-scaleofseveral109yr.However，theeffectofthepossibleinstabilityoftheorbitofMercurymaynotfatallyaffecttheglobalstabilityofthewholeplanetarysystemowingtothesmallmassofMercury.Wewillmentionbrieflythelong-termorbitalevolutionofMercurylaterinSection4usinglow-passfilteredorbitalelements.