在任何给定的时间，是否存在三个米兰科维奇因子的辐射强迫方程?还是价值数据库?
\ begingroup美元
\ endgroup美元
添加评论
|
2答案
\ begingroup美元
\ endgroup美元
0
有许多模型可以在长时间尺度上估计轨道参数。所有这些都依赖于相当复杂的公式，这些公式解释了其他行星、月球和许多其他因素产生的扰动。一些模型甚至考虑了相对论效应这个答案．
我经常看到在处理气候强迫的论文中引用和使用的一个轨道解是LA93，由Laskar等人(1993)计算，并在论文中进行了描述地球的轨道、进动量和日照量从-20 MYR到+10 MYR．然而，现在已经有了更精确的解决方案，比如LA2004和更著名的LA2010。我将在下面描述如何下载这些解决方案的输出数据和代码。
对于LA93，该模型的代码(fortran)和输出可在维齐尔目录，在这里您可以查询模型输出，或者直接进入这个FTP文件夹下载代码文件和每千年输出一次的模型长表。
如果“米兰科维奇的三个因素”你指的是斜度，偏心和进动指数(偏心$ \ * $Sin(近日点经度))。该数据可以下载并在以下两个文件中使用:
CLIVAR0N。ASC (-20Myr至0Myr)
CLIVAR0P。ASC(0密珥至+10密珥)
方法中描述的这些文件的逐字节描述自述是
-------------------------------------------------------------------------------- 字节格式单位标签的解释 -------------------------------------------------------------------------------- 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 * sin(经度的近日点移动equinox ) --------------------------------------------------------------------------------使用该代码，您还可以计算常见的气候代理，例如在65°N或任何其他纬度的隔离。
对于LA2004年(-50 Myr到+20 Myr)，同一研究小组实施了非常有用的在线计算器这样就可以计算出任何纬度的轨道和进动量的日晒量。它以表格形式输出结果。计算器界面看起来像这样，因此您可以了解可用选项。
本文介绍了LA2010解决方案La2010:地球长期运动的新轨道解决方案，并回到250亿美元以前。此解决方案的类似数据文件和代码可以在地球古气候的天文解决方案页面．根据我的理解，LA2010代表了最先进的轨道参数解决方案，它也包括相对论修正。
\ begingroup美元
\ endgroup美元
如果你想认真对待数据，我认为卡米洛·拉达的答案是值得遵循的。
然而，为了更直接地回答这个问题，如果你想使用实际的方程，那么如果你只研究全新世(也许是所有的第四纪?)，另一种方法应该也很好，那就是使用André Berger的原始表格
https://journals.ametsoc.org/view/journals/atsc/35/12/1520-0469_1978_035_2362_ltvodi_2_0_co_2.xml
在开头描述了这些方程，并将轨道力学“分解”为正弦曲线的和:
ε=ε* +∑Aicos(我)适合+δe = e0 +∑Ei cos(λ+φi) e sin(ω)=∑π罪(α+ζ我)
这些值是从论文中提供的表格中读取的。如果你只是随便玩玩，把每个轨道参数的前三项相加会得到很好的结果。
然而，正如另一个答案所说明的那样，即使最终结果看起来像一个相当简单的正弦曲线，实际上它有很多次级振荡项。你看一下表格，前5+项具有“显著”振幅。所以运动不只是几个有噪声叠加的主频。
再说一次，这些都很有趣，也很有教育意义，但如果你在做严肃的工作，那么你可能想要得到完整的Laskar表
