
Последнее десятилетие принципиально изменило наши представления о строении, динамической эволюции и устойчивости Солнечной системы. Привычными стали сообщения об открытии новых объектов, выявлении новых динамических структур, проявлении свойств неустойчивости движения или хаотического поведения у тех или иных групп объектов.
Это вызвано несколькими причинами: появление новых инструментов и модернизация старых, применение высокочувствительных ПЗС–матриц и новых методов математической обработки результатов наблюдений. Все это позволяет наблюдать новые объекты, имеющие очень малую яркость и существенное собственное движение.
Новые аналитические и численные методы небесной механики в совокупности с современными вычислительными системами дают возможность моделировать движение тел Солнечной системы на интервалах времени, сравнимых с ее возрастом и даже многократно превышающих его.
На наших глазах происходит смена представлений о динамике Солнечной системы: от регулярной и устойчивой к хаотической и неустойчивой. Все это напоминает ситуацию в физике начала XX века, когда совершался переход от классической к релятивистской картине Мира. Нам предстоит разобраться где, когда и при каких условиях мы можем рассматривать Солнечную систему регулярной и устойчивой, а в каких случаях проявляются признаки хаоса и неустойчивости.
Начнем рассмотрение с современных представлений о структуре Солнечной системы. Затем обсудим понятия устойчивости и неустойчивости движения, условия возникновения резонансов и хаотического поведения. После этого проанализируем динамику малых тел Солнечной системы и обратимся к большим планетам. В заключение рассмотрим динамику Солнечной системы как целого на временах, сравнимых с ее возрастом.
| Начало | "Движение тел Солнечной системы" | Закономерности движения | Структура, динамика и устойчивость |
Все объекты Солнечной системы можно разделить на четыре группы: Солнце, большие планеты, спутники планет и малые тела. Мы пока ничего не говорим о спутниках малых тел, поскольку к настоящему времени таких объектов открыто всего два, а наблюдательной информации недостаточно, чтобы детально исследовать их динамику.
Солнце — динамический центр системы. Его гравитационное влияние является доминирующим в Солнечной системе за исключением малых областей в окрестности других объектов.
Большие
планеты — визитная карточка Солнечной системы.
Пять ближайших к Земле больших планет были
известны с ранней истории человечества. Это — Меркурий, Венера, Марс, Юпитер и Сатурн.
История открытия трех других больших планет
показывает как менялось отношение астрономов к
вопросу о размерах и населении Солнечной
системы.
Открытие Урана явилось сюрпризом. Весной 1781 г. Вильям Гершель на своем 7-футовом (2.1 м) телескопе проводил наблюдения по программе определения параллаксов звезд. 13 марта 1781 г. он сделал запись об обнаружении туманной звезды или кометы. Спор о природе открытого объекта продолжался до 1787 г., когда Гершель открыл два спутника Урана: Оберон и Титанию.
Открытие Нептуна стало триумфом теории тяготения Ньютона. Анализируя неравенства в движении Урана, Бессель в Кенигсберге в 1840 г., Адамс в Кембридже в 1841 г. и Леверье во Франции в 1845 г. независимо друг от друга рассчитали орбиту планеты, ответственной за эти возмущения. 23 сентября 1846 г. Галле и д’Аррест из Берлинской обсерватории по эфемеридам Леверье открыли Нептун.
Открытие Плутона можно назвать запрограммированным. В 1896 г. Персиваль Ловелл обнаружил остаточные невязки в движении Урана после учета возмущений от Нептуна и высказал гипотезу, что эти возмущения производятся неизвестной занептунной планетой. В середине 90-х годов XIX века в Аризоне Ловелл построил обсерваторию, которая стала центром поиска новой планеты. В течение почти 30 лет было проведено несколько компаний по поиску Плутона. Но безрезультатно. В 1916 г. умер Ловелл. В 1929 г. Клод Томбо на 13-дюймовом (0.33 м) рефракторе начал новую атаку на Плутон. Открытие пришло 18 февраля 1930 г., когда Томбо сравнивал фотопластинки, полученные 23 и 29 января 1930 г. Директор Ловелловской обсерватории сообщил об открытии 13 марта 1930 г. в 149-ю годовщину открытия Урана Гершелем и 75-ю годовщину со дня рождения Персиваля Ловелла. За время поиска Плутона было проведено сравнение около 90 млн. изображений звезд в течение 7000 часов на блинк-компараторе.
Существуют
ли большие планеты за орбитой Плутона? Анализ
траекторий движения тел Солнечной системы и
космических аппаратов Пионер–10,
Пионер–11, Вояджер–1,
Воджер–2 позволяют утверждать, что объектов,
сравнимых с Плутоном, и более крупных во внешней
области Солнечной системы не существует.
История открытия спутников планет не менее драматична, но мы не будем на ней останавливаться. Отметим только, что спутниковые системы планет-гигантов сложностью своего устройства зачастую превосходят Солнечную систему. Не до конца решен вопрос о происхождении двойных планет Земля–Луна и Плутон–Харон.
Малые тела Солнечной системы — пробный камень и золотая жила небесной механики, кладезь новых открытий. Самые известные малые тела — кометы. Упоминания о кометах можно найти в легендах и летописях практически всех народов Земли. По динамическим признакам кометы разделяются на долгопериодические и короткопериодические.
Долгопериодические
кометы движутся по орбитам, большие полуоси
которых достигают десятков тысяч
астрономических единиц, а периоды обращения —
десятков миллионов лет. Орбиты сильно вытянуты,
их эксцентриситеты близки к единице. Ориентация
орбит и их наклоны к плоскости эклиптики
распределены случайным образом. В настоящее
время имеются сведения более, чем о 700 таких
комет.
Короткопериодические кометы имеют периоды менее 200 лет, умеренные эксцентриситеты, для большинства из них наклон орбит к плоскости эклиптики не превышает 35°
. Короткопериодические кометы делятся на семейства по признаку планеты-гиганта, определяющей динамику кометы. В настоящее время известно около 180 короткопериодических комет. Большинство из них принадлежит семейству Юпитера.Самая многочисленная популяцию малых тел Солнечной системы — астероиды. Первый астероид — Церера — был открыт в первый день XIX века сицилийским астрономом Пиацци. Хотя открытие и носило случайный характер, оно послужило толчком к разработке Гауссом классического метода определения орбит по трем наблюдениям и метода наименьших квадратов, благодаря которым удалось вычислить орбиту и переоткрыть Цереру спустя почти год после первых наблюдений. В настоящее время известно несколько десятков тысяч астероидов. И это число стремительно растет.
Популяция астероидов неоднородна. Большинство астероидов движутся по орбитам близким к круговым в поясе астероидов между орбитами Марса и Юпитера. В 1866 г. Кирквуд исследовал зависимость числа астероидов от больших полуосей их орбит и обнаружил, что полученное распределение имеет несколько глубоких минимумов. Позднее выяснилось, что эти минимумы соответствуют соизмеримости средних движений Юпитера и астероида. Они получили название люков Кирквуда.
На рисунке показана зависимость числа астероидов N от значения большой полуоси. Пунктирной линией отмечены резонансные значения больших полуосей. Соответствующие соизмеримости средних движений астероида и Юпитера даются над верхней рамкой. Видно, что большинству резонансных значений большой полуоси соответствует существенное уменьшение числа астероидов. (Yoshikawa M. A survey of the motions of asteroids in the commensurabilities with Jupiter. Astronomy and Astrophysics. 1989. 213. 436–458.)
Две группы астероидов расположены в окрестности треугольных точек либрации системы Солнце–Юпитер–астероид — это греки и троянцы. Они движутся в окрестности орбиты Юпитера: греки опережают Юпитер примерно на 60°, а троянцы на столько же отстают. Хотя астероиды движутся по эллиптическим орбитам, треугольник Солнце–Юпитер–астероид всегда остается близким к равностороннему. Иногда обе группы астероидов называют троянцами. По состоянию на 1 апреля 1999 г. известно 476 астероидов-троянцев (474 у Юпитера и 2 у Марса).
Еще одна группа астероидов — астероиды, сближающиеся с Землей. Их перигелийные расстояния меньше 1.33 а.е. В настоящее время известно несколько тысяч таких астероидов. Около сотни из них представляют реальную угрозу для Земли: они пересекают ее орбиту и имеют размер более 1 км. Столкновение Земли с подобным астероидом вызовет глобальную катастрофу, подобную той, что привела к вымиранию динозавров. Имеется еще около тысячи астероидов размером от 30 до 50 м, также пересекающих орбиту Земли. Столкновение Земли с таким астероидом способно вызвать локальную катастрофу типа тунгусской. Однако, ни один из известных астероидов не столкнется с Землей в ближайшем будущем, в течение 33 лет, в 21 веке.
После
открытия Плутона неоднократно предпринимались
попытки поиска десятой большой планеты
Солнечной системы. Во время одного из таких
обзоров 18 октября 1977 г. Коваль открыл малую
планету 2060 Хирон, которая движется между
орбитами Юпитера и Урана, пересекая орбиту
Сатурна. Вблизи перигелия у этого “астероида”
проявляются признаки газоизвержения и комы.
Более 14 лет этот объект оставался единственной
малой планетой, наблюдаемой глубоко внутри
области движения планет-гигантов. 9 января 1992
на автоматическом телескопе Космический дозор
(Аризона, США) был открыт еще один астероид этой
группы — 5145 Фолус. К настоящему времени
известно 7 астероидов группы Кентавра,
движущихся среди планет-гигантов между орбитами
Юпитера и Нептуна. Название группы отражает тот
факт, что объекты одновременно имеют признаки и
астероидов и комет. В табл. 1 приводится список
астероидов группы Кентавра по состоянию на
1 августа 1997 г. В таблице даны: имя астероида,
его предварительное обозначение, перигелийное и
афелийное расстояния в астрономических
единицах, наклон орбиты в градусах,
эксцентриситет орбиты, большая полуось в
астрономических единицах и дата открытия. Полный
регулярно обновляемый вариант таблицы доступен
по адресу http://cfa-www.harvard.edu/cfa/ps/lists/Centaurs.html.
Таблица 1. Астероиды группы Кентавра (по состоянию на 1 августа 1997 г.)
Имя |
Обозначение |
q, а.е. |
Q, а.е. |
i |
e |
a, а.е. |
Дата открытия |
| 1997 CU26 | 12.913 |
18.724 |
23.4 |
0.184 |
15.819 |
15.02.1997 |
|
| 1995 GO | 6.842 |
29.354 |
17.6 |
0.622 |
18.098 |
5.04.1995 |
|
| 1995 DW2 | 18.863 |
31.036 |
4.2 |
0.244 |
24.950 |
27.02.1995 |
|
| 1994 TA | 11.717 |
21.939 |
5.4 |
0.304 |
16.828 |
2.10.1994 |
|
| 7066 Nessus | 1993 HA2 | 11.825 |
37.469 |
15.6 |
0.520 |
24.647 |
26.04.1993 |
| 5145 Pholus | 1992 AD | 8.675 |
31.812 |
24.7 |
0.571 |
20.244 |
9.01.1992 |
| 2060 Chiron | 1977 UB | 8.454 |
18.877 |
6.9 |
0.381 |
13.666 |
18.10.1977 |
В 1949 г. К.Эджеворт высказал предположение о существовании остаточного неизрасходованного при формировании Солнечной системы материала за орбитой Нептуна. Однако, эта работа была малоизвестна до последнего времени. В 1951 г. Койпер предположил, что кометы и астероиды формировались в существенно различных областях Солнечной системы и, что за орбитой Плутона должен существовать пояс комет. 30 августа 1992 г. Джевитт и Лю (Гавайский университет, США) открыли первый объект, принадлежащий поясу Койпера. Он получил обозначение 1992 QB1. Сейчас известно 53 объекта, движущихся за орбитой Нептуна. В табл. 2 приводится список объектов пояса Койпера по состоянию на 1 августа 1997 г.
Таблица 2. Объекты пояса Койпера (по состоянию на 1 августа 1997 г.)
Обозначение |
q, а.е. |
Q, а.е. |
i |
e |
a, а.е. |
Дата открытия |
| 1997 CW29 | 36.271 |
42.480 |
19.0 |
0.079 |
39.375 |
8.02.1997 |
| 1997 CV29 | 39.984 |
48.469 |
7.8 |
0.096 |
44.227 |
6.02.1997 |
| 1997 CU29 | 44.571 |
44.571 |
1.5 |
0.0 |
44.571 |
6.02.1997 |
| 1997 CT29 | 40.803 |
44.985 |
1.0 |
0.049 |
42.894 |
2.02.1997 |
| 1997 CS29 | 43.320 |
43.945 |
2.3 |
0.007 |
43.633 |
3.02.1997 |
| 1997 CR29 | 41.996 |
41.996 |
20.2 |
0.0 |
41.996 |
3.02.1997 |
| 1997 CQ29 | 41.153 |
47.671 |
2.9 |
0.073 |
44.412 |
4.02.1997 |
| 1996 TS66 | 37.831 |
50.571 |
7.4 |
0.144 |
44.201 |
12.10.1996 |
| 1996 TR66 | 35.955 |
42.727 |
12.9 |
0.086 |
39.341 |
8.10.1996 |
| 1996 TQ66 | 34.611 |
44.598 |
14.6 |
0.126 |
39.605 |
8.10.1996 |
| 1996 TP66 | 26.354 |
52.985 |
5.7 |
0.336 |
39.669 |
11.10.1996 |
| 1996 TO66 | 38.714 |
48.172 |
27.4 |
0.109 |
43.443 |
12.10.1996 |
| 1996 TL66 | 35.055 |
133 |
24.0 |
0.583 |
84.127 |
9.10.1996 |
| 1996 TK66 | 42.474 |
42.474 |
4.9 |
0.0 |
42.474 |
9.10.1996 |
| 1996 SZ4 | 30.487 |
48.125 |
4.7 |
0.224 |
39.306 |
16.09.1996 |
| 1996 RR20 | 33.654 |
45.223 |
5.3 |
0.147 |
39.439 |
15.09.1996 |
| 1996 RQ20 | 33.437 |
61.523 |
31.8 |
0.296 |
47.480 |
6.09.1996 |
| 1996 KY1 | 35.712 |
43.322 |
30.9 |
0.096 |
39.517 |
16.05.1996 |
| 1996 KX1 | 35.704 |
43.381 |
1.5 |
0.097 |
39.543 |
22.05.1996 |
| 1996 KW1 | 46.602 |
46.602 |
5.5 |
0.0 |
46.602 |
22.05.1996 |
| 1996 KV1 | 41.188 |
44.744 |
8.4 |
0.041 |
42.966 |
21.05.1996 |
| 1995 YY3 | 30.695 |
47.913 |
0.4 |
0.219 |
39.304 |
24.12.1995 |
| 1995 WY2 | 43.574 |
47.697 |
1.7 |
0.045 |
45.636 |
18.11.1995 |
| 1995 QZ9 | 33.384 |
46.368 |
19.5 |
0.163 |
39.876 |
29.08.1995 |
| 1995 QY9 | 28.978 |
51.787 |
4.8 |
0.282 |
40.382 |
31.08.1995 |
| 1995 KK1 | 31.981 |
46.969 |
9.3 |
0.190 |
39.475 |
30.05.1995 |
| 1995 KJ1 | 43.468 |
43.468 |
2.7 |
0.0 |
43.468 |
30.05.1995 |
| 1995 HM5 | 29.485 |
49.394 |
4.8 |
0.252 |
39.439 |
26.04.1995 |
| 1995 GY7 | 41.347 |
41.347 |
0.9 |
0.0 |
41.347 |
6.04.1995 |
| 1995 GA7 | 34.751 |
44.160 |
3.5 |
0.119 |
39.455 |
3.04.1995 |
| 1995 GJ | 39.006 |
46.808 |
22.9 |
0.091 |
42.907 |
3.04.1995 |
| 1995 FB21 | 42.426 |
42.426 |
0.7 |
0.0 |
42.426 |
29.03.1995 |
| 1995 DC2 | 40.756 |
46.988 |
2.3 |
0.071 |
43.872 |
24.02.1995 |
| 1995 DB2 | 40.083 |
52.451 |
4.1 |
0.134 |
46.267 |
24.02.1995 |
| 1995 DA2 | 33.704 |
38.628 |
6.6 |
0.068 |
36.166 |
24.02.1995 |
| 1994 VK8 | 41.500 |
44.027 |
1.5 |
0.030 |
42.763 |
8.11.1994 |
| 1994 TG2 | 42.448 |
42.448 |
2.2 |
0.0 |
42.448 |
8.10.1994 |
| 1994 TH | 40.940 |
40.940 |
16.1 |
0.0 |
40.940 |
3.10.1994 |
| 1994 TG | 42.254 |
42.254 |
6.8 |
0.0 |
42.254 |
3.10.1994 |
| 1994 TB | 27.052 |
52.505 |
12.1 |
0.320 |
39.779 |
2.10.1994 |
| 1994 JR1 | 34.756 |
44.265 |
3.8 |
0.120 |
39.510 |
12.05.1994 |
| 1994 JQ1 | 41.819 |
46.260 |
3.8 |
0.050 |
44.039 |
11.05.1994 |
| 1994 JV | 35.251 |
35.251 |
18.1 |
0.0 |
35.251 |
13.05.1994 |
| 1994 JS | 33.042 |
51.707 |
14.0 |
0.220 |
42.375 |
11.05.1994 |
| 1994 GV9 | 41.023 |
46.055 |
0.6 |
0.058 |
43.539 |
15.04.1994 |
| 1994 EV3 | 40.914 |
44.741 |
1.7 |
0.045 |
42.828 |
13.03.1994 |
| 1994 ES2 | 40.299 |
50.835 |
1.1 |
0.116 |
45.567 |
13.03.1994 |
| 1993 SC | 32.216 |
47.456 |
5.1 |
0.191 |
39.836 |
17.09.1993 |
| 1993 SB | 26.912 |
52.039 |
1.9 |
0.318 |
39.476 |
16.09.1993 |
| 1993 RP | 34.863 |
43.795 |
2.6 |
0.114 |
39.329 |
15.09.1993 |
| 1993 RO | 31.487 |
47.615 |
3.7 |
0.204 |
39.551 |
14.09.1993 |
| 1993 FW | 41.544 |
45.631 |
7.8 |
0.047 |
43.587 |
28.03.1993 |
| 1992 QB1 | 40.884 |
47.587 |
2.2 |
0.076 |
44.236 |
30.08.1992 |
Полный регулярно обновляемый вариант таблицы доступен по адресу http://cfa-www.harvard.edu/cfa/ps/lists/TNOs.html. Некоторые исследователи относят к объектам пояса Койпера и Плутон. Возможно, что пояс Койпера является внутренней областью облака Оорта — сферического образования радиусом от тысяч до сотен тысяч астрономических единиц, являющегося резервуаром долгопериодических комет.
Таким образом, по современным представлениям Солнечная система имеет следующую структуру: вокруг Солнца вращаются 9 больших планет, между орбитами Марса и Юпитера находится пояс астероидов, часть астероидов движется среди планет земной группы и в окрестности треугольных точек либрации Юпитера, среди планет-гигантов движутся объекты группы Кентавра и короткопериодические кометы, за орбитой Нептуна располагается пояс Койпера, а вся система окружена облаком Оорта.
Нам предстоит разобраться каким образом в Солнечной системе сформировалась подобная динамическая структура, как она эволюционирует и является ли устойчивой. Мы обсудим сценарии происхождения люков Кирквуда. Узнаем откуда берутся астероиды, угрожающие Земле. Выясним, почему Кентавров намного меньше, чем объектов пояса Койпера. Найдем источник короткопериодических комет. И, наконец, попытаемся ответить на главный вопрос: устойчива ли Солнечная система?
Но прежде, чем исследовать свойства движения тел Солнечной системы, мы определим такое важное понятие как устойчивость.
| Начало | "Движение тел Солнечной системы" | Закономерности движения | Структура, динамика и устойчивость |
До открытия закона всемирного тяготения вопрос об устойчивости системы Мира решался a priori. Сначала делалось предположение о фундаментальном свойстве Мира — устойчивость, хаотичность и т.п. — затем создавалась система. По свой сути эти системы были кинематическими.
Ньютон
первым построил динамическую модель Солнечной
системы и сразу же столкнулся с вопросом о ее
устойчивости. Он вышел из этого затруднения с
помощью Великого Часовщика, который время от
времени должен возвращать планеты на их орбиты.
В дальнейшем понятие устойчивости
развивалось параллельно с исследованиями
движения планет.
Лагранж
считал движение устойчивым, если оно происходит
в замкнутой области пространства. Согласно
теореме Лапласа–Лагранжа (1773, 1776 гг.) об
отсутствии вековых возмущений больших полуосей
планетных орбит, их изменение с точностью до
величин первого порядка малости относительно
возмущающих масс можно представить в виде суммы
тригонометрических слагаемых.
На основе этой теоремы в 1773 г. Лаплас
сформулировал теорему об устойчивости Солнечной
системы: если движение планет происходит в
одном направлении, их массы одного порядка,
эксцентриситеты и наклоны малы, а большие
полуоси испытывают лишь небольшие колебания
относительно среднего положения, то
эксцентриситеты и наклоны орбит будут
оставаться малыми на рассматриваемом интервале
Позднее, в 1809 г.
Пуассон показал, что во втором приближении
решение для большой полуоси может содержать
тригонометрические слагаемые с амплитудой,
пропорциональной времени. Соответствующее
определение устойчивости по Пуассону требует,
чтобы частица за бесконечное время бесконечное
число раз проходила через сколь угодно малую
окрестность начальной точки.
Наиболее удачное понятие устойчивости сформулировал в конце XIX века русский математик А.М.Ляпунов. Исследуемое движение считается устойчивым, если все возможные движения, мало отличающиеся от него в начальный момент, в последующем будут мало отклоняться от него на всем интересующем интервале времени
. Если же найдется хотя бы одно (!) движение, в начальный момент мало отличающееся от исследуемого, которое постепенно, пусть и через большой промежуток времени, заметно отклонится от него, то исследуемое движение — неустойчиво. Это определение считается основным по сей день.В задачах небесной механики рассматривается устойчивость по части переменных: большой полуоси (задает размер орбиты), эксцентриситету (определяет вытянутость орбиты) и наклону орбиты. Совместная эволюция этих трех элементов должна удовлетворять определению устойчивости по Ляпунову. Отсюда, Солнечная система устойчива по Ляпунову, если размеры, форма и наклоны орбит остаются близкими к начальным на всем рассматриваемом интервале времени.
Говоря об устойчивости Солнечной системы, как правило, имеют ввиду устойчивость движения больших планет на бесконечном или очень большом, сравнимом с ее возрастом, интервале времени. В этом случае крайними проявлениями неустойчивости являются уход из Солнечной системы, падение на Солнце или столкновение с другой планетой. Такое событие способно существенно изменить структуру и динамику Солнечной системы.
В последнее время появилось большое число исследований по динамике малых тел Солнечной системы. Для этих объектов неустойчивость чаще всего проявляется в существенном изменении размеров, формы и наклона орбиты, зачастую многократном, и только после этого происходит уход из Солнечной системы, падение на Солнце или столкновение с планетой. Благодаря тому, что массы этих объектов незначительны, а их число велико, неустойчивое движение одного объекта, или даже группы, почти не влияет на устойчивость Солнечной системы в целом, хотя и играет важную роль в формировании ее структуры.
Классический
метод исследования движения небесных тел
заключается в представлении решения
соответствующих уравнений возмущенного
движения в виде отрезков рядов. Однако, в конце XIX
века Анри Пуанкаре показал, что ряды, применяемые
для описания движения небесных тел, расходятся.
Следовательно их нельзя использовать для
анализа поведения Солнечной системы на
бесконечном интервале времени. Казалось бы это
не страшно, ведь Солнечная система имеет хотя и
большую, но ограниченную эволюционную шкалу. Но
полученные Пуанкаре интервалы применимости
классических рядов, оказались значительно
короче возраста Солнечной системы.
Эта трудность была преодолена лишь в начале 60-х гг. нашего столетия. Советские математики А.Н.Колмогоров и В.И.Арнольд и американский математик Ю.Мозер разработали теорию, получившую название КАМ–теория. Ее приложение к Солнечной системе дает следующую теорему: если массы планет достаточно малы, эксцентриситеты и наклоны орбит малы, то для большинства начальных условий (исключая резонансные и близкие к ним) движение будет условно-периодическим, эксцентриситеты и наклоны будут оставаться малыми, а большие полуоси будут вечно колебаться вблизи своих первоначальных значений
, то есть Солнечная система будет устойчивой по Ляпунову на бесконечном интервале времени.Результат замечательный! Однако, с очень существенной оговоркой: при условии отсутствия резонансов. К сожалению, в реальной Солнечной системе резонансы играют очень важную роль. Поэтому, выводы КАМ–теории не могут быть применены к Солнечной системе в целом на всем интервале ее существования.
Каковы условия возникновения резонансов? Главное — соизмеримость частот. Различают два основных типа резонансов: резонанс средних движений, возникающий вследствие соизмеримости средних движений объектов, и вековой резонанс между средними движениями долгот восходящих узлов и аргументов перицентров объектов. Наиболее сильно проявляются резонансы низких порядков, когда отношение частот можно представить в виде рациональной дроби с небольшими значениями числителя и знаменателя, например, 2/1, 3/1, 3/2, 4/3 и т.п.
При наличии резонансов эволюция динамической системы может идти следующим образом. 1) Система пройдет через резонанс, что приведет к резкому скачкообразному изменению элементов орбиты, например, эксцентриситета или наклона. 2) Система застрянет в резонансе и перейдет в новое состояние с либрационным режимом движения, в котором позиционные элементы (большая полуось, эксцентриситет, наклон) вместе или по отдельности будут испытывать колебания, иногда достаточно большой амплитуды. Любой из этих сценариев может привести к тому, что объект перейдет на новую орбиту. В этом случае движение будет неустойчивым по Ляпунову.
Как формируется хаотический режим поведения динамической системы? Как правило, хаотическое поведение является следствием взаимодействия резонансов в результате перекрытия соседних резонансных зон. В этом случае формируется зона хаоса, характерная тем, что основные динамические параметры системы, в частности позиционные элементы орбиты, очень быстро изменяют свои значения, причем последующее состояние системы практически не зависит от предыдущего.
Хаотическое движение может наблюдаться и при отсутствии резонансов, например, при гравитационных маневрах и сближениях малых тел с большими планетами. В этом случае направление движения тела после маневра или сближения сильно зависит от начальных данных. После двух-трех маневров, тела, стартовавшие практически из одной точки, могут оказаться на существенно различных орбитах. Такое движение можно сравнить с перемещением по крупному городу с большим числом пересадок. Небольшая разница в положении или моменте старта может привести к критически различным финалам.
Для анализа скорости разбегания соседних траекторий применяют специальную характеристику, называемую временем Ляпунова. Она определяет промежуток времени, в течение которого расстояние между соседними траекториями увеличивается в
e раз.Наиболее сильно хаотическое поведение проявляется у малых тел Солнечной системы. Именно с них мы и начнем обзор динамической эволюции объектов Солнечной системы.
| Начало | "Движение тел Солнечной системы" | Закономерности движения | Структура, динамика и устойчивость |
Обзор динамики малых тел Солнечной системы мы начнем с астероидов и попытаемся разобраться как образовались люки Кирквуда и почему они существуют по сей день.
Как мы уже упоминали, люки соответствуют резонансным значениям больших полуосей орбит астероидов. Мы рассмотрим три резонанса средних движений между астероидом и Юпитером: 3/1 — люк Гестии, 2/1 — люк Гекубы и 3/2 — группа Гильды.
В 1982 г. Ж.Висдом (Калифорнийский технологический институт, США) исследовал движение астероидов, находящихся в резонансе 3/1 в рамках ограниченной задачи трех тел Солнце–Юпитер–астероид на интервале времени 20 млн. лет. В отличие от своих предшественников он учел эксцентриситет орбиты Юпитера. Висдом обнаружил несколько астероидов, которые на рассматриваемом интервале скачком увеличивали эксцентриситет своей орбиты от 0 до 0.35. Такого эксцентриситета достаточно, чтобы астероид стал двигаться по орбите, пересекающей орбиту Марса. Эффективное время сближения с Марсом порядка 100–200 млн. лет. После сближения с Марсом астероид уйдет из главного пояса в область планет земной группы. Более поздние исследования показали, что в окрестностях резонанса 3/1 эксцентриситет может скачкообразно увеличиваться до 0.9, что должно приводить к сближению астероидов с Землей и Венерой. Эффективное время сближения с этими планетами значительно меньше, чем с Марсом. Таким образом, на космогонически коротком интервале времени почти все астероиды, двигавшиеся вблизи резонанса 3/1 с Юпитером, были выброшены в область планет земной группы.
Многочисленные попытки использовать механизм, предложенный Висдомом, для объяснения происхождения люка Гекубы (резонанс 2/1) не дали желаемого увеличения эксцентриситета до 0.49, когда становятся возможны сближения астероидов с Марсом. В 1994 г. С.Ферраз-Мелло (Университет Сан-Пауло, Бразилия) опубликовал результаты исследования движения астероидов в резонансе 2/1 с Юпитером в рамках задачи четырех тел Солнце–Юпитер–Сатурн–астероид. Оказалось, что эволюция эксцентриситета идет очень быстро. За 5 млн. лет он увеличивается от 0.05–0.1 до 0.4—0.5, а за 12 млн. лет — до 0.7–0.9. В последнем случае астероид может сблизиться Землей.
Если в зоне резонансов 3/1 и 2/1 астероиды практически не наблюдаются, то в резонансе 3/2 обнаружена целая группа из 53 астероидов. Она носит название группы Гильды. Все астероиды этой группы имеют эксцентриситеты, превышающие 0.1, и малые наклоны орбит. В 1993 г. Ф.Франклин, М.Лекар и М.Морисон из Гарвардского астрофизического центра исследовали причины устойчивости движения астероидов в резонансе 3/2 с Юпитером. Они рассмотрели задачу четырех тел Солнце–Юпитер–Сатурн–астероид. Результаты моделирования показали, что время, в течение которого влияние резонансов ведет к распаду группы, превышает возраст Солнечной системы и составляет около 40 млрд. лет. Для астероидов с эксцентриситетами орбит меньшими 0.1 интервал времени, на котором происходит уход из группы, значительно короче и не превышает возраста Солнечной системы. Поэтому сейчас в группе Гильды нет астероидов с эксцентриситетами меньшими 0.1. В случае резонанса 3/2 мы встретились со случаем, когда движение объектов остается устойчивым на очень больших временах, но влияние резонансов в конце концов приведет к хаотизации их движения и последующему распаду группы.
Рассмотренные примеры показывают, что движение объектов по резонансным орбитам в поясе астероидов приводит к выбросу этих объектов в область планет земной группы. Динамическая эволюция астероидов, движущихся среди внутренних планет была рассмотрена в 1994 г. Ш.Фрошле (Обсерватория Лазурный берег, Франция). По результатам численного моделирования были сделаны следующие выводы. Под действием вековых резонансов с планетами земной группы эксцентриситеты орбит астероидов будут значительно увеличиваться. Если в процессе динамической эволюции астероид не столкнется с планетой, то после того, как перигелийное расстояние его орбиты станет меньше радиуса Солнца, он упадет на Солнце. Динамическое время жизни астероида среди планет земной группы очень короткое и для того, чтобы популяция астероидов, сближающихся с внутренними планетами, не исчезла, необходим постоянный приток новых объектов, выбрасываемых из главного пояса за счет резонансных возмущений.
В июне 1997 г. П.Вигерт (Университет Йорка, Канада) сообщил о том, что астероид 3753 Cruithne является компаньоном Земли. Этот объект относится к группе астероидов, сближающихся с Землей. Его движение относительно Земли можно представить в виде подковообразной орбиты, складывающейся из годичных витков спирали, по форме напоминающих фасоль. Путь астероида по подкове занимает 385 лет. Каждые 385 лет астероид сближается с Землей на минимальное расстояние и ее тяготение изменяет большую полуось его орбиты от 0.997 до 1.003 а.е. или обратно. Из-за большого эксцентриситета (0.5) и прецессии орбиты такое движение является неустойчивым. С большой вероятностью астероид 3753 столкнется с Венерой около 8000 г. Это является еще одним подтверждением того, что астероиды, сближающиеся с Землей, имеют короткую шкалу динамической эволюции.
В последнее время проявляется повышенный интерес к динамике астероидов типа троянцев, движущихся в окрестности треугольных точек либрации планет-гигантов. Треугольные точки либрации расположены на орбитах планет по обе стороны от текущего положения планеты и удалены от него на 60°
. В задаче трех тел Солнце–планета–астероид движение в окрестности треугольной точки либрации является устойчивым, если масса астероида достаточно мала по сравнению с массой планеты. В действительности же на тела, находящиеся в окрестности треугольных точек либрации планет, действуют возмущения и от других объектов. Как это влияет на устойчивость?В январе 1997 г. Г.Левисон (Юго-западный исследовательский институт, США), Е.Шумейкер (Ловелловская обсерватория, США) и К.Шумейкер (Северный университет Аризоны, США) опубликовали результаты исследований по динамической эволюции орбит астероидов-троянцев Юпитера. Оказалось, что время устойчивого существования астероидов в окрестности треугольной точки либрации существенно зависит от собственного эксцентриситета орбиты астероида и амплитуды либрационного движения. Чем меньше собственный эксцентриситет и амплитуда либрации, тем дольше астероид будет находиться в зоне устойчивости. Для 178 известных на момент исследования астероидов-троянцев Юпитера граница интервала устойчивого движения находится на уровне 100 млн. — 1 млрд. лет. Было проведено численное моделирование движения 36 астероидов-троянцев, находящихся вблизи указанной границы, на интервале 4 млрд. лет. В итоге 21 астероид с временем устойчивой эволюции менее 1 млрд. лет покинул группу астероидов-троянцев. Ушедшие астероиды перешли на орбиты, подобные орбитам комет семейства Юпитера. Дополнительное статистическое исследование показало, что среди наблюдаемых 150 комет семейства Юпитера, только одна могла быть ранее астероидом-троянцем. По результатам сравнения орбиты кометы Шумейкеров–Леви 9 с элементами орбит выброшенных троянцев сделан вывод, что эта комета не могла быть ранее астероидом-троянцем.
С.Миккола (Университет Турку, Финляндия) и К.Иннанен (Университет Йорка, Канада) в 1992 г., а М.Холман и Ж.Висдом (Массачусетский технологический институт, США) в 1993 г. исследовали движение астероидов типа троянцев в окрестностях треугольных точек либрации Сатурна, Урана и Нептуна. Троянцы Сатурна ушли из окрестностей точек либрации за время около 10 млн. лет. Это следствие сильных возмущений со стороны Юпитера. Орбиты троянцев Урана и Нептуна эволюционируют более медленно. Можно выделить области устойчивого движения (на рассмотренном интервале 20 млн. лет), в которых можно ожидать обнаружение новых объектов.
Таким образом, движение большинства астероидов-троянцев планет-гигантов можно считать устойчивым только на ограниченном интервале времени. Этот интервал изменяется от нескольких миллионов лет для троянцев Сатурна до нескольких миллиардов лет для троянцев Юпитера.
Для тел, движущихся за внешней границей пояса астероидов, разделение на кометы и астероиды достаточно условно. Один и тот же объект в афелии может наблюдаться как астероид, а в перигелии — как комета. Малые тела, движущиеся во внешней области Солнечной системы, вне зависимости от их наблюдательного статуса, связаны общей динамической эволюцией, в процессе которой возможен переход между различными группами объектов.
Большинство известных короткопериодических комет принадлежат семейству Юпитера. Их периоды не превышают 20 лет, а перигелии орбит лежат внутри орбиты Юпитера. Поэтому динамическая эволюция этих комет определяется их взаимодействием с Юпитером. В 1995 г. Танкреди (Уругвай) исследовал динамику комет семейства Юпитера. Было проведено численное интегрирование уравнений движения 145 комет этого семейства (всех известных на момент исследования) на интервале 200 тыс. лет. В итоге 5 комет были выброшены на гиперболические орбиты. Из оставшихся 140 комет для 137 время Ляпунова концентрировалось между 50 и 150 годами при среднем значении 60 лет. Короткие времена Ляпунова свидетельствуют, что движение комет происходит в зоне сильного хаоса, которая является следствием частых сближений комет с Юпитером. Результаты моделирования показывают, что после двух тесных сближений с Юпитером комета “забывает” свою прошлую траекторию.
Результаты, полученные М.Холманом и Ж.Висдомом в 1993 г. позволяют объяснить малое число наблюдаемых объектов группы Кентавра. Исследование эволюции орбит 3000 пробных частиц на первоначально круговых орбитах в области движения планет-гигантов показало, что большинство частиц покинуло область между орбитами Юпитера и Сатурна на интервале 10–100 тыс. лет, через 1 млн. лет все частицы покинули эту область. Причина такого поведения в частых сближениях частиц с Юпитером и Сатурном. Уход частиц и областей между орбитами Сатурна и Урана и орбитами Урана и Нептуна происходит значительно медленнее. Большинство частиц покинули эти области через 10 млн. лет. После 800 млн. лет все частицы между орбитами Сатурна и Урана были выброшены. Между орбитами Урана и Нептуна осталось всего 6 частиц в области треугольных точек либрации Нептуна. Итак, динамическое время жизни малых тел в области планет-гигантов намного меньше возраста Солнечной системы. Это объясняет почему группа Кентавра столь малочисленна.
М.Холман и Ж.Висдом рассмотрели также вопрос о происхождении короткопериодических комет и динамике объектов пояса Койпера. Было показано, что объекты с малыми наклонами и эксцентриситетами, имеющие большие полуоси меньшие 42 а.е., могут существенно увеличить эксцентриситет своей орбиты в результате сближений с Нептуном на интервале 10–100 млн. лет. После чего эти объекты могут вести себя либо как короткопериодические кометы, либо как астероиды группы Кентавра. Если начальные значения наклона и эксцентриситета орбиты объекта в поясе Койпера не являются малыми, то распределение областей устойчивого и неустойчивого движения становится весьма сложным. Но можно сделать вывод, что источниками короткопериодических комет являются тела, имеющие большие полуоси от 35 до 45 а.е.
В целом
происхождение и динамическую эволюцию
короткопериодических комет можно представить
следующим образом. Объекты из внутренней области
пояса Койпера под действием резонансных
возмущений от Нептуна увеличивают
эксцентриситеты своих орбит. Когда перицентры
орбит попадают в область движения
планет-гигантов, происходит сближение с одной из
планет и объект переходит на орбиту, целиком
лежащую в этой области, то есть становится
Кентавром. Но подобные орбиты неустойчивы и
объект за сравнительно короткое время покидает
эту область. Может оказаться, что в процессе
эволюции объект подойдет слишком близко к Солнцу
и будет наблюдаться как короткопериодическая
комета. Однако, как мы видели на примере комет
семейства Юпитера, движение объекта из-за частых
сближений с планетами-гигантами будет
хаотическим. То, что такие сближения происходят
достаточно часто, подтверждается наблюдениями.
Наиболее эффектное событие такого рода
произошло в июле 1992 г., когда после сближения с
Юпитером комета Шумейкеров–Леви 9
стала его спутником и через два года в июле
1994 г. вошла в атмосферу Юпитера. Менее
зрелищные, но не менее важные для небесной
механики, события происходили с кометами Вольфа,
Отерма 3 и многими другими.
В 1978 г. М.Харт, анализируя влияние астрономических факторов на возникновение жизни на Земле, ввел понятие зоны обитаемости в зависимости от удаления от центральной звезды. По его оценкам приближение Земли на 5% к Солнцу создаст условия, аналогичные венерианским, а удаление от Солнца на 5% приведет к глобальному оледенению планеты.
Через 15 лет
Ж.Ласкар (Бюро долгот, Франция)
показал, что к этим условиям необходимо добавить
требование стабилизации инсоляции — потока
солнечного тепла на единицу поверхности.
Количество тепла, передаваемого высоким широтам,
является важнейшим климатическим фактором,
который управляет ростом полярных шапок и
наступлением ледниковых периодов.
Численное моделирование поведения системы Земля–Луна на интервале 20 млн. лет в прошлое и 10 млн. лет в будущее показало, что наклон оси вращения Земли колеблется около среднего значения 23.3° с амплитудой 1.3°
. Это вызывает колебания инсоляции на 20% на широте 65° . Такие изменения потока солнечного тепла могут вызывать существенные климатические изменения, например, оледенения. Перераспределение воздушных и водных масс, рост или уменьшение ледникового покрова в высоких широтах ведут к изменению динамического сжатия Земли. А это ведет к изменению скорости прецессии оси вращения Земли. Напомним, что прецессия вызывается притяжением экваториального избытка массы Земли Луной и Солнцем и возмущениями гелиоцентрического движения Земли со стороны планет. В современную эпоху скорость прецессии составляет 50.4712'' в год. Частота возмущающей силы, вызываемой движением перицентров Юпитера и Сатурна и долготы восходящего узла Юпитера равна 50.3021'' в год.
Томсон в 1990 г. оценил изменение скорости прецессии за счет изменения динамического сжатия Земли при переходе от межледниковья к ледниковому периоду. Скорость прецессии уменьшается на 0.35'' в год. Это значительно больше разности между скоростью прецессии и возмущающей частотой. Следовательно, при смене климатических ситуаций Земля должна проходить через резонанс, вызываемый возмущениями от Юпитера и Сатурна. Для того чтобы наблюдались неоднократные переходы между оледенениями и межледниковьями, необходимо иметь устойчивые состояния, соответствующие этим климатическим ситуациям.
Ж.Ласкар, Ф.Жутел (Бюро долгот, Франция) и Ф.Будин (Обсерватория Лазурный берег, Франция) (Laskar J., Joutel F. and Boudin F. Orbital, precessional, and insolation quantities for the Earth from -20 Myr to +10 Myr. Astronomy and Astrophysics. 1993. 270. 522–533.) исследовали поведение частоты прецессии и наклона оси вращения Земли на интервале 18 млн. лет при значениях динамического сжатия Земли от 0.9940 до 1.0015 (современное значение принято за единицу). Результаты расчетов представлены на рисунках. Для наклона оси вращения Земли определялись среднее, минимальное и максимальное значения на рассматриваемом интервале. Точный резонанс достигается при динамическом сжатии 0.9977 от современного значения. Переход через резонанс в сторону уменьшения динамического сжатия приводит к уменьшению частоты прецессии с 50.5'' до 50.15'' в год. Максимальный наклон оси вращения увеличивается на 0.5° и достигает 25° .
На основании результатов моделирования можно предложить следующий сценарий наступления ледниковых периодов. В межледниковую эпоху за счет уменьшения ледяного покрова в высоких широтах уменьшается динамическое сжатие Земли. После прохождения через резонанс, резко увеличивается наклон эклиптики к экватору. Это ведет к уменьшению инсоляции в высоких широтах. Начинается оледенение. Рост ледяного покрова в высоких широтах ведет к увеличению динамического сжатия Земли. После прохождения через резонанс наклон эклиптики к экватору резко уменьшается. Увеличивается инсоляция в высоких широтах. Наступает межледниковая эпоха. После чего цикл повторяется снова. Таким образом, повторяемости ледниковых эпох на Земле мы обязаны совместному влиянию Луны, Солнца, Юпитера и Сатурна.
Как сильно влияет Луна на устойчивость оси вращения Земли? Оказывается, Луна — мощный стабилизатор скорости прецессии и угла наклона оси вращения Земли. Ж.Ласкар, Ф.Жутел и Ф.Будин рассмотрели эволюцию оси вращения Земли в отсутствии Луны. Частота прецессии уменьшилась до 15.6'' в год, а амплитуда изменения наклона оси вращения увеличилась: переход от 15° до 30° совершается менее, чем за 200 тыс. лет. Это должно вести к существенным изменениям инсоляции Земли и, как