Преамбула.

          В истории Земли прослеживается чередование холодных и теплых периодов. Несмотря на множество мнений, остаются неясными два вопроса: 1) когда эти периоды происходили? 2) какая причина таких долгопериодических изменений климата? В настоящей статье представлены ответы на эти вопросы.

          В результате гравитационного взаимодействия тел Солнечной системы орбита Земли изменяется. Она вращается в своей плоскости со средним периодом 147 тыс. лет. Эксцентриситет орбиты, который определяет отличие орбиты от окружности, колеблется с периодами 94.5, 413, 2310 тыс. лет. Плоскость орбиты изменяется в пространстве: перпендикуляр к ней вращается со средним периодом 68.7 тыс. лет. Вращение или прецессия происходит вокруг вектора кинетического момента Солнечной системы. Кроме того, этот перпендикуляр колеблется с периодами 97.4, 1160 и 2320 тыс. лет.

          Кроме орбитального движения вокруг Солнца, Земля совершает вращательное движение вокруг своей оси, которая проходит через Северный и Южный полюса. За счет вращения Земля растягивается в экваториальной плоскости. Поэтому окружающие тела: Солнце, Луна и планеты создают момент силы, в результате которого ось Земли вращается (прецессирует) и колеблется подобно оси юлы. Прецессия оси Земли происходит со средним периодом 25.74 тыс. лет. Прецессия происходит с колебаниями. Ось Земли также колеблется в перпендикулярном прецессии направлении. То есть, угол наклона оси Земли в пространстве изменяется. Он колеблется с разными периодами: полумесячными, полугодовыми, 18.6 тыс. лет и с периодами в десятки тысяч лет. Чем больше период колебания оси, тем больше амплитуда колебаний (для других параметров зависимость может быть другой). В итоге угол наклона ε плоскости экватора к плоскости орбиты Земли колеблется от 15° до 32° при современном его значении ε = 23.4°.

          Эти колебания оси Земли, а также колебания эксцентриситета орбиты и ее вращение в своей плоскости (вращение перигелия орбиты) приводят к колебаниям количества тепла, поступающего от Солнца на разные широты Земли. Колебания таковы, что могут быть эпохи, когда на широте 80° тепла летом больше, чем сейчас на экваторе. А в другие эпохи на широте 60° летом тепла может быть меньше, чем сейчас на полюсе. В статье показано, что эти колебания тепла (по-другому, колебания инсоляция) совпадают с колебаниями палеоклимата в Западной Сибири.

          Таким образом, день сменяется ночью за счет вращения Земли вокруг своей оси. На смену зимы приходит лето, потому что Земля движется по орбите, а ее ось вращения наклонена на угол ε к перпендикуляру ее орбиты. Из-за этого летом Северное полушарие лучше освещается Солнцем. А Южное полушарие, наоборот, в это время получает меньше Солнечного тепла и там – зима.

          Если угол ε приблизится к 32°, то высокие широты будут получать, как уже отмечалось выше, больше тепла, чем его получают экваториальные широты. А, если угол ε приблизится к 15°, то из-за малого количества тепла ледниковые щиты покроют высокие и средние широты.

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

          А в какие моменты наступают эти изменения? На интервале 200 тыс. лет, как показано в статье, существует 14 периодов изменения климата. Они названы инсоляционными периодами. Моменты наступления изменений климата не одинаковы и нерегулярны. Но они строго детерминированы, и мы их определяем с точностью до 6 сек к концу интервала 200 тыс. лет и до четырех дней к концу промежутка 100 млн. лет.

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

          В статье показано, что инсоляционные периоды совпадают с палеоклиматическими изменениями в Западной Сибири. Читателю и исследователю предоставляется возможность убедиться в этом по другим территориям: по Европейской части России, по Западной Европе, по Северной Америке? Весь арсенал средств, необходимых для этого, я предоставил в свободный доступ: в статье об этом говорится.

          Обоснование полученных результатов, их дальнейшее развитие, изменение инсоляции и палеоклимата за 1 млн. лет, а также эволюция инсоляции за 20 млн. лет представлено в монографии: Смульский И.И. Эволюция оси Земли и палеоклимата за 200 тысяч лет. Saarbrucken, Germany: "LAP Lambert Academic Publishing", 2016. 228 с. ISBN 978-3-659-95633-1.

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

          А изменение за более долгие периоды, в сотни миллионов и миллиарды лет, возможно, обусловлены другими причинами. Однако, если из этих изменений вычесть инсоляционные изменения, то достоверность установления причин таких сверхдолгих изменений повысится.

          Настоящая статья опубликована в журнале “Геология и Геофизика”, который известен и в нас и за рубежом. Рецензенты предпочли не скрывать свои имена. Им: академику РАН Николаю Леонтьевичу Добрецову и чл.-корр. РАН Михаилу Всеволодовичу Кабанову в статье я выразил свою благодарность за тщательный анализ статьи, полезные рекомендации и замечания. Рецензенты – известные ученые, а кто их не знает, то в Интернете найдет подробные сведения о них. Ссылаться на статью так:

Смульский И.И. Новые результаты по инсоляции Земли и их корреляция с палеоклиматом Западной Сибири в позднем плейстоцене // Геология и Геофизика, 2016, т. 57, № 7, с. 1393-1407. http://dx.doi.org/10.15372/GiG20160709.

-------------------------------------------------------------------------------------------------------------------

УДК 523.2 + 551.34 +551.324

НОВЫЕ РЕЗУЛЬТАТЫ ПО ИНСОЛЯЦИИ ЗЕМЛИ И ИХ КОРРЕЛЯЦИЯ С ПАЛЕОКЛИМАТОМ ЗАПАДНОЙ СИБИРИ В ПОЗДНЕМ ПЛЕЙСТОЦЕНЕ

Смульский И.И.

Институт криосферы Земли СО РАН, 625000, Тюмень, а/я 1230, Россия

---------------------------------------------------------------------------------------------------------------------

Аннотация

          Три проблемы, составляющие Астрономическую теорию палеоклимата, решены по-новому. Две из них: изменение орбитального движения Земли и ее инсоляция подтвердили результаты предшественников. По третьей проблеме, изменении вращательного движения Земли, получены колебания оси вращения Земли, амплитуда которых в 7-8 раз больше чем в прежних решениях. Они приводят к изменениям инсоляции, которые объясняют колебание палеоклимата. Рассмотрены изменение инсоляции и ее структура за 200 тыс. лет. Показано, что ключевые события позднего плейстоцена Западной Сибири, например, последние оледенения и потепление между ними, совпадают с экстремумами инсоляции. Приведены инсоляционные периоды изменения палеоклимата и дана их характеристика.

Ключевые слова: угол наклона, инсоляция, палеоклимат, плейстоцен, Западная Сибирь.

Key words: obliquity, insolation, paleoclimate, Pleistocene, Western Siberia.

1. Введение

           В начале 20ого века Милутин Миланкович [1939] создал Астрономическую теорию изменения климата. В ней на основании трех параметров: эксцентриситета e орбиты Земли, углового положения перигелия φ и угла наклона ε плоскости орбиты к плоскости экватора Земли рассчитывается инсоляция Земли на разных ее широтах. В астрономическую теорию климата Земли составными частями входят задачи об орбитальном движении тел, о вращательном движении Земли и задача об инсоляции Земли в зависимости от параметров ее орбитального и вращательного движения.

           Решения М. Миланковича последовательно повторяли несколько поколений исследователей [Вурком, 1958], [Шараф и Будникова, 1969], [Berger and Loutre, 1991], [Edvardsson и др., 2002], [Laskar и др., 2004]. Однако все они шли одним и тем же путем, который за столетия сложился в небесной механике. Уравнения орбитального и вращательного движений, начиная с их вывода, приспосабливались для решения приближенными аналитическими методами. Мы пошли другим путем. Мы не копировали уравнения предшественников, а сами их выводили из первоначальных основ [Smulsky J. и Smulsky Ya., 2012]. Во-вторых, при выводе стремились ввести минимальные упрощения [Smulsky, 2011], [Smulsky, 2012а]. И, в-третьих, задачи решались численными методами, при этом использовались самые высокоточные их варианты [Smulsky and Krotov, 2014], [Смульский, 2014] или создавались новые [Smulsky, 2012б]. По первой и третьей задачам наши независимые исследования подтвердили исследования предшественников [Мельников и Смульский, 2009], [Smulsky и Krotov, 2014]. А результаты вращательного движения – другие: амплитуда колебаний угла наклона ε экватора Земли к ее орбите в семь-восемь раз больше [Смульский, 2013а], [Смульский, 2014]. Эти колебания дают такие изменения инсоляции, которые могут объяснить имевшие место колебания климата в прошлом. Вначале рассмотрим изменения инсоляции по времени на северной широте 65°, а затем на других широтах.

