С учетом усложнений, которые обсуждались выше, для начального представления статистики лунных кратеров удобно следовать логике Хартманна, который избрал процессом, достаточно быстро обновившим большие участки лунной поверхности, излияние морских базальтов на видимой стороне Луны. Согласно имеющимся данным, базальты излились на поверхность довольно быстро (в геологическом смысле) — в интервале времени от 3,5 до 2,8 млрд лет назад образовалось почти 60 % общей площади лунных базальтовых морей [Иванов, 2005a; Hiesinger et al., 2003; Hiesinger et al., 2000; Shoemaker and Wolfe, 1982]. Распределение по размером кратеров, наложенных на поверхность лунных морей, показано на рис. 9.1 в инкрементальном виде (а) и в R-представлении (б). Для удобства использования эти данные можно выразить в аналитическом виде как зависимости числа NH кратеров на 1 км2 площади в интервалах диаметров с отношением D2/D1 = √2 в виде [Hartmann, 2005; Ivanov et al., 2001]:
lgNH = −2,61 − 3,82 lgDL, 0, 3 <DL < 1,41 км, (9.1)
lgNH = −2,920 − 1,80 lgDL, 1,41 <DL < 64 км, (9.2)
lgNH = −2,198 − 2,20 lgDL, DL > 64 км, (9.3)
где для определения интервала диаметров использовано левое (меньшее) значение граничного диаметра DL.
Для кратеров менее 300 м в диаметре для получения производящей функции нужно использовать более молодые участки поверхности, еще не успевшие насытиться постоянно образующимися новыми кратерами. Обычно на Луне такие участки находятся на покровах выбросов и днищах больших кратеров. При использовании этих данных производящее распределение по размерам можно продлить для диаметров < 300 м. Чтобы формально это продолжение соответствовало распределению на лунных морях, его можно записать в виде
lgNH = −2,0 − 2,90 lgDL, 0,01 < DL < 0,125 км. (9.4)
С таким добавлением кривая N(D) может двигаться вверх и вниз вдоль вертикальной оси для более молодых и более древних участков лунной поверхности.
Как видно из рис. 9.1, в диапазоне от метровых до километровых кратеров зависимость числа кратеров от их диаметра имеет сложную форму, лишь кусочно соответствующую простым степенным соотношениям.
Рис. 9.1. Распределение по размерам кратеров на поверхности лунных базальтовых морей в инкрементальном виде (а) и в R-представлении (б). Штриховая прямая 1 на рисунке б показывает примерный уровень «эмпирического» насыщения поверхности кратерами [Hartmann, 1984]. Рисунок первоначально был опубликован автором в работе [Neukum et al., 2001], позднее опубликован в русском варианте [Иванов, 2005a]. Сегменты степенных зависимостей (прямые участки на линии 2) соответствуют формулам (9.1) — (9.3). На рисунке б показана также штриховая линия 3, соответствующая аналитической кривой Нойкума [Ivanov, 2001; Neukum et al., 2001]
Возраст поверхности, на которой накопились эти кратеры, можно в среднем считать равным 3,3–3,4 млрд лет. Чтобы описать накопление кратеров со временем, или, наоборот, по известному количеству кратеров оценить возраст поверхности, применяется следующий прием: в каком бы диапазоне диаметров ни были произведены измерения, через измеренные точки проводится кривая, соответствующая производящей функции (показанной на рис. 9.1) и находится значение, на котором пересекается вертикальная линия, соответствующая D = 1 км (даже если истинное число кратеров с D∼ 1 км находится в зоне насыщения). Это значение N(D > 1 км) в кумулятивном представлении или соответствующее ему значение NH(DL = 1 км) позволяют ввести временну́ю шкалу, показывающую, с какой скоростью накапливаются кратеры. Такая шкала была построена с помощью измерения радиологического возраста образцов, доставленных с Луны КА «Аполлон» и автоматическими станциями серии «Луна». Параллельно по снимкам высокого разрешения измерялось распределение по размерам кратеров вокруг точек отбора образцов. С помощью стандартного (производящего) распределения по размерам все измерения были приведены к значению N(D > 1 км), которое сравнивается со значениями радиологического возраста образцов. Если поток кратерообразующих тел постоянен, значение N(D > 1 км) должно линейно расти с увеличением возраста пород. На рис. 9.2 показаны результаты такой обработки данных. Данные по возрасту кратера Коперник признаются весьма сомнительными — на лунной поверхности был собран материал, лежащий в пределах яркого луча, исходящего из кратера. Однако выживание вещества на поверхности Луны в течение 1 млрд лет является маловероятным из-за потока микрометеоров [Stoffler et al., 2006].
Рис. 9.2. Накопленное количество кратеров на различных участках лунной поверхности в зависимости от ее возраста, определенного по возвращенным лунным образцам: а) логарифмический масштаб по оси ординат, б) линейный масштаб по оси ординат.
Сплошная кривая — аналитическое представление [Stoffler et al., 2006]
На рис. 9.2. видны следующие особенности, важные для анализа АКО (Луна и Земля находятся в едином потоке бомбардирующих тел): в течение последних 3 млрд лет поток кратерообразующих тел был примерно постоянным; ∼ 4 млрд лет назад поток кратерообразующих тел был в 100–500 раз выше современного. Потоки в течение первых 0,5 млрд лет (не оставившие известных кратерных записей) являются в настоящее время предметом активных исследований [Hartmann et al., 2002]. К настоящему времени показано, что первоначальная идея о постепенном уменьшении темпа бомбардировки, начиная со стадии аккумуляции планет, не выдерживает проверки [Bottke et al., 2007a]. Причины резкой активизации бомбардировки примерно 4 млрд лет назад не ясны. Предложена интересная гипотеза (так называемая модель Ниццы — по названию города, где в обсерватории местного университета работают основные авторы модели). Согласно этой модели, миграция орбит планет-гигантов привела примерно 4 млрд лет назад к попаданию в резонанс периодов обращения Юпитера и Сатурна вокруг Солнца. Резонансные явления «встряхнули» всю Солнечную систему, вызвав поток «новых» кратерообразующих тел. Затухание этого потока и наблюдается на временно́м ходе кривой N(D) [Gomes et al., 2005; Tsiganis et al., 2005].
Несмотря на относительно низкую точность данных, показанных на рис. 9.2, Нойкум [Neukum, 1983; Neukum et al., 2001] предложил часто используемое аналитическое выражение для изменения параметра N(> 1 км) от времени накопления кратеров T (в данном случае N измеряется как число кратеров диаметром более 1 км на площади 1 км2, а время — в млрд лет):
N(> 1 км) = 5,44 10-14[exp(6,93T) — 1] + 8,38 10-4 Т. (9.5)
Сплошные кривые на рис. 9.2 построены именно по этому уравнению.
Проверка общей картины для последних 100 млн лет может быть проведена по данным о малых (D < 100 м) кратерах Луны. Сравнивая их с современным потоком метеороидов на границе земной атмосферы [Brown et al., 2002; Halliday et al., 1996], можно проверить постоянство потока малых тел метровых размеров [Ivanov, 2006]. На рис. 9.3 показаны данные измерения распределения по размерам на наиболее молодых датированных участках лунной поверхности. Эти данные показывают, что немногочисленные точки измерений позволяют предположить (пока не доказано обратное), что поток кратерообразующих тел был примерно постоянен в последние 100 млн лет. Более того, этот поток примерно соответствует современному потоку болидов, фиксируемых в земной атмосфере [Ivanov, 2006].
Примерное постоянство потока малых тел на Землю и Луну при времени усреднения 0,1–1 млрд лет не означает, конечно, что не существует кратковременных (≪ 100 млн лет) вариаций потока.
Рис. 9.3. а) Кумулятивные зависимости N(> D) для малых (D < 100 м) кратеров, наложенных на покровы выбросов 4 лунных кратеров, возраст которых был определен по длительности экспозиции в потоке космических лучей (cosmic ray exposure = CRE) для возвращенных на Землю образцов: 1 — кратер Южный Лучевой; 2 — кратер Конус; 3 — кратер Северный Лучевой [Moore et al., 1980]; 4 — кратер Тихо [Koenig et al., 1977]. Для сравнения показаны данные для лунных морей, иллюстрирующие в данном случае положение кривой насыщения 5. Пунктирная линия показывает средний уровень N(D) для насыщенной кратерами поверхности N(> D) = 0,047D-1,83 [Hartmann, 1984], хорошо соответствующий самым маленьким из измеренных кратеров на площадке у кратера Конус. б) Cравнение возраста площадок измерений по хронологии Нойкума [Neukum et al., 2001] для кратеров с D > 10 м и зависимости N(> D) ∼ D-2,9 для меньших кратеров. Возраст космической экспозиции (CRE) определен в работах [Arvidson et al., 1975; Drozd et al., 1974]