Циклы Миланковича — Википедия

Ци́клы Мила́нковича описывают влияние изменений в движении Земли на её климат. Названы так в честь сербского астрофизика Милутина Миланковича, который в 1941 году сформулировал предположение, что регулярные и закономерные изменения формы орбиты Земли и ориентации её оси вращения приводят к циклическим изменениям количества солнечной радиации, попадающей на Землю, и эти изменения определяют наступление ледниковых эпох^[1].

Движение Земли[править | править код]

В своём суточном вращении вокруг своей оси и годичном обращении вокруг Солнца, Земля испытывает небольшие возмущения, вызванные гравитационным взаимодействием с другими телами Солнечной системы. Отклонения, вызванные этими возмущениями, сложны; среди них доминируют несколько циклов^[2]. Поскольку периоды этих циклов не кратны друг другу, регулярно возникают достаточно продолжительные эпохи, в которых они оказывают кумулятивное влияние, усиливая друг друга.

Орбита Земли, её эксцентристет, изменяется от почти круговой до слегка эллиптической. При более вытянутой орбите, мы имеем бо́льший разброс в расстоянии между Землёй и Солнцем, как и в количестве солнечной радиации, в разные времена года^[⇨]. Кроме того, наклон оси вращения Земли тоже слегка изменяется. Бо́льший наклон оси делает времена года более выраженными^[⇨]. Наконец, направление, в котором указывает ось Земли, изменяется относительно неподвижных звёзд (прецессия земной оси), а эллиптическая орбита Земли вокруг Солнца тоже поворачивается (Прецессия линии апсид). Совместный эффект этих движений заключается в том, что Земля приближается к Солнцу на наменьшее расстояние (перигелий) в разные астрономические сезоны^{[3][комм. 1]}.

Миланкович изучал изменения в тех движениях Земли, которые изменяют количество (и распределение) солнечной радиации, достигающей Земли. Он особенно выделял изменения на широте 65°N, на которой в настоящее время нижняя граница ледников опускается до уровня моря^[4], из-за большого количества суши на этой широте. Суша нагревается и остывает быстрее чем океан благодаря тому, что грунт имеет меньшую объёмную теплоёмкость чем вода^{[5][комм. 2]}, а в океане происходит перемешивание поверхностной и глубинной воды.

Эксцентричность орбиты Земли[править | править код]

Круговая орбита, эксцентриситет равен 0

Орбита с эксцентриситетом, равным 0.5; значение преувеличено для целей иллюстрирования; орбита Земли лишь слегка эксцентрична

Орбита Земли аппроксимируется эллипсом; мера отличия эллипса от окружности — это эксцентриситет. Форма земной орбиты изменяется от почти круговой (теоретически эксцентриситет может стать и нулём) до слегка эллиптической (за последние 250 миллионов лет наибольшим значением эксцентриситета было 0.0679)^[6]. Самый большой компонент этих изменений вносит свой вклад с периодом 405 тысяч лет (изменение эксцентриситета ±0.012). Другие компоненты цикличны с периодом 95 тысяч and 124 тысячи лет (с периодом биения 400 000 лет). Они свободно объединяются в цикл с периодом 100 тысяч лет (изменение эксцентриситета от −0.03 до +0.02). В настоящее время эксцентриситет равен 0.017 (0.0167) и уменьшается^[7].

Эксцентриситет изменяется в первую очередь из-за гравитационного приятяжения Юпитера и Сатурна. Большая полуось орбитального эллипса, однако, остаётся неизменной: в соответствии с теорией возмущений, вычисляющей эволюцию орбиты, большая полуось — инвариант. Орбитальный период (продолжительность сидерического года) тоже инвариант, так как определяется большой полуосью орбиты (третий закон Кеплера). Более долгосрочные вариации вызваны взаимодействиями, включающими перигелии и узлы Меркурия, Венеры, Земли, Марса и Юпитера^[6].

Влияние на температуру[править | править код]

Если большая полуось — константа, то, когда земная орбита становится более эксцентричной, малая полуось сокращается. Это увеличивает размах сезонных изменений^[8].

Влияние на продолжительность сезонов[править | править код]

