The astronomical theory of the paleoclimate serves for interpreting the paleoclimate flactuations and,especially,the Ice Age problem basing on the variations of astronomical orbital elements.It initiated in nineteenth century, but was formally astablished in 1941 when Milankovitch’s Kanon was published.Milankovitch emplasized the importance of caloric summer insolation on glaciation,i.e.,the cool and short summer is in favour of the development of glaciers.His computation comfirmed his idea.Unfortunately,his theory was argued for many years because of the insufficient geological records dated with poor accurant techniques.It was not until the implementation of CLIMAP(Climate Mapping Project) that found the paleoclimate oscillations really occurred in Milankovitch frequency belts.Meanwhile,A.Berger,a Belgian astr geoscientist,re calculated long term variation of the insolation with simpler method and more accuracy than ever before.The insolation he computed has even been dated up to past 10 000 000 years or more.His spectral analysis for the insolation clearly shows the 400 and 100 41 23 and 19 ka cycles corresponding to the astronomical elements of the eccentricity,obliquity and astronomical precession,respectively.In addition,his calculation scheme could be easily applied in numerical climate model as an orbital forcing which drives the general circulation of atmosphere. On the other hand,numerical simulation of paleoclimate is rapidly developing.A sectorially averaged two dimensional climate model like LLN 2D model consisting of coupled atmospheric,ocean and ice sheet models has been built,it can be run from remote geological past up to now,and some of which results may be compared directly with geological records or other proxy data,for instance,the δ 18 O,so as to re exame the Milankovitch’s theory as well as the impacts of other factors in paleoclimate.