Астероидно-кометная опасность: вчера, сегодня, завтра

В методе ТГС исследуемая траектория аппроксимируется последовательностью кеплеровых гелиоцентрических орбит соударения с планетами. При «соударении» происходит мгновенное преобразование скорости — поворот вектора планетоцентрической скорости на угол между асимптотами планетоцентрической гиперболы. Модуль планетоцентрической скорости остается при этом постоянным. Элементы кеплеровой орбиты преобразуются в этом случае по формулам классической задачи двух тел. Как известно, по координатам и скоростям в данный момент времени элементы орбиты определяются однозначно. Угол поворота вектора планетоцентрической скорости зависит от прицельного расстояния и в рамках рассматриваемой аппроксимации (сжатия гравитационной сферы в точку, обнуления прицельного расстояния) не определен, т. е. может быть произвольным в некоторых пределах.

Для исследования траекторий, содержащих опасные сближения и соударения астероида Апофис с Землей сразу после возможного тесного сближения в 2036 г., в работе [Соколов и др., 2008] используются следующие методы. В предположении, что движение в промежутке между 2029 и 2036 гг. происходит по резонансной орбите соударения с соизмеримостью средних движений 6:7, методом ТГС строится множество возможных резонансных кеплеровых орбит соударения после тесного сближения в 2036 г. Затем с использованием численного интегрирования уравнений движения с учетом всех возмущений ищутся траектории тесных сближений и соударений с Землей, близкие к построенным кеплеровым резонансным орбитам соударения. Для преодоления основной трудности — потери точности из-за тесных сближений, область допустимых начальных данных можно транспортировать вдоль траекторий в 2035 г., в результате чего ее размеры увеличиваются на несколько порядков. После этого численное нахождение начальных данных, соответствующих опасным траекториям, не представляет принципиальных трудностей. Полученные методом ТГС приближенные значения минимальных геоцентрических расстояний используются при нахождении интересующих нас опасных сближений в качестве первых приближений, начальные данные для численного интегрирования берутся из новой области в 2035 г. Аналогичный подход использовал Ж. Ласкар при исследовании динамического хаоса в Солнечной системе и поиске «уходящих» из нее траекторий Меркурия.

В результате с использованием методов, указанных в работе [Соколов и др., 2008], были построены порождающие эллиптические резонансные траектории соударений в 2037–2052 гг. (95 траекторий). Некоторые из них не могут быть реализованы, так как проходят в 2036 г. на геоцентрическом расстоянии менее радиуса Земли, однако большинство соответствует реальным траекториям. Подробнее были рассмотрены порождающие траектории с соударениями до 2041 г. включительно: численно получены траектории, содержащие соударения Апофиса с Землей в 2040 и 2041 гг. (резонансы 3:4 и 6:5 после 2036 г.), а также тесные сближения в 2037, 2038 и 2039 гг. Численное интегрирование подтвердило корректность использования метода ТГС для построения приближенных решений в случае Апофиса. Для численного интегрирования использовался интегратор Эверхарта. Движение планет и Луны описывалось известными динамическими моделями DE403 и DE405. Минимальные геоцентрические расстояния после 2036 г., полученные по DE403 и DE405, различаются незначительно, разница обычно менее 1 км. Области начальных значений координат в 2035 г., соответствующие указанным соударениям в 2040 и 2041 гг., имеют размеры порядка сотен и десятков метров соответственно. Отметим, что размеры областей начальных данных, соответствующие соударениям, по DE403 и DE405 практически совпадают, в то время как сами эти области отстоят друг от друга на десятки километров.

Представляет также интерес история сближений Земли и Апофиса. Проведенное численное интегрирование от 2006 г. до 1700 г. показало отсутствие тесных сближений на этом интервале. Все сближения происходят 12–14 апреля. Самое тесное сближение имело место в 1819 г. до расстояния 0,84 млн км. Остальные шесть сближений происходили до расстояний 3–4 млн км.