2. Эволюция угла наклона и инсоляции на широте 65° северного полушария

          На рис. 1 линией 1 показана эволюция угла наклона ε за прошедшие 200 тыс. лет. Изменения угла ε по нашим решениям вначале совпадают с аппроксимацией данных наблюдения, а также до 2000 лет совпадают с результатами решений других авторов [Laskar и др., 2004], [Шараф Ш. Г. и Будникова, 1969]. Далее рассчитанный нами угол наклона ε начинает отличаться от результатов решений [Laskar и др., 2004], [Шараф Ш. Г. и Будникова, 1969]. По сравнению с расчетами по прежним теориям (линия 2) наши решения дают в 7-8 раз большие амплитуды колебаний. В начальный период несколько тысяч лет в прошлое, начиная с T = 0, угол наклона 1, как и угол наклона 2, увеличивается. Затем он начинает уменьшаться и достигает минимума, в то время как угол наклона 2 по прежним теориям находится в максимуме. На остальном интервале времени также максимумы и минимумы колебаний углов наклона 1 и 2 не совпадают. Однако более существенным являются величины этих экстремумов. По прежним теориям на этом интервале времени угол наклона изменяется от 22.21° до 24.43°. А по нашим решениям наклон плоскости экватора Земли к плоскости ее орбиты изменяется от 14.8° до 32.1°. Примерно такой же диапазон изменений угла ε получен при решении задачи за будущие 200 тыс. лет [Смульский, 2013а].

           Средний период колебаний угла наклона по прежним теориям равен 41.1 тыс. лет. Этот же период имеется в колебании инсоляции QS65N (линия 2). Как видно из новой зависимости для угла наклона ε (линия 1) характерный период колебаний имеет в 1.5 - 2 раза меньшую длительность.

           В астрономических теориях палеоклимата рассматривается инсоляция не за астрономические, а за равные калорические полугодия. Начало и конец летнего калорического полугодия определяется так, чтобы инсоляция за любой его день была больше инсоляции за любой день зимнего полугодия. Далее будем рассматривать инсоляцию на широте 65° северного полушария, которое обозначается буквой N. Мы рассчитали изменение инсоляции QS65N за 200 тыс. лет в прошлое как по определенным нами параметрам: e, ε и φ (линия 1 на рис. 1), так и по тем же параметрам, рассчитанных Дж. Ляскаром и др. [Laskar и др., 2004] (линия 2). Как видно из графиков, инсоляция QS65N за летнее калорическое полугодие на северной широте 65° по нашим решениям имеет амплитуду колебаний также в 7-8 раз большую, чем по прежним теориям. Моменты потепления и похолодания по нашим расчетам 1 и по прежним теориям 2 также не совпадают. Вначале, начиная с T = 0, как видно из QS65N на рис. 1, летняя инсоляция растет в течение 4-х - 5-и тысяч лет, затем начинает падать и достигает минимума 16 тысяч лет назад. После этого минимума идет потепление, которое завершается большим максимумом инсоляции 31 тыс. лет назад.

Рис. 1. Эволюция угла наклона ε, а также летних инсоляций QS65N и I за 200 тыс. лет в прошлое. Сравнение новых результатов 1 с результатами прежних теорий 2 на примере работы Дж. Ляскара и др. [2004]: ε – угол наклона в радианах экватора Земли к плоскости ее орбиты; QS65N – инсоляция в ГДж/м2 за летнее калорическое полугодие на северной широте 65°; I – инсоляция в эквивалентных широтах за летнее калорическое полугодие на северной широте 65°. В градусах приведены максимальные и минимальные значения угла ε. T – время в тыс. лет (kyr) от 30.12. 1949 г.

           Итак, по нашим решениям получены в 7-8 раз большие колебания инсоляции. Насколько они существенны? На этот вопрос может дать ответ представление инсоляции в эквивалентных широтах I, которая рассчитывается следующим образом. Если в эпоху T летняя инсоляция на широте φ была такая, как в современную эпоху на широте φ0, то инсоляция в эквивалентных широтах будет I = φ0. На рис. 1 представлена инсоляция I в эквивалентных широтах, рассчитанная нами для широты 65° северного полушария как по нашим данным (линия 1), так и по данным Дж. Ляскара и др. [Laskar и др., 2004] (линия 2). Начиная с T = 0, инсоляция I по нашим данным 1 с широты 65° уменьшается на несколько градусов, т.е. на широте 65° становится теплее. Затем после максимума величина I начинает увеличиваться до широт 80° и 90°. При T = -15.88 тыс. лет летняя инсоляция на широте 65° меньше, чем современная летняя инсоляция на полюсе, поэтому на графике I проведена горизонтальная линия. Таким образом, горизонтальная линия 12 - 19 тыс. лет назад означает, что инсоляция на широте 65° была меньше, чем сейчас на полюсе. Такое малое количество тепла могло привести к оледенению территорий на широте 65°.

           При дальнейшем изменении времени до T = -30 тыс. лет летняя инсоляция I в эквивалентных широтах достигает широт 50°, 40°, 30°, т.е. на широте 65° становится значительно больше солнечного тепла. Горизонтальная линия в эпоху 28 - 34 тыс. лет назад означает, что на широте 65° поступает больше тепла, чем сейчас на экваторе.

          Линией 2 приведена инсоляция I в эквивалентных широтах по прежним теориям. Как видно, на широте 65° летняя инсоляция I на рассмотренном отрезке времени в 50 тыс. лет изменяется в пределах от 60° до 70°. Маловероятно, что изменения количества тепла на широте 65° до таких величин, которые сейчас имеются на широтах 60° и 70°, могут привести к существенным потеплениям климата или его похолоданиям. Такие малые изменения инсоляции всегда вызывали сомнения [Большаков и Капица, 2011].

           Полученные в наших расчетах колебания инсоляции могут вызвать наблюдавшиеся изменения климата. Уменьшение летней инсоляции I 19 - 12 тыс. лет назад до значений меньших, чем на полюсе (линия 1 на рис. 1), могло привести к оледенению. А увеличение летней инсоляции I 34 - 28 тыс. лет назад до значений больших, чем в экваториальной области, могло привести к благоприятным условиям существования мамонтовой фауны.

3. Изменение инсоляции по широте

          Выше мы исследовали изменение инсоляции во времени на широте 65° северного полушария. Теперь рассмотрим изменение инсоляции по широте в отдельные моменты времени. На рис. 2 показано изменение по широте φ° летней Qs, зимней Qw и уменьшенной в два раза годовой инсоляции QT в эпоху T = -31.28 тыс. лет. В это время на широте 65° северного полушария была наибольшая за 200 тыс. лет летняя инсоляция Qs (см. рис. 1).

          Инсоляция за год QT (рис. 2) изменяется симметрично относительно экватора (φ = 0°). На полюсах она наименьшая QT = 7.43 ГДж/м2 – а на экваторе наибольшая QT = 12.98 ГДж/м2 Летняя инсоляция Qs в своем изменении по широте φ имеет несколько экстремумов. Вблизи экватора (φ = 0°) она имеет минимум, а на тропиках – максимальные значения. В эту эпоху тропики находятся на больших широтах, чем в современную эпоху. Широта тропиков φ = 32.1° равна углу наклона ε на рис.1 в эпоху T = -31.28 тыс. лет. Как видно из рис. 2, летняя инсоляция Qs в высоких широтах от 90° до 60° практически остается постоянной и она существенно больше, чем на экваторе. Зимняя инсоляция Qw имеет максимальное значение вблизи экватора, а на полюсах стремится к нулю.

