И. И. Смульский Институт криосферы Земли со ран, г. Тюмень
Скачать 128.49 Kb.
|
Некоторые Астрономические Проблемы Эволюции КриосферыИ.И. Смульский Институт криосферы Земли СО РАН, г. Тюмень Реферат Рассмотрены проблемы Астрономической теории ледниковых периодов с позиций небесной механики. При интегрировании упрощенных уравнений движения за большие промежутки времени прежние исследователи получали неустойчивые движения. Потому они пришли к выводу о невозможности рассчитать инсоляцию Земли за периоды большие 20 млн. лет. В настоящей работе решается неупрощенная задача об орбитальном движении новым численным методом и уравнения проинтегрированы за 100 млн. лет. Получены все периоды и амплитуды колебаний орбит планет и Луны и установлена их устойчивость. Дифференциальные уравнения вращательного движения также решаются численным методом без упрощения. Определены результаты воздействия планет и Солнца по отдельности на ось Земли. Эволюция оси Земли рассмотрена также на основе составной модели вращения Земли. Получены периоды ее колебаний, которые совпадают с наблюдаемыми. Установлено, что только решение неупрощенных уравнений вращения Земли позволит достоверно рассчитать эволюцию инсоляции и определить все периоды ее изменения. Ключевые Слова: Вычисления; уравнения; инсоляция; наклон; перигелий; эксцентриситет. Введение При изучении проблемы потепления климата, произошедшего во второй половине 20-ого века, и при разработке его моделей (Изменение климата 2007) исследователи учитывают также Астрономическую теорию ледниковых периодов, которую создал М. Миланкович (1939). Неопределенности моделей потепления климата в какой-то мере зависят от неопределенностей этой теории. Ряд исследователей, например, Большаков и Капица (2001) считает, что со времен М. Миланковича накопились палеоклиматические данные, которые свидетельствуют о необходимости дальнейшего ее развития. В настоящей работе рассматриваются результаты по уточнению Астрономической теории палеоклимата с позиций небесной механики или астрономии.  Рис. 1. Сравнение инсоляции и параметров орбиты Земли: I - инсоляция в эквивалентных широтах на широте 65° северного полушария по данным Шараф и Будниковой (1969) и угол наклона орбиты Земли к подвижной плоскости экватора (ε) согласно (Berger & Loutre 1991); эксцентриситет e и углы (в радианах) положения перигелия от неподвижного φp и от подвижного φpγ восходящего узла γ – по нашим расчетам (Мельников и Смульский 2009). T – время в kyr (тыс. лет), отсчитывается в тысячах лет в прошлое от современной эпохи 30.12 1949 г. Зависимость инсоляции от орбитального и вращательного движений Земли В астрономической теории (Миланкович 1939) установлено, что инсоляция земной поверхности определяются эксцентриситетом орбиты Земли e и двумя параметрами орбиты: углом положения перигелия φpγ от подвижного восходящего угла γ и углом наклона орбиты ε к подвижному экватору. Типичное изменение инсоляции I за прошедшие 500 тыс. лет представлено на рис. 1. Здесь же приведены эволюции параметров: угла наклона подвижного экватора Земли к подвижной плоскости её орбиты, эксцентриситета e и угла перигелия p. Угол перигелия p показан в виде циклов от 0 до и в непрерывном виде. На рис. 1 видно, что за 500 тыс. прошедших лет инсоляция I в северном полушарии изменялась колебательно: максимумы свидетельствуют о возможных потеплениях, а глубокие минимумы при T = -120 kyr и T = -230 kyr – о возможных оледенениях. Например, максимум при T = -10 kyr последнего потепления произошёл 10 тыс. лет назад. Он совпадает с моментом исчезновения ледникового щита в Евразии и в Северной Америке. Как видно из рис 1, период изменения угла наклона , равный T = 41.1 тыс. лет, являются основным периодом колебания инсоляции I. Период изменения эксцентриситета e, равный Te1 = 95 тыс. лет, является вторым существенным периодом изменения инсоляции. Угол перигелия p изменяется неравномерно. Перигелий орбиты Земли вращается в направлении орбитального её движения со средним периодом Tp = 147 тыс. лет относительно неподвижного пространства. Плоскость экватора Земли прецессирует в направлении противоположном орбитальному движению Земли с периодом Tpr = 25.7 тыс. лет. Поэтому, относительно подвижной плоскости экватора период обращения перигелия (угол p) в среднем составляет Tp = 21.9 тыс. лет. Это минимальный период колебания инсоляции I, который имеется на рис. 1. Недостатки прежних теорий и программа исследований Итак, эволюция инсоляции Земли определяется изменениями орбиты Земли и оси её вращения. Со времён И. Ньютона орбитальная задача упрощалась и сводилась к взаимодействию двух тел: Солнца и планеты или планеты и спутника. А все остальные тела рассматривались как факторы малых возмущений, и их воздействие разлагалось в ряды. Таким образом, удавалось решить эту задачу приближёнными аналитическими методами. Этот подход сохранился и в 20-м веке, но только точность представления решений рядами непрерывно улучшалась. Например, сейчас эти ряды содержат несколько сот членов. Дифференциальные уравнения вращательного движения упрощались более радикально. В этих уравнениях пренебрегали вторыми производными и произведениями первых производных. Учитывалось воздействие только Луны и Солнца на вращательное движение Земли, а их движение описывалось приближённо. Таким образом, были получены аналитические выражения для прецессии экватора Земли. Эти выражения не содержали короткопериодических колебаний параметров вращательного движения Земли. Скорость прецессии оси Земли по существу определялась наблюдаемой скоростью прецессии, т.е. её эволюция оставалась неизвестной. А полученный единственный период колебаний в 41.1 тыс. лет оси Земли нельзя было подтвердить каким-либо другим способом. Приближённый характер решения орбитальной и вращательной задач привёл к тому, что при решении их на большие промежутки времени результаты начинали расходиться. Поэтому ряд авторов, например Ж. Ляскар и др. (2004) пришли к выводу о неустойчивости Солнечной системы и о невозможности определения инсоляции Земли за период времени больший 20 млн. лет. В связи с этим возникает необходимость решения этих двух задач с минимальными упрощениями. Это можно сделать, используя современные численные методы и суперкомпьютеры. Вначале должна быть решена первая задача об эволюции орбитального движения. Используя её результаты, можно будет решить вторую задачу об эволюции оси Земли. Решение двух задач позволит рассчитать эволюцию инсоляции Земли. Затем при сопоставлении изменения инсоляции с эволюцией природных процессов можно будет установить зависимости, которые позволят прогнозировать развитие криосферных процессов на Земле. Решение орбитальной задачи Уравнения орбитального движения представляют систему 3n нелинейных дифференциальных уравнений (Смульский 1999)  i = 1,2,…,n, (1)где  – радиус-вектор относительно центра масс Солнечной системы; G – гравитационная постоянная; – радиус-вектор от тела с массой mk до тела с массой mi; п = 11 (девять планет, Солнце и Луна).Для решения уравнений (1) мы разработали новый метод численного интегрирования (Смульский 1999), который позволил решить её за 100 млн. лет (Мельников и Смульский 2009). На рис. 2 показана эволюция параметров орбиты Земли за 3 млн. лет в прошлое. Эксцентриситет е испытывает короткопериодические изменения с главным периодом Те1 = 95 kyr (тыс. лет) вокруг среднего значения еm = 0.028. Кроме того, наблюдаются более долгие колебания с периодами Те2 = 413 kyr и Те3 = 2.31 Myr, которые приводят к достижению крайних значений эксцентриситета е = 0.0003 и е = 0.065. На рис. 3 показано изменение скользящих средних эксцентриситета es на интервале 2·Tel, при этом видны колебания эксцентриситета с наибольшим периодом Te3.  Рис. 2. Эволюция орбиты Земли за 3 млн. лет в прошлое: е - эксцентриситет; φΩ - угловое положение восходящего узла плоскости орбиты; i - угол ее наклона к неподвижной плоскости экватора; φр - угловое положение перигелия; Т – время в миллионах лет от 1950 г. Углы даны в радианах. При решении орбитальной задачи углы положения орбиты рассматриваются относительно неподвижной плоскости экватора. Долгота восходящего узла φΩ орбиты и угол ее наклона і испытывает колебания с периодом ТS = 68.7 тыс. лет вокруг среднего значения іm = 0.402 радиан. Диапазон колебаний составляет 5.64º. Мы установили, что эти изменения обусловлены вращением оси орбиты  с периодом ТS = 68.7 kyr вокруг неподвижного вектора момента  всей Солнечной системы в направлении обратном обращению планеты вокруг Солнца. Это вращение или прецессия оси орбиты показано на рис. 3 в виде изменения угла прецессии s. Кроме прецессионного движения ось орбиты претерпевает нутационные колебания угла наклона qS к вектору. Максимальное отклонение оси от момента составляет qSmax = 2.94º. Главный период нутационных колебаний Т1 = 97.35 тыс. лет. На графике приведены скользящие средние угла нутации qSs, которые усреднены на интервале 2·Т1. На рис. 3 виден второй период нутационных колебаний Т2 = 1.164 млн. лет.Перигелий орбиты (см. φp на рис. 2) перемещается в направлении обращения Земли вокруг Солнца, совершая в среднем один оборот за Tp = 147 kyr. При этом угол перигелия φр немонотонно увеличивается со временем и наряду с вращением перигелия против стрелки часов, наблюдаются возвратные движения по часовой стрелке. А периоды полного оборота перигелия, как уже отмечалось ранее, изменяются в несколько раз. Было установлено, что эволюция орбит планет и Луны происходит в результате четырех движений: 1) вращения или, другими словами, прецессии оси орбиты; 2) нутационных колебаний оси орбиты; 3) колебаний эксцентриситета орбиты; 4) вращения орбиты в своей плоскости (вращение перигелия). Исследования показали, что оси орбит Земли, планет и Луны, как и ось вращения Земли, ведут себя идентично. Вместе с тем, в отличие от осей планетных орбит, которые вращаются вокруг неподвижного вектора момента количества движения Солнечной системы , ось орбиты Луны, как и ось вращения Земли, прецессирует вокруг подвижной оси орбиты Земли.Вращение перицентриев планет и Луны происходит в направлении орбитального движения, за исключением Плутона, перигелий которого вращается в обратном направлении. Аналогично рис. 3 выглядит изменение параметров орбиты Земли на интервале -50 Myr ≥ T ≥ -100 Myr, т.е. параметры орбиты колеблются в неизменном режиме. Такие же исследования выполнены по всем планетам и получены неизменные режимы колебания параметров их орбит. Таким образом, эволюция орбит планет на исследованном интервале в 100 млн. лет неизменна и устойчива. В работах других авторов (Laskar et al 2004) хаотичность решений упрощенных уравнений движения объясняют наличием резонансов и неустойчивостей. Они не возникают при интегрировании неупрощенных уравнений орбитального движения планет.  Рис. 3. Эволюция параметров орбиты Земли за -50 млн. лет в прошлое: es – скользящие средние эксцентриситета; φp – угол перигелия; s – угол прецессии; qSs – скользящие средние угла нутации оси орбиты. Углы jр и s даны в радианах, qSs – в градусах. Myr – млн. лет. Исследование вращательного движения Исследования вращательного движения мы проводим одновременно двумя способами. При первом способе интегрируются дифференциальные уравнения вращательного движения, а при втором – вращательное движение Земли моделируется совокупностью материальных точек. Мы по-новому вывели дифференциальные уравнения вращательного движения (Smulsky 2011), взяв за основу теорему изменения момента количества движения. В результате анализа следствий этой теоремы было установлено, что ось вращения Земли будет испытывать прецессионные и нутационные колебания: 1) с периодами равными полупериодам обращения планет, Солнца и Луны относительно подвижной оси вращения Земли; 2) с периодами равными интервалам сближения ближних планет с Землёй в наиболее отдалённых от экватора Земли точках и 3) с периодами равными периодам прохождения воздействующих тел в этих точках в моменты наибольшего наклона плоскости орбиты тела к плоскости экватора. Последние периоды будут иметь продолжительность от десятков до сотен тысяч лет. При решении задачи о вращении Земли углы положения оси ее вращения: угол наклона к оси орбиты (угол нутации) и угол прецессии рассматриваются относительно неподвижной плоскости орбиты Земли. Дифференциальные уравнения вращения Земли мы получили в следующем виде  (2)  (3) Рис. 4. Нутационные колебания и тренды оси Земли при одиночном воздействии планет: а – на интервале 100 лет = 1 cyr, б – на интервале 1000 лет;  – разность углов нутации, где 0 – угол нутации в начальную эпоху;  – средняя скорость нутации в радианах в столетие; на графиках периоды основных колебаний оси Земли равны от Меркурия до Луны, соответственно: 6.6, 8.1, 15.8, 5.9, 14.7, 42, 82.4, 248, 0.5, 18.6 лет. Углы даны в радианах. , (4)где Jx, Jy и Jz – моменты инерции Земли на оси вращающейся системы координат xyz;  – динамическая эллиптичность Земли; – комплекс параметров; = const – проекция абсолютной скорости вращения Земли на ее ось z – относительная скорость вращения Землиn – количество действующих на Землю тел; mi и x1i, y1i, z1i – масса и координаты i-того воздействующего тела. Вышеупомянутые выводы из теоремы о периодах колебаний были подтверждены при интегрировании уравнений вращательного движения (2)-(4) при воздействии планет, Луны и Солнца на Землю по отдельности (см. рис. 4). Как видно из графиков, под воздействием Солнца, Меркурия и Венеры угол нутации увеличивается. Воздействие других планет приводит к его уменьшению. На графиках приведены средние скорости изменения угла нутации . Они обусловлены движениями плоскостей орбит воздействующих тел.Прослеживаются разные типы нутационных колебаний. Большинство тел создает колебания с полупериодом своего орбитального движения. Для ближних планет наибольшие колебания обусловлены периодом их сближения с Землей в наиболее отдаленной от экватора точке. Например, периоды (см. рис. 4): 5.9, 14.7, 42, 82.4 лет и 0.5 года равны полупериодам обращения планет и Солнца; периоды: 6.6, 8.1 и 15.8 лет – периодам сближения планет с Землёй в наиболее отдалённых от экватора точках; период 18.6 лет равен периоду колебания плоскости орбиты Луны относительно плоскости экватора. Наибольшее воздействие на вращательное движение Земли оказывает Луна, затем Солнце, а из планет – Венера. Приведенные на рис. 4 периоды основных колебаний оси Земли не учитываются в Астрономической теории эволюции криосферы. Дальнейшие наши работы связаны с доработкой нашего подхода, с целью интегрирования уравнений вращательного движения при совместном воздействии на Землю всех тел за большие промежутки времени.  Рис. 5. Нутационные колебания оси вращения составной модели Земли No. 3 на разных интервалах времени: Δε – отклонение угла нутации от скользящего среднего; Δε1, Δε2, Δε3, и Δε4 – амплитуды колебаний, 1 cyr =100 лет. Углы даны в радианах. Составная модель вращательного движения Земли Вращательное движение моделируется совокупностью тел, которые осесимметрично расположены в экваториальной плоскости вокруг центрального тела. Под воздействием притяжения Луны, планет и Солнца эта система эволюционирует, и ось орбиты периферийного тела имитирует эволюцию оси вращения Земли. Было исследовано несколько моделей с разными параметрами. На рис. 5 показаны нутационные колебания оси вращения, полученные для модели 3. Рассматривается угол ε положения оси Земли относительно подвижной оси ее орбиты. Во всех моделях ось орбиты периферийного тела прецессирует относительно подвижной оси орбиты Земли. Этот результат является принципиально важным. Он позволяет на больших интервалах времени контролировать точность интегрирования дифференциальных уравнений вращательного движения. Получены (см. рис. 5) периоды колебаний оси Земли полумесячные (13.9 дней), полугодовые (0.5 г.), 18.6 года (обусловленные прецессией орбиты Луны) и 2580 лет. Первые три периода следуют из теоремы моментов и подтверждаются данными наблюдений. В составных моделях отсутствует период колебаний 41.1 тыс. лет. Дальнейшие решения задачи о вращательном движении позволят прояснить вопрос об этом периоде и определить другие долгие периоды. Заключение Мы остановились на проблемах астрономической теории эволюции криосферы, которые связаны с точностью решения орбитальной и вращательной задач. Это – сложнейшие задачи, в которых имеется ещё ряд нерешённых вопросов. В отличие от предшественников мы не упрощаем уравнения, а решаем их численно. Поэтому, возможно, на часть нерешённых проблем удастся пролить свет. Например, теперь мы можем с уверенностью сказать об устойчивости Солнечной системы. Установлены многие особенности орбитального движения. В результате решения задачи о вращении Земли двумя способами получены дополняющие друг друга представления о ее эволюции. Мы надеемся, что этот подход позволит уточнить положения астрономической теории эволюции криосферы. Благодарности В работе на разных этапах принимали участие: Л.И. Смульский, Я.И. Смульский, О.И. Кротов, К.Е. Сеченов, А.А. Павлова, И.А. Шаболина, М.Л. Панова, Е.Ф. Сафина и др. Основные вычисления были выполнены на суперкомпьютерах Сибирского Суперкомпьютерного Центра СО РАН (г. Новосибирск) и Института Прикладной Математики им. М.В. Келдыша. Данная работа выполнена при поддержке грантов губернатора Тюменской области и интеграционной программы Президиума РАН № 13. Литература Berger A. & Loutre M. F. 1991. Insolation values for the climate of the last 10 million years. Quaternary Science Reviews 10: 297 - 317. Большаков В.А. и Капица А.П. 2011. Уроки развития орбитальной теории палеоклимата. Вестник Российской Академии Наук 81, No. 7: 603-612. (на русском языке) Изменение климата, 2007 г.: физическая научная основа. Вклад Рабочей группы I в Четвертый доклад МГЭИК об оценках. Из-во кембриджского ун-та., 163 с. Laskar J., Robutel P., Joutel F., Gastineau M., Correia A. C. M., and Levrard B. 2004. A Long-term numerical solution for the Earth. Icarus 170. Iss. 2: 343-364. Мельников В.П. и Смульский И.И. 2009. Астрономическая теория ледниковых периодов: Новые приближения. Решенные и нерешенные проблемы. – Новосибирск: Академическое изд-во «Гео», 192 с. Книга на двух языках. С обратной стороны: Melnikov V.P. & Smulsky J.J. Astronomical theory of ice ages: New approximations. Solutions and challenges. http://www.ikz.ru/~smulski/Papers/AsThAnR.pdf. Миланкович М. 1939. Математическая климатология и астрономическая теория колебаний климата. – М.-Л. –ГОНТИ, 207 с. (на русском языке) Смульский И. И. 1999. Теория взаимодействия. - Новосибирск: Из-во Новосибирского ун-та, ННЦ ОИГГМ СО РАН. 294 с. (на русском языке) http://www.ikz.ru/~smulski/TVfulA5_2.pdf. Smulsky J.J. 2011. The Influence of the Planets, Sun and Moon on the Evolution of the Earth’s Axis. International Journal of Astronomy and Astrophysic.  1: 117-134.Шараф Ш. Г. и Будникова Н. А. 1969. Вековые изменения элементов орбиты Земли и астрономическая теория колебаний климата. Тр. Инст. теоретич. Астрономии. Вып. XIV. - Л.: Наука. 48 - 109. (на русском языке) |
Добавить документ в свой блог или на сайт
Похожие:
 | И. И. Смульский Институт криосферы Земли со ран, г. Тюмень Дэвид Карузо Ph. D., один из авторов наиболее авторитетной методики оценки эмоционального интеллекта msceit, коллега Питера Саловея... | | Кризисные явления, возможности и пути их преодоления в социально-экономических... Биробиджан. Институт комплексного анализа региональных проблем дво ран. Амурский государственный университет. Тихоокеанский институт... |
