方程米兰柯维奇因素?

Question

方程有三米兰柯维奇辐射营力的因素在任何给定的时间?或数据库的值?

Answer 1

有许多模型估计轨道参数耗时过长。所有的依赖,而复杂的公式,占其他行星产生的扰动,月亮,和许多其他因素。甚至一些模型考虑相对论效应中描述这个答案。

一个轨道的解决方案,我经常看到在论文中引用和使用处理气候迫使LA93,计算拉斯卡尔et al .(1993)和描述岁差的轨道,地球日晒量从-20年最高产量研究+ 10最高产量研究。然而,更准确的解决方案是可用的现在,更喜欢LA2004和特别是LA2010。下面我将描述如何下载输出数据和编码这些解决方案。

LA93的代码(在fortran)和输出的模型是可用的维齐尔目录,你可以查询模型输出,或进入这个FTP文件夹和下载的文件代码和长表与模型输出每一千年。

如果“三米兰柯维奇因素”你的意思是倾斜,离心率和旋进指数(偏心$ \ * $罪(近日点经度))。数据已经准备好这两个文件下载和使用:

CLIVAR0N。ASC(-20最高产量研究0最高产量研究)
CLIVAR0P。ASC(0 + 10最高产量研究最高产量研究)

的byte-per-byte描述这些文件中描述自述是

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -单位标签字节格式解释- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 4 - 11 F8.0 1000年t时间从1000年来J2000 13 - 31 D19.12——e偏心33 - 51 D19.12 rad eps倾斜(弧度)53 - 71 D19.12 - - - CP e *罪(经度的近日点移动equinox) - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

与代码还可以计算常见的气候代理如隔离在65°N或任何其他纬度。

LA2004(-50 + 20最高产量研究最高产量研究)的研究团队实现了一个非常有用的在线计算器,允许计算日晒量来源于轨道和岁差的数量对任何感兴趣的纬度。以表格的形式输出结果。计算器的界面是这样的所以你了解可用的选项。

LA2010解决方案描述La2010:新地球的轨道运动的长期解决方案250年,最高产量研究。你可以找到类似的数据文件和代码为这个解决方案天文地球古气候解决方案页面。我的理解,LA2010轨道参数代表了最先进的解决方案,它包括相对论修正。

Answer 2

我认为Camilo Rada的回答是遵循如果你想做任何严重的数据。

然而,直接回答这个问题更多,如果你想用实际方程的另一种选择,应该很好如果你看着全新世(或许所有的第四纪?)从安德烈·伯杰使用原始的表

https://journals.ametsoc.org/view/journals/atsc/35/12/1520 - 0469 _1978_035_2362_ltvodi_2_0_co_2.xml

一开始描述的方程和“分解”的轨道力学的正弦曲线:

ε=ε* +∑Aicos(适合+δi) e = e0 +∑Ei cos(λit +φi) e sin(ω)=∑π罪(αit +ζi)

从表中读取这些值提供。如果你只是玩,甚至只是总结前3项为每个轨道参数会给你好看的结果。

然而,随着其他回答说明了,即使最后的结果看起来像一个非常简单的正弦信号,在现实中是有很多二次振动的条件。你看看表,第一个5 +的术语是“重大”振幅。因此,运动不仅仅是几个主导频率噪声叠加。

这是所有乐趣和教育,但如果你做严肃的工作,那么你可能想要完整的拉斯卡尔表