【米兰科维奇循环】
VSOP想要非常准确的预测过去和未来的轨道参数：ε是轨道倾角（转轴倾角），e是离心率，ω是近日点黄经，esin(ω)是进动指数，它与轨道倾角一起控制著日照的季节周期。  是夏至日那一天在北纬65度的计算所得的大气层顶日平均日照量。有孔虫和南极冰蕊显示出过去的全球海平面和温度的特征，这两个替代性指标分别来自海洋沉积物和南极冰芯。垂直的灰线是目前的西元2000年。
米兰科维奇循环（Milankovitch cycles）是塞尔维亚的地球物理学家兼天文学家米卢廷·米兰科维奇的气候变化理论。米兰科维奇计算了过去数百万年地球的离心率、转轴倾角和轨道的进动的变化，发现了这些参数与地球上气候模式，尤其是冰川期的关系。
地球轨道倾角的大约每26,000年完成绕行一周的完整进动周期。在这同时，椭圆轨道旋转也以缓慢的21,000年引导著季节和轨道之间的变化。另一方面，地球的自转轴和轨道平面之间的倾角以41,000年的周期在22.1度到24.5度之间摇摆著，现在的角度是23.44度，并且还在减少中。此运动称为章动。
其它还有Joseph Adhémar、James Croll和其他人提出先进的天文理论，但是仍有所疑虑，由于要和过去的时间完全确切结合是很重要的证据，因此很难得到验证。直到深海岩蕊和James Hays、John Imbrie和Nicholas Shackleton的沉积层报告：《地球轨道的变动：冰河期的定标》，发表在1976年的《科学》期刊，理论才呈现目前的状态。
== 地球的变动 ==
地球围绕自转轴自转和在轨道上绕著太阳公转的过程中，会有几个准周期的变化发生。运动曲线上虽然有大量的正弦成分，但仅有几个成分在主导。米兰科维奇研究轨道离心率、倾角、和进动的变化。这些在运动和方向上的变化改变了太阳辐射抵达地球的方向和数量，这是所知的太阳强迫作用（一个辐射强迫作用的例子）。北极附近地区由于有大量的陆地，变动被认为比较重要，这是因为陆地的反应比海洋快速。
=== 轨道形状（离心率） ===
地球的轨道是椭圆形，而离心率是测量椭圆与圆形的偏差。地球轨道的形状在接近圆形（低离心率的0.005）到轻度的椭圆（高离心率的0.058）之间变化，平均的离心率是0.028，此变化的主要周期是413,000年（离心率改变±0.012）。其它较主要的周期是是95,000年和125,000年（联结的周期是400,000年），而分异松散的周期是100,000年)。目前的离心率是0.0167。
如果地球是唯一环绕著太阳的行星，它的离心率即使经过数百万年也不会有感觉得到的轻微变化，地球轨道离心率的改变主要是受到木星和土星不同引力的交互作用影响。椭圆轨道的离心率虽然会改变，但椭圆轨道的半长轴不会改变。从摄动理论的观点，使用天体力学计算轨道的演化，半长轴是绝热不变量。依据克卜勒第三定律，轨道周期是由轨道半长轴测定的。因此，地球的轨道周期，恒星年，当轨道逐渐变化时，长度也仍然保持著不变。当半短轴随著离心率的增加缩短时，季节的变化会加剧，但根据克卜勒第二定律，行星的平均太阳辐射变化在低离心率时只有微量的变化。
同样的平均辐射与平均温度并不会有相对应的关联性（由于史蒂芬-波兹曼定律是非线性的）。一个与温度20°相对应的辐射，可以有±50 %的对称变化（例如，来自季节的变化）。我们观测到的温度变化对应于平均16℃（也就是说偏差有n-4℃）。并且在一天之中的辐射变化（仍然对应于平均温度20℃），我们观察到的平均温度是（对零热容量）-113℃。
在最靠近太阳时（近日点）相对增加的太阳辐射大约是在距离太阳最远时离心率的4倍。以目前的轨道离心率，这相当于增加6.8%入射的太阳辐射，而目前近日点和远日点的差异只有3.4%（510万公里）。现在通过近日点的日期大约在每年的1月3日，而经过远日点的日期大约是7月4日。当地球轨道最椭圆时，在近日点的太阳辐射量将比远日点时大23%。
轨道力学要求季节的长度与季节的象限领域成正比，因此在离心率的极端值时，在轨道远心点端的季节持续的时间也会大幅的增长。当秋天和冬天是在出现在近心端时，如同目前北半球的状态，地球在轨道上的移动速度是最快的，因此秋天和冬天会比春天和夏天稍短一些。像这样，夏天比冬天长4.66天，春天比秋天长2.9天。
=== 转轴倾角（倾角/黄道赤道交角）===
地轴倾斜的范围在22.1-24.5°。
地球的转轴倾角（倾斜）是地球的转轴相对于轨道平面的角度。角度变化的范围是2.4°，在大约41,000年的周期内从倾斜22.1°缓慢的变化至24.5°并且再复原。当倾角增加时，日照（进入的太阳辐射）在季节周期上的振幅也增加；在两个半球的夏季都会接收到更多的太阳辐射通量，而冬季的辐射通量减少。但是，这种冬夏两季的反相变化在地表各处幅度不尽相同。当倾角增加时，高纬度的全年日照量会增加，而低纬度接收的全年日照量会减少。
现在地轴倾角处在逐渐变小的趋势，使四季变化较和缓（暖冬及凉夏），因此有气候冷却的倾向。地球上大部分的冰雪位在高纬度，倾角减小时，两个半球在夏季接收更少的日照量，并且高纬度地区的日照量亦缩减，凉爽的夏天造成先前在冬天的冰雪融化量减少，促使冰河期的开始。但是，没有气候变化与轴倾角极端值变化之间有意义的关联。
科学家利用电脑模型研究更极端的地轴倾斜，比那些实际发生过更高的轴倾角，引发高纬度的极端气候，以个别的研究能否威胁目前已经存在地球上的高等生物型态的生命。它们注意到高倾角并不能完全的消除一颗行星的生物，但会使它更难以茁状和成长，而且也会让现今存在于土地上的温血生物衰弱.。
目前地球相对于轨道平面的倾角是23.44度，大约是在各个极端值的中间。倾角是在周期的减少阶段中，大约西元10,000年会达到其最小值。这一趋势的本身往往会造成温暖的冬季和凉爽的夏季，但是增强的温室气体可能会超越这种变化的影响。
=== 轴向进动 ===
进动的运动。
地球自转轴的方向相对于恒星的变化称为进动，周期大约是26,000年。这种陀螺的运动是由于太阳和月球对固体的地球，加上地球的形状是扁椭球而不是理想的球，所施加的潮汐力，而太阳和月球有著大约一致的效果。
当自转轴的方向在轨道的近日点朝向太阳时，一个极半球的季节有著较大的变化而另一个极半球的季节变化较为温和。在近日点时是夏季的半球，接收到的太阳辐射会相对应的增加，而这个半球在冬季也会相对的较为寒冷。另一个半球则会有较温暖的冬季和较为凉爽的夏季。
当地球的近日点和远日点是朝向分点时，北半球和南半球有著相似的季节分布状态。
在目前，当南半球的夏季时地球位于近日点，并且在远日点时是南半球的冬季。因此，当其它的因素都相同时，南半球的季节会比北半球的较为极端。
=== 拱线进动 ===
以椭圆轨道（卵型）环绕太阳的行星，其轨道会随着时间改变（拱线进动），图中的离心率在视觉上被夸大了。在太阳系内大多数天体的离心率都很小，几乎接近圆形。
拱线进动对季节的影响。
另一方面，椭圆轨道本身在空间中的行进，主要是受到木星和土星交互作用的结果。这种轨道进动和陀螺旋转轴的运动是有著相同的意义，会使分点岁差相对于近日点进动的周期从25,771.5年缩短至大约为21,636年。
=== 轨道倾角 ===
地球轨道的倾角会在现在的轨道平面上以70,000年的周期向上和向下漂移。米兰科维奇没有研究这种三维空间的运动，这种运动被称为黄道进动或是行星进动。
近期的研究人员注意到这种轨道的漂移也会相对于其他的行星轨道移动。Invariable plane，代表太阳系角动量的平面，大约就是木星的轨道平面。地球的轨道以100,000年的周期相对于不变的平面倾斜；很偶然的机会，这非常类似于100,000年的离心率周期。这个100,000年的周期与100,000年的冰河期周期模式相符合。
不变的平面中曾被认有尘埃和其它碎片形成的盘面，并且会通过几种可能的手段影响到地球的气候。地球目前在每年的1月9日和7月9日通过不变平面，雷达的侦测显示这时流星和相关的夜光云的数量都有增加
。
使用南极冰核中困住的气泡中气体的氢氮比率进行的年代学研究，这似乎可以直接反应当地的日照，从冰核中记录的气候对应出北半球日照的结论，一如米兰科维奇的假说（Kawamura et al., Nature, 23 August 2007, vol 448, p912-917）。这是验证米兰科维奇假说追加的一个新方法，并且与100,000年周期的"倾斜"理论是不吻合的。
== 问题 ==
因为观测到的气候周期正好符合了轨道周期，使轨道理论获得了压倒性的支持。然而，在理论的满足和观测上还是有几个困难。
沉积物的性质以循环的方式在改变，从沉积的纪录中可以显示这些周期。此处，循环可以从不同的沉积层的颜色和强度观测出来。
===十万年的问题 ===
100000年的问题是，按照理论，离心率的变化，相比于与太阳强迫有关的转轴倾角变化或轴向进动变化，对气候的影响较小。但是同位素分析表明，气候变化的主导周期为约10万多年。
== 参考文献 ==
=== 引用 ===
=== 来源 ===
* 王干盈、王作台 等.《基础地球科学》. 台北：全华图书股份有限公司，1996年：p. 145. 2008年7月4日. ISBN 978-957-21-5808-1.
* English Wikipedia， Milankovitch cycles date=20210418105807 ，2008年6月30日
== 延伸读物 ==
* cite journal | author=Roe G | title=In defense of Milankovitch | journal=Geophysical Research Letters | volume=33 | year=2006 | pages=L24703 | doi=10.1029/2006GL027817  This shows that Milankovitch theory fits the data extremely well, over the past million years, provided that we consider derivatives.
*  pmid=11326091  This review article discusses cycles and large-scale changes in the global climate during the Cenozoic Era.
== 外部连结 ==
*  Milankovitch Cycles and Glaciation
*  The Milankovitch band
*  Some history of the adoption of the Milankovitch hypothesis (and an alternative) date=20080627174136
*  More detail on orbital obliquity also matching climate patterns date=20140827150228
*  Graph of variation in insolation date=20120112043906   Note 20,000 year, 100,000 year, and 400,000 year cycles are clearly visible.
*  dead-url=no
*  Potential Problems with Milankovitch Theory date=20080521053155   by Sean Pitman date=20080701191459
*  The Seasons date=20130802165648
*  The NOAA page on Climate Forcing Data date=20080430171840  includes (calculated) data on orbital variations over the last 50 million years and for the coming 20 million years.
*  The orbital simulations by Varadi, Ghil and Runnegar (2003) provide another, slightly different series for orbital eccentricity, and also a series for orbital inclination
* ABC: Earth wobbles linked to extinctions date=20081118050429
分类:地理术语