Рис. 2. Распределение по широте φ° Земли удельного количества тепла в ГДж/м2 в самую жаркую в высоких широтах эпоху T = -31.28 тыс. лет: Qs – за летнее калорическое полугодие; Qw – за зимнее калорическое полугодие; QT – за весь год: на графике величина QT уменьшена в два раза; φ > 0 – северное полушарие; φ < 0 – южное полушарие.

          На рис. 3 представлено изменение этих же составляющих инсоляции в наиболее холодную за 200 тыс. лет эпоху T = -46.44 тыс. лет. Годовая инсоляция QT, как и в предыдущую теплую эпоху, изменяется симметрично относительно экватора: на полюсах: наименьшая – QT = 3.58 ГДж/м2, а на экваторе – наибольшая QT = 13.79 ГДж/м2. То есть, по сравнению с теплой эпохой, в рассматриваемую холодную эпоху годовая инсоляция на экваторе немного увеличилась, а в высоких широтах уменьшилась более чем в два раза.

          Летняя инсоляция Qs (рис. 3) также имеет максимумы на тропиках, широта которых согласно рис. 1 равна ε = 14.8°. Однако в высоких широтах летняя инсоляция становится значительно меньше, чем на экваторе, почти в два раза. Характер изменения зимней инсоляции Qw не изменился: она имеет максимальное значение вблизи экватора и нулевые значения на полюсах. Максимум зимней инсоляции Qw немного увеличился по сравнению с теплой эпохой. В эпоху T = -46.44 тыс. лет сменилось также положение максимума Qw и минимума Qs в экваториальной зоне: они находятся в северном полушарии. В эпоху T = -31.28 тыс. лет эти экстремумы находятся в южном полушарии.

Рис. 3. Распределение по широте φ° Земли удельного количества тепла в ГДж/м2 в самую холодную в высоких широтах эпоху T = -46.44 тыс. лет. Обозн. см рис. 2.

          На рис. 4 сопоставлены распределения летних инсоляций Qs по широте за три эпохи: современную (1), самую теплую (2) и самую холодную (3). Самые большие изменения летней инсоляции происходят на высоких широтах, например, на полюсах Qs изменяется больше чем в два раза. В экваториальной зоне летняя инсоляция изменяется мало. При этом на экваторе в самую жаркую эпоху 2 инсоляция меньше, нежели в самую холодную 3, т.е. лето в жаркую эпоху может быть немного холоднее, нежели в холодную эпоху.

Рис. 4. Сравнение летних инсоляций Qs за три эпохи: 1 – современная; 2 – самая жаркая; 3 – самая холодная.

          Из графиков на рис. 4 следует, что принятая М. Миланковичем широта 65° в качестве опорной для характеристики климата является оправданной. На этой широте летняя инсоляция при переходе от холодной эпохи Qs65 = 4.72 ГДж/м2 к теплой Qs65 = 7.43 ГДж/м2 изменяется в 1.57 раза.

Рис. 5. Сравнение зимних инсоляций Qw за три эпохи: 1 – современная; 2 – самая жаркая; 3 – самая холодная.

          Изменения зимних инсоляций Qw, как видно из рис. 5, происходит однотипно. На полюсах величины Qw = 0. В наибольшей мере зимние инсоляции изменяются в средних широтах, и в меньшей мере – на экваторе. При этом в более теплые эпохи зимняя инсоляция меньше, чем в холодные. Другими словами, в теплые эпохи зимы холоднее, а в холодные – теплее.

          Изменение годовой инсоляции QT, как видно из рис. 6, четко разбивается на две зоны: высокоширотную (φ > 45°), и низкоширотную (φ < 45°). В высокоширотной зоне в более теплые эпохи тепла QT поступает больше, чем в холодные эпохи. Например, на полюсах – в два раза больше. В низкоширотной области зависимость обратная: в более теплые эпохи тепла QT поступает меньше, чем в холодные. Однако величина этих изменений небольшая: она в четыре раза меньше, чем на полюсах.

          На рис. 6, как уже отмечалось, приведены графики для двух эпох: самой теплой и самой холодной за 200 тыс. прошедших лет. В будущие 200 тыс. лет [Смульский, 2013а] годовая инсоляция QT в самую теплую T = 93.6 тыс. лет и самую холодную T = 109.1 тыс. лет эпохи ведут себя аналогичным образом. Если две зависимости QT за эти эпохи добавить к трем, уже имеющимся на рис. 6, то все пять кривых QT будут пересекаются при φ = 45°. Поэтому отмеченные закономерности справедливы на всем исследованном интервале ±200 тыс. лет.

Рис. 6. Сравнение инсоляций QT за год при трех разных эпохах: 1 – современная; 2 – самая жаркая; 3 – самая холодная.

          Как видно из графиков на рис. 4 – рис. 6, отличие инсоляции южного и северного полушарий Земли в наибольшей мере проявляется для летней инсоляции (рис. 4). При этом в экваториальных широтах это различие более выражено, чем в высоких широтах. Графики по годовой инсоляции QT (рис. 6) симметричны по широте, т.е. одинаковое количество тепла за год поступает в северное и южное полушария. Поэтому вариации зимней инсоляции по полушариям будут обратные летним вариациям.

          Итак, от теплой эпохи к холодной уменьшение летней инсоляции Qs на широте 65° происходит в 1.57 раза. А на полюсе годовая инсоляция QT уменьшается еще сильнее, а именно – в 2 раза. В то же время в экваториальной зоне происходят обратные по знаку изменения, но величина их в 4 раза меньшая. При этом на широте φ = 45° годовая инсоляция QT практически не изменяется, т.е. изменения летней инсоляции Qs компенсируются изменениями зимней инсоляции Qw.

4. Эволюция инсоляции на других широтах

          На рис. 1 была рассмотрена эволюция летней инсоляции Qs на широте 65°. Теперь рассмотрим эволюцию годовой, летней и зимней инсоляций на других широтах. На рис. 7 сопоставлены изменения за 200 тыс. прошедших лет инсоляций QT, Qs и Qw на пяти широтах 80°, 65°, 45°, 25° и 0° северного полушария. На широте 65° летняя инсоляция Qs хорошо отражает изменение годовой QT инсоляции. А зимняя инсоляция Qw изменяется в противофазе с летней: при теплом лете зима немного холоднее. Таким образом, для широты 65° годовая QT и летняя Qs инсоляции изменяются софазно, а зимняя Qw – в противофазе первым двум.

Рис. 7. Эволюция инсоляции за прошедшие 200 тыс. лет на пяти разных широтах северного полушария: Qs – за летнее калорическое полугодие; Qw – за зимнее калорическое полугодие; QT – за весь год.

          Такой же характер изменения инсоляций наблюдается для широт 80° и 45°. Для широты 25° почти софазно годовой инсоляции QT изменяется зимняя Qw. Отличия заключаются в том, что короткопериодические колебания Qw более ярко выражены. Летняя инсоляция Qs изменяется почти в противофазе годовой QT на начальном отрезке времени. Затем противофазность нарушается. На широте 0° снова софазно годовой инсоляции QT изменяется летняя Qs и в некоторой степени зимняя Qw.

          Софазность и противофазность изменения свидетельствует, что на широте 45° и выше в теплую эпоху лета – более теплые, а зимы – более холодные. А на широтах φ < 45 в теплую эпоху зимы более теплые.

          Амплитуды колебаний инсоляций с широты 80° до широты 45° уменьшаются. А после широты 45° амплитуды снова увеличиваются. Количественный порядок изменения амплитуд отражается ценой деления на рис. 7 годовой инсоляции QTg по широте φ: QTg = 2.5; 1.5; 0.05; 0.4; 0.5 Гдж/м2 при φ = 80°; 65°; 45°; 25°; 0°, соответственно.

          Таким образом, отсюда следует, что в высоких широтах происходят наиболее существенные изменения инсоляции. Кроме того подтверждается прежний вывод о том, что на широте 45° годовая инсоляция практически не изменяется.

          На рис. 8 показаны изменения летней инсоляции I в эквивалентных широтах на восьми разных широтах: 80°; 65°; 45°; 25°;-25°; -45°; -65°; -8° земного шара. Напомним, что тонкой горизонтальной линией отмечена широта местности, к которой относится инсоляция I. Верхние горизонтальные участки графиков отмечают отрезок времени, когда инсоляция в эпоху T летом больше чем в экваториальной зоне в современную эпоху T = 0. А нижние горизонтальные, как уже отмечалось, участки отмечают эпохи, когда летняя инсоляция меньше чем в области полюсов.