| Российской Академии Наук Институт проблем нефти и газа со ран министерство... Председатель – Александр Федотович Сафронов, чл корр. Ран, председатель Президиума Якутского научного центра со ран, директор ИПНГ... | | Методы дистанционного зондирования земли радарами с синтезированной апертурой Работа выполнена во Фрязинском филиале Федерального государственного бюджетного учреждения науки Института радиотехники и электроники... |
| Российской академии наук институт европы ран промышленная политика европейских стран Н. В. Говоровой.]. – М. Ин-т Европы ран : Рус сувенир, 2010. – 214 с. – (Доклады Института Европы = Reports of the Institute of Europe... | | Численное моделирование ионосферных предвестников сильных землетрясений... Западное отделение Учреждения ран «Институт земного магнетизма, ионосферы и распространения радиоволн им. Н. В. Пушкова» ран |
| Тезисы докладов VIII международной конференции. Москва, 4-6 октября... Российская академия наук институт биохимической физики им. Н. М. Эмануэля ран институт химической физики им. Н. Н. Семенова ран | | Кошкарев А. В., Медведев А. А., Атаева О. М., Серебряков В. А., Теймуразов... Муниципальное бюджетное общеобразовательное учреждение «Хочашевская основная общеобразовательная школа» Ядринского района Чувашской... |
| Учреждение Российской Академии Наук Институт проблем передачи информации... Алексей Леонтьевич Бызов родился 2 января 1926 года в Москве. В 1950 году окончил биофак мгу, кафедру физиологии животных | | А. С. Шишикин Зам директора, д Х. н., Ил со ран Ответственный исполнитель, зам директора Институт леса им. В. Н. Сукачева со ран, с н с., д б н |
 | Доклад о состоянии и использовании Земли промышленности, энергетики, транспорта, связи, радиовещания, телевидения, информатики, земли для обеспечения космической деятельности,... | | Доклад о состоянии и использовании Земли промышленности, энергетики, транспорта, связи, радиовещания, телевидения, информатики, земли для обеспечения космической деятельности,... |
| Фгбниу «Государственный институт искусствознания» Институт экономики ран Механизмы поддержки музейной деятельности в условиях мирового экономического кризиса (на примере музеев Европы и сша) | | Учреждение Российской Академии наук Институт народнохозяйственного... Разработка единой системы анализа и прогнозирования баланса денежных доходов и расходов населения, уровня и структуры бедности с... |
| Министерство экономического развития российской федерации Земли промышленности, энергетики, транспорта, связи, радиовещания, телевидения, информатики, земли для обеспечения космической деятельности,... | | Развития российской федерации федеральная служба Земли промышленности, энергетики, транспорта, связи, радиовещания, телевидения, информатики, земли для обеспечения космической деятельности,... |