Времена года это квадранты земной орбиты определяемые двумя солнцестояниями и двумя равноденствиями. Второй закон Кеплера гласит: за равные промежутки времени радиус-вектор, соединяющий Солнце и планету, описывает собой равные площади. Это означает, что орбитальная скорость планеты наибольшая при прохождении перигелия и наименьшая в афелии, и что Земля проводит меньше времени возле перигелии и больше времени возле афелия, а следовательно — продолжительность сезонов разная. В настоящее время перигелий происходит возле 3 Января, поэтому бо́льшая скорость Земли сокращает зиму и осень в северном полушарии. Лето в северном полушарии на 4,66 дня длиннее зимы, а весна на 2,9 дня длиннее осени. Бо́льший эксцентриситет увеличивает размах изменений орбитальной скорости Земли. Однако в настоящее время орбита Земли становится менее эксцентричной (более близкой к круговой). Это сделает времена года в ближайшем будущем более схожими по продолжительности^[10].

Наклон оси вращения Земли[править | править код]

Угол наклона оси Земли к плоскости орбиты (наклон эклиптики) варьируется от 22.1° до 24.5° в ходе долгопериодических (так называемых вековых) колебаний с протяжённостью цикла около 41 000 лет. Текущий наклон оси Земли составляет 23.44°, что примерно на полпути между крайними значениями. Последний раз наклон достигал своего максимума в 8700 году до нашей эры. Сейчас он находится в фазе убывания своего цикла и достигнет своего минимума примерно в 11 800 году нашей эры. Увеличение наклона оси увеличивает амплитуду сезонного цикла инсоляции, обеспечивая больше солнечной радиации летом в каждом полушарии и меньше зимой. Однако эти эффекты не одинаковы повсюду на поверхности Земли: увеличение наклона увеличивает общую годовую солнечную радиацию в более высоких широтах и уменьшает ее ближе к экватору^{[10] [комм. 3]}

Текущая тенденция уменьшения наклона сама по себе будет способствовать более мягким сезонам (более теплая зима и более холодное лето), а также общей тенденции к похолоданию. Поскольку бо́льшая часть снега и льда на планете лежит в высоких широтах, «уменьшающийся» наклон может способствовать окончанию межледниковья и наступлению ледникового периода по двум причинам: 1) общая летняя инсоляция меньше, и 2) меньше инсоляции в более высокие широты (где тает меньше снега и льда предыдущей зимы)^[10].

Прецессия земной оси[править | править код]

Прецессия земной оси — это изменение направления оси вращения Земли относительно неподвижных звезд с периодом около 25 700 лет. Это движение, также известное как предварение равноденствий, означает, что в конечном итоге Полярная звезда больше не будет северной полярной звездой. Прецессия вызвана приливными силами, оказываемая на вращающуюся Землю телами Солнечной системы, в первую очередь Солнцем и Луной, которые вносят примерно одинаковый вклад в этот эффект^[10].

Осевая прецессия делает сезонные контрасты более резкими в одном полушарии и менее резкими в другом. В настоящее время перигелий возникает зимой в Северном полушарии и летом в Южном полушарии. Это делает лето в Южном полушарии более жарким и смягчает сезонные колебания в Северном полушарии. Но примерно через 13 000 лет осевая прецессия заставит эти условия измениться: в Северном полушарии будут наблюдаться более экстремальные значения солнечной радиации, а в Южном полушарии будут наблюдаться более умеренные сезонные колебания^[10].

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

Прецессия земной орбиты[править | править код]

Сам орбитальный эллипс прецессирует в пространстве нерегулярным образом, совершая полный цикл относительно неподвижных звезд примерно за 112 000 лет.^{[11] [комм. 4]}. Апсидальная прецессия происходит в плоскости эклиптики и изменяет ориентацию орбиты Земли относительно эклиптики. Происходит это прежде всего в результате взаимодействия с Юпитером и Сатурном. Меньший вклад также вносят сжатие Солнца и эффекты общей теории относительности, которые хорошо известны для Меркурия^[12].

Прецессия земной орбиты в сочетании с 25700-летним циклом осевой прецессии (см. ^[⇨]) смещает момент прохождения перигелия. Апсидальная прецессия сокращает этот период примерно до 21 000 лет в настоящее время. Согласно относительно старому источнику (1965 г.), среднее значение за последние 300 000 лет составило 23 000 лет, варьируясь от 20 800 до 29 000 лет^[11].

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

Наклонение орбиты[править | править код]