Рис. 8. Эволюция инсоляции I в эквивалентных широтах северного полушария за 200 тыс. лет на восьми разных широтах Земли: горизонтальной прямой линией приведена широта, к которой относится изменение инсоляции I.

          На широте 65° (рис. 8) по горизонтальным участкам видно, что имеются две теплые эпохи при T = -31.28 и T = -110.8 тыс. лет и 5 холодных. На широтах 80° и 45° инсоляция I изменяется подобным образом. Отличие заключается в том, что на широте 80° длительность холодных периодов увеличивается и добавляется еще один при T = -75 тыс. лет. На широте 45° нет холодных периодов с инсоляцией меньшей, чем на полюсе. А теплые периоды с инсоляцией больше, чем на экваторе, становятся более продолжительными, и добавляется еще четыре теплых периода.

          На широте 25° продолжительность теплых периодов еще больше увеличивается, а по сравнению с широтой 45° добавляется еще три теплых периода. А в холодные периоды значение инсоляции I не опускается ниже широты 45°. То есть, на широте 25° не бывает лета холодней, чем сейчас лето на широте 45°.

          В южном полушарии на широте φ = -25° (рис. 8) амплитуды похолоданий имеют такой же порядок, как и на широте φ = 25° в северном полушарии. Однако моменты их наступлений, как правило, другие. Инсоляция в южном полушарии на широте φ = -65° примерно такая же, как в северном полушарии. Немного лишь отличаются величины в экстремальных точках. Соотношение между инсоляцией на широте φ = -65° с инсоляцией на широтах φ = -80° и φ = -45° остается такое же, как и в северном полушарии: на широте φ = -80° удлиняются холодные периоды, а на широте φ = -45° удлиняются теплые периоды.

          Инсоляция I на южной широте φ = -80° аналогична инсоляции в северном полушарии на той же широте. А на широте φ = -45° инсоляция в ряде эпох может существенно отличаться от инсоляции в северном полушарии, на той же широте.

          Итак, летняя инсоляция Qs65N на широте φ = 65 северного полушария адекватно представляет периоды похолодания и потепления на высоких широтах φ > 45° как в северном, так и в южном полушарии. При этом годовая инсоляция QT изменяется софазно летней Qs, а зимняя Qw – противофазно. Амплитуда колебаний зимней инсоляции в 3 с лишним раза меньше колебаний летней Qs. На широте 45° колебания инсоляции несущественны. А на широтах φ – 45° амплитуда колебаний инсоляции меньше, чем в высоких широтах. Кроме того в этой области имеются различия между изменениями инсоляции в северном и южном полушарии.

5. Инсоляция и последние оледенения в Западной Сибири

5.1. Самое значительное последнее оледенение

          В табл.1 приведены отдельные экстремумы инсоляции и их параметры. Потепления отмечены как максимумы (max), а похолодания – как минимумы (min). В результате изучения разными специалистами палеоклимата Западной Сибири выработано согласованное представление о его эволюции в позднем плейстоцене [Гросвальд, 2009; Svendsen et al, 1999]. Особенно совпадает их мнение о двух последних оледенениях. Хорошо согласуются с их мнением наибольшие похолодания, приведенные в табл.1, 15.88 тысяч лет назад (т.л.н.) и 46.44 т.л.н.

          Большая часть исследователей считает, что на севере Западной Сибири на шельфе Баренцева и Карского морей и на островах неоднократно образовывался ледниковый щит. По Гросвальду [2009], в ледниковые эпохи, которые многократно повторялись в последние 2.3 - 2.5 м.л., Баренцево-Карский шельф превращался в сплошную область эрозии (экзарации) и сноса. Этот щит объединялся со скандинавским щитом, а их ледосбросы на Печорский бассейн взаимодействовали между собой, периодически отклоняя друг друга. Имеются только различия относительно границ распространения ледниковых щитов и их положения в разные эпохи.

Таблица 1. Отдельные экстремумы инсоляции: max – потепления отмечены как максимумы; min – похолодания.

Tт.л.н.4.1615.8831.2846.44110.8
Типmaxminmaxminmax
Qs65N ГДж/м25.975.367.434.727.38

          В поздний Вейхелий, эквивалентный сартанскому оледенению [Svendsen et al, 1999], совместный ледниковый щит на востоке захватывал Новую Землю, на севере шел по шельфам северных морей, на западе достигал нижней трети Исландии и Британии, в Европе опускался до широт 56°, а в России поднимался на Север по Западной части Белого моря. Авторы работы [Svendsen et al, 1999] считают, что в Карском море могли быть островные оледенения с опорой на северную часть Таймыра. А в средне-раннем Вейхелии, который эквивалентен Ермаковскому оледенению, восточная часть совместного ледникового щита распространялась по северной части шельфа Карского моря и захватывала почти весь Таймырский полуостров и плато Путорана в Восточной Сибири, а по югу опускалась чуть ниже полярного круга и соединялась с южной границей ледника в Европе. В отличии от Svendsen et al [1999], у Гросвальда [2009] восточная граница этого ледника захватывает весь Таймыр и западную часть Среднесибирского плоскогорья, на юге опускается до широты 50°, а на западе захватывает всю Исландию и почти всю Британию. В связи с этим Гросвальд считает, что Баренцево-Карский ледниковый щит существовал в обе последние ледниковые эпохи.

          На определенном этапе образования ледникового щита в Западной Сибири сток рек бассейнов Оби и Енисея в ледовитый океан прекращается [Гросвальд, 2009]. Низменные части Западной Сибири начинают заполняться водой. Исследователи отмечают образование Хантыйского, Енисейского и Пуровского озер [Волков, Архипов, 1978; Волков и др., 1969], которые в последующем соединялись между собой и образовывали озеро-море Западной Сибири. С ростом высоты ледникового щита и сменой похолодания на потепление начинался сток льда с ледникового щита. Сток льда происходил как на север, в Ледовитый океан, так и на юг. Сток на север прослеживается по подводным желобам: Медвежинскому, Франц-Виктории, Св. Анны и Воронина. А сток на юг отмечен Сибирскими увалами в Западной Сибири и их продолжениям в Восточной Сибири. Эту границу Гросвальд назвал линией Волкова.

          С дальнейшим потеплением начинается таяние ледникового щита и происходит дальнейшее увеличение пресноводного моря до такой степени, что начинается его сток через Тоболо-тургайский желоб в Туранскую низменность Приаралья и, возможно, с проходами на юге – в Каспийскую низменность.

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

          В северной части Западной и Восточной Сибири обнаруживаются следы морских отложений. Они дали основания для утверждений о существовании периодов наступления моря (трансгрессий). Морские отложения на современной суше могут быть обусловлены разными причинами. Уровень суши мог опускаться ниже уровня моря в результате опускания всего массива, на котором находится ледниковый щит. Морские отложения могут быть выпаханы ледником со дна шельфа и перенесены им на материк. В теплый период возможно также таянием всех покровных льдов, включая ледники Антарктиды и Гренландии, в результате чего уровень океана может превысить в отдельных местах уровень суши. На них в последствие образуются морские отложения. Каждая из этих причин имеет свои особенности. Например, максимальное опускание суши будет происходить при наибольшей толщине ледника, а поднятие уровня моря – в эпоху его отсутствия, т.е. в период максимального потепления.

          Моренный пояс вдоль возвышенности Сибирские Увалы является самым южным [Архипов и др., 1980; Архипов, 2000]. Севернее, в зоне 65.5° и 67° с.ш., существует еще несколько моренных поясов. Вследствие неоднократности похолоданий моренные пояса могли возникать в разные периоды. Пути перемещения ледников разных эпох могли не совпадать. Поэтому морены могут принадлежать к разным ледниковым эпохам. Если последний ледниковый максимум (ПЛМ) был менее мощный чем предыдущие, то самый южный мореный пояс является результатом деятельности предыдущих ледниковых эпох. Самый южный пояс морен, доходящий до подножия Сибирских Увалов, относят к раннезырянскому времени [Архипов и др. 1980; Архипов, 2000], к Ермаковскому горизонту. К этому ледниковому периоду относится максимальная стадия Баренцево-Карского щита около 50 тыс. лет назад [Svendsen et al, 1999]. После максимальной стадии в Печорском Бассейне этот щит растаял к 40 т.л.н. Как видно из табл. 1, этому Ермаковскому ледниковому периоду соответствует минимум инсоляции 46.44 т.л.н. Инсоляция летом была Qs65 = 4.72 ГДж/м2 – самая низкая за весь период в 200 тыс. лет.

