In nature, the planetary albedo mainly depends on the ice extent, cloud cover and land surface properties.
In SimClimat, only the effect of ice sheet extent is taken into account. The albedo is calculated as a function of ice sheet extent by a piece-wise linear function. The albedo is bounded between the albedo of the ice (taken at 0.9) and the albedo of the Earth without ice, taken at 0.25. This formulation of the albedo as a function of the latitude of ice sheets, which itself depends on the temperature (section 5.4.2, explains the shape of the curve (the solar energy absorbed by the Earth) as a function of temperature in figure 9.
The latitude of the ice sheets is in degrees of latitude. It is calculated as a function of global temperature and of the summer insolation at 65°N, (in order to take into account the variations of orbital parameters).
We calculate the ice sheet extent at equilibrium :
is calculated as a function of the solar constant, eccentricity, obliquity and precession (section 5.4.3).
The parameters , and are tuned to satisfy the constraints summarized in section 2.1: =0.73, =49.53 and =0.2.
The ice sheets respond to climate forcing with a time scale = 3000 years. To represent this effect, the ice sheet latitude is calculated as a function of assuming that tends towards with the time constant :
The summer insolation at 65°N, , is calculated as a function of the solar constant , eccentricity , obliquity and precession following this formula:
where and the present-day eccentricity and precession. Angles and are given in °.