На примере астероида Апофис хорошо видны некоторые важные аспекты проблемы астероидно-кометной опасности, на которые раньше не обращали должного внимания. Так, соударению астероида с планетой могут предшествовать сближения с ней, хотя бы потому, что сближения на умеренные расстояния более вероятны, чем соударения. Этот вопрос обсуждается, в частности, в работе [Елькин, Соколов, 1995]. Отмеченное обстоятельство позволяет заблаговременно обнаружить потенциально опасный объект. С другой стороны, рассеяние возможных траекторий при тесных сближениях ведет к потере точности, затрудняет прогнозирование движения и требует применения специальных методов. Среди возможных движений астероида после тесных сближений появляются опасные траектории, имеющие резонансные возвраты к Земле. Соответствующие «замочные скважины», или зоны резонансного возврата, имеют очень малые размеры и, следовательно, мала вероятность столкновения с Землей. Точное положение «скважин» зависит от величины возмущающих сил: малые недостаточно известные эффекты могут изменить место «скважин» в пространстве начальных данных. Точное определение их положения является сложной актуальной задачей. В то же время само наличие «скважин» и их размеры мало зависят от возмущений. Следует отметить сложную структуру соответствующего множества, аналогичную фрактальной. Небольшое искусственное изменение траектории астероида (удар по нему и т. п.) с целью предотвратить соударение в апреле 2036 г. не гарантирует отсутствия соударений в последующие несколько лет. Отклонив траекторию, можно попасть в близлежащую «замочную скважину», ведущую к соударению. Недостаточная точность знания орбиты Апофиса не позволяет пока исключить возможность попадания в эти зоны резонансного возврата. Поэтому наряду с совершенствованием методов прогнозирования движения астероидов особую ценность имели бы наблюдения АСЗ из космоса, а также использование для уточнения орбиты Апофиса сигналов радиопередатчика, доставленного на орбиту искусственного спутника этого астероида.

7.7.6. Влияние эффекта Ярковского на движение Апофиса. Во всех рассмотренных случаях модель движения Апофиса не включала влияние эффекта Ярковского (см. раздел 3.6). Он вызывается неравномерным нагревом тела солнечными лучами в результате осевого вращения тела и его движения по орбите. Эффект Ярковского зависит от положения оси вращения, орбиты и массы тела, от теплопроводности его поверхностных слоев. Поскольку большинство этих параметров неизвестны, явным образом учесть эффект невозможно. Однако неявным образом этот эффект частично учитывается в результате подгонки элементов орбиты, прежде всего большой полуоси, к наблюдениям, которые, естественно, отражают влияние эффекта, если он достаточно велик. Тем не менее, значительная часть эффекта остается неучтенной, что может вызвать заметные ошибки при прогнозировании движения тела на основе найденной орбиты. В работе [Giorgini et al., 2008] выполнена оценка максимального смещения Апофиса по орбите за период от эпохи оскуляции (2006 г.) до 2029 г. под действием не учитываемой части эффекта. При предполагаемой массе в зависимости от направления вращения и гипотетически заданной теплопроводности поверхностных слоев максимальное смещение вдоль орбиты может составить к 2029 г. от –720 км до +780 км. Эти значения находятся в качественном согласии с оценками влияния эффекта Ярковского на движение ряда АСЗ, найденными Ю. А. Чернетенко иным путем по сравнению с работой [Giorgini et al., 2008].

В работе [Чернетенко, 2007] учет эффекта Ярковского в движении астероидов производился без каких-либо предположений о физических характеристиках этих тел. Предполагалось лишь, что зависимость этого эффекта от гелиоцентрического расстояния r имеет вид 1/r², а его величина характеризуется в общем случае тремя составляющими ускорения: радиальной, трансверсальной и нормальной, которые определяются из наблюдений (радарных и оптических) совместно с параметрами орбиты. При этом для астероида (6489) Голевка включение в число определяемых по наблюдениям запаздывания параметров трансверсальной составляющей A₂ позволило уменьшить ошибку единицы веса с 2,3 мкс до 0,5 мкс. Величина A₂ оказалась равной (-2,00 ± 0,14) 10^{– 14} а.е./сут².

