Структура, динамика и устойчивость Солнечной системы

 

Введение
Современные представления о строении Солнечной системы
Эволюция взглядов на устойчивость Солнечной системы
Особенности динамической эволюции малых тел Солнечной системы
Роль Луны в эволюции жизни на Земле
Хаос в наклонах осей вращения планет
Эволюция орбит планет на очень больших интервалах времени

 

Введение

Последнее десятилетие принципиально изменило наши представления о строении, динамической эволюции и устойчивости Солнечной системы. Привычными стали сообщения об открытии новых объектов, выявлении новых динамических структур, проявлении свойств неустойчивости движения или хаотического поведения у тех или иных групп объектов.

Это вызвано несколькими причинами: появление новых инструментов и модернизация старых, применение высокочувствительных ПЗС–матриц и новых методов математической обработки результатов наблюдений. Все это позволяет наблюдать новые объекты, имеющие очень малую яркость и существенное собственное движение.

Новые аналитические и численные методы небесной механики в совокупности с современными вычислительными системами дают возможность моделировать движение тел Солнечной системы на интервалах времени, сравнимых с ее возрастом и даже многократно превышающих его.

На наших глазах происходит смена представлений о динамике Солнечной системы: от регулярной и устойчивой к хаотической и неустойчивой. Все это напоминает ситуацию в физике начала XX века, когда совершался переход от классической к релятивистской картине Мира. Нам предстоит разобраться где, когда и при каких условиях мы можем рассматривать Солнечную систему регулярной и устойчивой, а в каких случаях проявляются признаки хаоса и неустойчивости.

Начнем рассмотрение с современных представлений о структуре Солнечной системы. Затем обсудим понятия устойчивости и неустойчивости движения, условия возникновения резонансов и хаотического поведения. После этого проанализируем динамику малых тел Солнечной системы и обратимся к большим планетам. В заключение рассмотрим динамику Солнечной системы как целого на временах, сравнимых с ее возрастом.

ВведениеСтруктура, динамика и устойчивость Солнечной системыСовременные представления о строении Солнечной системы

| Начало | "Движение тел Солнечной системы" | Закономерности движения | Структура, динамика и устойчивость |

 

Современные представления о строении Солнечной системы

Все объекты Солнечной системы можно разделить на четыре группы: Солнце, большие планеты, спутники планет и малые тела. Мы пока ничего не говорим о спутниках малых тел, поскольку к настоящему времени таких объектов открыто всего два, а наблюдательной информации недостаточно, чтобы детально исследовать их динамику.

Солнце — динамический центр системы. Его гравитационное влияние является доминирующим в Солнечной системе за исключением малых областей в окрестности других объектов.

nineplan.jpg (30000 bytes)Большие планеты — визитная карточка Солнечной системы. Пять ближайших к Земле больших планет были известны с ранней истории человечества. Это — Меркурий, Венера, Марс, Юпитер и Сатурн. История открытия трех других больших планет показывает как менялось отношение астрономов к вопросу о размерах и населении Солнечной системы.

Открытие Урана явилось сюрпризом. Весной 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 часов на блинк-компараторе.

tour_v12.gif (11833 bytes)Существуют ли большие планеты за орбитой Плутона? Анализ траекторий движения тел Солнечной системы и космических аппаратов Пионер–10, Пионер–11, Вояджер–1, Воджер–2 позволяют утверждать, что объектов, сравнимых с Плутоном, и более крупных во внешней области Солнечной системы не существует.

История открытия спутников планет не менее драматична, но мы не будем на ней останавливаться. Отметим только, что спутниковые системы планет-гигантов сложностью своего устройства зачастую превосходят Солнечную систему. Не до конца решен вопрос о происхождении двойных планет Земля–Луна и Плутон–Харон.

Малые тела Солнечной системы — пробный камень и золотая жила небесной механики, кладезь новых открытий. Самые известные малые тела — кометы. Упоминания о кометах можно найти в легендах и летописях практически всех народов Земли. По динамическим признакам кометы разделяются на долгопериодические и короткопериодические.


west2.gif (5632 bytes)Долгопериодические кометы движутся по орбитам, большие полуоси которых достигают десятков тысяч астрономических единиц, а периоды обращения — десятков миллионов лет. Орбиты сильно вытянуты, их эксцентриситеты близки к единице. Ориентация орбит и их наклоны к плоскости эклиптики распределены случайным образом. В настоящее время имеются сведения более, чем о 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 веке.

sc93_mov.gif (4608 bytes)После открытия Плутона неоднократно предпринимались попытки поиска десятой большой планеты Солнечной системы. Во время одного из таких обзоров 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. Сначала делалось предположение о фундаментальном свойстве Мира — устойчивость, хаотичность и т.п. — затем создавалась система. По свой сути эти системы были кинематическими.

newton.jpg (12208 bytes)Ньютон первым построил динамическую модель Солнечной системы и сразу же столкнулся с вопросом о ее устойчивости. Он вышел из этого затруднения с помощью Великого Часовщика, который время от времени должен возвращать планеты на их орбиты.

