Ви зробили ту саму принципову помилку, яку зробив Антон Громов у своєму питанні про сестринський веб-сайт космічної розвідки StackExchange: Ви використовували барицентр Сонячної системи, а не Сонце як джерело кадру. Якби ви використовували Сонце, явна розбіжність зменшилась би майже на два порядки.
Набагато меншою вадою є те, що ви, мабуть, використовували опівночі, а не полудень. Епоха J2000.0 - це 12 полудня земного часу 1 січня 2000 року. Це не сильно змінить речі, але все-таки важливо. Та компенсація за півдня перекусила мене не раз.
Координата z положення Землі не впала б до нуля, якби ви внесли вищевказані виправлення (барицентр, орієнтований на Сонце, а не на сонячну систему, і 12 полудня, а не півночі). Це не впало б до нуля, навіть якби ви використовували барицентр Земля-Місяць, а не центр Землі.
Однією з причин того, що вона не опустилася б до нуля, є те, що орбіта Землі навколо Сонця не є еліпсом. Земна орбіта навіть не лежить на площині! Основна причина цього непланарного руху полягає в тому, що орбіта Місяця навколо Землі нахилена приблизно на 5,15 ° по відношенню до екліптики, завдяки чому Земля піднімається вгору і вниз по відношенню до екліптики 1 . Оскільки орбіта Землі не є справді еліптичною або навіть плоскою, площина екліптики в певну епоху є середньозваженою за часом площиною, орієнтованою на час епохи, що робить положення Землі відносно Сонця (або положення Сонця щодо Землі) в середньому до нуля.
Меншою причиною є те, що середній рівень екліптики та середнього рівня рівнодення J2000.0 був визначений понад 36 років тому. HORIZONS використовує DE431, який був випущений 6 років тому. Здатність JPL моделювати сонячну систему значно покращилася за три десятиліття.
1 Строго кажучи, це змушує Сонце битись вгору і вниз по відношенню до екліптики. Поняття площини екліптики датується ще древніми греками, і воно зберігає рудицію цього визначення і донині.