5.2. Последний ледниковый максимум

          Последнюю ледниковую эпоху связывают с сартанским горизонтом [Архипов и др. 1980; Архипов, 2000]. В этот период сформировался ледниковый рельеф Западной Сибири, в том числе краевые морены в зоне 65.5° - 67° с.ш. по южному подножию Салехардских увалов и Хадатейским грядам Тазовского полуострова. Севернее их расположены Ямало-Гыданские пояса морен на широте около 68° с.ш. от озёр Ярро-то на Ямале к востоку вдоль Гыданской гряды. Севернее расположены морены, которые являются самыми молодыми. Последние два гряда морен связаны с деградацией сартанского ледника. Как видно из табл. 1, эта эпоха 15.88 т.л.н. с инсоляцией Qs65 = 5.36 ГДж/м2 была теплее предыдущей на 12%. Поэтому моренами последней ледниковой эпохи не был достигнут южный моренный пояс.

          Одновременным сартанскому оледенению в Западной Сибири является последнее оледенение в Скандинавии, которое характеризуют горизонтом позднего Вейхелия [Svendsen et al, 1999]. В Архангельском регионе максимум оледенения относят к 17 т.л.н. и дегляциацию – к 16 т.л.н. Восточнее озера Онега дегляциация относится к 14.4 - 12.9 т.л.н. На Таймыре возраст ледника позднего Вейхелия или сартанского времени относится к 18 - 7.5 т.л.н. Пространство между Норвегией и Новой Землёй покрыто ледником – 10.7 т.л.н. По Архипову [1997] радиометрический возраст сартанского ледникового горизонта находится в пределах 23 - 10 т.л.н.

          После прохождения максимума похолодания крупнейшие северные подводные желоба: Франц-Виктория, Св. Анны и Воронина отводили ледяные потоки с Баренцево-Карского ледника в ледовитый океан [Гросвальд, 2009]. Колонки, поднятые с глубины 470 м желоба Франц-Виктории содержали материал морены, покрытый ледниково-морскими илами с возрастом 12.9 т.л.н. и голоценовыми отложениями –10 т.л.н. Деградация льда в желобе Св. Анны началась около 13 т.л.н. и завершилась к 10 т.л.н.

          Последствия образовавшегося в сартанское время Мансийского озера имеют примерно такой же возраст 20 - 10 т.л.н. [Архипов, 1997; Пятосина, 2005]. Причем озеро сартанского времени занимало меньшие площади, чем более древнее. Отложения более древнего озера перекрыты отложениями с остатками мамонтов, возраст которых примерно 18 т.л.н. [Пятосина, 2005].

          Возраст Колпашевской террасы на средней Оби на высоте 55 м., образованной Мансийским морем, находится в пределах 12.8 - 10.6 т.л.н. [Архипов, 1997]. К ней примыкает транзитная терраса-равнина [Арсланов и др., 1983], которая прослеживается через всю зону сартанского оледенения до устья Оби. Её возраст 12.26 – 0.17 т.л.н. На среднем Енисее терраса высотой 60 - 70 м у Фарково имеет возраст от 16.4 до 11.7 т.л.н.

          По Архипову [1997] грандиозный конус выноса в бассейн верхней Оби выходит из Алтая по долинам рек Бии, Катуни и др., проходит по Бийско-Барнаульской впадине до Новосибирска и далее протягивается в акваторию Мансийского палеоозера. Он формирует ряд террас, высота которых на Алтае достигает 220 - 240 м, у Новосибирска – 140 - 120 м., в Томском Приобье – 100 - 80 м. Далее конус выноса уходит по ложбинам Приобского плато в Кулундийскую равнину и в долину Иртыша. Возраст конуса выноса и его террас находится в пределах 17.6 - 10.4 т.л.н. [Бутвиловский, 1993; Панычев, 1979; Архипов, 1997].

          Как видим, приведенные датировки сартанского оледенения и его последствий совпадают с минимумом инсоляции 15.88 т.л.н.

5.3. Теплый период между оледенениями

          Между Ермаковским оледенением с минимумом инсоляции 46.44 т.л.н. и сартанским с минимумом инсоляции 15.88 т.л.н. находится Каргинский межледниковый горизонт с максимумом инсоляции 31.28 т.л.н. Во многих долинах Печорской низменности: Шапкиной, Хвостовой Созьвы и Соймы древесина и торф из-под рельефообразующих морен имеет возраст от 25 до 40 т.л.н. [Гросвальд, 2009]. Под распространенной к северу от Сибирских увалов мореной расположены озерно-болотные отложения, возрастом от 25 до 40 - 50 т.л.н. [Арсланов и др. 1983]. Раковины на восточном Баренцево-Карском побережье и берегах Таймыра и Северной Земли имеют возраст от 24 до 38 т.л.н. [Гросвальд, 2009].

          Как уже отмечалось [Svendsen et al, 1999], существовавший ранее Баренцево-Карский ледниковый щит к 40 т.л.н. полностью исчез. По Архипову [1997] казымская пачка Каргинского горизонта (дер. Казым-Мыс, правый берег нижней Оби) простирается по долине Оби до г. Колпашево и бассейна р. Васюган, а также до дер. Липовка на Тоболе. Ее возраст 33 - 31 т.л.н. На нижнем Енисее от г. Игарка до устья Бахты прослеживаются конощельские слои, аналог Казымскому слою с возрастом 33 - 32 т.л.н. Как правило, это озерно-аллювиальные отложения с прослойками торфа.

          Третья терраса Иртыша и Тобола с абсолютными отметками 70 - 75 м. у дер. Липовка сложена озерно-аллювиальными отложениями [Илларионов, 2013]. Возраст древесных и растительных остатков, имеющихся в ней, составляет от 31.78 до 32.77 т.л.н. Такой же возраст имеют костные остатки бизонов, шерстистых носорогов и лошадей.

          Аналогом Липовской террасы на р. Оби является кирьянская терраса (левый берег Оби в 15 км выше д. Покур). Она сложена озерно-аллювиальными отложениями с линзами торфа, возрастом от 27.5 т.л.н. до 36.3 т.л.н. [Лаухин и др, 2006]. В Тургайском желобе из скважины с абсолютным уровнем 37 м (25 км южнее оз. Кушмурун) обломок древесины имел возраст 27.8 т.л.н. [Гроссвальд, 1983]. В обнажении Красный Яр, в 15 - 20 км южнее Новосибирска, находятся укоренившиеся пни под озерной толщей слоем 8 - 10 м. Возраст пней 28 - 29 т.л.н. [Архипов и др., 1980]. В Кас-Кетском «проливе», между Енисейским и Мансийским древними озерами, под слоем озерно-аллювиальных отложений остатки древесины и торфа имеют возраст от 27.3 до 29.5 т.л.н. [Гросвальд, 2009].

          Таким образом, максимум инсоляции T = 31.28 т.л.н. в табл. 1 соответствует периоду Каргинского межледнековья.

6. Оптимум инсоляции в голоцене

          Как видно из табл. 1, в эпоху T = 4.16 т.л.н. имеется небольшой максимум инсоляции в голоцене (10 т.л.н. – 0). Ему предшествовал Сартанский ледниковый период, поэтому от него до этого максимума шло непрерывное потепление. Оптимум голоцена однозначно проявляется в интервале 9 – 3.3 т.л.н. [Васильчук, 1982; Ершов, 1989]. Исходя из полинологических данных потепление в нем было выражено более слабо, чем в предыдущее межледнековье [Ершов, 1989]. С конца ледникового периода до оптимума голоцена, т.е. 8 т.л.н. – 5 т.л.н., принято считать, что происходила трансгрессия моря, с 5 т.л.н. и по настоящее время море отступает [Ломанченков, 1966]. В связи с этим, формирование новой современной террасы и пойменных террас в долинах рек происходят после оптимума голоцена [Ломанченков, 1966; Сакс, 1953]. Возраст древесины из пойменных отложений Индигирки – 4.125 т.л.н., а с отложений ее притока, р. Большой Эрги – 4.770 т.л.н.

          В процессе анализа данных по Западной Сибири было установлено [Баулин, 1959; Некрасов, 1990], что верхняя толща многолетнемерзлых пород (ММП) сформировалась за последние 5 - 6 т.л.н. Формирование бугристых торфяников в южной зоне распространения ММП Западной Сибири началось примерно 3 т.л.н. [Шполянская и Евсеев, 1970; Некрасов и др. 1990].

          Кроме приведенных свидетельств существует множество других о наличии потепления, которое совпадает по времени T = 4.16 т.л.н. с максимумом летней инсоляции на широте 65°N (см. табл. 1). Чтобы облегчить специалистам интерпретацию палеоданных, рассмотрим более детально изменение инсоляции за последние 5 тыс. лет. На рис. 9 представлены изменения инсоляций QT, Qs и Qw на пяти широтах 80°, 65°, 45°, 25° и 0° Северного полушария. Эти данные получены с интервалов 1 год. Для широты 65°N видно, что годовая QT и летняя Qs инсоляции изменяются однотипно с наличием максимума. Однако максимумы их наступают в разное время: летней Qs – 4.16 т.л.н., а годовой QT – 3 т.л.н. В то же время зимняя инсоляция Qw на всем интервале 5 тысяч лет монотонно увеличивается.

          На широте 80° максимумы QT и Qs смещены ближе к современной эпохе, а в зимней инсоляции Qw в эпоху T = 0.499 т.л.н. появляется максимум. На широте 45° по сравнению с широтой 65° летняя инсоляция Qs имеет минимум в современную эпоху T = 0, а зимняя Qw монотонно увеличивается и в недалекое время в будущем достигнет максимума.

          На широте 25° инсоляция за год QT имеет ярко выраженный минимум при T = 2.8 т.л.н., а летняя достигает минимума при T = 0.7 т.л.н. Зимняя инсоляция Qw имеет максимум примерно при T = 0.5 т.л.н. На экваторе годовая QT и зимняя Qw инсоляции минимумы имеют примерно при T = 2.5 – 2.8 т.л.н. А летняя инсоляция Qs имеет минимум (T = 3 т.л.н.) и максимум (T = 1.5 т.л.н.).

          Амплитуды колебаний инсоляций на этом интервале, также как на интервале 200 т.л.н., с широты 80° до широты 45° уменьшаются. Качественный порядок величин изменения ΔQT годовой инсоляции по широтам φ = 80°; 65°; 45°; 25° и 0° имеют такой вид: Δ QT = 0.06; 0.023; 0.0015; 0.008 и 0.012 ГДж/м2, соответственно.