Наклон орбиты Земли периодически изменяется относительно ее нынешней орбиты. Это трехмерное движение известно как «прецессия эклиптики» или «планетарная прецессия»; Миланкович не изучал его. Оно было обнаружено совсем недавно, период колебаний относительно земной орбиты составляет около 70 000 лет. Однако относительно неизменной плоскости, прецессия эклиптики имеет период около 100 000 лет. Этот период очень похож на 100 000-летний период эксцентриситета. Оба периода близко соответствуют 100 000-летней схеме ледниковых событий^[13].

Подтверждения и проблемы теории[править | править код]

Керны антарктического льда содержат захваченные пузырьки воздуха, соотношение изотопов кислорода ¹⁶O и ¹⁸O в которых является надежным индикатором (прокси) глобальных температур в период формирования льда. Изучение этих данных привело к выводу, что климатические изменения, зафиксированные в ледяных кернах, были обусловлены инсоляцией Северного полушария, как и предполагает гипотеза Миланковича^[14].

Проблема 100 тысяч лет[править | править код]

Миланкович считал, что из всех орбитальных циклов наибольшее влияние на климат оказывает наклонение земной оси, и что это происходит за счет изменения летней инсоляции в высоких широтах Северного полушария. Отсюда он вывел ледниковые периоды продолжительностью 41 тысяча лет^[15][16]. Однако последующие исследования^[1][17][18] показали, что циклы ледникового периода четвертичного оледенения за последний миллион лет имели период 100 тысяч лет, что соответствует основному циклу эксцентриситета. Были предложены различные объяснения этого несоответствия, включая частотную модуляцию^[19] или различные влияния (от углекислого газа или динамики ледникового покрова). Некоторые модели могут воспроизводить 100 000-летние циклы в результате нелинейного взаимодействия между небольшими изменениями орбиты Земли и внутренними колебаниями климатической системы^[20][21]. В частности, для описания этого взаимодействия изначально был предложен механизм стохастического резонанса^[22][23].

Некоторые утверждают, что прецессия имеет значение только тогда, когда эксцентриситет велик. Вот почему мы видим более сильный темп за 100 тысяч лет, чем за 21 тысячу лет^[24]. Другие считают, что история известных климатических данных недостаточно длинна, чтобы установить статистически значимую связь между климатом и изменениями эксцентриситета^[25].

Проблема перехода от 41 000-летнего к 100 000-летнему циклу[править | править код]

1–3 миллиона лет назад климатические циклы соответствовали 41 000-летнему циклу вариаций наклона оси вращения Земли (см. ^[⇨]). Спустя миллион лет назад произошел среднеплейстоценовый переход на 100 000-летний цикл, соответствующий основному циклу изменения эксцентриситета (см. ^[⇨]). «Проблема перехода» означает необходимость объяснить, что изменилось миллион лет назад^[26].

Среднеплейстоценовый переход теперь можно воспроизвести с помощью численного моделирования, которое учитывает удаление реголита ледниками и плейстоценовый тренд на уменьшение содержания углекислого газа^[27].

Проблема неразделённых циклов эксцентриситета[править | править код]

Даже хорошо датированные климатические данные за последний миллион лет не совсем соответствуют форме кривой эксцентриситета. Эксцентриситет имеет циклы продолжительностью 95 тысяч и 125 тысяч лет. Некоторые исследователи, однако, сообщают, что данные не показывают этих пиков, а указывают только на один цикл продолжительностью 100 000 лет^[28]. Однако такое разделение на два компонента эксцентриситета наблюдалось по меньшей мере один раз в керне, взятом в скандинавской формации Alum Shale (глиноземный сланец) возрастом 500 миллионов лет.

Проблема причинности[править | править код]

Это образное название описывает проблему с началом морской изотопной стадии под номером 5 (МИС 5 или MIS 5). Глубоководные керны показывают, что этот межледниковый интервал начался 130 тысяч лет назад. Однако это на 10 тысяч лет раньше, чем подъём уровня инсоляции, предсказываемый гипотезой Миланковича (проблемой причинности это названо из-за того, что отклик опережает причину)^[29].

Следствия теории Миланковича[править | править код]

Эпохи, способствующие возникновению оледенения[править | править код]

Это эпохи, когда происходит сочетание следующих факторов:

• Эксцентриситет орбиты Земли достигает умеренных и высоких значений;
• Дата прохождения Землёй перигелия близка к дате зимнего солнцестояния в северном полушарии.