В дальнейшем понятие устойчивости развивалось параллельно с исследованиями движения планет. lagrange.jpg (32862 bytes)Лагранж считал движение устойчивым, если оно происходит в замкнутой области пространства. Согласно теореме Лапласа–Лагранжа (1773, 1776 гг.) об отсутствии вековых возмущений больших полуосей планетных орбит, их изменение с точностью до величин первого порядка малости относительно возмущающих масс можно представить в виде суммы тригонометрических слагаемых. laplace.jpg (39225 bytes)На основе этой теоремы в 1773 г. Лаплас сформулировал теорему об устойчивости Солнечной системы: если движение планет происходит в одном направлении, их массы одного порядка, эксцентриситеты и наклоны малы, а большие полуоси испытывают лишь небольшие колебания относительно среднего положения, то эксцентриситеты и наклоны орбит будут оставаться малыми на рассматриваемом интервале. Вывод обнадеживающий. Однако, в настоящее время эта теорема имеет только исторический интерес. Она не применима на интервалах времени, сравнимых с возрастом Солнечной системы, поскольку учитывает возмущения только первого порядка. Кроме того, массы тел Солнечной системы различаются существенно.

poisson.jpg (19091 bytes)Позднее, в 1809 г. Пуассон показал, что во втором приближении решение для большой полуоси может содержать тригонометрические слагаемые с амплитудой, пропорциональной времени. Соответствующее определение устойчивости по Пуассону требует, чтобы частица за бесконечное время бесконечное число раз проходила через сколь угодно малую окрестность начальной точки.

Наиболее удачное понятие устойчивости сформулировал в конце XIX века русский математик А.М.Ляпунов. Исследуемое движение считается устойчивым, если все возможные движения, мало отличающиеся от него в начальный момент, в последующем будут мало отклоняться от него на всем интересующем интервале времени. Если же найдется хотя бы одно (!) движение, в начальный момент мало отличающееся от исследуемого, которое постепенно, пусть и через большой промежуток времени, заметно отклонится от него, то исследуемое движение — неустойчиво. Это определение считается основным по сей день.

В задачах небесной механики рассматривается устойчивость по части переменных: большой полуоси (задает размер орбиты), эксцентриситету (определяет вытянутость орбиты) и наклону орбиты. Совместная эволюция этих трех элементов должна удовлетворять определению устойчивости по Ляпунову. Отсюда, Солнечная система устойчива по Ляпунову, если размеры, форма и наклоны орбит остаются близкими к начальным на всем рассматриваемом интервале времени.

Говоря об устойчивости Солнечной системы, как правило, имеют ввиду устойчивость движения больших планет на бесконечном или очень большом, сравнимом с ее возрастом, интервале времени. В этом случае крайними проявлениями неустойчивости являются уход из Солнечной системы, падение на Солнце или столкновение с другой планетой. Такое событие способно существенно изменить структуру и динамику Солнечной системы.

В последнее время появилось большое число исследований по динамике малых тел Солнечной системы. Для этих объектов неустойчивость чаще всего проявляется в существенном изменении размеров, формы и наклона орбиты, зачастую многократном, и только после этого происходит уход из Солнечной системы, падение на Солнце или столкновение с планетой. Благодаря тому, что массы этих объектов незначительны, а их число велико, неустойчивое движение одного объекта, или даже группы, почти не влияет на устойчивость Солнечной системы в целом, хотя и играет важную роль в формировании ее структуры.

poincare.jpg (30612 bytes)Классический метод исследования движения небесных тел заключается в представлении решения соответствующих уравнений возмущенного движения в виде отрезков рядов. Однако, в конце 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 а.е.

komet_2.jpg (2698 bytes)В целом происхождение и динамическую эволюцию короткопериодических комет можно представить следующим образом. Объекты из внутренней области пояса Койпера под действием резонансных возмущений от Нептуна увеличивают эксцентриситеты своих орбит. Когда перицентры орбит попадают в область движения планет-гигантов, происходит сближение с одной из планет и объект переходит на орбиту, целиком лежащую в этой области, то есть становится Кентавром. Но подобные орбиты неустойчивы и объект за сравнительно короткое время покидает эту область. Может оказаться, что в процессе эволюции объект подойдет слишком близко к Солнцу и будет наблюдаться как короткопериодическая комета. Однако, как мы видели на примере комет семейства Юпитера, движение объекта из-за частых сближений с планетами-гигантами будет хаотическим. То, что такие сближения происходят достаточно часто, подтверждается наблюдениями. Наиболее эффектное событие такого рода произошло в июле 1992 г., когда после сближения с Юпитером комета Шумейкеров–Леви 9 стала его спутником и через два года в июле 1994 г. вошла в атмосферу Юпитера. Менее зрелищные, но не менее важные для небесной механики, события происходили с кометами Вольфа, Отерма 3 и многими другими.

особенности динамической эволюции малых тел Солнечной системыСтруктура, динамика и устойчивость Солнечной системыРоль Луны в эволюции жизни на Земле

| Начало | "Движение тел Солнечной системы" | Закономерности движения | Структура, динамика и устойчивость |

 

Роль Луны в эволюции жизни на Земле

В 1978 г. М.Харт, анализируя влияние астрономических факторов на возникновение жизни на Земле, ввел понятие зоны обитаемости в зависимости от удаления от центральной звезды. По его оценкам приближение Земли на 5% к Солнцу создаст условия, аналогичные венерианским, а удаление от Солнца на 5% приведет к глобальному оледенению планеты.

earth5.jpg (9090 bytes)Через 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 тыс. лет. Это должно вести к существенным изменениям инсоляции Земли и, как