Рис. 9. Эволюция инсоляции за прошедшие 5 тыс. лет на пяти разных широтах северного полушария: Qs – за летнее калорическое полугодие; Qw – за зимнее калорическое полугодие; QT – за весь год; интервал между точками по времени – один год.

          Как было показано в работе [Smulsky и Krotov, 2014], существуют колебания инсоляции с периодами: полумесячным, полугодовым, 2.75 лет, 3.58 лет, 11.86 лет и 18.6 лет. Амплитуды колебаний возрастают с увеличением их периода. Наибольший период 18.6 года обусловлен воздействием прецессионного движения орбиты Луны на вращательное движение Земли. Для широты φ = 80°N амплитуда колебаний годовой инсоляции QT с периодом 18.6 лет составляет 532 кДж/м2, а с периодом 0.5 года – 31 кДж/м2. Амплитуды этих колебаний по сравнению с величиной годовой инсоляции QT порядка 5.8•106 кДж/м2 – малы и на графиках не видны. На рис. 9 наличие колебаний с периодом 18.6 лет наблюдается на графиках Qs и Qw для φ = 0°. Как уже отмечалось, данные на рис. 9 получены с интервалом 1 год, поэтому колебания инсоляции с периодом 18.6 лет содержатся в них с хорошим разрешением по времени.

          Изменение инсоляции на других широтах, в Южном полушарии, а также другие ее составляющие могут быть получены любым исследователем с помощью программы Insl2bd.mcd, которая имеется в свободном доступе [Смульский, 2013б]. Теория расчета инсоляции содержится в работе [Smulsky и Krotov, 2014], а параметры орбитального и вращательного движений Земли за 5 т.л.н. представлены в файле OrAl-5kyr.prn, а за 200 т.л.н. – в файле OrAl-200ky.prn.

7. Инсоляционные периоды изменения климата

          Экстремумы инсоляции за 50 т.л.н., как показано выше, хорошо согласуются с колебаниями палеоклимата. На рис. 10 пронумерованы все экстремумы инсоляции, за исключением первого, числами от 1I до 12I. Первый экстремум в эпоху T = 4.16 т.л.н., связанный с оптимумом голоцена, обозначен буквой OI.

          Как было показано выше эволюция летней инсоляция отличается от эволюции зимней и инсоляции за год. По-разному также изменяется инсоляция по широтам Земли. Тем не менее, летняя инсоляция Qs65N на широте 65° северного полушария адекватно отражает колебания инсоляции на высоких широтах φ > 45°, как в северном, так и в южном полушарии. Именно на этих широтах наблюдаются существенные изменения палеоклимата. Поэтому колебания инсоляции Qs65N были приняты за основу для установления инсоляционных периодов и их границ. При этом следует отметить, что при выборе другой характеристики инсоляции, например QT65N или Qs80N моменты наступления экстремумов и их границы будут немного изменены. Об этом необходимо помнить при интерпретации палеособытий, которые в большей мере определяются не летней инсоляцией, а годовой QT или зимней Qw, или тех событий, для которых их широта является определяющим фактором. В этих случаях возможна ситуация, что событие обусловлено потеплением, а его необходимо относить к более холодному периоду на широте 65° северного полушария, или наоборот.

          Итак, инсоляционные периоды будут устанавливаться по характеристике инсоляции Qs65N. Чтобы установить границы периодов, необходимо выбрать среднее значение Qs65N. Тогда по отклонению от среднего можно определять начало и конец периода. Для качественной характеристики инсоляционного периода необходимо выбрать количество градаций величины инсоляции. Для этого воспользуемся градацией климата в работе Фотиева [2009] из трёх уровней для холодного и трёх уровней для тёплого периодов, например, для холодного климата: умеренно холодный (у.х.), холодный (х.) и очень холодный (о.х.). При выборе среднего значения Qs65N мы руководствовались не математическим средним, а максимальным и минимальным значениями за 20 млн. лет: Qsmax = 7.581 ГДж/м2 при T = -9.2475 млн. лет и Qsmin = 4.505 ГДж/м2 при T = -18.95 млн. лет. Тогда среднее значение инсоляции можно записать так:

Qsm = 0.5·(Qsmax + Qsmin) = 6.043 ГДж/м2.
(1)

          При шести уровнях инсоляции величина ступени между уровнями будет:

ΔQs = (QsmaxQsmin)/6 = 0.5127 ГДж/м2.
(2)

          Тогда границы уровней холодного климата запишутся так:

Qsc,k = Qsmk·ΔQs, где k = 1, 2, 3,
(3)

          а границы уровней теплого климата запишутся аналогично:

Qst,k = Qsm + k·ΔQs, где k = 1, 2, 3.
(4)

          Следует отметить, что при выборе Qsm по экстремальным инсоляциям на интервале 200 тыс. лет, его значение равно 6.075 ГДж/м2, т.е. близкое к величине (1). Для современной эпохи T = 0 значение инсоляции Qs,065N = 5.922 ГДж/м2, т.е. меньше среднего значения 6.043 ГДж/м2. Таким образом, современная эпоха является немного холодней по сравнению со средним уровнем инсоляции. В связи с этим для оптимума голоцена начальная эпоха Tb определена по пересечению инсоляционной кривой со значением Qs,065N в современную эпоху. Для остальных инсоляционных периодов их границы начала Tb и конца Tf определены по пересечению инсоляционной кривой со средней инсоляцией Qsm. На рис. 10 средний уровень инсоляции Qsm указан позицией 1, границы первого и второго теплых уровней отмечены 1t и 2t, а границы холодных уровней – 1c и 2c. В таблице 2 приведены границы периодов: Tb и Tf , а также величины инсоляции Qs,ex65N в экстремальные эпохи Tex.

          Следует отметить, что по аналогичному алгоритму были введены уровни загрязнения атмосферы города [Смульский, 1987]. Они хорошо согласовались с уровнями, которые опытно-чувственным путем специалисты разработали при ежедневной оценке качества атмосферы на протяжении многих лет.

8. Корреляция инсоляционных периодов с существующими классификациями палеоклимата

          Как уже отмечалось, главные экстремумы инсоляции 1I, 2I и 3I (см. рис. 10 и табл. 2) за 50 т.л.н. хорошо согласуются с последними двумя ледниковыми периодами: Сартанским и Ермаковским и Каргинским межледниковьем между ними. Каждый ледниковый период перекрывает следы предыдущего, поэтому обнаружение и датировка более ранних колебаний палеоклимата затруднена. Перед Ермаковским похолоданием отмечают Казанцевское межледнековье. По словам Архипова [1997] оно было одним из самых теплых в плейстоцене. В современной южно-таежной зоне наряду с хвойными породами произрастали липа, вяз, дуб и орешник. В лесостепной зоне формировались черноземные почвы [Волкова, 1991; Архипов и др.,1995]. Как видно из рис. 10, в течение 50 тысяч лет, начиная с экстремума 3I похолоданий не было. Были два небольших потепления 60.8 т.л.н. (4I) и 83.16 т.л.н. (6I). Это, по-видимому, подтверждается двумя слоями торфа, возрастом 65 т.л.н. и 80 т.л.н. в отложениях Белогорской возвышенности по правому берегу нижней Оби [Архипов, 1997]. Эти слои с включением торфа перекрывают морену возрастом 100 т.л.н. Как видно из рис. 10, в то время 96 т.л.н. был пик похолодания 7I. Это похолодание в литературе также относят к нижним слоям Ермаковского горизонта. В долине средней Оби, в разрезе Кирьяс, нижние озерно-ледниковые отложения Ермаковского горизонта подстилаются палеопочвами 120 т.л.н. [Архипов, 1997]. Эти палеопочвы согласуются с пиком потепления 8I, который наступил 110.8 т.л.н.