При таком сочетании Земля движется по удалённой части своей орбиты тогда, когда в северном полушарии лето. В результате лето северного полушария становится более длительным (интервал между датами весеннего и осеннего равноденствия становится больше полугода, так как орбитальная скорость Земли при движении по удалённой части эллиптической орбиты становится меньше средней) и прохладным (расстояние от Земли до Солнца больше среднего), что является фактором, способствующим росту оледенения. Миланкович писал: «Не суровая зима, но прохладное лето способствует надвиганию ледников»^[30].

Эпохи, способствующие потеплению[править | править код]

Спустя примерно 11 тысяч лет с перигелием совпадает момент летнего солнцестояния, а эксцентриситет не успевает существенно измениться. Теперь лето в северном полушарии становится коротким и жарким, что ведёт к уменьшению ледникового покрова. При этом в южном полушарии устанавливаются условия, способствующие оледенению. Но там почти нет суши в умеренных и субантарктических широтах, где могли бы увеличиваться ледники. В целом по Земле площадь ледников сокращается, альбедо планеты сокращается, среднегодовая температура растёт^[4].

Ситуация в настоящее время[править | править код]

В нынешнюю эпоху разница между зимним солнцестоянием (21 декабря) и прохождением перигелия (3 января) составляет всего 13 дней, но эксцентриситет сейчас равен 0,0167, что существенно меньше среднего (максимальное значение 0,0658), и продолжает уменьшаться. В связи с этим сезонные колебания орбитальной скорости Земли и расстояния до Солнца невелики, и вносимые ими сезонные изменения приходящей к Земле солнечной энергии незначительны.

Прогноз[править | править код]

В настоящее время Земля переживает пик температуры — межледниковье — один из самых тёплых за последний миллион лет. Подобный пик с аналогичными значениями можно видеть 400 тыс. лет назад.

Хотя периоды межледниковья продолжаются от 10 до 30 тыс. лет, климатический оптимум держится всего несколько столетий. Вполне возможно, он уже закончился вместе со Средневековым климатическим оптимумом.

Приложения[править | править код]

Астрохронология[править | править код]

Климатические изменения в существенной мере определяются циклами Миланковича. Взаимно складываясь и вычитаясь, они регулярно и на продолжительное время влияют на изменения количества солнечной радиации, попадающей на Землю. При изменениях температуры изменяется содержание стабильных изотопов кислорода ¹⁶O/¹⁸O в морской воде. Поэтому изменение величины δ¹⁸O (см. δ18O^[en], Оценка палеоклиматических условий), извлекаемое, к примеру, из ракушек бентосных фораминифер, захороненных в последовательных слоях эвапоритов и нормальных морских отложений^[31], отражает циклы Миланковича, а они могут быть рассчитаны из моделей движения Земли, что позволяет получать абсолютный возраст осадочных слоёв. Астрономическая шкала, к которой привязэываются измерения δ¹⁸O, получается корреляцией характеристических паттернов седиментарных циклов в нормальных морских отложениях с вариациями инсоляции на разных широтах (часто вслед за Миланковичем выбирается широта 65°N, на которой в настоящее время нижняя граница ледников опускается до уровня моря^[4]), рассчитанным по колебаниям параметров земной орбиты в используемой модели движения Земли^[32]. Для независимого подтверждения такой шкалы используются методы радиоизотопного датирования. См. также Морская изотопная стадия и Циклостратиграфия.

Комментарии[править | править код]

Примечания[править | править код]

См. также[править | править код]

Ссылки[править | править код]

• Циклы Миланковича  (неопр.). Элементы. Дата обращения: 14 декабря 2012. Архивировано 30 мая 2012 года.

Литература[править | править код]