При оценке влияния эффекта Ярковского на движение астероида Апофис было принято, что зависимость ускорения от гелиоцентрического расстояния имеет вид 1/r², а для A₂ принимались некоторые возможные значения (+2, –2, +6, –6) 10^{– 14} а.е./сут². Значения минимального расстояния от Апофиса до Земли и смещения астероида вдоль орбиты на 13 апреля 2029 г., полученные в результате уточнения параметров движения и последующего интегрирования, приведены в табл. 7.6.

Таблица 7.6. Влияние принятых значений ускорения, вызываемого эффектом Ярковского, на величину минимального расстояния от Апофиса до Земли и смещение астероида вдоль орбиты на 13 апреля 2029 г.

Изменения смещения вдоль орбиты при A₂ = (+6,-6) 10^{– 14} а.е./сут² близки к оценкам, полученным в работе [Giorgini et al., 2008]. Конечно, подобные расчеты имеют только качественный характер, поскольку реальное значение A₂ из имеющихся наблюдений найти невозможно. Но метод может быть использован в будущем, если будет получено необходимое число достаточно точных наблюдений.

Оценим теперь, как влияет эффект Ярковского на положение эллипса рассеяния на плоскости цели и на вероятность столкновения Апофиса с Землей в 2036 г. Будем исходить из предположения, что величина A₂ для астероида равна -2 10^{– 14} а.е./сут², что приводит к сокращению его полуоси и к увеличению среднего движения. В результате к апрелю 2029 г. астероид сместится вдоль орбиты на величину +272 км (табл. 7.6). Таким образом, астероид, соответствующий номинальному решению, найденному с учетом дополнительного ускорения, пересечет плоскость цели несколько раньше, и его минимальное расстояние от Земли составит, как показано в табл. 7.6, 37 964 км по сравнению с расстоянием в 38 220,5 км, на котором астероид должен проследовать мимо Земли в 2029 г. согласно решению, полученному ИПА.

7.7.7. Перспективы уточнения орбиты Апофиса. Для планирования возможного изменения орбиты Апофиса в случае, если это изменение будет вести к неизбежному или весьма вероятному столкновению, важно иметь представление о том, каковы перспективы проведения дальнейших астрометрических и физических наблюдений данного астероида и какого уточнения модели движения астероида можно ожидать в предстоящий период. При рассмотрении данного вопроса важно иметь в виду, что активные действия по изменению орбиты Апофиса должны быть предприняты до 2029 г. Величина корректирующего одиночного импульса, прикладываемого к астероиду до апреля 2029 г., чтобы избегнуть столкновения в 2036 г., примерно на три порядка меньше необходимого импульса после сближения в 2029 г. [Ивашкин, Стихно, 2008]. Причина такого различия вполне понятна. В случае, если коррекция выполняется до 2029 г., необходимо лишь несколько изменить минимальное геоцентрическое расстояние астероида (изменить орбиту так, чтобы ее пересечение с плоскостью цели произошло вне ближайшей окрестности «замочной скважины» размером около 600 м). Остальную работу по изменению траектории выполнит гравитационное поле Земли. После 2029 г. гелиоцентрическую орбиту астероида придется смещать уже на величину порядка нескольких радиусов Земли.

Размеры астероида и особенности его орбиты таковы, что благоприятные условия для его наблюдения складываются лишь в определенные весьма непродолжительные периоды времени. На рисунках 7.6–7.8 показаны элонгация астероида Апофис от Солнца, его звездная величина и расстояние от Земли в период 2006–2024 гг. На протяжении 2009–2011 гг. видимая звездная величина астероида будет оставаться в пределах 19,5^m–22,1^m. Но относительно большой блеск имеет место при малых элонгациях от Солнца. Сомнительно, что при звездной величине, превышающей 21^m, астероид можно будет наблюдать в элонгациях, меньших 50°. Таким образом, ближайшее «окно» для оптических наблюдений откроется не ранее декабря 2011 г., когда элонгация превысит 50°, а звездная величина будет меньше 21^m.