Рис. 10. Периоды инсоляции OI, 1I, 2I,... 12I за 200 т.л.н. и границы их смены: 1 – средняя инсоляция Qsm; 1t и 2t – первая и вторая границы теплых уровней; 1c и 2c – первая и вторая границы холодных уровней; у.т., т., о.т. – теплые уровни; у.х., х., о.х. – холодные уровни.

          Все периоды инсоляции 4I - 12I перед самым значительным похолодание 3I в настоящее время не могут быть надежно связаны с палеособытиями. С одной стороны, это обусловлено тем, что первыми двумя ледниковыми периодами следы более древних палеособытий во многих местах уничтожены. А с другой стороны, исследователи не имели перед собой руководящей линии, которая бы позволяла им привязать разнообразные палеособытия в отдаленных друг от друга местах к одной и той же причине их возникновения. Мы рассчитываем, что инсоляционная кривая, представленной на рис. 10, станет для палеоклиматологов таким инструментом, позволяющим более определенно согласовывать разрозненные палеособытия.

          В табл. 2 приведена корреляция инсоляционных периодов с некоторыми горизонтами верхнего и среднего плейстоцена Западной Сибири. Период OI соответствует оптимуму голоцена. Период 1I согласуется с Сартанским оледенением. Ему, согласно Svendsen et al (1999], Никифоровой и Певзнеру [1982] также соответствуют горизонты других регионов: верхний Вюрм, верхний Вейхелий, верхний Висконсин и Осташковский ледниковый период.

          Периоду 2I соответствует Каргинское потепление. С ним согласуется средний Вюрм, средний Вейхелий и средний Висконсин. Периоду 3I следует сопоставить позднюю холодную фазу Ермаковского похолодания. С ней согласуется нижний Вюрм, нижний Вейхелий, нижний Висконсин и Днепровский ледниковый период.

          С периодами 4I, 5I, 6I, по-видимому, необходимо связать позднюю фазу Казанцевского потепления. Возможно, ей соответствует Эем, Рисс-Вюрм и Сангамон. Период 7I можно сопоставить с ранней фазой Ермаковского похолодания, а фазы 8I, 9I, 10I – с ранней фазой Казанцевского потепления, которой, по-видимому, соответствуют Микулинское и Эемское межледниковья.

Таблица 2. Периоды инсоляции и их характеристика за последние 200 тыс. лет. T – время в тыс. лет от 30.12.1949 г. Tb – начало периода, Tf – конец периода, Tex – эпоха экстремума, Qs65N – инсоляция за летнее полугодие на широте 65° северного полушария; градации климата: х. – холодный; т. – теплый; у. – умеренно; о. – очень.

№ периодаИнтервал Tb – Tf т.л.н.Экстремумы инсоляцииГрадации климатаКорелляция с горизонтами Западной Сибири
Texт.л.н.Qs65N ГДж/м2
OI6.86-04.165.973у.хОптимум голоцена
1I22.08-6.8615.885.364х.Сартанский (I ледниковый)
2I39.5-22.0831.287.4316о.т.Каргинский
3I53.8-39.546.444.7174о.х.Поздняя фаза Ермаковского оледенения (II ледниковый)
4I69.1-53.860.86.929т.Поздняя фаза Казанцевского потепления
5I76.96-69.172.85.7946у.х.
6I88.52-76.9683.46.5197у.т.
7I102.56-88.5295.924.9187о.х.Ранняя стадия Ермаковского оледенения (III ледниковый)
8I120.08-102.56110.87.3757о.т.Ранняя стадия Казанцевского потепления
9I137-120.08127.565.1832х.
10I161.08-137144.86.4903у.т.
11I180.24-161.08171.085.4419х.?
12I200.6-180.24190.366.6781т.?

          При рассмотрении периодов 4I - 10I мы разбили Ермаковское похолодание и Казанцевский теплый период на две фазы: верхнюю и нижнюю, которые между собой пересекаются. Но это пересечение обусловлено теми свойствами, которые были наделены исследователями Ермаковский и Казанцевский горизонты при их введении. По-видимому, в дальнейшем с детализацией наших знаний о палеоклимате, идентификация горизонтов и палеоклиматических периодов будет уточняться.

          Изменение инсоляции Земли является причиной изменения палеоклимата. Поэтому выше причина сравнивалась со следствиями, которые происходили в Западной Сибири. Изотопные стадии изменения отложений в океане также являются следствием. С ними тоже необходимо сравнить изменение инсоляции. Как известно, набор бентосных записей изотопа кислорода δ18O за 5.3 млн. лет, например, LR4 в работе Lisiecki and Raymo [2005] составлен из данных, полученных в 57 местах мирового океана. С помощью статистических методов он был подогнан к изменениям инсоляции по прежним теориям за один день, а именно 21 июня на 65° северной широты. Экстремумы набора LR4 на интервале 500 т.л.н. обозначены как морские изотопные стадии МИС-1 ÷ МИС-13. Сейчас они широко используются при анализе палеоклимата. Например, по ним определяют объемы льда на Земле [Imbrie et al, 2011].

          Однако, как было показано выше, инсоляция имеет широкий спектр изменения, как по времени, так и по пространству. В холодную эпоху для высоких широт, в экваториальной зоне – теплее. При теплом лете зима – холоднее. Имеется отличие в инсоляции Северного и Южного полушарий. Поэтому необходимо устанавливать обусловленные инсоляцией механизмы изменения изотопного состава морских отложений. Полученная таким образом зависимость его от инсоляции позволит более достоверно изучать изменение климата и понимать его. В статье мы рассматриваем изменение инсоляции за 200 тыс. лет. В 2016 году планируется получить изменение инсоляции за период больший 5 млн. лет. Это позволит произвести сопоставление инсоляции, как причины изменения климата с реакцией донных отложений не только по изотопу кислорода δ18O, но и по другим параметрам.

Благодарности

          Выражаю признательность рецензентам: академику РАН Н.Л. Добрецову и чл.-корр. РАН М.В. Кабанову за тщательный анализ статьи, полезные рекомендации и замечания.

          Задачи об орбитальном и вращательном движении Земли решались на суперкомпьютерах Сибирского суперкомпьютерного центра СО РАН (г. Новосибирск).