• Бялко А.В. Наша планета — Земля. — 2-е изд., перераб. и доп.. — М.: Наука. Гл. ред. физ.-мат. лит., 1989. — С. 222-228. — 240 с. — (Б-чка «Квант». Вып. 29). — ISBN 5-02-014079-1.
• Миланкович М. Математическая климатология и астрономическая теория колебаний климата. — М., Л., 1939.
• Н.Г. Нургалиева. Изотопная стратиграфия  (рус.) (PDF). Казанский ун-т (2017).
• Abe-Ouchi A, Saito F, Kawamura K, Raymo ME, Okuno J, Takahashi K, Blatter H (August 2013). "Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume". Nature. 500 (7461): 190—3. Bibcode:2013Natur.500..190A. doi:10.1038/nature12374. PMID 23925242. S2CID 4408240.
• Abu-Hamdeh (2020). "Thermal Properties of Soils as affected by Density and Water Content". Biosystems Engineering. 86 (1): 97—102. doi:10.1016/S1537-5110(03)00112-0. Дата обращения: 16 мая 2021. Объёмная теплоёмкость варьируется от 1.48 до 3.54 МДж/m³/°C для глины и от 1.09 до 3.04 МДж/м³/°C для песка при содержании влаги от 0 до 0·25 (кг/кг) [...]
• Barbieri, L.; Talamucci, F. (20 February 2018). "Calculation of Apsidal Precession via Perturbation Theory". Advances in Astrophysics. 4 (3). arXiv:1802.07115. doi:10.22606/adap.2019.43003. S2CID 67784452.
• Benzi, R; Sutera, A; Vulpiani, A (1 November 1981). "The mechanism of stochastic resonance". Journal of Physics A: Mathematical and General. 14 (11): L453—L457. Bibcode:1981JPhA...14L.453B. doi:10.1088/0305-4470/14/11/006. ISSN 0305-4470. S2CID 123005407.
• Benzi, Roberto; Parisi, Giorgio; Sutera, Alfonso; Vulpiani, Angelo (February 1982). "Stochastic resonance in climatic change". Tellus (англ.). 34 (1): 10—16. doi:10.1111/j.2153-3490.1982.tb01787.x.
• Berger A, Loutre MF, Mélice JL (2006). "Equatorial insolation: from precession harmonics to eccentricity frequencies" (PDF). Climate of the Past Discussions. 2 (4): 519—533. doi:10.5194/cpd-2-519-2006.
• Buis, Alan Milankovitch (Orbital) Cycles and Their Role in Earth's Climate  (неопр.). climate.nasa.gov. NASA (27 февраля 2020). — «За последний миллион лет он менялся от 22.1 до 24.5 градусов. [...]Чем больше угол наклона оси Земли, тем экстремальнее наши времена года[...] Большие углы наклона благоприятствуют периодам дегляциации (таяния и отступления ледников и ледниковых щитов). Эти эффекты не являются одинаковыми в глобальном масштабе: в более высоких широтах изменение общего количества солнечной радиации больше, чем в районах, расположенных ближе к экватору.[...]Ось Земли в настоящее время наклонена на 23.4 градуса[...] При увеличении ледяного покрова он отражает бо́льшую часть солнечной энергии, которая возвращается в космос, способствуя еще большему охлаждению.» Дата обращения: 10 мая 2021.
• Ghil M (1994). "Cryothermodynamics: the chaotic dynamics of paleoclimate". Physica D. 77 (1—3): 130—159. Bibcode:1994PhyD...77..130G. doi:10.1016/0167-2789(94)90131-7.
• G. K. Gilbert (February–March 1895). "Sedimentary Measurement of Cretaceous Time". The Journal of Geology. University of Chicago Press. 3 (2): 121—127. Bibcode:1895JG......3..121G. doi:10.1086/607150. JSTOR 30054556. S2CID 129629329. По мере того, как земная ось медленно описывает свою окружность на набесной сфере, наступление времен года по отношению к перигелию постоянно смещается.
• Gildor H, Tziperman E (2000). "Sea ice as the glacial cycles' climate switch: Role of seasonal and orbital forcing". Paleoceanography. 15 (6): 605—615. Bibcode:2000PalOc..15..605G. doi:10.1029/1999PA000461.
• Girkin AM (2005). A Computational Study on the Evolution of the Dynamics of the Obliquity of the Earth (Master of Science thesis). Miami University. Архивировано из оригинала (PDF) 30 сентября 2014.
• James D. Hays, John Imbrie, Nicholas Shackleton. Variations in the Earth's Orbit: Pacemaker of the Ice Ages (англ.) // Science : журнал. — 1976. — Vol. 194, no. 4270. — P. 1121–1132 pmid = 17790893. — doi:10.1126/science.194.4270.1121. — Bibcode: 1976Sci...194.1121H.
• F.J. Hilgen, W. Krijgsman, C.G. Langereis, L.J. Lourens, A. Santarelli, W.J. Zachariasse. Extending the astronomical (polarity) time scale into the Miocene. (англ.) // Earth and Planetary Science Letters. — 1995. — Vol. 136, no. 3—4. — P. 495—510. — ISSN 0012-821X. — doi:10.1016/0012-821X(95)00207-S.
• van den Heuvel EP (1966). "On the Precession as a Cause of Pleistocene Variations of the Atlantic Ocean Water Temperatures". Geophysical Journal International. 11 (3): 323—336. Bibcode:1966GeoJ...11..323V. doi:10.1111/j.1365-246X.1966.tb03086.x.
• Imbrie, John. Ice Ages: Solving the Mystery / John Imbrie, Katherine Palmer Imbrie. — Harvard University Press. — P. 158. — ISBN 978-0-674-44075-3.
• Karner DB, Muller RA (June 2000). "PALEOCLIMATE: A Causality Problem for Milankovitch". Science. 288 (5474): 2143—4. doi:10.1126/science.288.5474.2143. PMID 17758906. S2CID 9873679.
• Laskar J, Fienga A, Gastineau M, Manche H (2011). "La2010: A New Orbital Solution for the Long-term Motion of the Earth" (PDF). Astronomy & Astrophysics. 532 (A889): A89. arXiv:1103.1084. Bibcode:2011A&A...532A..89L. doi:10.1051/0004-6361/201116836. S2CID 10990456..
• Laskar, J. Astrochronology // Geologic Time Scale 2020. — Elsevier, 2020. — ISBN 978-0-12-824360-2. — doi:10.1016/b978-0-12-824360-2.00004-8.
• Lee JE, Shen A, Fox-Kemper B, Ming Y (1 January 2017). "Hemispheric sea ice distribution sets the glacial tempo". Geophys. Res. Lett. 44 (2): 1008—1014. Bibcode:2017GeoRL..44.1008L. doi:10.1002/2016GL071307.
• Kawamura K, Parrenin F, Lisiecki L, Uemura R, Vimeux F, Severinghaus JP, et al. (August 2007). "Northern Hemisphere forcing of climatic cycles in Antarctica over the past 360,000 years". Nature. 448 (7156): 912—6. Bibcode:2007Natur.448..912K. doi:10.1038/nature06015. PMID 17713531. S2CID 1784780.
• Milankovitch, M. Canon of Insolation and the Ice Age Problem. — Belgrade : Zavod za Udz̆benike i Nastavna Sredstva. — ISBN 978-86-17-06619-0.
• Muller RA, MacDonald GJ (August 1997). "Spectrum of 100-kyr glacial cycle: orbital inclination, not eccentricity". Proceedings of the National Academy of Sciences of the United States of America. 94 (16): 8329—34. Bibcode:1997PNAS...94.8329M. doi:10.1073/pnas.94.16.8329. PMC 33747. PMID 11607741.
• Rial JA (October 2003). "Earth's orbital Eccentricity and the rhythm of the Pleistocene ice ages: the concealed pacemaker" (PDF). Global and Planetary Change. 41 (2): 81—93. Bibcode:2004GPC....41...81R. doi:10.1016/j.gloplacha.2003.10.003. Архивировано (PDF) из оригинала 20 июля 2011.
• Shackleton NJ, Berger A, Peltier WR (3 November 2011). "An alternative astronomical calibration of the lower Pleistocene timescale based on ODP Site 677". Transactions of the Royal Society of Edinburgh: Earth Sciences. 81 (4): 251—261. doi:10.1017/S0263593300020782. S2CID 129842704.
• Willeit M, Ganopolski A, Calov R, Brovkin V (April 2019). "Mid-Pleistocene transition in glacial cycles explained by declining CO₂ and regolith removal". Science Advances. Наклон оси вращения Земли5 (4): eaav7337. Bibcode:2019SciA....5.7337W. doi:10.1126/sciadv.aav7337. PMC 6447376. PMID 30949580.
• Wunsch C (2004). "Quantitative estimate of the Milankovitch-forced contribution to observed Quaternary climate change". Quaternary Science Reviews. 23 (9—10): 1001—12. Bibcode:2004QSRv...23.1001W. doi:10.1016/j.quascirev.2004.02.014.
• Yu, Zhiwei; Ding, Zhongli Nonlinear coupling between 100 ka periodicity of the paleoclimate records in loess and periodicities of precession and semi-precession  (неопр.).
• Zachos JC, Shackleton NJ, Revenaugh JS, Pälike H, Flower BP (April 2001). "Climate response to orbital forcing across the Oligocene-Miocene boundary". Science. 292 (5515): 274—8. Bibcode:2001Sci...292..274Z. doi:10.1126/science.1058288. PMID 11303100. S2CID 38231747. Архивировано из оригинала 3 декабря 2017. Дата обращения: 24 октября 2010.