Литература

  1. Арсланов Х.А., Лавров А.С., Потапенко Л.М. Новые данные о позднеплейстоценовом оледенении севера Западной Сибири. – Оледенения и палеоклиматы Сибири в плейстоцене. Новосибирск, «Наука», 1983, с. 27-35.
  2. Архипов С.А. Хронология геологических событий позднего плейстоцена Западной Сибири. Геология и геофизика, 1997, т. 38, № 12, с. 1863-1884.
  3. Архипов С.A. Главные геологические события позднего плейстоцена (Западная Сибирь). – Геология и геофизика, т. 41, № 6, 2000, с. 792-799.
  4. Архипов С.А., Астахов В.И., Волков И.В., Волкова В.С., Панычев В.А. Палеогеография Западно-Сибирской равнины в максимум позднезырянского оледенения. – Новосибирск : Наука, 1980, 109 с.
  5. Архипов С.А, Волкова В. С., Зыкина В. С., Бахарева В.А., Гуськов С.А., Левчук Л.К. Природно-климатические изменения в Западной Сибири в первой трети будущего столетия // Геология и геофизика, 1995, т. 36, № 8, с. 51-71.
  6. Большаков В.А., Капица А.П. Уроки развития орбитальной теории палеоклимата // Вестник Российской Академии Наук, 2011, т. 81, № 7, с. 603-612.
  7. Баулин В.В. История развития многолетнемерзлых пород в районе Нижней Оби в четвертичный период: Автореф. дис. канд. геол.-мин. наук. – М., МГУ, 1959. – 19 с.
  8. Бутвиловский В.В. Палеогеография последнего оледенения и голоцена Алтая. Событийно-катастрофическая модель. Томск, Изд-во Томск, ун-та, 1993, 218 с.
  9. Васильчук Ю.К. Закономерности развития инженерно-геологических условий севера Западной Сибири в голоцене. Автореф. канд. геол.-мин. наук. – М., МГУ, 1982. 27 с.
  10. Волков И.А., Архипов С.А. Четвертичные отложения района Новосибирска: оперативно-информационный материал. Новосибирск, «Наука», 1978, 90 с.
  11. Волков И.А., Волкова В.С., Задкова И.И. Покровные лессовидные отложения и палеогеография юго-запада Западной Сибири в плиоцен-четвертичное время. Новосибирск, «Наука», 1969, 332 с.
  12. Волкова В. С. Колебания Климата в Западной Сибири в позднеплиоценовое и четвертичное время // Эволюция климата, биоты и среды обитания человека в позднем кайнозое Сибири. Новосибирск, ОИГГМ СО АН СССР, 1991, с. 30-39.
  13. Вурком А. Астрономическая теория изменения климата/В кн. Изменение климата. – М.: ИЛ. 1958, с. 168-178.
  14. Гросвальд М.Г. Покровные ледники континентальных шельфов. M., «Наука», 1983, 216 с.
  15. Гросвальд М.Г. Оледенение Русского Севера и Северо-Востока в эпоху последнего великого похолодания// Материалы гляциологических исследований. Вып. 106. М.: «Наука», 2009. 152 с.
  16. Ершов Э.Д. Геокриология СССР. Западная Сибирь – М.: Недра, 1989. – 454 с.
  17. Илларионов А.Г. Тургайский спиллвей. К истории становления и развития речной сети Арало-Иртышского региона. Ижевск, 2013. 103 с.
  18. Лаухин С.А., Арсланов X.А., Шилова Г.Н., Величкевич Ф.Ю., Максимов Ф.Е., Кузнецов В.Ю., Чернов С.Б., Тертычная Т.В. Палеоклиматы и хронология средневюрмского мегаинтерстадиала на Западно-Сибирской равнине // Доклады РАН. 2006. Т. 411. № 4. С. 540-544.
  19. Ломанченков В.С. Об основных этапах геологического развития Лено-Колымской Приморской низменности в позднечетвертичную и современную эпохи. – Четвертичный период Сибири. Материалы докладов Всесоюзного совещания по изучению четвертичного периода, т. II. М., «Наука», 1966, с. 283-288.
  20. Мельников В.П., Смульский И.И. Астрономическая теория ледниковых периодов: Новые приближения. Решенные и нерешенные проблемы. – Новосибирск: Академическое изд-во “Гео”, 2009. – 98 с. Книга на двух языках. С обратной стороны: Melnikov V.P., Smulsky J.J. Astronomical theory of ice ages: New approximations. Solutions and challenges. – Novosibirsk: Academic Publishing House “GEO”, 2009. – 84 p. http://www.ikz.ru/~smulski/Papers/AsThAnR.pdf.
  21. Миланкович М. Математическая климатология и астрономическая теория колебаний климата. – М.-Л.: ГОНТИ, 1939. – 207 с.
  22. Некрасов И.А., Коновальчик Н.Г., Семенова Г.В., Скорбилин Н.А. История геокриологического исследования Западной Сибири. – Новосибирск: Наука. Сиб. отд-ние, 1990. – 270 с.
  23. Никифорова К.В., Певзнер М.А. Проблемы стратиграфии, геохронологии и корреляции четвертичной системы. – В кн.: Проблемы геологии и истории четвертичного периода (антропогена). М.: Наука, 1982, с. 5-99.
  24. Панычев В. А. Радиоуглеродная хронология аллювиальных отложений Предалтайской равнины. Новосибирск, Наука, 1979; 132 с.
  25. Пятосина Н. Мансийское озеро - море оказалось древнее, чем думали // Наука и жизнь. - 2005. - N 2. - с. 138-139
  26. Сакс В.Н. Четвертичный период в Советской Арктике. М.; Л., Морское и речное изд-во, 1953, 627 с.
  27. Смульский И.И. Комплексный показатель и уровни загрязнения атмосферы в периоды неблагоприятных метеорологических условий // Метеорология и Гидрология, – 1987. – N8. – С. 48 – 56.
  28. Смульский И.И. Анализ уроков развития астрономической теории палеоклимата // Вестник Российской Академии Наук, 2013 а. Т. 83. № 1. С. 31-39.http://elibrary.ru/item.asp?id=18448265. http://dx.doi.org/10.7868/S0869587313010118.
  29. Смульский И.И. Программа расчета инсоляции Земли в среде MathCad. / Институт криосферы Земли СО РАН. – Тюмень, 2013 б.http://www.ikz.ru/~smulski//Data/Insol/.
  30. Смульский И.И. Основные положения и новые результаты астрономической теории изменения климата / Институт криосферы Земли СО РАН. – Тюмень, 2014. – 30 с.: ил: 16.- Библиогр.: 44 назв. - Рус. Деп . в ВИНИТИ РАН 30.09.2014, № 258-В2014. http://www.ikz.ru/~smulski/Papers/OsPoATLP3.pdf.
  31. Фотиев С.М. Геокрилогические летописи Сибири // Криосфра Земли, 2009, т. XIII, № 3, с. 3-16.
  32. Шараф Ш. Г. и Будникова Н. А. Вековые изменения элементов орбиты Земли и астрономическая теория колебаний климата // Тр. Инст. теоретич. астрономии. – Вып. XIV. – Л.: Наука, 1969 г. – С. 48-109.
  33. Шполянская Н.А., Евсеев В.П. Выпуклобугристые торфяники северной тайги Западной Сибири // Тезисы докладов Всесоюзного совещания по мерзлотоведению. – М.: Изд-во МГУ, 1970. – С. 125-126.
  34. Berger A. and Loutre M. F. Insolation values for the climate of the last 10 million years // Quaternary Science Reviews. 1991. № 10. Р. 297-317.
  35. Edvardsson S., Karlsson K.G. and Engholm M. Accurate Spin Axes and Solar System Dynamics: Climatic Variations for the Earth and Mars // Astronomy &Astrophysics, 384, 2002, 689-701.http://dx.doi.org/10.1051/0004-6361:20020029.
  36. Imbrie J.Z., Imbrie-Moore A., Lisiecki L.E. A phase-space model for Pleistocene ice volume // Earth and Planetary Science Letters, 2011, 307 pp. 94–102.
  37. Laskar J., Robutel P., Joutel F., Gastineau M., Correia A.C.M., and Levrard B. A Long-term numerical solution for the Earth // Icarus 170, 2004. Iss. 2: 343-364.
  38. Lisiecki L.E., Raymo M.E. A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18O records // Paleoceanography, 2005, Vol. 20, PA1003, pp. 1-17.http://dx.doi.org/10.1029/2004PA001071.
  39. Smulsky J.J. The Influence of the Planets, Sun and Moon on the Evolution of the Earth's Axis // International Journal of Astronomy and Astrophysics, 2011, 1, 117-134.http://dx.doi.org/10.4236/ijaa.2011.13017. http://www.SciRP.org/journal/ijaa.
  40. Smulsky J.J. The System of Free Access Galactica to Compute Interactions of N-Bodies. I.J.Modern Education and Computer Science, 2012, 11, 1-20. Published Online December 2012 in MECS (http://www.mecs-press.org/) http://dx.doi.org/10.5815/ijmecs.2012.11.01.
  41. Smulsky J.J. Galactica Software for Solving Gravitational Interaction Problems // Applied Physics Research, 2012, Vol. 4, No. 2. P. 110-123.http://dx.doi.org/10.5539/apr.v4n2p110.
  42. Smulsky J.J. and Krotov O.I. New Computing Algorithm of the Earth's Insolation // Applied Physics Research, Vol. 6, No. 4; 2014, p. 56-82.http://dx.doi.org/10.5539/apr.v6n4p56.
  43. Smulsky J.J. and Smulsky Ya.J. Dynamic Problems of the Planets and Asteroids, and Their Discussion // International Journal of Astronomy and Astrophysics, 2012, 2, 129-155. http://dx.doi.org/10.4236/ijaa.2012.23018.
  44. Svendsen J.I., Astakhov V.I., Bolshiyanov D.Yu., Demidov I., Dowdeswell J.A., Gataflin V., Hjort C., Hubberten H.W., Larsen E., Mangerud J., Melles M., Moller P., Saarnisto M., Siegert M.J. Maximum extent of the Eurasian ice sheets in the Barents and Kara Sea region during the Weichselian. — Boreas, v. 28, № 1, 1999, p. 234-242.

Рекомендована к печати 30 октября 2015 г.Поступила в редакцию 29 апреля 2015 г.,
Н.Л. Добрецовым после доработки — 19 октября 2015 г.